similar to: Bivariate lme

Dear RHelpers, Does anyone know of any R code to perform bivariate Poisson regression (including random effects)? Best wishes Simon Simon D.W. Frost, M.A., D.Phil. Department of Pathology University of California, San Diego Antiviral Research Center (Formerly: UCSD Treatment Center) 150 W. Washington St., Suite 100 San Diego, CA 92103 USA Tel: +1 619 543 8080 x275 Fax: +1 619 298 0177 Email:

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

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 -

Hello. I've tried to implement dovecot's deliver LDA to my home mailserver. I've run into troubles :) dovecot-deliver.conf: -- snip -- mail = maildir:%h/Maildir:INBOX=maildir:%h/Maildir mbox_read_locks = flock mbox_write_locks = flock -- snip -- when i run: deliver -d bfg < email.txt dovecot delivers message to ~/Maildir/maildir/new It acctually creates "maildir"

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,

Full_Name: Simon Frost Version: 1.0.0 OS: Windows Submission from: (NULL) (129.215.60.79) Hi everyone, I've just installed R for Windows - I've previously been using it on Linux - and I've run into a problem when using lme. The lme routine throws up an error in the getGroups.data.frame command, complaining about an error in the formula: Error in getGroups.data.frame(dataMix,

hello! i'm using redhat 7.2 with ext3 as my primary fs on kernel 2.4.17 + grsecurity + acl after 2-3 days of uptime i'm expiriencing problems... i attached below excert from my system logs. machine stops responing for a few seconds and after then it looks, like it's in normal operation again. the only problem is load, which is incrementing constantly, but cpu is 99% idle...

hello! i'm using redhat 7.2 with ext3 as my primary fs on kernel 2.4.17 + grsecurity + acl after 2-3 days of uptime i'm expiriencing problems... i attached below excert from my system logs. machine stops responing for a few seconds and after then it looks, like it's in normal operation again. the only problem is load, which is incrementing constantly, but cpu is 99% idle... after

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

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, I am using sm.density() to find the bivariate joint PDFof events: For eg, x<-cbind(rnorm(30),rnorm(30)) den<-sm.density(x) Then I get the joint pdf from den$estimate in order to constructthe joint cumulative PDF. However, summing up all the values from den$estimateisnot equal to 1(have multipliedby the grid size). Anyone could help? Thanks. mc [[alternative HTML version

Dear all, I have a bivariate dataset from a preliminary study. I want to do two things: (1) estimate the probability density of this bivariate distribution using some nonparametric method (kernel, spline etc); (2) sample a big dataset from this bivariate distribution for a simulation study. Is there any good method or package I can use in R for my work? I don?t want parametric models like

Am trying to do some study on bivariate binomial distribution. Anyone knows if there is package in R that I can use to calculate the density function of bivariate binomial distribution and to generate random samples of it. Thanks, -- View this message in context: http://n4.nabble.com/How-to-calculate-density-function-of-Bivariate-binomial-distribution-tp992002p992002.html Sent from the R help

Any one know is there any package or function to generate bivariate exponential distribution? I gusee there should be three parameters, two rate parameters and one correlation parameter. I just did not find any function available on R. Any suggestion is appreciated. -- View this message in context:

Hello all!!! 1) I am wondering is there a way to generate random numbers in R for Bivariate Uniform distribution? 2) Does R have  built-in function for generating random numbers for any given bivariate distribution. Any help would be greatly appreciated !! Good day! Haneef Anver [[alternative HTML version deleted]]

x <- iris$Sepal.Length[1:50]/iris$Sepal.Width[1:50] y <- iris$Petal.Length[1:50]/iris$Petal.Width[1:50] I want to check whether (x,y) follows a bivariate normal distribution or not, using density plot or scatter plot. Is it possible to plot a bivariate density in R. I cant find any. Arindam Fadikar M.Stat Indian Statistical Institute. New Delhi, India [[alternative HTML version

Hello, I'm quite new with R and so I would like to know if there is a command to calculate an integral. In particular I simulated a bivariate normal distribution using these simple lines: rbivnorm <- function(n, # sample size mux, # expected value of x muy, # expected value of Y sigmax, # standard deviation of

Hello, I have used the univariate Ripley K function in R, but does anyone know if there is a bivariate function built in? I have two species that I am dealing with. Also, how might I add error bars into the graphs (univariate and/or bivariate)? Thank you, Karin Leiderman k_leiderman at hotmail.com Graduate Student/Research Assistant Department of Mathematics Univesity of New Mexico

Hi, I woul like to know if it is possible to have a "R code" to generate EM Algorithm for a normal bivariate mixture. Best regard, S.F.