Displaying 20 results from an estimated 900 matches similar to: "[R-sig-ME] interpretation of main effect when interaction term being significant (ex. lme)"
2012 Jun 06
2
Main effects and interactions in mixed linear models
Dear all,
This question may be too basic quesition for this list, but if someone has
time to answer I will be happy. I have tried to find out, but haven't found
a consice answer.
As an example I use "Pinheiro, J. C. & Bates, D. M. 2000. Mixed-effects
models in S and S-PLUS. Springer, New York." page 225, where rats are fed
by 3 different diets over time, which body mass has
2006 Mar 13
2
Error Message from Variogram.lme Example
When I try to run the example from Variogram with an lme object, I get
an error (although summary works):
R : Copyright 2005, The R Foundation for Statistical Computing
Version 2.2.1 (2005-12-20 r36812)
ISBN 3-900051-07-0
...
> fm1 <- lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat)
Error: couldn't find function "lme"
> Variogram(fm1, form = ~ Time | Rat, nint =
2002 Feb 13
1
nlme package dependencies?
Dear R-help,
The following nlme example, from help(Variogram.gls), does not work on my
system.
> data(BodyWeight)
> fm1 <- gls(weight ~ Time * Diet, BodyWeight)
> Variogram(fm1, form = ~ Time | Rat)[1:10,]
Error in FUN(X[[1]], ...) : unused argument(s) (method ...)
It looks like I might not have loaded a necessary package. I did load nls
and lattice (the latter not actually being
2012 May 02
3
Consulta gráfica
Hola,
Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?
http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5
Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.
Muchas gracias.
Eva
[[alternative HTML
2015 Oct 27
2
pregunta
Otras variantes con y sin paquetes adicionales...
> sapply(split(datIn$Gain, as.factor(datIn$Diet)), mean)
d1 d2 d3
280 278 312
> by(datIn$Gain, datIn$Diet, mean)
datIn$Diet: d1
[1] 280
--------------------------------------------------------------
datIn$Diet: d2
[1] 278
--------------------------------------------------------------
datIn$Diet: d3
[1] 312
>
> library(dplyr)
>
2011 Mar 01
2
Entering table with multiple columns & rows
Hi,
I'm having difficulty with getting a table to show with
multiple rows and columns. Below is the commands that I've typed in and
errors that I am getting. Thank you.
Laura
Table trying to enter:
Diet: Binger-yes: Binger-No: Total:
None 24 134 158
Healthy 9 52 61
Unhealthy 23 72 95
Dangerous 12 15 27
>
2015 Oct 28
2
pregunta
Me gusta la respuesta uqe has dado, pero si por ejemplo, alguno de los datos tiene datos faltantes, entonces devuelve NA.
He probado con:
sapply(split(datos$uno, as.factor(datos$dos)), mean(na.rm=TRUE))
pero da fallo.
¿Cómo se podría hacer para que devolviera además la media obviando los NA y que contara el numero de NA por categoria?
> Date: Wed, 28 Oct 2015 00:13:45 +0100
> From: cof
2013 Jan 12
1
panel failure in xyplot
Hello
I ran the code below but it said:
no object "'panel.xyplot.intermediate.hh'"
Please kindly advise how to modify the code.
thank you.
(It works with panel.bwplot.intermediate.hh)
Elaine
code
library(HH)
# data input
dataN <-read.csv("H:/R_data/Mig_bird_586.csv",header=T, row.names=1)
dim(dataN)
dataN[1,]
str(dataN)
diet.code <-
2013 Jan 12
1
color in xyplot
Hello
I want to draw a xyplot.
Its dots will have three colors: red for meat, green for vegetable, and
blue for both.
I used the code below but could not make the dot in the same group show the
same color.
Please kindly advise how to modify it.
Thank you.
code
library (lattice)
diet.code <- c("Herbivore", "Omnivore", "Carnivore")
Diet.colors <-
2008 Jan 22
1
anova function to test the difference between two coefficients in nlme package
Dear Dr. Bates, and R-help,
I've tried the anova function to test the difference between two
coefficients, as shown on page 225 of your book "Mixed Effects Models
in S and S-Plus (Statistics and Computing)".
When I type: anova( fm2BW.lme, L = c(TimeDiet2 = 1, TimeDiet3 = -1) )
I got the following error message:
Error: unexpected '=' in "anova( fm2BW.lme, L =
2008 Jul 30
1
bug in 'margins' behavior in reshape - cast
according to the documentation of the cast function in the reshape function,
I would expect this bit of code from the examples to calculate marginal
means over only the 'diet' variable.
#Chick weight example
names(ChickWeight) <- tolower(names(ChickWeight))
chick_m <- melt(ChickWeight, id=2:4, na.rm=TRUE)
cast(chick_m, diet + chick ~ time, mean, margins="diet")
But,
2010 Mar 16
1
simple line graphics, labels and legend
Dear users,
I think my questions are pretty simple, but I got lost in the hundreds
of par() and plot() arguments and plot functions, so I don't know in
which direction I should go.
Here is my sample dataset:
test <- structure(list(DIET = structure(c(1L, 1L, 1L, 1L, 1L, 3L, 3L,
3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L), .Label = c("G",
"GG", "L",
2009 Aug 11
1
nested repeated measures MANOVA using adonis
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
2007 Feb 09
0
cluster "non-diet", "diet" example.
Suppose, we have 3 people called: Francis, Cedric and Nina. Base on what
they have eaten, we want to cluster people by "diet", "non-diet".
# original data file, named as filename "food.csv".
Francis|potato
Francis|chocolate
Francis|chocolate
Francis|milk
Cedric|vegetable
Cedric|vegetable
Cedric|potato
Nina|potato
Nina|chocolate
Nina|chocolate
Nina|potato
# Step 1: I
2018 May 04
2
Regression model fitting
Hi all ,
I have a dataframe (Hypertension) with following headers :-
> Hypertension
ID Hypertension(before drug A) Hypertension(On drug A) On drug B? Healthy diet?
1 160 90 True True
2 190
2004 Aug 26
1
Error TukeyHSD
I am running the following code on the coagulation data and I am getting
an error. Please let me know
if I am missing anything from my code.
coag<- matrix( scan("//Samba3/nair/R/blood.dat", sep=","), 24, 3,
byrow=TRUE)
colnames(coag) <- c("time","diet","order")
coag <- as.data.frame(coag)
oneway.test(time ~ diet, data=coag,
2009 Jul 06
2
ReShape chicks example - line plots
Hi,
In the examples from the ReShape package there is a simple example
of using melt followed by cast that produces a smallish amount of
output about the chicks database. Here's the code:
library(reshape)
names(ChickWeight) <- tolower(names(ChickWeight))
chick_m <- melt(ChickWeight, id=2:4, na.rm=TRUE)
DietResults <- cast(chick_m, diet + chick ~ time)
DietResults
My challenge
2007 Aug 09
1
strange warning
Hi all,
I am using an asterisk as a client to connect to another asterisk server by
registering with the register string. Registration is done without any
hassel, but after sometime my asterisk loses the registration with the
server and the server starts displaying the following msgs repeatedly:
[Aug 9 06:37:59] NOTICE[8380]: chan_sip.c:8151 check_auth: Correct auth,
but based on stale nonce
2010 Sep 10
1
lmer output
Hi
I have a question regarding an output of a binomial lmer-model.
The model is as follows:
lmer(y~diet * day * female + (day|female),family=binomial)
The corresponding output is:
Generalized linear mixed model fit by the Laplace approximation
Formula: y ~ diet * day * female + (day | female)
AIC BIC logLik deviance
1084 1136 -531.1 1062
Random effects:
Groups Name Variance
2005 Apr 23
2
I already used --recursive
In /usr/share/doc/rsync/README.gz:
BUG REPORTS
-----------
If you have web access then please look at...
Well what if you don't have web access? You don't mention what to do
then. "Send bugs to rsync@lists.samba.org" probably. OK. Will do.
Here's one now:
$ rsync -avz debian.linux.org.tw:var/lib/apt/diet.lists var/lib/apt/diet.lists
receiving file list ...