similar to: glm package

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

Sys.setlocale(,"C") x1 <- as.character(date()) # I use date to record the time, and save it to sqlite database, to it converted to character x1_2 <- strptime(x1, "%a %b %d %H:%M:%S %Y") x2 <- as.character(date()) x2_2 <- strptime(x2, "%a %b %d %H:%M:%S %Y") X<-c(x1_2,x2_2) order(X) ## I want to get the permutation other than the sorted vector. ##

A colleague from China is visiting our lab is working with us on some R-based remote sensing projects, and he was curious if there is a r-help listserv in Chinese and, if not, if it would be possible to set one up? I didn't see one listed on the Mailing List page on the cran-r website. Thanks! --j

i have a loop which looks likes this: > data.info <- rbind(data.info, cbind(station.id, year, date, > max.discharge)) + y <- split(data.info, data.info[station.id]) + for (i in names(y)) {write.csv(y[[i]], file=paste(i, ".csv", sep=","))} i am wondering, where the file (which i am about to write in .csv format) will be saved? i looked at ?write.csv and

Hello helpers, This is probably quite simple, but I'm stuck. I want to create a summary statistics table with frequencies and summary statistics for a large number of variables. The problem here is that (1) there are two different classes of categories (sex, type of substance abuse and type of treatent) which overlap, (2) the data for different variables should be presented in different

Dear R-helpers, I have a data frame with 3 variables, each record is a unique combination of the three variables. I would like to count the number of unique values of v3 in each v1, and save it as a new variable v4 in the same data frame. e.g. df1 [v1] [v2] [v3] [1,] "a" "C" "1" [2,] "b" "C" "2" [3,] "c" "B"

i use the GUI version of R on OS X. I launch it by double-clicking on the R icon. The process always starts in my home directory which means that the only .RData file that is ever read in is the one in my home directory. I would instead like to have different R sessions saved in different directories, but I can't figure out how to do this. A workaround is to use the shell version of R and

I am trying to update the packages that I have installed but I get the following warning messages: package 'tseries' successfully unpacked and MD5 sums checked Warning: cannot remove prior installation of package 'tseries' bundle 'forecasting' successfully unpacked and MD5 sums checked Warning: cannot remove prior installation of package 'forecast' What does that

Hi, I have a simple linear regression for which I want to obtain HAC corrected standard errors, since I have significant serial/auto correlation in my residuals, and also potential heteroskedasticity. Would anyone be able to direct me to the function that implements this in R? It's a basic question and I'm sure I'm missing something obvious here. I looked up this post:

This is my first excursion into using reshape2 and I want to ensure that the melt() function call is syntactically correct. The unmodifed data frame is organized this way: head(tds.anal) site sampdate param quant 1 UDS-O 2006-12-06 TDS 10800 4 STC-FS 1996-06-14 Cond 280 7 UDS-O 2007-10-04 Mg 1620 9 UDS-O 2007-10-04 SO4 7580 19 JCM-10B 2007-06-21 Ca 79 20

Dear all, I am trying to install RSQLite package since I want to install "sqldf", and I used >> install.packages("RSQLite") first, which gave Error message as below: make: *** [RS-DBI.o] Error 1 chmod: cannot access `/usr/lib/R/library/RSQLite/libs/*': No such file or directory ERROR: compilation failed for package 'RSQLite' ** Removing

Dear All, I use R to conduct multilevel modeling. However, I have a problem about the interpretation of random effect. Unlike the variables in fixed effects, the variables in random effects have not shown the p-value, so I don't know whether they are significant or not? I want to obtain this figure to make the decision. Thanks a lot! Below is the syntax and output of my program:

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>

The incorrect result incurs when the file path contains Chinese character. It seems that dirname/basename action on unit of byte instead of char, so the result in the following example is half of what is expected. > g<-"d:\\$BG!2L4^M-CfJ8(B\\$BG!2L4^M-CfJ8(B.txt" > dirname(g) [1] "d:/$BG!2L4^(B" > basename(g) [1] "$BG!2L4^M-(B" -- HUANG Ronggui,

. I just copied all of my Paint Shop Pro 7 files from another partition into a "Paint Shop Pro 7" subdirectory within the "Program Files" directory on a new wine installation. When I run psp.exe, I got a very long list of complaints/errors similar to this one: err:module:import_dll Library MFC42.DLL (which is needed by L"C:\\Program Files\\Paint Shop Pro

Hi, I am trying to run a logistic regression to look at the risk of malaria infection in individuals. I want to account for intra household correlation and so want to include a household level random effect. I have been using the lmer command in lme4 package but am getting some strange results that are completely different to those I get using STATA. Can I just check that this is the correct

Working with 5 subset streams from my source data frame, three of them successfully call dcast(), but two fail: jerritt.cast <- dcast(jerritt.melt, site + sampdate ~ param) Aggregation function missing: defaulting to length and winters.cast <- dcast(winters.melt, site + sampdate ~ param) Aggregation function missing: defaulting to length Yet both data frames have the values in their

I want to write a function mimic the function of select.list(), here is my preliminary version. select <- function(x,multiple=TRUE,...){ ans<-new.env() g <- gwindow(title=title,wid=200,heigh=500) x1<-ggroup(FALSE,con=g) x2<-gtable(x,multiple=multiple,con=x1,expand=TRUE) gbutton("OK",con=x1,handler=function(h,...){ value <- svalue(x2) if (length(value)==0)