search for: wzep

Displaying 4 results from an estimated 4 matches for "wzep".

Did you mean: wep
2011 Aug 17
0
vglm regression with weibull distribution
...y to apply solution provide in help of the package that consisted to use lschape="logoff" and eschape=list(offset=-2) regWeib0<-vglm(eqn0,family=weibull(lshape='logoff',eshape=list(offset=-2)),subset(br, fdh<1)) Erreur dans if ((temp <- sum(wz[, 1:M, drop = FALSE] < wzepsilon))) warning(paste(temp,? : ? l'argument n'est pas interpr?table comme une valeur logique De plus : Message d'avis : In log(theta + offset) : production de NaN Finaly i use the above function to find initial value for shape (a) and scale (b)? weibullPar<-function(x){ ?f<-fu...
2007 Jan 06
2
Using VGAM's vglm function for ordinal logistic regression
...ntercept):1 (Intercept):2 log(exposure.time) -1.6306499 2.5630170 -0.1937956 Degrees of Freedom: 48 Total; 45 Residual Residual Deviance: 503.7251 Log-likelihood: -251.8626 Warning messages: 1: 4 elements replaced by 1.819e-12 in: checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) 2: 4 elements replaced by 1.819e-12 in: checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) 3: 4 elements replaced by 1.819e-12 in: checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) 4: 4 elements replaced by 1.819e-12 in: checkwz(wz, M = M, trace = trace...
2007 Jul 16
2
Error while fitting Partial Proportional Odds model using vglm
...number of items to replace is not a multiple of replacement length In addition: Warning messages: 1: NaNs produced in: log(x) 2: fitted values close to 0 or 1 in: tfun(mu = mu, y = y, w = w, res = FALSE, eta = eta, extra) 3: 19 elements replaced by 1.819e-12 in: checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) Ironically, it works ok with Proportional Odds (change parallel = T above), but fails with PPO (parallel = F). I will appreciate any help to fix this problem. Thanks Rizwan Younis Grad Student University of Waterloo Waterloo, ON Canada
2008 Oct 10
0
Problems and bugs in vgam()
...I frequently run into the error below. Any advice on how to overcome it? Applying Greenstadt modification to 357 matrices Error in ans[, index] <- tmp777 : number of items to replace is not a multiple of replacement length In addition: Warning message: In checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) : 357 elements replaced by 1.819e-12 In case someone has the interest or time in helping me get models to converge, I provide an example model and data. (Data are attached and the model code is below.) Thanks for your help. library(VGAM) Cnictdat=read.csv("Cnictd...