search for: magagnoli

All R experts, How do I fit a dynamic Random effects model with a binary dependent variable in R Thanks JCM [[alternative HTML version deleted]]

Morning All: Having problems installing this one into my R setup and wondered if anyone on the list had the same problem and found a workable solution. First to the system particulars: OS is the Hardy Heron edition of Ubuntu Linux R version is 2.13.1 Java installation is sun-java6 Next up is the result of javareconf: R CMD javareconf Java interpreter : /usr/bin/java Java version : 1.5.0

Hi all, I was wondering how to construct a seasonal differenced time series variable. I used the following code to construct a 12 span seasonal difference seasonal<-diff(V2, lag=12, differences=1) is this correct? thank you in advance 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: >

Dear All, I am practicing data manipulation and I would like to generarte a count variable. My data looks like this: Country MID 1 NA 1 0 1 0 1 1 1 0 2 0 2 1 2 0 2 0 2 0 I would like to to generate a variable

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

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

Hi all, I am getting an error in my code and I don't know what the problem is. I am using R 2.9 on ubuntu. my code is as follows: ## Libraries ## library(survival) library(foreign) ## reading data ## data<-read.dta("http://psfaculty.ucdavis.edu/bsjjones/cabinet.dta") head(data) attach(data) fit1<-survfit(Surv(durat,censor)) and I get the following 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>

...is really compatible with hpux or not . If compatible , then which version of hpux reqd ? what are other dependency to load "R" ? What will be the patch level and FORTRAN version ? Also is there any other s/w equivalent to "R" which can be tested on hpux ? " -- Joseph C. Magagnoli [[alternative HTML version deleted]]

Can someone please recommend to me a book on the programming language that R is based on? I'm looking for a foundational book that teaches the logic of the S language. It seems that knowing the underpinnings of the language can only make using R a bit easier. Any leads are greatly appreciated . . . Matt. [[alternative HTML version deleted]]

Hi everyone, I have the following: sqldf("select Premie,count(tpounds) N,avg(tpounds) Avg_Weight, stddev_samp(tpounds) StdDev from children group by Premie having !is.na(Premie)") sqldf() does not like the !is.na(Premie) specification. How does one exclude a "missing" group in an aggregated query using sqldf()? Thanks! Dan [[alternative HTML version deleted]]

Hello people, I am in the process of migrating from Stata to R and I would like to check if my results are similar under the two softwares: Here is my GLM command under R nurse.model<-glm(pQSfteHT~dQSvacrateHTQuali3_2 + dQSvacrateHTQuali3_3 + dQSvacrateHTQuali3_4 + dQSvacrateHTQuali3_5 + cluster_32 + cluster_33 + cluster_34 ,family=binomial(link = "logit")) and below the stata

...roduct of every two elements in two vectors (Henrique Dallazuanna) 18. Re: Book on R's Programming Language (Joshua Wiley) 19. Re: how to calculate the product of every two elements in two vectors (Gabor Grothendieck) 20. Re: Book on R's Programming Language (Joseph Magagnoli) 21. Re: Constrain density to 0 at 0? (Greg Snow) 22. matest function for multiple factors (shabnam k) 23. Re: How to deal with more than 6GB dataset using R? (Greg Snow) 24. Using R to fill ETM+ data gaps? (Abdi, Abdulhakim) 25. How to generate a sequence of dates without...

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

I estimated a negative binomial model using zelig. z.out<- zelig(NEWBHC~ PW80 + CHNGBLK + XBLK,data=data, model="negbin") How do I calculate predicted probabilities for this model? Is it the same process as a poisson regression? Thanks in advance Joe [[alternative HTML version deleted]]