similar to: converting BMDP 8V mixed model treatment to R ?

[Sorry, In case this is repeating a message already sent to the list]. I am trying to move a project to R (base or nlme), for which I have a partial solution in BMDP 8V. Here is the 8V control language: /input title='Augenbewegungen'. variables=4. file='latm.dat'. format='11x,f7.0,f7.0,f12.4,f10.0'. /variables names= llat,rlat,vg,diff. /design

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

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

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

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 wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

Data from Fisher's paper: Confidence Limits for a Cross-Product Ratio. > y col1 col2 [1,] 10 3 [2,] 2 15 fisher.test(y) Fisher's Exact Test for Count Data data: y p-value = 0.0005367 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 2.753438 300.682787 sample estimates: odds ratio 21.30533 The crude odds

Hi there, I need to do the same thing as cumsum but with the variance and skewness. I have tried to do a loop for like this: var.value <- vector(mode = "numeric", length = length(daily)) for (i in (1:length(daily))) { var.value[i] <- var(daily[1:i]) } But because my dataset is so huge, I run out of memory..... Any ideas?!?! Much appreciate

SAS has DUD (Does not Use Derivatives)/Secant Method for nonlinear regression, does R offer this option for nonlinear regression? I have read the helpfile for nls() and could not find such option, any suggestion? Thanks, Derek [[alternative HTML version deleted]]

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

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

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

I’m not sure if it will help you but we’re seeing something similar that I call the “10 hours problem” here. On any Samba - after 4.7.6 (we’re currently trying out 4.9.4 but still no go) - at 10 hours after it (smbd & winbindd) was restarted it new connections often (but not 100% all the time, but we think it is related to the number of connections, or active connections, or how many are

********************************************************************** Este email assim como os ficheiros que possa ter em anexo s?o confidenciais e para uso exclusivo da pessoa ou organiza??o para o qual foi enviado. Se recebeu este email por engano por favor notifique Redes@bnc.pt Esta nota confirma que esta mensagem foi verificada pelo MIMEsweeper n?o tendo sido encontrados virus.

All, using Samba as an AD (2k12) domain member in Stretch (2:4.5.12+dfsg-2+deb9u4) with tdb as default and rid as domain backend. No overlapping. Everything works fine. Setup was done as in the wiki [1]. If you're connecting from a Windows 10 client and do not add dedicated keytab file = /etc/krb5.keytab kerberos method = secrets and keytab winbind refresh tickets = Yes to

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