R语言 每列每十个数求一次平均值

Python011

R语言 每列每十个数求一次平均值,第1张

#把年份取出来，按10年分成一组

decade <- ceiling(Temp_m[1] / 10)

names(decade) <- "decade" 

#把分组信息合并成原数据框

Temp_m <- cbind.data.frame(Temp_m, decade) 

#加载plyr包。

#把Temp_m这个数据框中除year列的数据外，按照decade分组分别把每列求算术平均值。

#并放在Temp_mean10数据框中。

#Temp_mean10含有decade信息,进行索引时需要注意。

library(plyr)

Temp_mean10 <- ddply(Temp_m[-1], .(decade), colwise(mean))

mean_24 <- function(x){

k=NULL

for(i in 24:length(x))

{k[i-23]=sum(x[seq(i-23,i)])/24}

return(k)

}