similar to: lme for repeated measures, one within, one between factor

Hello and thanks in advance I am running a glm in R the code is as follows with residual diagnostic code below model4<-glm(Biomass~(Treatment+Time+Site)^2, data=bobB, family=quasi(link="log", variance="mu")) par(mfrow=c(2,2)) plot(model2) to test the effect of grazing exclusion of feral horses for a Phd with following factors: Treatment - 3 levels which are grazed

Dear R help group, I am teaching myself linear mixed models with missing data since I would like to analyze a stats design with these kind of models. The textbook example is for the procedure "proc MIXED" in SAS, but I would like to know if there is an equivalent in R. This example only includes two time-measurements across subjects (a t-test "with missing values"), but I

Hi Everyone, I’m trying to figure out how to get R to analyze this experiment properly. I have a series of subjects each with two legs. Within each leg there are two bones that I am interested in. There are also two treatments that I am interested in. That results in four different combinations of treatments. Obviously, since the subjects only have two legs, they can’t receive each

Sorry, I just posted the email below but realised I did not give a name or details, would it be possible to adjust before posting and send what is below, sorry again, first time user... From: Joanne Lenehan [mailto:jlenehan@une.edu.au] Sent: Tuesday, 1 December 2009 3:51 PM To: 'r-help@r-project.org' Subject: GLM Repeated measures test of assumptions: e.g. test for sphericity e.g.

Dear list, I have a dataset where we sampled multiple individuals either 1 or 9 times. Our measurement variable is 'Delta13C' (see below sample dataset). I cannot figure out how to efficiently use a vector command (preferably) or a loop to create a new vector of the names of the individuals sampled 9 times. Note that the 'FeatherPosition' variable will only be "P1" for

I suppose that multcomp in R and multicomp in S-Plus are related and it appears that it is possible to use multicomp with lme in S-Plus given the following correspondence on s-news sally.rodriguez at philips.com 12:57 p.m. 24/04/03 -0400 7 [S] LME summary and multicomp.default() Is it possible to use multicomp with lme in R and if so what is the syntax from a simple readily available

Dear list, I have 9 repeated measures (measurement variable == 'Delta13C') for individuals (ID variable == 'Individual_ID'. Each repeated measure is "indexed" (right term?) by the variable 'FeatherPosition' and given as c('P1', 'P2', 'P3', 'P4', 'P5', 'P6', 'P7', 'P8', 'P9'). I would like to

hello, i found it most convenient to use package contrast for planned comparisons on mixed models. for instance i have a model with 2 fixed factors, one with 4 levels (stage) and one with 2 levels (gap) and a nested random factor (site) and i tested gap within level A of factor stage: library(contrast) library(nlme) m1<-lme(rich ~ stage*gap, random=~1|site,data=richness) contrast(m1, a =

Hi, Although I set check.names to FALSE in read.table, the duplicate names get modified. What should be done in this case? the text file to be read by read.table AM2      AM2     AM2     FAL 2               3             4         5 1               -1           -3        -2 t = read.table ("my_file", check.names = FALSE, header = T) > t AM2      AM2.1     AM2.2     FAL

HI everyone, I try to get some bootstrap CIs for coefficients obtained by quantile regression. I have influencial values and thus switched to quantreg.. The data is clustered and within clusters the variance of my DV = 0.. Is this sensible for the below data? And what about the warnings? Thanks in advance for any guidance, Kay > dput(d) structure(list(Porenfläche = c(4990L, 7002L, 7558L,

Hello Jeremy, I think you missed a "/" in: mount --bind /share1/DVR /medium/DVR mount --bind /share1/Audio /medium/Audio mount --bind /share2/Video /medium/Video mount --bind /share2/Backups /medium/Backups BR, Matt On Thu, Feb 13, 2020 at 5:15 PM Jeremy Allison via samba < samba at lists.samba.org> wrote: > On Thu, Feb 13, 2020 at 04:45:16PM -0500, Hal Vaughan via samba

Hello list, I' d very much appreciate some help with a two sample repeated measures ANOVA. I did the analysis yielding sign. main effects (between subj.=site, within subj.=cover) and a sign. interaction: Univariate Type II Repeated-Measures ANOVA Assuming Sphericity SS num Df Error SS den Df F Pr(>F) site 18.7620 1 18.831 10 9.9631 0.010220 *

On Thu, Feb 13, 2020 at 04:45:16PM -0500, Hal Vaughan via samba wrote: > I haven?t used Samba for over a decade and, from what I remember, I don?t think this is possible. Whenever I search for "two mount points" or "multiple mount points," I get hits for "multiple users," but not for mount points. > > I have the following directories I want to share through

This question teeters on the brink of general statistics rather than an R-specific issue, but I'm hoping it boils down to how I coded V&R's EM algorithm for finite mixture modeling and not something worse (i.e. that my plan is impossible). Briefly, I'm trying to simulate an experiment where the hypothesis concerns bimodal vs. unimodal distributions, in order to see what

Dear list members, First of all thank you for your helpful advices. After your answeres to my firt mail I studied a lot (R-News n?5) and I tried to perform my analysis: First, to fit a GLM with a nested design I decided to use the function "lmer" in package "lme4" as suggested by Spencer Graves and Filippo Piro. I remember to you that my data were: land use classes, 3 levels

Hi Lists, sorry for reposting this, but since no response has arrived so far and my time on this project is passing by I need to make a conclusion. Do I better use my DVR cards in dom0 or moving it to a separate physical machine rather than using them in a domU? Could someone shed light on why it is trickier to pass-thru one of these cards than any other pci devices? Do they try to read-write

Hi Lists, sorry for reposting this, but since no response has arrived so far and my time on this project is passing by I need to make a conclusion. Do I better use my DVR cards in dom0 or moving it to a separate physical machine rather than using them in a domU? Could someone shed light on why it is trickier to pass-thru one of these cards than any other pci devices? Do they try to read-write

I am using R for a class which uses S. Often I can not find functions that in S but not in R, example, multicomp(). What is your advise. Thanks. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body",

Dear R people, I have heard rumors of a function called (I believe) multicomp in S+ for doing multiple comparisons for linear models. It doesn't seem to be in the version of R I am using, 0.90.0. Is there something equivalent to this in R, or does it simply not exist? I see there is are more recent versions of R. For my platform, SuSE Linux, the most recent rpm available appears to be

Does anyone have any "security appliance" DVR, and can you a) telnet or b) ssh into it? Also, I have read that on a lot of the low-end ones, like Swann, port 9000 is completely open, but it's not clear to me from the hack post <http://console-cowboys.blogspot.com/2013/01/swann-song-dvr-insecurity.html> whether I could ssh in there, or rsync from a server through that port.