0

嗨,我是 r 新手,我必须在下面解决这个问题

将每个组(因子水平)内的最大值和最小值与其各自的组均值进行比较。您的一个值与其组均值之间最大的绝对差异是多少?假设数据是正态分布的,并且以各自的组平均值为中心,标准差为 1,那么获得这样一个值的机会是多少?

数据集和框架由

  fact<-rep(c("E","F","G","H"),each=12)
  variable2=rnorm(48,10)*(rep(rpois(4,.2),each=12)/8+1)
  ds<-data.frame(fact,variable2)

任何帮助将不胜感激

这是我尝试过的

library(“plyr”)
ddply(ds,~fact,summarise,maximum=max(variable2),min=min(variable2),mean=mean(variable2))
4

1 回答 1

1

你快到了。该dnorm功能将在这里为您提供帮助

res <- ddply(ds, ~fact , 
                summarise , 
                maxi = max(variable2) - mean(variable2),
                mini = min(variable2) - mean(variable2) )

res$probmax <- dnorm( res$maxi )
res$probmin <- dnorm( res$mini )
#  fact      maxi      mini    probmax      probmin
#1    E 1.7736537 -1.622157 0.08275571 0.1070311818
#2    F 1.7733593 -2.269254 0.08279894 0.0303883803
#3    G 2.6621257 -3.708242 0.01153470 0.0004120085
#4    H 0.8461922 -1.749625 0.27888407 0.0863339664
于 2013-06-10T10:08:40.610 回答