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

2012 Oct 31

0

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

0

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<- read.gal("E:/splm/zj.GAL",override.id=TRUE)
ww<-nb2listw(w,zero.policy=TRUE)
#????
library(foreign)

2008 Oct 31

2

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

1

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

1

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

1

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

1

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

2

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

1

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

1

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

3

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

1

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

1

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

1

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

1

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

1

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

2

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

3

An 'orrible question: Outlook 365 under wine on CentOS?

On 09/26/2016 01:28 PM, m.roth at 5-cent.us 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

1

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

1

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