J'ai un résultat avec une censure à droite comme ceci:
y<-c(rep(2.83,3),
rep(3.17,4),
rep(3.83,4),
rep(4.17,5),
rep(4.83,8),
rep(5.5,3),
rep(7.17,5),
rep(8.17,7),
rep(8.83,12),
rep(9.5, 12),
rep(9.83,17),
rep(10.17,30),
rep(10.50,100))
où y=10.5
sont les bonnes valeurs de censure. Ensuite, j'essaierais d'utiliser quantreg::crq
pour ajuster un modèle de régression quantile censuré et commencer avec une variable d'intervention binaire:
set.seed(123)
require(quantreg)
yc<-rep(10.5, length(y))
treat<-rbinom(length(y), 1, 0.5)
age<-as.integer(rnorm(length(y), 50, 2))
modèle1
fit1<-crq(Curv(y, yc, "right")~treat, taus=(1:4)/5, , method="Powell")
Error in solve.default(x[h, ]) :
Lapack routine dgesv: system is exactly singular: U[2,2] = 0
Error in crq.fit.pow(X, y, cen, tau = taus[i], weights, left = left, ...) :
Singular basic solution generated by 'start'
modèle2
fit2<-crq(Curv(y, yc, "right")~treat+age, taus=(1:4)/5)
Error in solve.default(x[h, ]) :
Lapack routine dgesv: system is exactly singular: U[2,2] = 0
Error in crq.fit.pow(X, y, cen, tau = taus[i], weights, left = left, ...) :
Singular basic solution generated by 'start'
Quelqu'un ici peut-il savoir ce qui ne va pas dans les modèles? Est-ce à cause des liens en y? Y a-t-il une solution?
la source