search for: diet

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 > diet=matrix(c(24,134,9,52,23,72,12,15),ncol=4,byrow=TRUE) > rownames(diet)=c("none", "healthy", "unhealthy", "d...

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) > sum...

...Wed, 28 Oct 2015 00:13:45 +0100 > From: cof en qualityexcellence.es > To: jbetancourt en iscmc.cmw.sld.cu > CC: r-help-es en r-project.org > Subject: Re: [R-es] 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$...

...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 been measured. Response: Body mass, fixed effects Time*Diet, random effect ~Time|Rat. The main question is if the interaction term is significant (i.e. growth rate). My question is could I also look at the p-values of the main effects to say if body mass increase sig...

Estimados Cuando existia epicalc, había una manera muy fácil de determinar la media de una variable (en esta caso Gain) por grupos, en este caso (Diet). ?Como se puede hacer ahora? Diet Gain 1 d1 270 2 d1 300 3 d1 280 4 d1 280 5 d1 270 6 d2 290 7 d2 250 8 d2 280 9 d2 290 10 d2 280 11 d3 290 12 d3 340 13 d3 330 14 d3 300 15 d3 300 Saludos José [[alternative HTML version deleted]]

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 plan to reform the food.csv...

...diate.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 <- c("Herbivore", "Omnivore", "Carnivore") Diet.colors <- c("forestgreen", "darkblue","red2") levels(dataN$Diet) <- diet.code # graph par(mai=c(1.03,1.03,0.4,0.4)) par(new=T) xyplot(GE_distance~Diet, d...

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 <- c("forestgreen", "darkblue","red2") levels(dataN$Diet) <- diet.code # graph par(mai=c(1.03,1.03,0.4,0.4)) par(new=T) xyplot(GE_distance~Diet, d...

...lier shape. Based on the following code, the outliers are vertical lines. However, I want the outliers to be empty circles. Please kindly help how to modify the code to change the outlier shapes. Thank you. code package (lattice) dataN <- data.frame(GE_distance=rnorm(260), Diet_B=factor(rep(1:13, each=20))) Diet.colors <- c("forestgreen", "darkgreen","chocolate1","darkorange2", "sienna2","red2","firebrick3","saddlebrown","coral4", "chocola...

...les 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 is to extract an plot only a portion of this data. I would like to plot the data for each chick that participated in diet 1 only. Assume that the numbered column names (0,2,4, ...) represent time on the diet and will b...

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, summary(cast(subset(chick_m, time == 0), diet + chick ~ time, mean, margins=&quo...

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 140 False False 3...

...xample of a design such as this for adonis. However, I have designed my factors for this test which hopefully will work. This design is just to understand the principals of this test, my real data set is comprised of over a thousand specimens. Column 1 represents the ecomorph, column 2 respresents diet, column 3 represents family (nested in ecomorph). This example would have 9 individuals, but being that measures are repeated 1 1 5 1 2 5 1 2 6 2 1 7 2 2 7 2 2 8 3 1 9 3 2 9 3 1 10 1 1 5 1 2 5 1 2 6 2 1 7 2 2 7 2 2 8 3 1 9 3 2 9 3 1 10 my factor matrix is called Z, so the call to adonis i...

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", "LG"), class = "factor"), SILIKATG = c(0.0152, 0.0161, 0.0141, 0.0131, 0.0108, 0, 0, 0, 0, 1e-04, 0.0056, 0.0078, 0.0072, 0...

...ear model.?? We report the p-values of the main effects and intenraction from the anova table: anova(fm1BW.lmer).? Of course, it is possible to interpret the main effects and interaction from summary table as ANOVA is only a special case of regression.? Time estimates at individual levels of factor Diet (in this case, baseline Diet0) is okay for interpretation.? A.K. ----- Original Message ----- From: David Winsemius <dwinsemius at comcast.net> To: Ron Stone <ronstone1980 at gmail.com>; arun <smartpink111 at yahoo.com> Cc: r-sig-mixed-models at r-project.org Sent: Wednesday...

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, var.eq=TRUE) coag.mod <- aov(time ~ diet, data=coag) options(scipen=5) options(digits=5) TukeyHSD( coag.mod ) I get the following error when I run the above code Read 72 items Error in rep.int(n, lengt...

On Apr 24, 2008, at 9:03 AM, Gabor Greif wrote: > > As you can see, the use-diet changes actually lower the build time > of kimwitu++! (this is as of yesterday's r50182). > Parity is not only reached, but surpassed. Thanks for these numbers. Do you know how much of this increase is due to co-allocating Use arrays with their users, and how much is due to the actual...

...g performance data about the "class Use"-related changes. I have nice measurements on a 8Gig MacPro with kimwitu++. This is important to say, because this machine is in plenty of memory, so swapping is not likely, which means that in more constrained setups (when swapping occurs) the use-diet approach is probably producing even better results. So here are my values when doing a cd ~/test-suite/MultiSource/Applications/kimwitu++ time make --always in the regular trunk and with my changes merged in: grep "^real" < kimwituRegular.scatter.backup2 | sort real 1m34.002s...

...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 ... link_stat var/lib/apt/diet.lists : No such file or directory done client: nothing to do: perhaps you need to specify some filenames or the --recursive option? It turns out I spelled the input filename wrong. The bug? Telling me to tr...

...g on my data. I was wondering if someone could look at my code and tell me what I'm doing wrong. I was hoping somebody's found a way to do this (I'm sure they have) and can tell me how to fix my code. # Input example data growth <- c(6.6,7.2,6.9,8.3,7.9,9.2,8.3,8.7,8.1,8.5,9.1,9.0) diet <- as.factor(rep(c("A","B","C"),2,each=2)) coat <- as.factor(rep(c("light","dark"),each=6)) growth.means <- aggregate(growth,list(coat,diet),mean) library(plotrix) growth.errs <- aggregate(growth,list(coat,diet),std.error) # Try usin...