怎样用aggregate估算数据的斜度

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我写了个算斜度的function, 要求要用这个:

我写的：
不过我不知道对不对，n就直接可以写n吗？可是n是我数据的大小，是不是要把n写成dataset的名字？

1. myskewness<-function(x){return(sum((x-mean(x))^3/sqrt(var(x))^3))*(n/((n-1)*(n-2)))}

复制代码

现在我要用aggregate这个命令和我写的function估算所有weight和sex的组合的斜度，
我的data是这样的：

   Dosis weight sex
1      0    6.62          m
2      0    6.65          m
3      0    5.78          m
4      0    5.63          m
5      0    6.05          m
6      0    6.48          m
7      0    5.50          m
8      0    5.37          m
9      1    6.25          m
10     1    6.95          m
11     1    5.61          m


从来没用过aggregate，求请教

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Aggregate Gate GRE ATE REG function 我不知道

1. myskewness<-function(x){
   
2.   n=length(x)
   
3.   return(sum((x-mean(x))^3/sqrt(var(x))^3))*(n/((n-1)*(n-2)))}
   
4. aggregate(data[,2],by=list(data$weight),myskewness)

复制代码

水天一色DIY 发表于 2013-5-14 16:29

weight和sex的所有组合不需要写sex吗？

kiotoqq 发表于 2013-5-14 21:20
weight和sex的所有组合不需要写sex吗？

你如果想根据sex的话，就把上面代码改成sex即可

thanks for sharing