| 查看: 1559 | 回复: 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]] |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
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
200