install.packages('xslx')
library(xlsx)
Sys.setlocale("LC_ALL", "zh_cn.utf-8")
a=read.xlsx2('d:/1.xlsx',1,header=F)
2.计算平均值
install.packages('Rmisc')
library(Rmisc)
单因素:求y的平均值,以x为分类变量
a1=summarySE(a,measurevar='y', groupvars=c('x'))
双因素:求y的平均值,以x1和x2为分类变量
a1=summarySE(a,measurevar='y', groupvars=c('x1','x2'))
如果是简单数据,数据集a中有变量:只有因变量y,没有自变量x。求y的平均值代码如下:
mean(a$y)
sma<-function(x,n){ #x是数组,求n项平均sapply(1:length(x),function(i){
if(i>=n) mean(x[(i-n+1):i])
else mean(x[1:i])
})
}
用的时候就直接
sma(x,3)就行了
mean_24 <- function(x){k=NULL
for(i in 24:length(x))
{k[i-23]=sum(x[seq(i-23,i)])/24}
return(k)
}