楼主: 盛1234567
2071 0

[问答] 新手求问packge手册里的lazydata,linking to,RoxygenNote等是什么意思? [推广有奖]

  • 0关注
  • 0粉丝

高中生

40%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
0.0089
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
462 点
帖子
39
精华
0
在线时间
7 小时
注册时间
2012-8-19
最后登录
2018-2-2

楼主
盛1234567 发表于 2017-4-6 21:38:40 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
最近下载了这个package: https://cran.r-project.org/web/packages/MAVE/index.html

想问问大家Reference manual提到的这些描述都是什么意思?
在使用其中的MAVE函数的时候R(version 3.3.3)总是会崩溃,下面是code,麻烦大家帮忙看看时怎么回事。

  1. library(mvtnorm)
  2. library(Rcpp)
  3. library(RcppArmadillo)
  4. library(knitr)
  5. library(MAVE)
  6. n<-400; p<-10        
  7. sigma<-diag(p)
  8. x<-rmvnorm(n,rep(0,p),sigma)
  9. beta1<-c(rep(1,2),rep(0,p-2))
  10. beta2<-c(rep(0,2),rep(1,2),rep(0,p-4))
  11. beta3<-c(rep(0,4),rep(-1,2),rep(0,p-6))
  12. beta4<-c(rep(0,6),1,-1,rep(0,p-8))
  13. beta<-cbind(beta1,beta2,beta3,beta4)
  14. prob<-exp(x%*%beta)
  15. y<-apply(prob,1,which.max)
  16. c1<-which(y==1);n1<-sum(y==1);x1<-x[c1,]
  17. c2<-which(y==2);n2<-sum(y==2);x2<-x[c2,]
  18. c3<-which(y==3);n3<-sum(y==3);x3<-x[c3,]
  19. c4<-which(y==4);n4<-sum(y==4);x4<-x[c4,]
  20. y12<-c(rep(1,n1),rep(2,n2));x12<-rbind(x1,x2)
  21. y13<-c(rep(1,n1),rep(3,n3));x13<-rbind(x1,x3)
  22. y14<-c(rep(1,n1),rep(4,n4));x14<-rbind(x1,x4)
  23. dr12<-MAVE(x12,y12,method="csMAVE" )
  24. theta22<-dr12$dir[[1]]
  25. dr13<-MAVE(x13,y13,method="csMAVE" )
  26. theta33<-dr13$dir[[1]]
  27. dr14<-MAVE(x14,y14,method="csMAVE" )
  28. theta44<-dr14$dir[[1]]
复制代码


二维码

扫码加我 拉你入群

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

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


您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-31 08:04