| 查看: 1528 | 回复: 3 | ||
[求助]
新手求问packge手册里的lazydata，linking to，RoxygenNote等是什么意思？ 已有1人参与
|
|
最近下载了这个package： https://cran.r-project.org/web/packages/MAVE/index.html 想问问大家Reference manual提到的这些描述都是什么意思？ 在使用其中的MAVE函数的时候R（version 3.3.3）总是会崩溃，下面是code，麻烦大家帮忙看看时怎么回事。 library(mvtnorm) library(Rcpp) library(RcppArmadillo) library(knitr) library(MAVE) n<-400; p<-10 sigma<-diag(p) x<-rmvnorm(n,rep(0,p),sigma) beta1<-c(rep(1,2),rep(0,p-2)) beta2<-c(rep(0,2),rep(1,2),rep(0,p-4)) beta3<-c(rep(0,4),rep(-1,2),rep(0,p-6)) beta4<-c(rep(0,6),1,-1,rep(0,p-8)) beta<-cbind(beta1,beta2,beta3,beta4) prob<-exp(x%*%beta) y<-apply(prob,1,which.max) c1<-which(y==1);n1<-sum(y==1);x1<-x[c1,] c2<-which(y==2);n2<-sum(y==2);x2<-x[c2,] c3<-which(y==3);n3<-sum(y==3);x3<-x[c3,] c4<-which(y==4);n4<-sum(y==4);x4<-x[c4,] y12<-c(rep(1,n1),rep(2,n2));x12<-rbind(x1,x2) y13<-c(rep(1,n1),rep(3,n3));x13<-rbind(x1,x3) y14<-c(rep(1,n1),rep(4,n4));x14<-rbind(x1,x4) dr12<-MAVE(x12,y12,method="csMAVE" ) theta22<-dr12$dir[[1]] dr13<-MAVE(x13,y13,method="csMAVE" ) theta33<-dr13$dir[[1]] dr14<-MAVE(x14,y14,method="csMAVE" ) theta44<-dr14$dir[[1]] |
回复此楼» 猜你喜欢
朋友圈看到的
已经有4人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
2楼2017-04-07 07:36:51
3楼2017-04-07 19:29:09
秋犯
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 578.8
- 沙发: 7
- 帖子: 331
- 在线: 37.5小时
- 虫号: 3427738
- 注册: 2014-09-19
- 性别: GG
- 专业: 环境污染化学
4楼2017-04-07 21:08:27
