农业数据中如何计算小区多株的平均值（R语言分享）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

问题描述

本来想着小区内多株是可以分析的，发现要用平均值代替，数据就已经按列拍好了，怎么办？

用R语言

知道R语言可以做，可是不知道怎么做，刚才看到了如何做，就迫不及待记录。

1. ```
   
2. > dat <- data.frame(cul=rep(rep(1:5,each=3),4),block=rep(1:4,each=15),
   
3. +                   plot=rep(1:3,20),y=rnorm(60,30,10))
   
4. # 生成数据
   
5. > head(dat)
   
6.   cul block plot        y
   
7. 1   1     1    1 29.92395
   
8. 2   1     1    2 47.77084
   
9. 3   1     1    3 18.61392
   
10. 4   2     1    1 43.67827
    
11. 5   2     1    2 43.29565
    
12. 6   2     1    3 33.36473
    
13. 
    
14. # 错误的模型，这里残差是错误的
    
15. > # wrong model, the df of residuals is wrong!
    
16. > for (i in 1:3) dat[,i] <- as.factor(dat[,i])
    
17. > str(dat)
    
18. 'data.frame':        60 obs. of  4 variables:
    
19. $ cul  : Factor w/ 5 levels "1","2","3","4",..: 1 1 1 2 2 2 3 3 3 4 ...
    
20. $ block: Factor w/ 4 levels "1","2","3","4": 1 1 1 1 1 1 1 1 1 1 ...
    
21. $ plot : Factor w/ 3 levels "1","2","3": 1 2 3 1 2 3 1 2 3 1 ...
    
22. $ y    : num  29.9 47.8 18.6 43.7 43.3 ...
    
23. > mod <- aov(y ~ cul + block,data=dat)
    
24. > summary(mod)
    
25.             Df Sum Sq Mean Sq F value Pr(>F)
    
26. cul          4    335   83.80   1.054  0.389
    
27. block        3     30   10.16   0.128  0.943
    
28. Residuals   52   4134   79.50               
    
29. 
    
30. #查看数据
    
31. > head(dat)
    
32.   cul block plot        y
    
33. 1   1     1    1 29.92395
    
34. 2   1     1    2 47.77084
    
35. 3   1     1    3 18.61392
    
36. 4   2     1    1 43.67827
    
37. 5   2     1    2 43.29565
    
38. 6   2     1    3 33.36473
    
39. > library(reshape2)
    
40. > dd <- dcast(dat, cul+block~plot)
    
41. Using y as value column: use value.var to override.
    
42. 
    
43. # 转换后的数据
    
44. > head(dd)
    
45.   cul block        1        2        3
    
46. 1   1     1 29.92395 47.77084 18.61392
    
47. 2   1     2 26.82941 28.22210 28.30006
    
48. 3   1     3 25.04417 33.55550 18.65392
    
49. 4   1     4 29.47895 28.04065 23.50930
    
50. 5   2     1 43.67827 43.29565 33.36473
    
51. 6   2     2 16.27698 28.26213 38.50232
    
52. 
    
53. # 求3个观测值的平均数
    
54. > dd$mean <- (dd[b]问题描述[/b][align=left]本来想着小区内多株是可以分析的，发现要用平均值代替，数据就已经按列拍好了，怎么办？[/align][b]用R语言[/b][align=left]知道R语言可以做，可是不知道怎么做，刚才看到了如何做，就迫不及待记录。[/align]1` +dd[b]问题描述[/b][align=left]本来想着小区内多株是可以分析的，发现要用平均值代替，数据就已经按列拍好了，怎么办？[/align][b]用R语言[/b][align=left]知道R语言可以做，可是不知道怎么做，刚才看到了如何做，就迫不及待记录。[/align]2` + dd[b]问题描述[/b][align=left]本来想着小区内多株是可以分析的，发现要用平均值代替，数据就已经按列拍好了，怎么办？[/align][b]用R语言[/b][align=left]知道R语言可以做，可是不知道怎么做，刚才看到了如何做，就迫不及待记录。[/align]3`)/3
    
55. > dd <- dd[,-(3:5)]
    
56. > head(dd)
    
57.   cul block     mean
    
58. 1   1     1 32.10291
    
59. 2   1     2 27.78386
    
60. 3   1     3 25.75120
    
61. 4   1     4 27.00964
    
62. 5   2     1 40.11288
    
63. 6   2     2 27.68048
    
64. 
    
65. # 正确的模型
    
66. > mod1 <- aov(mean ~ cul + block,data=dd)
    
67. > summary(mod1)
    
68.             Df Sum Sq Mean Sq F value Pr(>F)
    
69. cul          4  111.7   27.93   0.692  0.612
    
70. block        3   10.2    3.39   0.084  0.968
    
71. Residuals   12  484.4   40.37               
    
72. ```

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：农业数据 R语言 小区

我用markdown语法，然后点击快捷键发现可以了，原来可以支持markdown啊，高兴中，然后提交，发现发表后成乱码了。。。。。