Guest User

Untitled

a guest
Jun 17th, 2018
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.72 KB | None | 0 0
  1. df <- data.frame(time = c(10,20,30,50,100,200,500,1000),
  2. A = c(1.5,2.5,2.8,4.5,5.4,6.2,7.4,8.1),
  3. B = c(1.6,2.6,2.7,4.4,4.9,6.1,7.5,8),
  4. C = c(1.4,2.4,2.6,4.2,4.9,5.9,7.8,7.9))
  5.  
  6. melted_df <- data.table::melt(df, id.vars = c("time"), variable.name = c("group"))
  7.  
  8. av <- aov(value~group+time+group:time, data = melted_df)
  9. av
  10.  
  11. Df Sum Sq Mean Sq F value Pr(>F)
  12. group 2 0.11 0.05 0.025 0.975
  13. time 1 83.26 83.26 39.506 6.32e-06 ***
  14. group:time 2 0.02 0.01 0.005 0.995
  15. Residuals 18 37.93 2.11
  16. ---
  17. Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Add Comment
Please, Sign In to add comment