similar to: SAMM package for mixed models

First, a big thanks to all of the developers and users that have worked to make R such useful software. It is only because I find the software so useful that I have the following opinions. A recent post to R-devel listed the 'Top 10 Features' for one person. I found it to be quite an interesting read. Over the past couple of years I have assembled my own lists. Retrospective. Some of

Years ago, when I used Splus, I used to do this sort of thing: assign("data", data, frame = 1) This was so that if the object data was used in, say a call to glm, and I wished to do a summary of the glm object, data wasn't otherwise visible. Since I moved to R, the lexical scoping has virtually eliminated the need to do that. However, I've been trying to use the package samm

Hi all, I need to run a program (asreml) thousands of times and each time I have to provide a slightly different dataset. Because I have to run asreml under Windows (DOS or scripts) environment, I have trouble to pass a command variable (or pointer variable or %counter in the following example) to R program so at each counter R can generate a different data set for asreml to run. Any

Dear R users,   I am trying to use "wgaim" package for QTL analysis using mixed model approach. But i am stuck with "asreml" function while using "wgaim" package. Do i need a separate package to activate "asreml" function beside "wgaim" package ?   If so, i tried to download "Asreml-R" package (i guess this is right package

Dear Sirs I would like to make you aware of the following training taking place in Chicago next month. Analysis of experiments using ASReml-R, Chicago, 24/25th April This workshop is aimed at scientist/practitioners that are interested in analyzing complex datasets by fitting linear mixed models, particularly users with experience in breeding programs and design and analysis of experiments.

How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Cue Direction (2) x Brain Hemisphere(2) Each of these has 2 levels, 'left' and 'right', so it's a simple 2x2 design matrix. We have 8 subjects in each cell (a balanced design) and we want to specify the interaction contrast so that: CueLeft>CueRght

Dear List, I have a richness data distributing across 20 N to 20 S latitude. (120 E-140 E longitude). I would like to draw the richness in the north hemisphere and a regression line in the plot (x-axis: latitude, y-axis: richness in the north hemisphere). The above demand is done using plot. Then, south hemisphere richness and regression are required to be generated using the same y-axis above

Hi Everyone, I am using ASReml to fit a spatial model. I do not have all the components of ASReml when I call; names(summary()) e.g. names(summary(fcov.asr2)) [1] "call" "loglik" "nedf" "sigma" "varcomp" I am trying to get the coefs but I get "NULL". Does anybody know the reason? Any help would be much appreciated. Regards

Hi all, I have a strange problem and rigth now I can't figure out a solution. Trying to calculate an ANOVA with one between subject factor (group) and one within (hemisphere). My dependent variable is source localization (data). My N = 25. My data.frame looks like this: > ML.dist.stack subj group hemisphere data 1 1 tin left 0.7460840 2 2 tin left

Why do I see this error? > library(stats) > require(stats) [1] TRUE > > tHSD <- TukeyHSD(aov) Error in TukeyHSD(aov) : no applicable method for "TukeyHSD" In case it helps: > aov Call: aov(formula = roi ~ (Cue * Hemisphere) + Error(Subject/(Cue * Hemisphere)), data = roiDataframe) Grand Mean: 8.195069 Stratum 1: Subject Terms: Residuals Sum

Estimados, Les tengo una pequeña consulta. Estoy trabajando con ASReml en Rstudio, me encuentro corriendo un modelo (CHL) para obtener el ANOVA de la interacción de los factores gen x env CHL<- asreml(fixed= MS~geno:env, random = ~rep, data = index) (interacción) Sin embargo, cuando ejecuto el modelo, R indica que siguiente mensaje: Current workspace: 128.000000Mb Warning message: Abnormal

Dear Experts, I am new to R and unfortunately cannot start with a simply statistical analysis: I manually determined the volume of the right and left hippocampus in a group of meditators and in a group of controls. My data-sheet looks as follows: observation subject group age gender hemisphere volume 1 am04 m 25 f left 3.637 2 am04 m 25 f right 3.713 3 ao08 m 47 m left 3.715 4 ao08 m 47

Hi, we need some help to define planned comparisons. I've based my understanding of the problem on reading Tabachnick and Fidell (2006), ie: http://www.ablongman.com/catalog/academic/product/0,1144,0205459382,00.html I don't understand how to specify planned comparisons in R. I've not found explanations for this in MASS or elsewhere. There is only discussion of the contrast

I'm trying to read in a CSV, with lines looking like: HEADER, Latitude DecDeg, Latitude Hemisphere, Longitude DecDeg, Longitude Hemisphere, Speed knots, Bearing Degrees, fixQualityGga, noOfSatGga, altGga, heightGga, selectionGsa, fixGsa, pdopGsa, hdopGsa, vdopGsa, noOfSatGsv, Time, *checkSum $GPS, 3747.0224, N, 12223.4522, W, 0.36, 348.21, 1, 3, 4.01, 175.5, -25.2, A, 2, 4.14,

Hello all. I've been running Samba for quite some time but am now trying to move it up from (some rather limited) file sharing into being a full PDC for our network. I figured out all the "hard stuff" but am having a heck of a time with the printing. What I need to make sure of is that certain people can pause/delete, etc print jobs for other users, as well as start or stop the

Hi I hope this is going to the right place. I am trying to write a program which uses KernSmooth library to estimate mutual information between two time series at various different lags. At the moment it’s producing negative values, which is supposed to be impossible (something is fishy). I am summing across one row of the matrix to get p(value is in bin x) and summing across the columns to get

Dear All: I want to do association study based on mixed linear model, My model not only includes serval fixed effects and random effects but also incorporates some covariates such as "birth weight". Otherwise, the size of the data are about 180 individuals and 12 variables and 60000 Fixed effect estimates As asreml-R is not free ,is there any packages for my study? I heard nlme or

With respect to calculating the epsilon index of sphericity for ANOVA, discussed on pp. 45-47 of: http://www.psych.upenn.edu/~baron/rpsych.pdf It notes that epsilon is not required for a repeated measures design with only k=2 levels, as the minimum value of epsilon (e) is given by: e = 1/(k-1) so for k=2, we have e = 1 (ie, no correction of the F test df; see p. 46). These notes apply to a

Full_Name: Arni Magnusson Version: 1.6.2 OS: Windows XP Submission from: (NULL) (210.48.49.68) Hi there. I'm experiencing unexpected behaviour from as.POSIXct: > as.POSIXct("1969-12-24") [1] "1969-12-23 23:00:00 New Zealand Standard Time" > as.POSIXlt("1969-12-24") [1] "1969-12-24" > as.POSIXlt("1969-12-24")+1 [1] "1969-12-23

Hi, I am looking to do some quantitative genetic analyses using animal models and was wondering if someone could suggest an appropriate package in R. It would help if it was similar to the ASReml genetic analyses software. Thanks, Deepa Senapathi Deepa Senapathi Centre for Agri-Environmental Research (CAER) School of Agriculture, Policy & Developement. University of Reading