| 查看: 1414 | 回复: 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]] |
回复此楼» 猜你喜欢
288求调剂 一志愿哈工大 材料与化工
已经有31人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有5人回复
-
085801 总分275 本科新能源 求调剂
已经有17人回复
-
324求调剂
已经有6人回复
-
学硕化学工程与技术，一志愿中国海洋大学320+求调剂
已经有7人回复
-
材料调剂
已经有5人回复
-
初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文
已经有6人回复
-
化学308分调剂
已经有15人回复
-
346求调剂
已经有3人回复
-
环境285分，过六级，求调剂
已经有6人回复
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
