| 查看: 1363 | 回复: 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]] |
回复此楼» 猜你喜欢
最失望的一年
已经有16人回复
-
存款400万可以在学校里躺平吗
已经有33人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
求助一下有机合成大神
已经有3人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
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