similar to: Simulation - Natrual Selection

-------- Original Message -------- Subject: Re: [R] Simulation - Natrual Selection Date: Wed, 05 Jan 2011 17:24:05 +0000 From: Ben Ward <benjamin.ward@bathspa.org> To: Bert Gunter <gunter.berton@gene.com> CC: Mike Marchywka <marchywka@hotmail.com> On 05/01/2011 17:08, Bert Gunter wrote: > Couple of brief comments inline below. -- Bert > > On Wed, Jan 5, 2011 at

Dear all, I'm trying reproduce an example of bioassays Yielding Concentration-Mortality Data particularly control - adjustment model from book Bioassay of Entomopathogenic Microbes and Nematodes chapter 7 with R. I used glm with family=binomial and link=probit, but I do not know how to implement parameter gamma (control mortality - mortality of the untreated control insect in this exaple)

> Can you file a bug with some example outputs that are unrelated to the search string? Here is the example (see attachment). This example does the following: 1)First, it indexes text from the "text.txt" file (see attachment) (actually, this is the text of the following book: "Abbas, Lichtman. Basic immunology"). 2)Next, it searches for the "extracellular

Hi, I'm doing some modelling (lm) for my 3rd year dissertation and I want to do some resampling, especially as I'm working with microbes, getting them to evolve resistance to antimicrobial compounds, and after each exposure I'm measuring the minimum concentration required to kill them (which I'm expecting to rise over time, or exposures), I have 5 lineages per cleaner, and

Hi, I have some values in a list format generated by the following: Path_Number <- 0010 ID.Path <- formatC(0001:Path_Number, width=4, flag=0) # Make vector of ID's. No_of_Effectors <- sample(1:550, length(ID.Path), replace=TRUE) # Define Number of Effectors each individual gets. Effectors <- split(sample(1:10000, sum(No_of_Effectors), replace=TRUE), rep(ID.Path, No_of_Effectors))

Hello, This is my first time addressing such a big audience so apologies in advance in case I fail to formulate this question. I am working with 13 species of trees, and the data I have are: 1 continuous (phenolic concentration in xylem and in phloem) and 2 categorical variables: lineage (3 subclades) and habitat (fire and non fire). I am trying to see how species can be splitted

Dear all: I am a new R-user and I have 2 questions about it. 1) I have a dataframe. Based on ?formulation? and ?subject?, a dataframe is split into 4 dataframes. The example is as follows. Moreover, I want to calculate ?test? value for these 4 dataframes. My question is that the ?test? values not correct and I do not know where the problem is. 2) There are 12 ?test? (y) values from

Hello R-help I have the following data for a standard curve concentration(nM),fluorescence 0,48.34 2,58.69 5,70.83 10,94.73 20,190.8 50,436.0 100, 957.9   (1)Is there function in R to plot a spline. (2)How can I interpolation,say 1000 point from 0nM-100nM and store this as a data frame of concentration,fluorescence (3)How can I modify the code below so that instead of retrieving a concentration

Hi, I am struggling to create a 2 by 2 multiple graphs in one page. I used par(mfrow=c(2,2)) to divide the screen into 4. In each screen I draw a pie chart (They are all same). For example, my data is like this Concentration value A1 69 A2 8 G1 51 G2 1

Hello-- I have a data for small population who took 1 drug at 3 different doses. I have the actual drug concentrations as well as predicted concentrations by my model. This is what I'm looking for: - Time vs Concentration by ID (individual plots), with each subject occupying 1 plot -- there is to be 9 plots per page (3x3) - Observed drug concentration is made up of points, and predicted drug

Dear R users, today I've got the following problem. Here you are a dataframe as example. There are some SAMPLES for which a CONCentration was recorded through TIME. The time during which the concentration was recorded is not always the same, 10 points for Sample A, 7 points for Sample B and 11 for sample C Also the initial concentration was not the same for the three samples. I would like

Dear all: I am a new R-user and I have 2 questions about it. 1) I need to find specific sub-dataframe, and then use specific column to calculate. For example, after splitting dataframe, I find specific the sub-dataframe, such as ?A.split [1]?. But, I don?t know how to find ?time? and ?concentration? columns of ?A.split [1]?. 2) The equation used to sub-dataframe is

Greetings, I have a question that I'd like to get input on. I have a classic toxicology study where I artificially fertilized and exposed embryos to a chemical and counted defects. In addition, I kept track of male-female pairs that I used to artificially fertilize and generate embryos with. I need to use logistic regression to model the response, but also check that the genetics of the

Greetings, I cannot find solution for this problem (I was searching on web, but without success): I want to plot dose-response models for one concentration and many responses (lets say 200) and I don?t want to do it manually. So I use loop for this: for (i in mydata[,2:201]){ #first column is concentration pdf(paste("plot_",i,".pdf",sep = ""))

Hey Lineage 2 lovers :D, does anybody know why I don't have any texts in Lineage ? Although Lineage works like a charm, I see no text, titles, names, etc. What should be the option ?

Hello, I am using R, through rsruby, to create a graph and best fit line for a set of data points, regarding data collected in a Chemistry class. The problem is that although the graph functions perfectly properly, the best fit line will not work. I initially used code I pretty much copied from a website with a tutorial on this, which was: graphData.png("/code/Beer's-Law

Hi, folks, ###### food=c('fruit','fruit','fruit','drink','drink','drink') type=c('apple','apple','orange','water','soda','soda') value=c(2,3,1,5,7,6) data=data.frame(food,type,value) share=c((2+3)/(2+3+1),5/6,1/6,5/(5+7+6),13/18,13/18)

A data frame has one factor, one date, and one numeric column. When I plot these using the default pch of the open circle (first attachment), xyplot(TDS ~ sampdate | she.s, data = sheep.cast, main = 'TDS in Sheep Creek', ylab = 'Concentration (mg/L)', xlab = 'Time') I see the higher concentration points toward the right on several panels. But, when I change to using a

Hi, I just started a new course this semester on R, I never used it in my life and i'm stuck on these questions from 3 days, it would be really nice if someone could explain me the answers with the relative commands. thanks a lot in advance The following 7 questions are based on the CO2 dataset of R. 1) How many of the plants in CO2 are Mc2 for Plant? 2) How many are either Mc2 or Mn2?

I am sorry to ask such question, but I can't find a solution... I have a dataframe 'd2004' and I want to remove two columns: 'd2004$concentration' and 'd2004$stade". I could do it just as follows: > names(d2004) [1] "Localite" "Date" "parcelle" "maille" "presence.plant"