- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 0 个
- 通用积分
- 0.0089
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 462 点
- 帖子
- 39
- 精华
- 0
- 在线时间
- 7 小时
- 注册时间
- 2012-8-19
- 最后登录
- 2018-2-2
高中生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
 - 0 个
- 通用积分
- 0.0089
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 462 点
- 帖子
- 39
- 精华
- 0
- 在线时间
- 7 小时
- 注册时间
- 2012-8-19
- 最后登录
- 2018-2-2
 | 开心 2014-4-25 22:33:59 |
|---|
签到天数: 1 天 连续签到: 1 天 [LV.1]初来乍到
|
经管之家送您一份
应届毕业生专属福利!
求职就业群
感谢您参与论坛问题回答
经管之家送您两个论坛币!
+2 论坛币
最近下载了这个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]]
复制代码
|
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
|
|
|