similar to: Plotting coxph survival curves

Here is a clean version. I did this with nls and it works (see below), but I need to do it with optim. Keun-Hyung # optim vol<-rep(c(0.03, 0.5, 2, 4, 8, 16, 32), 3) time<-rep(c(2,4,8),each=7) p.mated<-c(0.47, 0.48, 0.43, 0.43, 0.26, 0.23, NA, 0.68, 0.62, 0.64, 0.58, 0.53, 0.47, 0.24, 0.8, 0.79, 0.71, 0.56, 0.74, 0.8, 0.47) eury<-data.frame(vol=vol, time=time, p.mated=p.mated)

I saw a similar question but I still don't fully understand how to implement optim. Can someone help me out with this? Thanks. Keun-Hyung > vol<-rep(c(0.03, 0.5, 2, 4, 8, 16, 32), 3) > time<-rep(c(2,4,8),each=7) > p.mated<-c(0.47, 0.48, 0.43, 0.43, 0.26, 0.23, "null", 0.68, 0.62, 0.64, 0.58, 0.53, 0.47, + 0.24, 0.8, 0.79, 0.71, 0.56, 0.74, 0.8, 0.47) >

Given a data frame consisting of a pointID and 12 pairs of (lat, long) variables, with names latA, longA, latB, longB, ... latL, longL, I want to reshape it into a data frame with the structure point source lat long 1 A ... ... 1 B ... ... I've looked at reshape and plyr, but can't figure out how to do this. Details of my data frame are below

I want to implememt auto mated email generation for every transaction. And then we will sent any uplods through our own emails. i wnat grep those email with attachments. This is very urgent for me. if any one know this please let me know. Thankyou for advance. Actually in paperlesspipeline they are did same functionslity. if you want check this url http://paperlesspipeline.com/ I

Hello R list, I am hoping to conduct a logistic regression with repeated measures, and would love an actual "code run through" for such an analysis. I found only one related post on this list, but a full answer was never provided. I understand that the routine lmer (or lmer2) in the lme4 package is often recommended in such a case, but actually implementing it is where I've hit a

Hello, I am having some trouble coding a two-way anova due to replicated treatments. I have a factorial design with three male parents and three female parents. They were mated in all combinations and their babies were grown out and measured for size. 50 babies were measured for each of the 9 crosses. If I stopped here, I would have no troubles. But I also have 5 replicate tanks for each of

Hello, My advisor and I have been working on some survivorship analyses in R and we are hoping to get some feedback on a particular issue involving nesting. We are interested in patterns of food discovery by ant species. Our observations consist of time to discovery by an ant for three different food types, each of two different sizes. These data were collected at 6 plots located in each of

This was accidentally addressed to me, personally: ------- start of forwarded message ------- Message-id: <39F70D7C.287323D5 at yale.edu> MIME-version: 1.0 From: Diana Fischer <diana.fischer at yale.edu> To: maechler at stat.math.ethz.ch Subject: R-Project Date: Wed, 25 Oct 2000 11:42:36 -0500 I briefly looked at the R documentation of available statistical procedures. It appears

Hi Everyone, I posted a similar question about a week ago, but haven't gotten any replies -- I'm afraid that's because my previous question was too vague. Let me try again with a more specific question, and I hope someone can help. NOTE, I know I should be using the newer lme4 package, I just haven't had a chance to update my version of R yet, so the question below relates

Dear list, I have data on insect survival in different cages; these have the following structure: deathtime status id cage S F G L S 1.5 1 1 C1 8 2 1 1 1 1.5 1 2 C1 8 2 1 1 1 11.5 1 3 C1 8 2 1 1 1 11.5 1 4 C1 8 2 1 1 1 There are 81 cages and

Hi all, I run samba 3.6.6 from debian wheezy (version 3.6.6) and i am experiencing some troubles with file name mangling. If i try "dir /x" on a mapped folder it gives me unexpected mangled names: the name mangling matches only the first character and not the first 5 as i expect. For example: if the long file name is LONGFILENAME.TXT, i expect the mangled sholud be something like

All, I had a question about the JGR console and whether or not I can manipulate the location where line wrapping occurs. I have searched 'JGR' in the R listserve archives and attempted to find console manipulation on the JGR website to no avail and could use some direction. TIA, Bret As an example, the below output wraps every 4th value, leaving about 2/3 of the console empty. >

Hi R community, I am new bird to R and moved recently from SAS. I am no means expert on either but very curious learner. So your help crucial for me to learn R. I have already got positive expression. I was trying to fit a mixed model in animal experiment but stuck at simple point. The following similar example is from SAS mixed model pp 212. # data genetic_evaluation <-

Dear friends, I have asked last few days about cross-random effects using PQL, but I have not receive any answer because might my question was not clear. My question was about analysing the salamander mating data using PQL. This data contain cross-random effects for (male) and for (female). By opining MASS and lme library. I wrote this code sala.glmm <- glmmPQL(fixed=y~WSf*WSM,

Dear List, Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... Thank you in advance for any help. Best, Brian [[alternative HTML version deleted]]

Hello, I am trying to compute the mdeian of the survival time from the function survfit: > fit <- survfit(Surv(time, status) ~ 1) > fit Call: survfit(formula = Surv(time, status) ~ 1) records n.max n.start events median 0.95LCL 0.95UCL 111 111 111 20 NA NA NA The results is NA? the fit$surv gives values between 1 and 0.749! Am I doing this correct?

Hi All, I am trying to do a survivorship analysis with library(survival)from a data set that looks like this: I followed a bunch of naturally germinated seedlings of an annual plant from germination to death (none made it to reproduce, and died in a period of ~60 days after germination.) I also know the size of the seed of every individual censused. So I am trying to analyze seedling survival as

All, I am struggling to create a survival plot using LTRC data for each year of a 10 year period. I have a set of individuals (birds) where 'entry' is the day of the year (1-365) they are released (let out of pens) into the wild (2 year data snip below). 'Entry' (e.g., day of year the first bird is released for each year) is highly variable, ranging from 48 to >250. When I

Hello all, I have two questions regarding nested anova. I've attached a simple example, where a male is mated with 4 females, and the intensity of the colour of the eyes is measured in two offspring from each female. 1) As I understand it, residuals in a nested anova are the difference between an observed value and the predicted value, which is equal to the sample mean for the cell that

I have a pcwin3.11. I was mapping solaris 2.5 with samba2.0.5 fine. I added the solaris2.7 system with samba2.0.7 to this small lan. Actaully this solaris 2.7 box was another 2.5 box that I upgraded it. I ftped the smb.conf file from the solaris2.5 box to the solaris 2.7 box so that this way the two solaris systems match and that pcwin3.11 will work with both equally well. I even mated the