Displaying 1 result from an estimated 1 matches for "cfev".
Did you mean:
cdev
2010 Dec 22
1
Seeking feedback on my first attempt at R programming
...o")))
fev <- within(fev,chg <- fev6 - fev0)
summaryfn <- function(x)
data.frame(n=sum(complete.cases(x)),mean=mean(x),std.dev=sd(x),
t.value=mean(x)/(sd(x)/sqrt(sum(complete.cases(x)))),
p.value=2*pt(-abs(mean(x)/(sd(x)/sqrt(sum(complete.cases(x))))),df=sum(complete.cases(x))-1))
cfev <- na.omit(fev)
by(cfev[3:5],cfev[2],summaryfn)
fligner.test(chg ~ trtgrp, data=fev)
t.test(chg ~ trtgrp, data=fev, var.equal=TRUE)
#summaryfn <- function(x)
# data.frame(n=sum(complete.cases(x)),mean=mean(x),std.dev=sd(x),
# t.value=t.test(x)$statistic,
# p.value=t.test(x)$p.valu...