인쇄 중요한 결과

Question

I는 다음과 같습니다 데이터가 있습니다

 KAFFE SAFT FERSKVARE HERMETIKK ..... K 
     35 23.91  17.3   70.2   1 
     53 30.13  27.4   75.7   1 
     35 34.54  17.3   81.1   2 
     32 31.93  47.3   92.7   2 
     34 54.69  17.2   88.4   2 
     35 34.69  19.1   81.5   3 
     54 25.63  18.3   92.1   3 
     34 46.65  18.6   76.6   3 
     36 38.21  23.3   79.9   3 

내가 분산 분석을 만들고 각 데이터에 대한 특별 검사를 게시하는 루프를 사용을

rcluso is my data frame 
K is the categorical variable located at the end column 



sink("AnovaTest.doc") 
     for(i in 1:(ncol(rcluso)-1)){ 
     columns <- names(rcluso[i]) 
     anovaresult<- summary(aov(rcluso[,i]~K,data=rcluso)) 
     posthocresult <- TukeyHSD(aov(rcluso[,i]~K,data=rcluso)) 

    print(columns) 
    print(anovaresult) 
    print(posthocresult) 
    } 
    sink() 

문제는 인쇄 할 anova의 중요한 결과를 얻고 싶을 뿐이고 각 중요한 결과에 대한 막대 그림을 만들고 싶습니다.

여러분이 도와 드릴 수 있습니까? 미리 감사드립니다.

Answer 1

summary() 대신 anova() 출력을 얻으면 p 값에 액세스 할 수 있습니다.

대체 당신의 당신은 조건부 예는 IF() 문으로 결과를 인쇄 할 수 있습니다

anovaresult<- summary(aov(rcluso[,i]~K,data=rcluso)) 

다음

anovaresult<- anova(aov(rcluso[,i]~K,data=rcluso)) 

로

if(anovaresult$Pr[1] < 0.05){ print(anovaresult) }