| 查看: 1390 | 回复: 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]] |
回复此楼» 猜你喜欢
今年春晚有几个节目很不错，点赞！
已经有6人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有3人回复
-
球磨粉体时遇到了大的问题，请指教！
已经有15人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有5人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
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