similar to: Code for bivariate Poisson regression?

Dear R Help, Does anyone know off-hand the syntax for a bivariate (two-outcome) lme? It's straightforward to fit a standard random-intercept random-slope model to each of the two outcome variables: lme(Y~X,random=list(ID=~X)...) lme(Z~X,random=list(ID~X)...) But it would be nice to fit the data (Y,Z) against X so that I can obtain estimates of the correlation between the slopes and the

Hi Folks, I compiled and installed Assist-0.1.0; I'm using R-2.0.0 running on SuSe Linux, on a dual AMD64 machine. However, I (apparently randomly) get segmentation faults when I use cubic2. Any ideas? Thanks Simon -- Simon D.W. Frost, MA DPhil Adjunct Assistant Professor Department of Pathology University of California, San Diego UCSD Antiviral Research Center 150 W. Washington St. San

Dear R-Help list, I have a nonlinear least squares problem, which involves a changepoint; at the beginning, the outcome y is constant, and after a delay, t0, y follows a biexponential decay. I log-transform the data, to stabilize the error variance. At time t < t0, my model is log(y_i)=log(exp(a0)+exp(b0)) at time t >= t0, the model is log(y_i)=log(exp(a0-a1*(t_i - t0))+exp(b0=b1*(t_i -

Dear R-Help, I'm trying to learn the sspir package for state space modeling. Has anyone coded a cubic spline smoother (continuous time) in state space format in sspir? The syntax for setting up the various matrices would be really helpful. Best Simon -- Simon D.W. Frost, D.Phil. Assistant Adjunct Professor of Pathology University of California, San Diego Mailcode 8208 UCSD Antiviral

Hello, I would like to know if there is a function in R that will test for normality and handle censored data sets. Currently, I evaluate each censored data set by the extent to which a normal scores plot approximate a straight line. For complete data sets I use shapiro.test(). Below is an example of a censored data set. data1<-c(0.00, 0.00, 0.00, 5.86, 5.17, 8.17, 5.12, 4.92, 7.08,

Rhelpers: I'd like to modify this RSQLite statement: rs_stations<-dbSendQuery(con_stations, "select * from stations") so that stations is actually an R variable, e.g.: stations=c("stationA","stationB") How would I modify the above statement to query from stations[[1]] (aka "stationA")? --j

Dear Rhelpers, Is there a faster way than below to set a vector based on values from another vector? I'd like to call a pre-existing function for this, but one which can also handle an arbitrarily large number of categories. Any ideas? Cat=c('a','a','a','b','b','b','a','a','b') # Categorical variable

Hello all, I need a help. I am modeling a disease and a create a R function like that: Lambda<-function (x,date1,r,h,a){ ndate1 <- as.Date(date1, "%d/%m/%Y") t1 <- as.numeric(ndate1) x[order(x$i),] t <-x[,"t"] i <-x[,"i"] CONTAGIEUX <-x[,"CONTAGIEUX"] while ( t1 < min(t) ){ for (i in 1:length(i) ){ {for (j in

Hi Rhelpers, I would like to have some help with a plot which is beyond my capabilities. This plot that I am seeking involves an overlay of two different barcharts and one xyplot. The code that I have used is the following : #save(df1,file="M:\\KBR\\df1.RData") load(file="M:\\KBR\\df1.RData") # df1$Year.ord created to obtain the right order i.e. 2015M < 2015K

Rhelpers: I'm trying to make a barchart of a 2-group dataset (barchart(x~y,data=data,groups=z,horizontal=FALSE)). My problem is that I can't, for the life of me, seem to get rid of the inter-bar space -- box.ratio set to 10000 doesn't do much. Any ideas? I'd ideally want zero space between the bars. Thanks! --j -- Jonathan A. Greenberg, PhD Assistant Project Scientist

Rhelpers: I recently installed the 64-bit version of R on my Debian system, and afterwards was asked if it was compiled using ATLAS. Is there a way to test to see if R is using ATLAS? --j

Rhelpers: Having a problem solving this. I have an xyplot call that looks like this: print(xyplot(temp_species_EAM_Pred_Pop$x+temp_species_NULL_Pred_Pop$x~temp_species_EAM_Pred_Pop$Action,main=current_species, xlab="Action",ylab="Predicted Pop", xlim=c(xmin,xmax),ylim=c(ymin,ymax),

Hello All, Would like to be able to convert a matrix to a dataframe without losing an assigned dimname. Here is an example that should illustrate what I'm talking about. tableData <- state.x77[c(7, 38, 20, 46), c(7, 1, 8)] names(dimnames(tableData)) <- c("State", "") tableData State Frost Population Area Connecticut 139 3100 4862

Hi all, I have lots of samba servers and want them to all authenticate against a single password file which will be a LDAP database. Where I want to be: Login to a domain called "FROST" that passes the username:passwd to a domain called "ACR" which checks LDAP. Where I'm at: I can login (from a W2K machine) to a test server called LIBIT which uses LDAP

Hello all, I've become confused by the output produced by a call to anova(model1,model2). First a brief background. My model used to predict final tree height is summarised here: Df Sum Sq Mean Sq F value Pr(>F) Treatment 2 748.35 374.17 21.3096 7.123e-06 *** HeightInitial 1 0.31 0.31 0.0178 0.89519

My system: Suse 9.3 Prof. I downloaded Wine: wine-0.9.7-SuSELinux100.i586.rpm. I get the message libstdc++.so.6 not available. What can I do? Please note: I am not a linux expert. What little I know, I acquired by myself. Regards, suren -- The woods are lovely dark and deep But I have promises to keep And miles to go before I sleep And miles to go before I sleep (Robert Frost)

Hello, I am having problems installing RSQLite under FC8. The error message seems to indicate Rversion.h and Rconfig.h are missing. To illustrate the problem, I have copied below the output shown when I first install R (binary distribution) and then try to install RSQLite. Because the problem does not occur when installing RSQLite with R under Windows I suspect the problem is perhaps due to my

Can anyone recommend a database modelling tool for the Mac? I''m building a new Rails app and want an interface to draw out the database''s entity relationship diagram, setup foreign keys and then output DDL at the end to MySQL. I''ve previously used DBDesigner, CASE Studio and some pricier enterprise software on Windows, but can''t find anything suitable for my

Hi, i have a problem with xen. the output of uname is:Linux xen 2.6.12-gentoo-r6 #9 SMP Mon Jan 23 18:13:25 CET 2006 i686 AMD Opteron(tm) Processor 844 AuthenticAMD GNU/Linux i user the SAME config of the kernel, for xen0 domain, and when i boot with xen0 kernel it reboots after this output... . . . ipt_recent v0.3.1: Stephen Frost <sfrost@snowman.net>.

I'm trying to plot a time series in ggplot, but a date column in my data frame is causing errors. Rather than provide my own data, I'll just refer to the scale_date example at: http://had.co.nz/ggplot2/scale_date.html , which reproduces the error. > df <- data.frame( date = seq(Sys.Date(), len=100, by="1 day")[sample(100, 50)], price = runif(50) ) > dt <- qplot(date,