search for: irl

...the weights used in the IWLS process wresid a working residual, weighted for "inv.var" weights only. How to use these input arguments? Anybody please shed some light? Also, is my understanding below correct? The input argument "w" is used for the initial values of the rlm IRLS weighting and the output value "w" is the converged "w". The "weights" input argument is actually what I want to apply - if I have my own set of weights. And the real/actual weights are the product of "weights"(I supplied) and the converged output "w&q...

Hi, In the help files in the ?mgcv package for the gam.control() function, there is an option irls.reg. The help files describe this option as: For most models this should be 0. The iteratively re-weighted least squares method by which GAMs are fitted can fail to converge in some circumstances. For example, data with many zeroes can cause problems in a model with a log link, because a mean of...

...2X86. Yet I am unable to print even a test page. I get the following error: "Test page failed to print. Would you like to view the print trouble shooter for assistance? Access is denied." And I get the following errors in my samba log file: tdb(/usr/local/samba/var/locks/printing/irl.tdb): tdb_oob len 1600941951 beyond eof at 67552 [2003/11/18 11:05:20, 0] tdb/tdbutil.c:tdb_log(724) tdb(/usr/local/samba/var/locks/printing/irl.tdb): tdb_oob len 858730825 beyond eof at 67552 [2003/11/18 11:05:20, 0] tdb/tdbutil.c:tdb_log(724) tdb(/usr/local/samba/var/locks/printing/irl.tdb)...

..."properties box to pop up, I cannot get the Test Page to print. I get the error: Test page failed to print. Would you like to view the print trouble shooter for assistance? Access is denied. Also, these strange errors show up in /var/log/samba/ tdb(/usr/local/samba/var/locks/printing/irl.tdb): tdb_oob len 976566084 beyond eof at 67552 [2003/11/14 14:05:21, 0] tdb/tdbutil.c:tdb_log(724) tdb(/usr/local/samba/var/locks/printing/irl.tdb): tdb_oob len 1530015816 beyond eof at 67552 [2003/11/14 14:05:21, 0] tdb/tdbutil.c:tdb_log(724) tdb(/usr/local/samba/var/locks/printing/irl.tdb)...

...M VSS) with 500 ntlm_auth subprocesses. Lesser number of ntlm_auth sometimes results in squid crash as it sometimes has several hundreds requests per second to authorize and is intolerant to exhaustion of free ntlm_auth. "squid -k rotatelog" at midnight results in crash: Feb 14 00:03:00 irl savecore: reboot after panic: get_pv_entry: increase vm.pmap.shpgperproc Feb 14 00:03:00 irl savecore: writing core to vmcore.1 Btw, I have coredump. vm.pmap.shpgperproc has default value (200) here, as well as m.v_free_min, vm.v_free_reserved, and vm.v_free_target and KVA_PAGES. These crashes a...

Here in Sweden, it's supposed to be springtime. A wonderful time of the year, with sunny skies and wonderful weather. Almost summer. Today, it's not. It's winter all over again with rain and only 3 degrees celsius outside. Better to stay inside and write a weekly Asterisk newsletter :-) This week's topics: ------------------- * Looking beyond Asterisk 1.0/1.1 - what's up? *

Dear R developers, I have noticed a discrepancy between the coefficients returned by R's glm() for logistic regression and SAS's PROC LOGISTIC. I am using dist = binomial and link = logit for both R and SAS. I believe R uses IRLS whereas SAS uses Fisher's scoring, but the difference is something like 100 SE on the intercept. What accounts for such a huge difference? Thank you for your time. Ben Godlove [[alternative HTML version deleted]]

...981 1371706.000 0 16513.0 1.070802e+03 10 1992 AUS GBR 0 1 17.4950008 321708.281 57.9023476 1071537.375 0 16602.3 2.279130e+03 11 1992 AUS GRC 0 0 17.4950008 321708.281 10.3699999 102022.352 0 14845.6 4.164985e+01 12 1992 AUS IRL 0 1 17.4950008 321708.281 3.5490999 54272.410 0 16895.0 1.076323e+02 13 1992 AUS ISL 0 0 17.4950008 321708.281 0.2611000 6976.168 0 16443.6 2.190602e+01 14 1992 AUS ITA 0 0 17.4950008 321708.281 56.7976494 1265800.125...

Hello, I would like to know if there is any package that allow us to fit Generalized Linear Models via Maximum Likelihood and Linear Models using Generalized Least Squarse in R as the functions glim and gls, respectively, from S-Plus. Also, anybody know if there is any package that fit Log-Linear Models using Generalized Least Squares? Any help will be very useful. Thanks, -- Frederico

...resence has been made for next year : http://www.fosdem.org/index/press/devrooms_call_for_presence_2006 to allow projects and communities to book a room. Ruby and Rails are becoming more and more popular in Europe so I think it could be great to have a developer room to allow lot of Rubyist to see irl and talk/work/hack/whatever. If nobody is already planning something (I didn''t see anything in the archives I can manage the thing receive emails and do the communication with the Fosdem''s guys. For a start a wiki with few pages could do the thing and allow everyone to tell when...

...stoph Lehman had problems with seperated data in two-class logistic regression. One useful little trick is to penalize the logistic regression using a quadratic penalty on the coefficients. I am sure there are functions in the R contributed libraries to do this; otherwise it is easy to achieve via IRLS using ridge regressions. Then even though the data are separated, the penalized log-likelihood has a unique maximum. One intriguing feature is that as the penalty parameter goes to zero, the solution converges to the SVM solution - i.e. the optimal separating hyperplane see http://www-stat.stan...

...l using family=poisson(link = "identity"). The problem is that the glm function fits a model with a negative intercept, which sounds like a nonsense to me, being the response a Poisson variable. >From a previous discussion on this list I've understood that the glm function uses IRLS for the fitting without any constraint so it is possible for it to end up in a region of values which are not admissible given the model, and in fact sometimes it fails asking for valid starting values. In this case I expected it to fail asking to supply starting values, and instead it fits the...

Hi, I am trying to build GAM with linear constraints, for a general link function, not only identity. If I understand it correctly, the function pcls() can solve the problem, if the smoothness penalties are given. What I need is to incorporate the constraints before calculating the penalties. Can this be done in R? Any help would be greately appreciated. -- View this message in context:

Is anyone else concerned with the way in which the R function parser relocates comments that occur after condional expressions in functions to before, i.e. fred <- function(x) { # wonder what x is like if (x>0) stop("Sorry non-positive x only") # that showed x big-time! x } but then fred is parsed and stored as "fred" <- function (x) { # wonder what x is

I have written a program to solve a particular logistic regression problem using IRLS. In one step, I need to integrate something out of the linear predictor. The way I'm doing it now is within a loop and it is as you would expect slow to process, especially inside an iterative algorithm. I'm hoping there is a way this can be vectorized, but I have not found it so far. The...

...ssage: " There were 12 warnings (use warnings() to see them)" So I checked them: Warnings 1 to 11 said: 1: optim returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH in: "LMEoptimize<-"(`*tmp*`, value = structure(list(maxIter = 50, ... and Warning 12 said: 12: IRLS iterations for glmm did not converge in: lmer(Enfermo ~ Bloque + Trat * Var * dia + (1 | IDfruto), data = desinf, ... There was no error message in the call of lmer function, so: What do these warnings mean? Do these errors occur at earlier iterations but finally the model converges? Can t...

Dear All, I have the following error message when I fitted lmer to a binary data with the "AGQ" option: Error in family$mu.eta(eta) : NAs are not allowed in subscripted assignments In addition: Warning message: IRLS iterations for PQL did not converge Any help? Thanks in advance, Abderrahim [[alternative HTML version deleted]]

Hi all! My question is why, and what I can do about that I sometimes, but not always, get warning-messages like nlminb returned message singular convergence (7) in: LMEopt(x = mer, value = cv) or IRLS iterations for PQL did not converge when trying to fit a model looking like this: lmer<-(cbind(Diseased,Healthy)~Fungus+(1|Family)+(1|Fungus:Family), family="binomial") to four similar data sets? All four data sets consists of four columns; Fungus (1 or 2), Family (1-30), Disease...

...e send all offers for help *off list* to us at info@astricon.net. Do not disturb the list with offers of your services, please. I repeat: Only the Developer's Meeting will be considered for broadcast at this time. In order to enjoy the conference, you will simply have to be there. It's an IRL experience - meeting all the other Asterisk user's from around the globe, listening to speakers, chatting, making business deals, meeting the exhibitors, developing code and submitting it to the bug tracker on the lunch break... If you have any questions around this - please take them off list...

Dear R users, I use the glm() function to fit a generalized linear model with gamma distribution function and log link. I have read in the help page that the default method used by R is "glm.fit" (iteratively reweighted least squares, IWLS). Is it possible to use maximum likelihood method? Thanks Silvia Narduzzi Dipartimento di Epidemiologia ASL RM E Via di S. Costanza, 53 00198