search for: indivs

Dear all, I'm having a few problems trying to reshape a data frame. I tried with reshape{stats} and melt{reshape} but I was missing something. Any help is very welcome. Please find details below: ################################# # data in its original shape: indiv <- rep(c("A","B"),c(10,10)) level.1 <- rpois(20, lambda=3) covar.1 <- rlnorm(20, 3, 1) level.2

Dear all, I have the following dataset with results from an experiment with individual bats that performed two tasks related to prey capture under different conditions: X variables: indiv - 5 individual bats used in the experiment; all of which performed both tasks task - 2 tasks that each individual bat had to perform dist - 5 repeated measures of individual bats at 5 different distances from

Dear R-Gurus I am a PhD student from South Africa working on chimpanzee behaviour. I am looking at patterns of shade utilization and am using generalized linear mixed models to examine the effects of various factors on whether chimpanzees choose to spend time in the sun or shade. I realise that the lme4 package and the outputs of the lmer functions have been discussed ad nauseum but I have been

Dear r-helpers, why do I get an output in the first iteration of the for-loop which contains the string values of the input vector, and how can I avoid that? Here's the output (only line 1 is wrong) latentVariable Indiv Group 1 rPlanning rIterat rTDD 2 rPlanning 0.79 0.84 3 rIterat 0.79 0.83 4 rTDD 0.9 0.96 5 rStandup 0.83 0.82 6

Dear List I am trying to fit a multinomial model using the mlogit package. Attempting to load the data into mlogit presents the following error. MLOG<-mlogit.data(Mult3,shape="long",choice="CHOICE",alt.var="mode.ids",indivs = "set3",chid.var = "obs") Error in `row.names<-.data.frame`(`*tmp*`, value = c("1.1", "1.2", "1.3", : duplicate 'row.names' are not allowed In addition: Warning message: non-unique values when setting 'row.names': '1.1'...

Hello world, I'm trying to do an anova on data in data.set, dependent variable is a column named "dep.var", grouping variable is in a column called "indep.var", and is.factor(indep.var) is TRUE... why can't I just do aov(dep.var ~ indep.var, data = data.set)? What have I done to deserve this?! What gives? Am I missing something totlly obvious? R-base-1.0.0-1,

...vations of individuals. I would like to make a list of individuals, with a data frame of observations for each individual. The following code usually works, but not always ---------------------------------------------------------------------- # Make a list of empty data frames animals = list() indivs = levels(Z$individual_id) donotprint <- sapply(indivs, function(i){ animals[[i]] = data.frame() }) # Add observations of each animal to the appropriate frame donotprint <- apply(Z, 1, function(r){ ind = r[["individual_id"]] bind = ind # Use different names to confirm that...

I'm trying to create a graph where different cells of a grid (a shapefile) will be painted with a color share scale, where the most easy way is to use gray(). Can I somehow get a vector (gradient) of colors, a vector of colors with other methods but gray()? I'm doing this until now quad_N_sp <-

Sorry to barge in on an old thread. I''m having the same trouble as the gent who started this thread. I''ve tried the options described and can''t seem to get the tunnel to pass packets through it. I''m using the Netscreen Remote VPN client (Safenet derivative) on a windows machine, trying to connect to a Netscreen 5xp at the other end. The connection fires

I am trying to run a GLMM on some binomial data. My fixed factors include 2 dichotomous variables, day, and distance. When I run the model: modelA<-glmmPQL(Leaving~Trial*Day*Dist,random=~1|Indiv,family="binomial") I get the error: iteration 1 Error in MEEM(object, conLin, control$niterEM) : Singularity in backsolve at level 0, block 1 >From looking at previous help

Hello all, The following line of code includes a right-to-left language text, yet the R graphics engine displays it from left to right. One problem this causes is when there are parenthesis in the test, here is a basic example? plot(1:10, main = "שלום (טקסט)") Is there a way to make sure the text is displayed from right to left? Many thanks for any suggestions, Tal

Estimados apañeRos: La duda o propuesta que os voy a plantear es a la vez metodológica y relacionada con R. Me encuentro trabajando con tres variables que son el resultado de un computo de porcentajes. Me explico, se toma una muestra de n casos (unos 6.500 aprox) pertenecientes a i individuos (unos 230 aprox) en la que se comprueba si un determinado evento ha ocurrido o no, anotándose 1 en

Is there a way to create a 'bubble plot' in R? For example, if we define the following data frame containing the level of y observed for 5 patients at three time points: time<-c(rep('time 1',5),rep('time 2',5),rep('time 3',5))

Dear R helpers, I would like to write a loop that makes 4 objects (called A, B, C, and D) each of which contains ten random numbers. This attempt: individuals<-c("A","B","C","D") for(i in 1:length(individuals)) { individuals[i]<-rnorm(10) } does not work because "individuals[i]" is not the proper way to extract each letter from the

Dear all, I would like to draw a 3-D horizontal cylinder preferably in RGL device (because this gives the look from different angles). Basic idea is from http://www.tau.ac.il/cc/pages/docs/sas8/insight/chap18/sect3.htm. Below is the description exactly what I want to do. Please see at figure 18.6, 1st plot. Here it is an confidence ellipsoid. Suppose now you put another 4 same ellipsoids on top

Dear all, i'm trying to label specific individuals (supplementary ones) after a PCA with the FactoMiner package. There is not much details (possibilities?) in the R-help of the plot.pca function. There is indeed a "label" parameter but i could only manage to label the supplementary individuals with there "row.names" (i.e. label="indiv.sup") and not with the

I do the these passages: library(Geneland) set.seed(1) data <- simdata(nindiv=200, coord.lim=c(0,1,0,1) , number.nuclei=5 , allele.numbers=rep(10,20), IBD=FALSE, npop=2, give.tess.grid=FALSE) geno <- data$genotypes coord <- t(data$coord.indiv) path.mcmc <-

Hola Isidro: También había sopesado esa posibilidad pero por una parte me parecía "complicar" el proceso y por otra tengo mis dudas acerca de que en que el análisis cluster pueda prescindir de la hipótesis de Normalidad en las dos variables "raras" que tengo. De cualquier forma muchas gracias por tu idea. Un saludo {In Archive} RE: [R-es] Clasificacion de

En cualquier caso, para nada necesitas normalidad en las variables. De las distribuciones que presentas, y abundando en lo que te indica Olivier, ¿sería interesante considerar sólo 2 grupos (?normal? y ?enhorabuena?) en las variables 2 y 3? Todo esto sin saber de lo que estamos tratando, claro? XD Un saludo, Isidro De: Olivier Nuñez [mailto:onunez en unex.es] Enviado el: martes, 30 de

> > Message: 5 > Date: Sat, 28 Apr 2007 11:37:51 -0700 > From: "J.H." <warthog19 at eaglescrag.net> > Subject: Re: [syslinux] USB stick meets ISO system > To: Geert Stappers <stappers at stappers.nl> > Cc: syslinux at zytor.com > Message-ID: <1177785471.25140.29.camel at localhost.localdomain> > Content-Type: text/plain > > On Sat,