similar to: glm with Newton Raphson

I contacted John Fox about this first, because parts of the file are attributed to him. He says that he didn't write rstandard.glm(), and suggests asking r-devel. As it stands, rstandard.glm() has summary(model)$dispersion outside the sqrt(), while in rstandard.lm(), the sd is already sqrt()ed. This seems to follow stdres() in VR/MASS/R/stdres.R. Of course for the c("poisson",

Hi, Is there available package on the optimization function using Newton-Raphson method (iterative quadratic approximation)? I have been using the 'optim' function in R and found it really unstable (it depends heavily on the initial values and functional forms). If I have to code it by myself, can I get some advice on how to start (any good reference or sample code)? I really

Hello, Does anyone know whether the gls function in the nlme library uses the Newton-Raphson or EM algorithm to find the restricted log-likelihood or maximum log-likelihood estimates? Brendan Klick bklick@jhsph.edu [[alternative HTML version deleted]]

Dear all, I want to solve a score function by using Newton-Raphson algorithm. Is there such a fucntion in R? I know there's one called optim, but it seems only doing minimizing or maximizing. Thanks, Jimmy

I have not included the previous postings because they came out very strangely on my mail reader. However, the question concerned the choice of minimizer for the zeroinfl() function, which apparently allows any of the current 6 methods of optim() for this purpose. The original poster wanted to use Newton-Raphson. Newton-Raphson (or just Newton for simplicity) is commonly thought to be the

Hi, Is there a routine available in R for the Newton-Raphson method for simulataneous equations in several unknowns? Thanks Robert -- Robert J. Chandran Department of Botany 3506 Miller Plant Sciences Building University of Georgia Athens, GA 30602 Phone: (706)-583-0943 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear all.. I have a respon variable y. Predictor variable are x1, x2, x3, x4, x5 (1) What is the syntax to get paramater estimation of ZINB Model by Newton Raphson (not BFGS) (2) What syntax to plot probability of observed & predicted of ZINB Thx. Regards Krist. [[alternative HTML version deleted]]

Hi I would like to use the newton raphson method to find the root if the equation x^3-0.165*x+0.0003993 without using any readliy available program in r but instead by writing my own code and loop. the problem is that i really cant understand how to write the loop so that it keeps using the last calcualted values. if anyone could help me or give me some tips i would deeply appriciate it thanks

Hola soy novata en el programa R, pero lo encuentro súper interesante, tengo un par de consultas... 1. necesito crear una nueva base de datos. 2. necesito saber como se codifica el sample, subset y el rbind. Por favor agradecería sus respuesta Saludos Cordiales Valeska Yaitul Yaitul. [[alternative HTML version deleted]]

I have the following function for which I need to find the root of a: f <- function(R,a,c,q) sum((1 - (1-R)^a)^(1/a)) - c * q To give context for the problem, this is a psychometric issue where R is a vector denoting the percentage of students scoring correct on test item i in class j, c is the proportion correct on the test by student k, and q is the number of items on the test in total. I

Dear Mr.Zeileis & all. (1)     Thx for your reply. Yes, I am talk about the function zeroinfl() from the package "pscl". I want to use Newton Raphson to get parameter             estimation ZINB, so I try this: ----------------------------------------------------------------------------------------------------------------------------------         > zinb <- zeroinfl(y

can anyone tell me how can i control the method of estimation (i.e. scoring method or Newton raphson method) in glm and compute deviance function ? -- View this message in context: http://r.789695.n4.nabble.com/changing-method-of-estimation-in-GLM-tp3170836p3170836.html Sent from the R help mailing list archive at Nabble.com.

Hi, I am thinking about IWLS vs ML estimation. When I use glm() for a 2-parameter distribution (e.g., Weibull), I can otain the MLE of scale parameter given shape parameter through IWLS. Because this scale parameter usually converges to the MLE. In this point, I am wondering: i) can you say that the direct MLE, which is obtained by maximizing a likelihood function, is equalvant to the indirect

I have a question about running an optimization function on an existing LOESS function defined in R. I have a very large dataset (1 million observations) and have run a LOESS regression. Now, I want to run a Newton-Raphson optimization to determine the point at which the slope change is the greatest. I am relatively new to R and have tried several permutations of the maxNR and nlm functions with

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

What kind of special magic does glm have? I'm working on a logistic regression solver (L-BFGS) in c and I've been using glm to check my results. I came across a data set that has a very high condition number (the data matrix transpose the data matrix) that when running my solver does not converge, but the same data set with glm was converging ( I love R :) ). I noticed that glm using

Dear R-Help List, I'm trying to implement Firth's (1993) bias correction for log-linear models. Firth (1993) states that such a correction can be implemented by supplementing the data with a function of h_i, the diagonals from the hat matrix, but doesn't provide further details. I can see that for a saturated log-linear model, h_i=1 for all i, hence one just adds 1/2 to each count,

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070312/5b30b712/attachment.pl

Full_Name: Carmen Fernandez Version: 1.3.1 OS: Submission from: (NULL) (138.251.202.115) I think there is a problem when computing Pearson residuals, in that they seem to be computed at the raw residuals divided by the square root of the corresponding diagonal element of the weight matrix W evaluated at the last step of the iterative model fitting procedure (IWLS), instead of dividing by the

Hi, is there a possibility to calculate the different "types" of residuals directly using the multinom function from MASS as it is possible for the functions gam, glm using type="deviance" or "working" or "pearson" or "response"? I tried it but got always the "response" type, I guess. thanx Matthias