similar to: translating sas proc mixed to lme()

Dear all, I would like to conduct a Pearson's correlation using cor.test separately for each individual in my data set. Each individual has nine observations of the measurement variable 'Delta13C'. I can figure our how to do it manually one 'Individual_ID' at a time, but I cannot figure out how to automate the process (a loop I am guessing). In addition, I would like to

Dear list, The below dataset and code creates a new dataset with the results from the function cor.test being performed on each individual ('Individual_ID') from my original dataset. How do I convert each variable from the cor.test results to a numeric data type, as it is passed into the new dataframe? For example, 'estimate', 'p.value', and 'conf.int' should be

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

hi, all I do a hotplug detach a network in a thread, because virDomainDetachDeviceFlags maybe asynchronous, so I do like follow: cb_para->cluster_id = info->cluster_id; cb_para->group_id = info->group_id; cb_para->vsys_id = info->vsysid; cb_para->vnf_id = info->vnf_id; cb_para->conn = conn; cb_para->time_out = 20*1000;//20s

I have a list of observations 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){

Hello R users, Doing My PhD I collected count data which I believe is zero inflated. I have run a statistical model with lmer and family=poisson and got summary(model)@sigma=1 so I believe there is no overdispertion. I would like to use the fmr function from the 'gnlm' library but I just cannot figure out from the examples in the help page and some forums out there how to convert the lmer

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

Dear all, I could not draw a forest plot for meta-analysis under ramdon models using the rmeta package. The rmeta has a default function for MH (fixed-effect) model. Has the rmeta package been updated for such a function? Or someone revised it and kept a private code? I would appreciate it if you could provide some information on this question. Thanks, Andrew This email is intended

Dear R experts, I have a labeled data set. Each data is assigned a binary label 0 or 1. Assume that I use some clustering algorithm to group the data by clusters (using some features of the data). Now I want to know how many data are labeled as 0/1 in each cluster. For example, assume that I have 9 labeled data grouped into three clusters. The ids of the clusters are 1, 2, and 3. The dataset is

Dear R colleagues, for a statistics tutorial I want to develop a nice 3d-graphic of the well known target comparison/analogy of accuracy and precision (see i.e. http://en.wikipedia.org/wiki/Accuracy_and_precision for a simple hand made 2d graphic). The code for a really beautiful graphic is already provided as demo(bivar) in the rgl package (for a picture see i.e

Hi everybody, I have tested CTDB+GFS2+CMAN under Debian. It works good but I do not understand some points. It is possible to run the CTDB defining it under services section in cluster.conf but running it on the second node shuts down the process at the first one. My CTDB configuration implies 2 active-active nodes. Does CTDB care if the node starts with clean_start="0" or

Hello friends of R list, Im a physician and Im not that good in statistics. I have posted similar email in the epi-sig list but no one aswered so far. Im cunducting a systematic review on two diagnostic test for a particular tropical disease. I found a free software that make most of the analysis called MetaDiSc. However there is a paticular analysis that I wuould like to try that it is not

Hi, I am working on calculating X^2 for some matrix (most of them have either two rows or 2 columns) by using chisq.test in R. However when there are 0s in the matrix, chisq.test does not work. For example: ==================== > elements <- matrix( c( 0, 0, 9, 5, 71, 168), nr = 2 ) > elements [,1] [,2] [,3] [1,] 0 9 71 [2,] 0 5 168 > chisq.test( elements )

==================================================================== "Harmony makes small things grow, lack of it makes great things decay." Sallust ==================================================================== Release Announcements --------------------- This is is the latest stable release of Samba 4.0. Major enhancements in Samba 4.0.3 include: o

==================================================================== "Harmony makes small things grow, lack of it makes great things decay." Sallust ==================================================================== Release Announcements --------------------- This is is the latest stable release of Samba 4.0. Major enhancements in Samba 4.0.3 include: o

Hi All, I'm trying to run a simulation of host-pathogen evolution based around individuals. What I need to have is a dataframe or table of some description - describing all the individuals of a pathogen population (so far I've implemented this as a matrix): ID No_of_Effectors Effectors (Sequences) [1,] 0001 3 ## 3

Let it slide. Brian Brunner brian.t.brunner at gai-tronics.com (610)796-5838 >>> craigwhite at azapple.com 11/07/05 06:39PM >>> On Mon, 2005-11-07 at 15:03 -0800, Bryan J. Smith wrote: > No sense in responding further. I think I've made my point. ******************************************************************* This email and any files transmitted with it are

Hi, I have a genotype data for both case and controls and would like to calculate the HW p-value. However, since the number of one genotype is 0, I got wired result. Would someone help me to figure it out? Or confirm it's right? Thanks a lot. ============ > library( "genetics" ) NOTE: THIS PACKAGE IS NOW OBSOLETE. The R-Genetics project has developed an set of enhanced

I have 10000 bivarate samples (x1, x2), and I want to estimate the marginal density of x2. I searched the R manual but couldn't find a function that can do this job. It seems "density" only works for single-variate samples. Can anybody help me with it? Thanks a lot! Best, J. Deng

On Tuesday, March 07, 2000 5:40 PM, Richard Bilonick wrote: >I need to merge several data.frames into one data.frame. In S-Plus I would >use >"merge" but I don't see a merge command in R. What is the best way to >accomplish >this? The easiest way to to this, I think, is as follows: if you have several data frames, dataframe1, dataframe2, . . . , dataframen, you