similar to: how to get loglik parameter from splm package?

Displaying 20 results from an estimated 300 matches similar to: "how to get loglik parameter from splm package?"

2012 Oct 31
pseudo R-squared for model generated with spgm (splm)
I am working with the splm package. I use the spgm function: general estimation of a panel data model. Based on this approach, I know it is possible to compute a R2, eg the ratio of variation explained by a given model. My model is : bivmod<-spgm(logIKA~NBLITRE0+NBLITRE1,data=mydatap,listw=comsKnn.nbW,spatial.error=TRUE) I know that we can calculate the R^2 as the variance of the fitted
2012 Sep 27
error while estimating spatial Durbin (mixed) model
Dear all, I am new here ,I attempted to use R to estimate the spatial Durbin (mixed) model,and mydata is a panel data form,and the matrix is generated by geoda software ,here is my Command and error,really hope your help ,thank you! #??gal library(spdep) w<-"E:/splm/zj.GAL", ww<-nb2listw(w,zero.policy=TRUE) #???? library(foreign)
2008 Oct 31
Problem of running R console from Java on linux box
Hi the list: I can not create R process from Java application with my linux box, the process never was created, or just terminated immediately after I called Runtime.getRuntime().exec(...) The command lines I have tried are(the paths are all correct): "/bin/sh /usr/bin/R" "/bin/sh -c /usr/bin/R"
2012 Apr 26
PLM package PGGLS strange behavior
When using the PLM package (version 1.2-8), I encounter the probem that calling the FGLS estimator evokes strange behavior, when choosing the "random" effects model. After calling the PGGLS function to estimate FGLS, PLM gives me a warning, stating that the "random" model has been replaced with the "pooling" model. I would, however, really like to estimate the random
2009 Dec 08
Serial Correlation in panel data regression
Dear R users, I have a question here library(AER) library(plm) library(sandwich) ## take the following data data("Gasoline", package="plm") Gasoline$f.year=as.factor(Gasoline$year) Now I run the following regression rhs <- "-1 + f.year + lincomep+lrpmg+lcarpcap" m1<- lm(as.formula(paste("lgaspcar ~", rhs)), data=Gasoline) ###Now I want to find the
2007 May 24
lme with corAR1 errors - can't find AR coefficient in output
Dear List, I am using the output of a ML estimation on a random effects model with first-order autocorrelation to make a further conditional test. My model is much like this (which reproduces the method on the famous Grunfeld data, for the econometricians out there it is Table 5.2 in Baltagi): library(Ecdat) library(nlme) data(Grunfeld)
2013 Nov 06
resdiuals of random model estimated by plm function
Hi all, I have estimated a random panel model using plm function. I have a question about the vector of resduals obtained with the object $residuals. example: data("Produc", package = "plm") zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, model="random", data = Produc, index = c("state","year")) res<-zz$residuals #
2010 Jul 05
Dovecot git mirror
I've been working on Evolution's 'imapx' back end to support QRESYNC, and was going to look at NOTIFY next. Then I realised I may have to look at implementing it on the server side first. Before I can even contemplate that, though, I had to mirror it into git. I've always resisted the one-version-control-system-per-child craze, and these days if I'm even contemplating a
2003 Nov 05
A real-life production scenario
Since it's all the craze, I might as well post our current Asterisk usage. :-) EQUIPMENT: - Beefyish box (dual Xeon 2.4GHz, gig of RAM, more-than-adequate disk space, etc) in a 1U chassis. - A second, slightly less beefyish box of specs I don't have handy right now, also in a 1U. - 2xTE410P CONNECTIONS: - 1 PRI to telco for local outbound/direct-dial inbound, 300 numbers
2008 Mar 20
logLik calculations
Does the ?logLik? function applied to a ?glm? and ?glm.nb? (from MASS package) calculate the complete log-likelihoods, or does it drop the constant terms of the equation? (It?s not clear from the associated help pages, and I?ve found no reference from searching the R help mailing list) Thank you, Kelly Young
2006 Jul 01
help me understand rest
I am having a heck of a time understanding the new rest craze in rails. what I think I understand so far: 1 rest is about the way we use http to access information on the internet. 2 http was created with nouns and verbs in mind, but the only verbs that are supported in browsers and server software today are ''post'' and ''get''. Other useful verbs include put
2003 Sep 23
loops in Sweave
Dear all, I was wondering whether there is a way to make loops in Sweave, i.e. for example to: 1) calculate a parameter, say, a=length(b) 2) according to that, add #a# chapters to the document, each including some repetitive analysis, each time done on a particular subset of the data indexed by the elements of 1:a. This would be of great help for repeating exploratory data analyses on, say,
2010 Feb 04
plm issues: error for "within" or "random", but not for "pooling"
Dear all I am working on unbalanced panel data and I can readily fit a "pooling" model using plm(), but not a "within" or "random" model. Reproducing the examples in vignette("plm") and in the AER package I encountered no such issues. ##unfortunately I cannot disclose the data, and it is too big anyway > dim(ibes.kld.exp.p[x.subs , ]) [1] 13189 34
2013 Jan 31
LogLik of nls
Hello there, Can anyone point me to the code for logLik of an nls object? I found the code for logLik of an lm but could not find exactly what function is used for calculating the logLik of nls function? I am using the nls to fit the following model to data - Model 1: y ~ Ae^(-mx) + Be^(-nx) +c and want to understand what is the likelihood function used by nls. Presumably it is using -
2009 Jun 12
Function for AIC or logLIK for nlsList object
Dear R users, Does anybody have a function to calculate logLik or AIC for nlsList objects? After receiving error messages, another user helped me ascertain that this function is not currently written into R. Many thanks Lindsay
2009 Jul 15
GLM Gamma Family logLik formula?
Hello all, I was wondering if someone can enlighten me as to the difference between the logLik in R vis-a-vis Stata for a GLM model with the gamma family. Stata calculates the loglikelihood of the model as (in R notation) some equivalent function of -1/scale * sum(Y/mu+log(mu)+(scale-1)*log(Y)+log(scale)+scale*lgamma(1/scale)) where scale (or dispersion) = 1, Y = the response
2006 Jan 10
Obtaining the adjusted r-square given the regression coefficients
Hi people, I want to obtain the adjusted r-square given a set of coefficients (without the intercept), and I don't know if there is a function that does it. Exist???????????????? I know that if you make a linear regression, you enter the dataset and have in "summary" the adjusted r-square. But this is calculated using the coefficients that R obtained,and I want other coefficients
2016 Sep 26
An 'orrible question: Outlook 365 under wine on CentOS?
On 09/26/2016 01:28 PM, m.roth at wrote: > John Jasen wrote: >> The only linux-based client that, if I recall, can speak native MS mail >> protocols, was Evolution. >> >> I don't know if it still does. >> > Yeah... and this is O365. > > Stupid question: if I check out evolution... will it munge my thunderbid > email inbox or
2010 Jul 20
Servreg $loglik
Dear R-experts: I am using survreg() to estimate the parameters of a Weibull density having right-censored observations. Some observations are weighted. To do that I regress the weighed observations against a column of ones. When I enter the data as 37 weighted observations, the parameter estimates are exactly the same as when I enter the data as the corresponding 70 unweighted observations.
2009 Feb 10
loglik and arima()
All - I am evaluating an arima(2,1,3) and arima(3,1,3) and notice the log-likelihood of the restricted model is higher than the log-likelihood of the unrestricted. Since these are nested models, I thought the unrestricted model would have a log-likelihood at least as large as that of the restricted model. Am I interpreting the "loglik" output incorrectly? Regards, Stephen