在摘要输出中更改“N”

时间:2018-01-15 22:11:39

标签: r summary weighted-average hmisc

我的summarySE输出如下所示:

summary <- summarySE(df, measurevar = "Cover", groupvars = c("Year","Loc")) 

  Year   Loc N     Cover      sd         se         ci
1 2010    1 104 0.10560533 0.2163059 0.02121054 0.04206610
2 2010    2  97 0.06513590 0.1505460 0.01528563 0.03034173
3 2010    3  48 0.12428977 0.2031627 0.02932401 0.05899230
4 2011    1  90 0.09785354 0.2184823 0.02303006 0.04576023
5 2011    2  88 0.04467975 0.1591636 0.01696690 0.03372355
6 2011    3  45 0.06313131 0.1494173 0.02227382 0.04488992

问题是,由于我的数据框的布局,我必须在采样期间为同一个样本包含多行,因为多个物种占据相同的样方。

E.g。

Quadrat_no. Year Site Cover
    1       2010   1    10
    2       2010   1     5
    2       2010   1     6
    3       2010   1     7  etc.

因此,R将Cover中的所有值相加,并将其除以(例如104)以计算加权平均值(例如0.105),而不是真实样本量。

是否可以更改我的数据帧,以便R识别多行代表单个样方?因此,更改summarySE中的N.

编辑:我后来发现:https://cran.r-project.org/web/packages/Hmisc/Hmisc.pdf我可以使用normwt,虽然我无法弄清楚如何实现它,因为没有例子。

0 个答案:

没有答案