similar to: repeated measures nonlinear regression

Dear group, I tried to run a Repeated Mesures Anova for Mixed effects model and I got a warnning after entering the model specification saying: "Reached total allocation of 254Mb: see help(memory.size)". here is part of the log: *********************************************************** > aphids<-read.table("aphid.txt",header=T) > attach(aphids) > names(aphids)

Dear R-users, I have 3 plant populations (fixed). Within each population there is the same number of “families” (random) – the seed progeny of the same plant. These families were exposed to 2 treatments (fixed) and their response was measured (mean values for 25 seedlings per family per treatment are presented in data table). I would like to know if there is a significant difference in the

Hi all - I have the latest samba up and running fine except for one small problem. The username I use to log into the NT domain is different than the username I log into samba with. I have no control over the former, and would *really* rather not change the latter. I've searched the net, and the archives, and got only saw bits and pieces of this issue, never quite fully answered. Is it

Hello, I read a book about statistics in psychology. The authors use SPSS. They talk about post hoc tests after ANOVA finds significant effects: - Gabriel's procedure (for equal or slightly different sample sizes) - Hochberg's GT2 (for different sample sizes) - Games-Howell procedure (for populations with unequal variances) I could not find them in R. Do they not exist in R

How to fit a nonlinear regression model with power exponentially distributed errors? I know gnlm has a function gnlr3 that could work, but I would be grateful if example R code is provided. Daniel [[alternative HTML version deleted]]

Hi all- I'm new to both these lists so if this isn't the right place, point me in the right direction. I've got a Sun that I want to use for file services for our desktops (Mac, W95). I was planning on using netatalk and samba until I thought about the lack of file locking b/n the two. So, I started looking at other options (like Syntax's Totalnet), but it's quite a bit of

Hi all - I'll be the first to admit I don't know a lot about Windows networking, but I've managed to get samba (latest version) up and running and even have a cool netlogon script :) My config has shares for [homes] and [netlogon] and [foobar] (a test). When I look around in Network Neighborhood it lists homes, netlogon, foobar, and philip (my username). Is there some way to make

Dear list members, I have a question concerning multiple comparisons after using glm. My response variable is days until emergence of an insect species. The explanatory variables are sex (two levels), parasitoids added (two levels) and populations (34 levels). I would like to know now which populations are different in days until insect emergence. For this I used multiple comparisons as

Hi - I'll be the first to say I don't completely understand Windows networking, but I did get the following to work when everybody (clients-w95/server-solaris) were on the same subnet. W95 boots up and prompts for username,password, and domain. I log in, and presto I'm connected and netlogon scripts run and everything is great. This is when I'm testing it on a completely

Dear R members, I have been trying to find out whether one can use multivariate regression trees (for example mvpart) to analyze repeated measures data. As a non-parametric technique, CART is insensitive to most of the assumptions of parametric regression, but repeated measures data raises the issue of the independence of several data points measured on the same subject, or from the same plot

Dear R Community, I am happy to introduce the next release of the *dclone* package (see R Journal, 2(2): 29-37, 2010) which now includes a suite of functions for parallel computations: - see functions parJagsModel, parUpdate, parCodaSamples using 'snow' type clusters; and the - jags.parfit function that can be used with 'snow' type clusters or multiple cores. - iterative model

Dear R Community, I am happy to introduce the next release of the *dclone* package (see R Journal, 2(2): 29-37, 2010) which now includes a suite of functions for parallel computations: - see functions parJagsModel, parUpdate, parCodaSamples using 'snow' type clusters; and the - jags.parfit function that can be used with 'snow' type clusters or multiple cores. - iterative model

Dear R Community, I am happy to introduce the latest version 1.3-0 of the 'dclone' R package. The package provides low level functions for implementing maximum likelihood estimating procedures for complex models using data cloning and Bayesian Markov chain Monte Carlo methods with support for JAGS, WinBUGS and OpenBUGS. Data cloning is a global optimization approach and a variant of

Dear R Community, I am happy to introduce the latest version 1.3-0 of the 'dclone' R package. The package provides low level functions for implementing maximum likelihood estimating procedures for complex models using data cloning and Bayesian Markov chain Monte Carlo methods with support for JAGS, WinBUGS and OpenBUGS. Data cloning is a global optimization approach and a variant of

Hi, I would like to obtain the dispersion parameter for a quasipoisson model for later use in calculating QAIC values for model comparison.Can anyone suggest a method of how to go about doing this? The idea I have now is that I could use the residual deviance divided by the residual degrees of freedom to obtain the dispersion parameter. The residual deviance is available in the summary

Hi all, I'm trying to calculate a nonlinear regression to get a simpler expression for a complicated formula. I expect that a function of the type 1 ---------------------------- - (a + b1*x + b2*y) 1 + e should be a good fit. Now I'm trying to calculate a regression with R and must admit that I'm slightly confused by the variety of possibilities to

I have a complicated nonlinear function, myfun(a,b,c), that I want to fit to data, allowing one or more of the parameters a, b, and c in turn to have linear dependence on other covariates. In other words, I'd like to specify something like nls(y~myfun(a,b,c),linear=list(a~f1,b~1,c~1)) I know would this work in nlme *if I wanted to specify random effects as well*, but I don't -- and

Hello, Does anyone know if there is any package in R to perform wombling (sensu Barbujani et al. 1990) or a related method? RSiteSearch only gives 2 results with "wombling" and I am not aware of any other name for the method. Some weeks ago there was a related question on the list but there was no response... I acknowledge that programming it wouldn't be so difficult (only

Dear colleagues, Von Bertalanffy model is commonly adjust to data on fish length (TL) and age (AGE) TL= Linf*(1-exp(-K*(AGE-t0)). Linf, K and t0 are parameters of the model. One main goal of the growth study is the comparison of growth parameter estimates between sexes of the same species, or estimates from different populations. The realibility statistical tests normally applied are highly

Dear All I am kind of stuck up with a code a part of which seems to be causing a problem, or at least I think so. May be the community can help me. It’s simple but I suppose I am missing something. I generate a data matrix X, say of order n*p, where n represents independent row-vectors and p correlated col vectors. Let the row representation be X = (X’_1, . . ., X’_n)’. I generate the