similar to: changing method of estimation in GLM

i am trying to fit a linear model with both continuous covariates and factors. When fitted with the intercept term the first level of the factor is treated by R as intercept and the estimate of the effects of remaining levels(say i th level) are given as true estimate of i th level - estimate of 1st level.can any please help me? thanks in advance..... -- View this message in context:

i am new to the quantmod package . so if the answer is trivial please excuse me. i want to study stock values within a day. i get current stock updates using getQuotes and then want to produce usual quantmod graphs with that values. also the graph should be able of adding technical indicators. please help. in addition it will be helpful if anyone suggests how to run that code continuously to get

i am trying to solve a partial differential equation analytically(PDE) in R . i have found some functions that do the stuff numerically. But that will not meet my purpose. is there any function to solve PDE analytically. please help. -- View this message in context: http://r.789695.n4.nabble.com/analytical-solution-of-partial-differential-equation-tp4291618p4291618.html Sent from the R help

i am using quantmod package.it get stock quotes from google finanace. but unfortunately i am not able to get the quotations of some stocks(e.g. NSE:TCS,NSE:SAIL ) through the "getSymbol" command of this package although they are available in the google finance website. anyone please help me. thanks in advance..... -- View this message in context:

dear community, I am a beginer in R , and can't predict with logistic model in package logistf, could anyone help me ? thanks ! the following is my command and result : >library(logistf) >data(sex2) >fit<-logistf(case ~ age+oc+vic+vicl+vis+dia, data=sex2) >predict(fit,newdata=sex2) Error in predict(fit, newdata = sex2) : no applicable method for "predict"

Hi, I am using Rstudio in Ubuntu . While installing R from the terminal I got the version R 2.13.1 . But the problem is that knitr 0.5 requires higher versions. So is there any way to get the latest version of R in ubuntu (or 2.13.1 is the maximum I can get)? In case the higher R version is not available can anybody please help me how do I use the "knitr Html" option in R studio.

Hi, Does anyone know if there is a function to find the maximum likelihood estimates of glm using Newton Raphson metodology instead of using IWLS. Thanks Valeska Andreozzi -------------------------------------------------------- Department of Epidemiology and Quantitative Methods FIOCRUZ - National School of Public Health Tel: (55) 21 2598 2872 Rio de Janeiro - Brazil

Sonu: I'm giving the whole history to the rsync list. If the tests i subscribed passed, yet you get that error, it's more in depth than i can go into. You've got 1 hour, let's hope somebody immediately recognizes the problem. I've got some ideas, but no time to explore any more for you. Gentlemen: does anybody see a quick resolution for Sonu? Looks like he's

I am using R 2.14.1. I am trying to configure Lyx with Sweave. I have read articles but I found none to be complete and I cannot import sweave document in lyx . Can anyone please help me with a stepwise procedure how to configure Lyx with Sweave so that I can run my R-code chunks from LYX ,the output being a pdf.? Thanks in advance. -Atanu -- View this message in context:

Dear all, I am encountering some issues with my data and need some help. I am trying to run glm analysis with a presence/absence variable as response variable and several explanatory variable (time, location, presence/absence data, abundance data). First I tried to use the glm() function, however I was having 2 warnings concerning glm.fit () : # 1: glm.fit: algorithm did not converge # 2:

Hi everyone, I'm wondering how to select the "best" model when using logistf? AIC does not work neither does anova. I tried fitting a glm model but got the separation warning message so I tried using the logistf package but as I stepwise simplify the model I don't know if the simplification is motivated or not... Can anyone explain to me how I should approach this problem? I

hi all, first i wann to explain my network & it''s problem, i am using a redhat el5 for multigateway routing, i have 3 adsl lines which is connected with eth0, eth2 and eth3. eth1 is my local network. those 3 lines comes with adsl modem from same isp with equal bandwidth (512 Kbps). In this adsl modems there is one lan and one wan port, adsl modem''s lan ip is gateway for my

Hi, I got two questions and would really appreciate any help from here. First, is the penalized maximum likelihood estimation(Firth Type Estimation) only fit for binary response (0,1 or TRUE, FALSE)? Can it be applied to multinomial logistic regression? If yes, what's the formula for LL and U(beta_i)? Can someone point me to the right reference? Second, when I used *logistf *on a dataset with

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,

Hi, we try to do a logistic regression with the function glm. But we notice that this function don't give the same results as the SAS proc catmod (differents estimate given). We try to change the contrast on R system with: > options(contrasts=c(unordered="contr.SAS",ordered="contr.poly")) We also try with brlr and logistf functions. Unfortunately, the estimate

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]]

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, Could anyone tell me what is the command used in R to do 1. Differentiation 2. Newton Raphson method (Numerical Analysis in general...) Are there any packages separately for this? Thanks for your help! BR, Shubha [[alternative HTML version deleted]]

Hi , I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likelihood but my problem after runing the code is the estimators were too high. with this email ,I provide