similar to: Combining Overlapping Data

dfhfsdhf at ghghgr.com I ran a simple lme model: modelrandom=lmer(y~ (1|Test) + (1|strain), data=tempsub) Extracted the BLUPs: blups=ranef(modelrandom)[1] Even wrote myself a nice .csv file....: write.csv(ranef(modelrandom)[1],paste(x,"BLUPs.CSV")) This all works great. I end up with a .csv file with the names of my strains in the first column and the BLUP in the second

So i have a simple function: bmass=function(y){ weight=y$WT*y$MSTR return(bio) } And want to apply to a whole bunch of rows in my data.frame: final1=apply(final,1,yldbu) BUT...recieve the following error: "Error in y$WT : $ operator is invalid for atomic vectors" However when i try: > final[1,]$WT*final[1,]$MSTR [1] 156.3 It gives me the correct answer....what is apply not

Dear R user, I have a question on using R to analyze data with repeated measurements. I have 2 species with several strains (12) per species, each of which has been measured twice with for a given trait. No particular covariance, just two measures. Now I want to analyze the data with an ANOVA (aov) considering these repeated measures to get the MSq and SSq for the species and strain level. I

Dear list, Has anyone created a version of SSmicmen that allows testing for group differences? The basic Michaelis-Menten equation is: (Bmax * X) / (Kd + X). The nlme package allows modeling of random effects for Bmax and Kd as needed, but I curious how I can build in group differences? I have receptor binding data for strains of mice, and following Pinheiro and Bates' lead in their

Dear List, I downloaded R for the first time yesterday, in the hopes that I might deal more effectively with a complex repeated measures experimental design involving inbred strains of laboratory mice. The design below, somewhat simplified, cannot be computed with standard ANOVA, because something called the X'X matrix is too large. The design has the following factors: Between-subject

I have assayed the concentrations of various metal elements in different anatomic regions of two strains of mice. Now, for each element, in each region, I want to do a t test to find whether there is any difference between the two strains. Here is what I did (using simulated data as an example): # create the data frame > elemconc = data.frame(expand.grid(id=1:3, geno=c('exp',

Having problems with the write.table function. I can write a tab delimited file just fine, but for each line in my matrix its inputs a carriage return when i dont want it to. For example my matrix might be: ID V1 V2 V3 FARY1004 1 2 3 FARY2067 2 3 1 FARY4587 2 2 2 And I want the written File to be: FARY1004 1 2 3FARY2067 2 3 1FARY4587 2 2 2 TIA -- View this

Hello. We have a webserver using FreeBSD, we read about tunning kern.ipc.somaxconn (http://www.freebsd.org/doc/en_US.ISO8859-1/books/handbook/configtuning-kernel-limits.html) so the OS can handle all the connections. Is there a way to know how many connections are established in a certain moment?. I know about netstat(1) but is there any other command that we can use to know the exact amount of

Hi, I have a data-frame, r (column names below), that needs subsetting: date, time, strain, gene, deltact When I try to subset r by applying selection criteria on two columns I get an empty data frame. For example I would like to extract all rows that have time == 0h and strain == ROC. So, t <- subset(r, (r$time == "0h" && r$strain == "ROC"), select= c(time,

Hello, I am new to R and currently have the following problem: I have successfully loaded my data in R which consists of two numeric columns (LI_F and female) and one character column (Strain). So far I can plot two different set of boxplots for each of the numeric columns plotted by the groups of the character column and the commands look like that: boxplot(LI_F~Strain, ylab="LI_F",

Hello Sir I am getting problem in plotting in CCA . Could you please help me? I wrote the below command but I don't know why it is taking only first 5 env data rather than all 9. > strain.data <- read.xlsx("Dee rhiz.xlsx", sheetName="strain", header = T, row.names = 1) > env.data <- read.xlsx("Dee rhiz.xlsx", sheetName="env", header = T,

In the vcdExtra package on R-Forge, I have functions and generic methods for calculating log odds ratios for R x C x strata tables. I'd like to define methods for fitting weighted lm()s to the resulting loddsratio objects, but I'm having problems figuring out how to do this generally. # install.packages("vcdExtra", repos="http://R-Forge.R-Project.org")

Hello, I seem unable to construct a legend which contains a substitution as well as math symbols. I'm trying to do the following: strain2 <- "YJG48" legend.txt <- c( substitute( strain * %==% * "YJG45, rpn10" * %Delta%, list(strain=strain2) ), "Verhulst/Logistic", "Malthus" ) legend( 100,2.5, legend.txt, cex=0.75,

Hi, I have some difficulty in figuring out whether I am doing correct or not. A brief introduction about the work: It is a light/dark choice test conducted in insect larvae.  The response is binary (0- present in dark area, 1-present in light area) and the experiment is run for 15 min, so there are 15 repeated measurements per individual larva at 1 min intervals.  The factors which affect

Dear R list, What is the proper way to specify a nested model so that the F values agree with the expected mean square errors? Specifically, suppose I have a design where "Heads" are nested within "Machines". I would like to model the following Y_ijk = Mu + Machine_i +Head_j(i) +Error_k(ij). Using the commands below, > summary(aov(Strain~Machine + Head%in%Machine ))

I was searching thru the internet and I found a wide variety of different web interfaces for asterisks I was curious which one is best suited for asterisks. Thanks

Hi, I get ssh connect attempts all the time, to my servers at home and at work. I've noticed lately they come from a certain ip address, hitting every 3 or 4 seconds, trying 50 or 100 different user names and passwords. And I get these sweeps from 2 or 3 ip addresses a day. I guess this is an automated attempt to guess a user/pass and break into a system. I tried to secure ssh better by

Hello, as some may remember, we're running very dense IMAP cluster here, in excess of 50k IMAP sessions per node (current record holder is 68k, design is for 200k+). The first issue we ran into was that the dovecot master process (which is single thread and thus a bottleneck) was approaching 100% CPU usage (aka using a full core) when trying to spawn off new IMAP processes. This was

Dear R, I have a question concerning quantile regression models. I am using the quantile regression model to test the relationship between abalone and the percentage cover of algae etc at different sites and depths. An example is fit<-rq(abalone~%coversessileinvertebrates+factor(Depth)+factor(Site),tau=0.7) In this model depth is fixed and site is random. My question is, is it possible

Dear Subscribers, I am using R for quite a while nowadays on Ubuntu 10.04 LTS. I a using R for doing my statistics. Furthermore I am using it as a tool to generate the graphics for my publications. I am currently working on a project which involves nls-fits of three dimensional surfaces such as ellipsoids or even more complex. I have been searching R help and the manuals for a possibility to