similar to: Tackling of convergence issues in gamlss vs glm2

Hello all: Please help me with this "can't find object" issue. I'm trying to get leave-one-out predicted values for Beta-binomial regression. It may be the gamlss issue because the code seems to work when %do% is used. I have searched for similar issues, but haven't managed to figure it out. This is on Windows 10 platform. Thanks in advance, Nik #

If the code you are running in parallel is complicated, maybe foreach is not sophisticated enough to find all the variables you refer to. Maybe use parallel::clusterExport yourself? But be a aware that passing parameters is much safer than directly accessing globals in parallel processing, so this might just be your warning to not do that anyway. -- Sent from my phone. Please excuse my brevity.

Hello: Could you please suggest the best way to become an "advanced" R programmer. I went through "R for dummies" by de Vries and Meys and I can see two ways to proceed: 1) Get a more advanced textbook. E.g. could you recommend Gentleman, "R for Bioinformatics"? 2) Because textbooks are limited and become obsolete fast, I can focus on learning state-of-the-art

Hello Mikis: Thanks a lot, it worked. Could you tell me what the problem was? Regards, Nik ----- Original Message ----- From: r-help-request at r-project.org To: "r-help" Sent: Sunday, March 11, 2018 6:00:02 AM Subject: R-help Digest, Vol 181, Issue 11 Send R-help mailing list submissions to r-help at r-project.org To subscribe or unsubscribe via the World Wide Web, visit

Hi, I'm a bit puzzled by the gamlss fitting of exponential and lognormal data. Gamlss seems to think that exponentially distributed data fits better with a lognormal distribution, and vice versa. For example, X <- rexp(1000) Gexp <- gamlss(X~1,family=EXP) # X~1 is X tilde 1 GAMLSS-RS iteration 1: Global Deviance = 2037.825 GAMLSS-RS iteration 2: Global Deviance = 2037.825 Glno

See Hadley's advanced R along Thomas Mailund's books. I haven't gone through them carefully but they both seem (from what I've looked at ) to be the best ones for that. Mentions of others are appreciated. On Tue, Mar 13, 2018 at 5:26 PM, Nik Tuzov <ntuzov at beacon.partek.com> wrote: > > Hello: > > Could you please suggest the best way to become an

Dear all, I am using gamlss (Package gamlss version 4.0-0, R version 2.10.1, Windows XP Service Pack 3 on a HP EliteBook) to relate bird counts to habit variables. However, most models fail because “the global deviance is increasing” and I am not sure what causes this behaviour. The dataset consists of counts of birds (duck) and 5 habit variables measured in the field (n= 182). The dependent

Hola a tod en s, tengo una duda que la comunidad R me puede ayudar. Estoy trabajando con gamlss, porque tengo una variable respuesta con valores entre 0 y 1 e incluidos estos. La distribución que utilizo com gamlss para este caso es "beta inflated" (Stasinopulos and Rigby 2007. Journal of Statistical Software 23(7)). El modelo que intento correr es: m1<-gamlss(Teleosteos ~

Dear R helpers, I have following loss data and I need to fit LEFT truncated Log Normal distribution to this data which is Truncated at 1000000. dat = c(1333834,5710254,9987567,7809469,6940935,3473671,1270209,1102523,1124002, 5830159,4302300,3925242,2638409,2324421,7238436,9088709,7439250,4976551,4864319, 8741334,1863770,7098310,4942288,4971829,4986372) library(gamlss.tr) gen.trun(5, LOGNO)

Dear all and Mikis I have the opportunity to compare fits with the 'classical' glm and gamlss and no smoother of any kind just the same model formula (both with the binomial family). I get exactly the same coefficients but very different residuals, gamlss giving residuals which are extremely close to 'normal' and glm very far... How can this be ? Thanks in advance for

Dear Nik Try the following code loo_predict.mu <- function(model.obj, input.data) { yhat <- foreach(i = 1 : nrow(input.data), .packages="gamlss", .combine = rbind) %dopar% { updated.model.obj <- update(model.obj, data = input.data[-i, ]) predict(updated.model.obj, what = "mu", data = input.data[-i, ], newdata = input.data[i,], type =

Dear All: I have succeeded in fitting a GAMLSS.dist model to growth data I am working with it. My aim is to create a matrix of predicted percentiles and the corresponding the fitted model's sigma mu nu by agebins. Q: How do it generate these parameters as in L M S per Cole and Green 1992? Here are my working codes. Name of fitted model is gamlssfit > Agebin<-seq(6,36,6)

Hi, I'm playing around with gamlss and don't entirely understand the sigma result from an attempted lognormal fit. In the example below, I've created lognormal data with mu=10 and sigma=2. When I try a gamlss fit, I get an estimated mu=9.947 and sigma=0.69 The mu estimate seems in the ballpark, but sigma is very low. I get similar results on repeated trials and with Normal and

Hi Folks! Just started using the gamlss package and I tried a simple code example (see below). Why the negative sigma? John > y <- rt(100, df=1)> m1<-fitDist(y, type="realline")Warning messages:1: In MLE(ll3, start = list(eta.mu = eta.mu, eta.sigma = eta.sigma, : possible convergence problem: optim gave code=1 false convergence (8)2: In MLE(ll4, start = list(eta.mu =

Dear all I am trying to build a package in Windows. I use the following command (which it used to work with previous versions ) and I am getting the following error #-------------------------------------------------------------------------------------------------------------- C:\PROGRA~1\R\rw2000\bin>Rcmd build --binary --use-zip C:\PROGRA~1\R\rw2000\src\library\gamlss * checking for

I am attempting to run a glm with a binomial model to analyze proportion data. I have been following Crawley's book closely and am wondering if there is an accepted standard for how much is too much overdispersion? (e.g. change in AIC has an accepted standard of 2). In the example, he fits several models, binomial and quasibinomial and then accepts the quasibinomial. The output for residual

Hi everyone, I''m trying to use function vcov to extract the covariance matrix from a GAMLSS object. But I''m getting some strange errors and I was hoping someone could help me out? Vcov works with the same model for lm and glm objects, but not gamlss objects. I''ve searched various help sites to no avail. Its very possible the reason is that vcov failed though,

Hi, I am running a negative binomial model using Gamlss and when I try to include random effect, I get the following message: Warning messages: 1: In vcov.gamlss(object, "all") : addive terms exists in the mu formula standard errors for the linear terms maybe are not appropriate 2: In vcov.gamlss(object, "all") : addive terms exists in the sigma formula standard

Dear all I am trying to create a package in R 1.9.0 and I a getting an error message which I do not understand. (I am using R in Windows XP and 2000) For example the following works well in 1.8.1 C:\Program Files\R\rw1081\src\gnuwin32>make pkg-gamlss ---------- Making package gamlss ------------ adding build stamp to DESCRIPTION installing inst files installing indices not zipping

This is a question regarding the 'regTermTest' function in the 'survey' package. Imagine Z as a three level factor variable, and code ZB and ZC as the two corresponding dummy variables. X is a continuous variable. In a 'glm' of Y on Z and X, say, how do the two test specifications test.terms = c("ZB:X","ZC:X") # and test.terms = ~ ZB:X + ZC:X in