similar to: 2 way repeated measures ANOVA using R: syntax and reporting question

Dear Gideon Unless you are willing to accept uncorrelated measures in time and homogeneity of variance I would go for a linear mixed model instead. You might be surprised on how diffrent can be your results. Try; >?lme I hope that this helps Francisco >From: GIDEON WASSERBERG <wasserberg at wisc.edu> >To: "R-help at lists.R-project.org" <R-help at

I am designing an interlaboratory validation study for a presence/absence alternative method test kit vs. a presence/absence reference method test kit. There will be 10 laboratories conducting tests using both methods. In each laboratory, there will be 5 specimens tested, each of the 5 specimens twice by both methods (alternative, standard). The total number of data are 10 x 5 x 4 = 200.

[Corrected the model formula to include "method".] I am designing an interlaboratory validation study for a presence/absence alternative method test kit vs. a presence/absence reference method test kit. There will be 10 laboratories conducting tests using both methods. In each laboratory, there will be 5 specimens tested, each of the 5 specimens twice by both methods (alternative,

Hi, I wonder if anyone can help me with specifying a right model for my analysis. I am a beginner to lme methods and though have spent already many hours studying from various books an on-line helps, I was unfortunately not able to find a solution to my problem on my own. Data structure: I studied escape behavior of three species of a prey to a predator. The prey specimens (many) were in a

Hi all, i have a dataset of 3D coordinates and can't figure out how to' tell R which One Are the individuals: I have 3 columns which i named x,y and z And then i have 2607 Lines, but each specimen is 33 Lines (79 specimens) How can i tell R to' categorize individuals every 33 Lines? Thanks in advance

Dear all, I'm building binomial mixed-model using lme4 package. I'm able to obtain outputs properly except when I include two particular variables: date (from 23 to 34; 1 being to first sampling day) and Latitude (UTM/100 000, from 55.42 to 56.53). No "NA" is any of those variables. In those cases, I get the warning message: "nlminb failed to converge" I tried to

I'm trying to create maps of reptile abundance in different states & counties using data from Herp.net, which provides lists of specimens with the places that they were found. First I would like to parse the list by state using 2-letter abbreviations, since I'm focusing on certain regions. To do this, I've been trying to create a vector (state2) that gives all state names as

*Thanks* all those who took the time to help me (even if the "question" was not related to - the use of - R). Now I think I can soundly make my point w/ the referee (can I use your replies? If so I intend to properly cite its use?!?). Regards, Eduardo Esteves ps - Sorry for not explaining the "biological details" of my posting: RNA/DNA is the ratio of RNA content to

Hello, Darik Horn sent a pull request adding support for LZ4. LZ4 is supposed to be a very fast compression algorithm, especially when it comes to decompression. I did a quick benchmark with zlib, LZO and LZ4. Now, there are many benchmarks you can find online, but most of them deal with compressing large files. Tinc on the other hand has to compress small packets individually. So I did the

Hi, I am encountering a problem in reading a file, the file looks like that: T 3 0 -- -- -- T -- -- -- 18.98 3 1 6.75 4.39 39 -- -- -- 18.58 3 2 6.90 4.90 43 -- -- -- 18.63 3 3 7.07 5.39 48 -- -- -- 18.78 G 4 0 7.41 5.54 47 G

I've enabled the GSS code in dovecot, but our Kerberos nerds are complaining that it doesn't work :) I probably have the thing totally misconfigured (so don't worry about that part for now), but I do have a crash: Info: imap-login: Disconnected: rip=XXX.YYY.229.8, lip= XXX.YYY.17.59, TLS handshake Error: auth(default): gssapi(?,XXX.YYY.229.8): While acquiring service credentials: No

Dear R-Helpers, I am working on a dataset containing ca. 500 specimens of 17 different species. The data is three text columns, followed by ca 70 columns of binary data. I'm trying to write a function that will allow me to interactively exclude species using the inofrmation in the third text column. My function, matrix.slicer, appears to work as I intend. However, when I pass the

On 06/10/2014 11:41 AM, Thomas Quinot wrote: > Dear all, > > I'm trying to boot a recent PXELINUX on a brand new Soekris net4501 > (comBIOS 1.33). This combination appears to not work out of the box... > > I see much advice floating around from a few years ago, hinting that > earlier versions of PXELINUX might have worked. However before resorting > to archaeology

* H. Peter Anvin, 2014-06-10 : > You probably want to test both pxelinux.0 and lpxelinux.0... if it is > sticking in pxe_init_isr you're using the latter. Correct. And indeed pxelinux.0 goes much further than lpxelinux.0! I have a few issues yet to fix: * Undef symbol FAIL: __syslinux_debug_enabled (I must have done something wrong while building) * find proper fix for

Dear Group Members, Forgive me if I am a little bit out of subject. I am looking for a good way to test the homogeneity of two variance-covariance matrices using R, prior to a Hotelling T test. Youll probably tell me that it is better to use a robust version of T, but I have no precise idea of the statistical behaviour of my variables, because they are parameters from the harmonics of

Hi, I got the following example from MASS (10.2 Classic Nested Designs) "A cooperative trial in analytical chemistry taken from Analytical Methods Committe (1987). Seven specimens were sent to six laboratories, each three times a month apart for duplicate analysis. The response is the concentration of (unspecified) analyte in g/kg". I understand what is written in the book about the

I try to fit a repeated measures model, in which: FA: The difference in the number of scales between the left and right sides for a given trait trait: the each of the measured features are 6 for each individual Condition: is the state of condition of each individual Individual: Each of the specimens studied By adjusting the model happens to me the following >

I am trying to apply a permuation-based MANOVA (Anderson 2001) to a set of morphological data from three ecomorphs of fish reared under two different conditions and measured at two points during ontogeny. I will supply a distance matrix based on Procrustes distances calculated outside of vegan. I have not found an example of a design such as this for adonis. However, I have designed my factors

Hi! I want to choose part of the data. But it has a question. > attach(PIKO) > PIKO[1:5,] code species dbh x y tag status branch 92 10142 PIKO 38.9 6.05 12.81 10165 alive 0 109 10213 PIKO 41.0 6.71 26.21 10202 alive 0 135 10222 PIKO 48.3 18.98 21.28 10214 alive 0 146 10223 PIKO 47.0 19.45 22.83 10224 alive 0 147 10223 PIKO 20.0 19.74

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =