r语言非线性分位数回归出现问题，有奖解答

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

r语言做非线性分位数回归出现 Error in rq.fit.br(x, y, tau = tau, ...) : Singular design matrix 此外: Warning messages:1: In lsfit(x, y, d^2, intercept = FALSE) : 'X' matrix was collinear
2: In lsfit(x, y, d^2, intercept = FALSE) : 'X' matrix was collinear
请问该如何改进？这个错误有时候出现，有时候不出现。若能解答，我会给奖励

1. n<-50
   
2. dat<-rsn(n,0,1,-1.8)
   
3. dat<-sort(dat)
   
4. f1<-psn(dat,0,1,-1)
   
5. sncdf<-function(x,alpha){
   
6.   n<-length(x)
   
7.   re<-rep(0,n)
   
8.   for(i in 1:n){
   
9.     snpdf<-function(x,alpha){
   
10.       p1<-dnorm(x,0,1)
    
11.       p2<-pnorm(alpha*x,0,1)
    
12.       2*p1*p2
    
13.     }
    
14.     re[i]<-integrate(snpdf,alpha,lower = -Inf,upper = x[i])$value
    
15.   }
    
16.   re
    
17. }
    
18. reg<-function(n,alpha){
    
19.   m=100
    
20.   res<-matrix(0,m,3)
    
21.   for(i in 1:m){
    
22.     ef<-(1:n-0.5)/n
    
23.     da<-rsn(n,0,1,alpha)
    
24.     da<-sort(da)
    
25.     fit1<-nlrq(ef~sncdf(da,alpha),start=list(alpha=0),tau = 0.5)
    
26.     fit1<-summary(fit1)$coefficients
    
27.     fit2<-nlrq(ef~sncdf(da,alpha),start=list(alpha=0),tau = 0.75)
    
28.     fit2<-summary(fit2)$coefficients
    
29.     fit3<-nlrq(ef~sncdf(da,alpha),start=list(alpha=0),tau = 0.85)
    
30.     fit3<-summary(fit3)$coefficients
    
31.     res[i,]<-c(fit1[1],fit2[1],fit3[1])
    
32.   }
    
33.   res
    
34. }  
    
35. 
    
36. paras<-seq(-5,5,by=0.5)
    
37. paras<-paras[paras!=0]
    
38. nn<-length(paras)
    
39. nres<-matrix(0,nn,6)
    
40. for(i in 1:nn){
    
41.   a<-reg(50,paras[i])
    
42.   jz<-apply(a,2,mean)
    
43.   fc<-apply(a,2,var)
    
44.   nres[i,1:3]=abs((jz-paras[i])/paras[i])
    
45.   nres[i,4:6]=fc+(jz-paras[i])^2
    
46.   print(i)
    
47. }

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：分位数回归 分位数 非线性 R语言 Intercept 非线性回归

其中rsn（）在sn包中，nlrq（）在[size=13.333333015441895px]quantreg包

请问问题解决了吗？可以求一份吗

您好，请问问题解决了嘛，我也遇到了相同的问题，能否帮忙解答一下