咨询一个R的小白问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

新人刚开始学习，教科书上有个地方不是很明白
> mystats<-function(x,na.omit=FALSE){
if(na.omit)
x<-x[!is.na(x)]
m<-mean(x)
n<-length(x)
s<-sd(x)
skew<-sum((x-m)^3/s^3)/n
kurt<-sum((x-m)^4/s^4)/n-3
return(c(n=n,mean=m,stdev=s,skew=skew,kurtosis=kurt))}

> sapply(mtcars[vars],mystats)

或者> sapply(mtcars[vars],mystats，na.omit=TURE)
前3行代码看的自己的逻辑比较混乱，只知道大概是处理缺失值用的，大神能帮着分析下思路不？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：小白问题 白问题 stats 教科书 教科书

1. mystats<-function(x,na.omit=FALSE){
   
2. if(na.omit) #进行if判定，na.omit=TRUE,代表你知道x变量你有缺失值，如果为FALSE,就不执行。
   
3. x<-x[!is.na(x)] # 有这一句，其它不需要上一句的，is.na是判定是否为缺失值，NA显示为TRUE,！是取非，!is.na(x)是显示的非缺失值的TRUE or FALSE状态，这一句就是把缺失值去掉了
   
4. m<-mean(x) # 均值
   
5. n<-length(x) # 样本量
   
6. s<-sd(x) # 标准差
   
7. skew<-sum((x-m)^3/s^3)/n #偏度
   
8. kurt<-sum((x-m)^4/s^4)/n-3 # 峰度
   
9. return(c(n=n,mean=m,stdev=s,skew=skew,kurtosis=kurt))} # 返回data.frame

复制代码

jiangbeilu 发表于 2016-4-12 08:37

明白了，非常感谢