https://blog.csdn.net/tanzuozhev/article/details/51106089
tg <- ToothGrowth
head(tg)
## len supp dose
## 1 4.2 VC 0.5
## 2 11.5 VC 0.5
## 3 7.3 VC 0.5
## 4 5.8 VC 0.5
## 5 6.4 VC 0.5
## 6 10.0 VC 0.5
# summarySE 计算标准差和标准误差以及95%的置信区间.
tgc <- summarySE(tg, measurevar="len", groupvars=c("supp","dose"))
tgc
## supp dose N len sd se ci
## 1 OJ 0.5 10 13.23 4.459709 1.4102837 3.190283
## 2 OJ 1.0 10 22.70 3.910953 1.2367520 2.797727
## 3 OJ 2.0 10 26.06 2.655058 0.8396031 1.899314
## 4 VC 0.5 10 7.98 2.746634 0.8685620 1.964824
## 5 VC 1.0 10 16.77 2.515309 0.7954104 1.799343
## 6 VC 2.0 10 26.14 4.797731 1.5171757 3.432090
以上是例子。
下面是我得部分数据
head(qq)
pl wz dl
1 44.44 SGPD 0.1
2 44.44 SGPD 0.1
3 25.00 SGPD 0.1
4 33.33 SGPD 0.1
5 24.15 SGPD 0.1
6 44.44 SGPD 0.5
但是当我用summarySE求我数据的sd和se时发现,我的均值都没有求出来,这是为什么。
QQ <- summarySE(qq, measurevar="pl", groupvars=c("wz","dl"))
错误:There were 12 warnings (use warnings() to see them)
我的代码都是照着例子输入的,出现了这种到底是怎么回事?
有人知道吗


雷达卡





京公网安备 11010802022788号







