Displaying 1 result from an estimated 1 matches for "otpt".
Did you mean:
opt
2012 Feb 09
1
Tukey HSD
...ODE: ( which occurs after TEST <- read.table("test_text.txt") )
i <- 2;
sink (file = "test_output.txt", append = FALSE)
mydf <- data.frame(TEST)
for (j in 1:ncol(mydf)-1) {
var1 <- mydf[,1]
var2 <- mydf[,i]
fm1 <- aov(var1 ~ var2)
tky <- TukeyHSD(fm1)
otpt <- capture.output(summary(fm1))
i <- i+1;
lines <- as.vector(unlist(strsplit(otpt[2]," ")),mode="list") # gets the
p-value
if (grepl("[1234567890]",lines[14],perl = TRUE)) { #make sure that slot for
p-value has a number
number <- as.numeric(lines[14])...