Displaying 1 result from an estimated 1 matches for "f_prob".
2008 Jan 31
1
difficulties computing a simple anova
...an(spiderdata[j])) ^ 2)
}
sqdata
N = 0
for(j in 1:J) {
N = N + length(sqdata[[j]])
}
SSW = sum(sqdata)
MSW = SSW / (N - J)
SSB = 0
for(j in 1:(length(spiderdata))) {
SSB = SSB + length(spiderdata[[j]]) * ((mean(spiderdata[j])[[1]] -
(sum(spiderdata) / N)) ^ 2)
}
MSB = SSB / (J - 1)
F = MSB / MSW
f_prob = pf(F, J - 1, N - J)
reject_point = qf(.95, J - 1, N - J)
cat("SSW:", SSW, ", MSW:", MSW, ", SSB:", SSB, ", MSB:", MSB, ", F:", F, ",
P(F):", f_prob, ", P(", reject_point, ") = .95\n", sep = "")
anova(lm(Zo...