similar to: converting a BMDP 8V mixed model to R / nlme

I am trying to translate a given BMDP 8V problem treatment to something approximately equivalent in nlme package (or base R, if sufficient). (Of course, it is not my goal to end there. I want to get access to the advanced flexibility of a more modern treatment.) Here is the problem statement in BMDP control language: /input title='Augenbewegungen'. variables=4.

I used both BMDP and SAS in my earlier years, side by side. At that time the BMDP statistical methods were much more mature and comprehensive: we treated them as the standard when the two packages disagreed. (It was a BMDP manual that clearly explained to me what the hypothesis of "Yate's weighted mean test" is, something SAS decided to call "type III" and eternally

I have a problem moving from multistratum aov analysis to lmer. My dataset has observations of ampl at 4 levels of gapf and 2 levels of bl on 6 subjects levels VP, with 2 replicates wg each, and is balanced. Here is the summary of this set with aov: >> summary(aov(ampl~gapf*bl+Error(VP/(bl*gapf)),hframe2)) > >Error: VP > Df Sum Sq Mean Sq F value Pr(>F) >Residuals

-------------------------------------------------- Eltemette az elef?nt-trutyi Paderborn, N?metorsz?g Friedrich Riesfeldt ?llatkerti gondoz? sz?kreked?ses elef?ntj?nak beadott 22 adag ?llati hashajt?t ?s megetette egy v?ka szederrel, f?g?vel ?s szilv?val. A bedugult vastagb?r? v?g?l els?lt - ?s megfojtotta a gondoz?t 200 font trutyival. A vizsg?lat szerint a 46 ?ves Friedrich megpr?b?lt a beteg

I'm looking for a solution to provide a live stream of the debate on Sofware Patents in Monday's EU Council of Ministers [1]. Here's what I've in mind: * Council --> mobile phone --> software answering machine --> Internet server * Council --> mobile phone --> VOIP --> Internet server Has anyone ever streamed something via a cell phone? How would you realize

I have design with two repeated-measures factor, and no grouping factor. I can analyze the dataset successfully in other software, including my legacy DOS version BMDP, and R's 'aov' function. I would like to use 'Anova' in 'car' in order to obtain the sphericity tests and the H-F corrected p-values. I do not believe the data are truly deficient in rank. I

Scot, Thanks for the info. I will try your code out to verify the result, but before I do that, will your code (SAS and R) work with variable names that are longer than 8 characters long without truncating the variable name in R? Also, I wonder about using your method or the PROC EXPORT method with larger data sets. The data sets I will be working with for the most part will not be that large,

Hi All, I've been fiddling around with various ways to estimate the popularity of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS, R-PLUS and Revolution R. It's not an easy task. You can see what I've come up with so far at http://r4stats.com/popularity . I'm sure people will have plenty of ideas on how to improve this, so please let me know what you think.

<headline>Short question</headline>: Do people have advice on debugging R programs running after CGIwithR inputting of data from forms? Is there a way of setting up fast local versions if your local machine has to be a windoze (2k) machine (R 2.1.0) and your server is a Debian, ssh shell only set up running R 1.8.0? Are there simple guides to ways of not having to invoke R each

Well, you don't give much of an example.... I'm replying CC to the R mailing list. Please ask questions there, rather than adressing individuals for basic help. Here is one; does it answer your question ? data(agriculture) ag.ag <- agnes(agriculture) class(ag.ag) pltree(ag.ag) # the dendrogram, if you want to see it ## cut the dendrogram -> get cluster assignments: (ck3 <-

Under linux/bash, it would be something like.... cd /MusicDirectory find . -type f -name "*.flac" -exec sh -c 'flac -t {} && flac -8V {}' \; Or to simplify and not test each file before re-encode.... find . -type f -name "*.flac" -exec flac -8 {} \; -print Can someone translate this into NTBatch? Basically, recursive "dir" *.flac files as a

I am including my config.log in the hopes that someone can tell me what is going on with dovecot installing on Mac OSX 10.2.6 Sincerely, ? Roger Cates, CCNA Vice President & Chief Technical Officer Xpower Internet, LLC Xpowerhosting.com | Xpoweronline.com P 888.245.7501 | F 270.338.4602 Internet to the power of X. ? begin 666 dovecot_config.log

Hi I've a Sangoma A200D with 2FXO and 2FXS. When using it with only the FXO module, it's all good. But when I put in the FXS module and connect the power, logs tells me not enough power. > Mar 31 14:11:54 phone kernel: [ 4761.246931] wanpipe1: Module 1: > Failed to powerup within 600 ms (8V : 72V)! > Mar 31 14:11:54 phone kernel: [ 4761.246937] wanpipe1: Module 1: Did

Hello, I am confused by a warning from mkswap : When making a swap on a LVM volume, I see the following warning : mkswap: /dev/vg_SDB1/swap_test6_64: warning: don't erase bootbits sectors on whole disk. Use -f to force. Setting up swapspace version 1, size = 1048572 KiB Is it safe to use mkswap -f in this case ? (If that matters, the swap is intended for a virtual machine, and the

Dear R users, I'm applying a correlation structure in a mixed model (nmle function) to control for spatial correlation between land parcels that are adjacent to each other. I generated X,Y coordinates in ArcGIS for each land parcel and used them in the correlation form like this: test.exp<-corExp(1, form = ~ X + Y) test.exp<- Initialize(test.exp,dataset) However, the correlation

I found this: >I would add -V and log all output to a file. also, someone has >reported that if the input FLAC file has errors and you use --force, >the transcoded file will only have data up to the error, and I still >have to look into this. so safer still would be>> > >... -exec sh -c 'flac -t {} && flac -8V {}' \; ... > >Josh > I wish I

Hi I am using R 1.7.1 on RH linux 9.0 > sum(unlist(lapply(ls(),function(x)object.size(get(x)))))/1024^2 [1] 2.424263 so I am not using much memory (I have a gig of ram on my machine) now in nlme > gtest<-groupedData(log(X8)~Time|sub,all[,c(names(all)[1:9],"X8")],outer=~A*B) > object.size(gtest)/1024 [1] 59.98438 > plot(gtest,outer=~Dose*chem,key=FALSE,asp=.5) Plotting

Following the Indomethicin example in Pinheiro & Bates, chapter 6, page 277 etc, coming to the following comand: fm2Indom.nlme <- update( fm1Indom.nlme, random = pdDiag(A1 + lrc1 + A2 ~ 1) ) debugging nlme gives the following output: Browse[1]> n debug: modelResid <- ~eval(model, data.frame(data, getParsNlme(plist, fmap, rmapRel, bmap, groups, beta, bvec, b, level,

Hello. I have 16 subjects with 1-4 obs per subject. I am using the package "nlme" to fit a simple random effects (variance components model) with 3 parameters: overall mean (fixed effect), between subject variance (random) and within subject variance (random). I need a 3x3 variance-covariance matrix that includes all 3 parameters in order to compute the variance of a specific

Hi, All: Does 'nlme' require the nonlinear function to be differentiable? I'm fitting structural change models to many related time series, and I'd like to improve the estimates of the change point times through some type of pooling. Unfortunately, I've so far been unable to get 'nlme' to work for me. The following toy example is the closest I've come