similar to: Dynamic random effects model

Hi all, can the glm package be used to estimate a logit model to panel data? I am asking this because stata has a standard logit model and then an xtlogit for longitudinal data. Is there something similar in R? Thank you jcm [[alternative HTML version deleted]]

Dear I am trying to find a 90% credible interval. I am using the following code. fit<-bugs( model.file=BUGScode, data=data, inits = list(geninits1,geninits2), parameters.to.save=keepers, n.chains=nchains, n.iter=runs, n.burnin=burn, n.thin=nthin, DIC= TRUE, bugs.directory="C:/Program Files/WINBUGS.14", \ ) But this is only giving 95%

the version 2 parser thinks I have unnecessary braces, but I can't find any. False positive or am I missing something? If a false positive, is there any way to work around the warning? * checking Rd files against version 2 parser ... WARNING Warning: ./man/dbetabinom.Rd:32-34: Unnecessary braces at ?{p(x) = % (C(N,x)*Beta(N-x+theta*(1-p),x+theta*p))/% Beta(theta*(1-p),theta*p)}?

with the following simple data frame dd = structure(list(z = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L ), .Label = c("a", "b"), class = "factor"), x = c(0.3, 0.2, 0.1, 0, 0, 0, 0.2, 0.3)), .Names = c("z", "x"), row.names = c(NA, -8L), class = "data.frame") I would like know if it's possible to use model.matrix() to construct the

I have an optimization question that I was hoping to get some suggestions on how best to go about sovling it. I would think there is probably a package that addresses this problem. This is an ordering optimzation problem. Best to describe it with a simple example. Say I have 100 "bins" each with a ball in it numbered from 1 to 100. Each bin can only hold one ball. This optimization is

I am having difficulty with evaluation/environment construction for a formula to be evaluated by model.matrix(). Basically, I want to construct a model matrix that first looks in "newdata" for the values of the model parameters, then in "object at data". Here's what I've tried: 1. model.matrix(~f,data=c(newdata,object at data)) -- fails because something (terms()?)

Is there a reason that Reduce() doesn't take a "..." argument that would allow arbitrary extra arguments to be passed through to the function? Here is a little example of how this would be convenient: z <- list( data.frame(state=c("California"), cases=0), data.frame(state=c("California","Massachusetts"),

From: Vitalie S. <vitosmail <at> rambler.ru> Subject: Bug in cut.POSIXt Newsgroups: gmane.comp.lang.r.general Date: 2009-10-15 15:47:48 GMT (1 hour and 29 minutes ago) Hello Everyone, Before reporting decided to post here first: tt <- structure(c(1254238817, 1254238859, 1254238969, 1254239080), class = c("POSIXt",

Consider the following situation: we have quantified algal concentrations for a variety of species using many samples at each of three years. It seems to make sense to generate a line plot (matplot-like), with each species plotted as a separate line, with the points connected to emphasize the temporal pattern. The problem: lots of overlapping error bars. The question: from both a

Hi, I'm trying to install rjags and R2jags with R version 2.14.1 (2011-12-22) on Linux Mint 12 machine using software repository: deb http://cran.fhcrc.org/bin/linux/ubuntu oneiric/ jags seems to install OK, but rjags gives a weird error about not finding fprec symbol. Here is a history of the problem including several attempts to work around it. (Because of Ben Bolker's recent

I have found the "simulate" method (incorporated in some packages) very handy. As far as I can tell the only class for which simulate is actually implemented in base R is lm ... this is actually a little dangerous for a naive user who might be tempted to try simulate(X) where X is a glm fit instead, because it defaults to simulate.lm (since glm inherits from the lm class), and the

Hi, I havw some porblem in R and i need your help. I have this vector and I want to change it to date. for example G=[05 12 2008] or g=[2] f=[3] y=[2208]

There is a thread on this topic already: http://finzi.psych.upenn.edu/Rhelp10/2010-August/250934.html I'm rather mystified by a similar problem and wondering whether I've overlooked something obvious. I'm running with latest versions of R and all packages, and latest version of JAGS running under Windows 7. Here's the problem. I have some source code. It's given below -

Hi all, I ran a Weibull model, and I am wondering if there is any way to extract the log likelihood. I tried loglik(model) but it does not seem to work any help would be greatly appreciated joe [[alternative HTML version deleted]]

Dear R experts, I am trying to estimate an ARMA 2,2 model with garch errors. I used the following code on R 2.9. #library library(fGarch) #data data1<-ts(read.table("C:/Users/falcon/Desktop/Time Series/exports/goods1.csv"), start=c(1992,1), frequency=12) head(data1) #garch garchFit(formula.mean= ~arma(2,2),formula.var=~garch(1,1), data=data1) but get this error: >

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091125/5ce76d9b/attachment-0001.pl>

Hello, I asked a question about what the most likely process to follow if after a time-series fit is performed the residuals are found to be non-normal. One peron responded and offered to help if I supplied a sample data set. Unfortunately now that I have a sample I have lost the emai addressl. If you are that person or have some ideas please email me back at rkevinburton at charter.net. Thank

Anyone used jags.parallel in latest version of R2jags? Ubuntu Oneiric, just ran updates in R to get latest R2jags package which supposedly has added jags.parallel command to send chains to multiple cores ... but when I submit, I get: test <- jags.parallel(d, inits, param, model.file="rats.bug", n.chains=3, n.iter=2000) Error: could not find function

Hello ALL! I am running Linux, Fedora 15 64-bits, and R on it. I need to use WinBUGS and R2WinBUGS, but as far as I read, WinBUGS is closed project, to be continued with/as OpenBUGS. Thus, I have found R2OpenBUGS on OpenBUGS Contributed Code (http://openbugs.info/w/UserContributedCode), not on CRAN. Author(s) states that it is equivalent for R2WinBUGS. I tried briefly, and realized few minor

Frank (or anyone else), can you offer any comments comparing runjags, R2jags, rjags ? I couldn't find any vignettes, nothing except a brief mention on Task Views. Kevin On Wed, Oct 7, 2009 at 7:48 AM, Frank E Harrell Jr <f.harrell@vanderbilt.edu > wrote: > Have you tried the rjags package which uses the jags system? It is much > more integrated into R and works quite well.