similar to: problem with glm(family=binomial) when some levels have only 0 proportion values

Dear all! I tried to fit Gompertz growth models to describe cummulative germination rates using nls. I used the following code: germ.model<-nls(percent.germ~a*exp(-b*exp(-k*day)),data=tab,start=list(a=100,b=10,k=0.5)) My problem is that I want that the fitted model goes through the origin, since germination cannot start before the experiment was started, and y-max should be 100. Does anyone

Hi all, I am trying to use nonlinear regression (nls) to analyze some seed germination data, but am having problems with error codes. The data that I have closely matches the germination dataset included in the drc package. Here is the head of the data temp species start end germinated TotSeeds TotGerminated Prop 1 10 wheat 0 1 0 20 0 0.0 2 10 wheat

Hi Keep your messages in the list, you increase your chance to get some answer. I changed your data to groupedData object (see below), but I did not find any problem in it. plot(wlg) gives reasonable picture and I am not such expert to see any problem with data. Seems to me, that something has to be wrong with nlsList function. > wheat.list <- nlsList(Prop ~ SSlogis(end,Asym, xmid,

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: https://stat.ethz.ch/pipermail/r-help/attachments/20061120/73240e63/attachment.pl

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

Does someone out there have an example of R-code for a split-split plot ANOVA using aov or another function? The design is not balanced. I never set up one in R before and it would be nice to see an example before I tackle a very complex design I have to model. Thanks, Mike Mike Saunders Research Assistant Forest Ecosystem Research Program Department of Forest Ecosystem Sciences University of

Hello everybody, I'm using the lattice package and the xyplot to make several graphs like below. However, I can just print the three grouped plots onto one page as I'm putting them into a pdf-file, which gives me a huge amount of pages... Is it possible to put them all, or at least more than one on the same page, for instance put three groups beside each other like columns? ...

Hi Everyone, I am trying to analyze a split plot experiment in the field that was arranged like this: I am trying to measure the fitness consequences of seed size. Factors (X): *Seed size*: a continuous variable, normally distributed. *Water*: Categorical Levels- wet and dry. *Density*: Categorical Levels- high, medium and solo *Plot*: Counts from 1 to 20 The *response variable *(Y) was the

I have a function "eval.delta" which does what I want but isn't very elegant. I have consulted the R documents, MASS, and S Programming. Is there a practical way to optimize the code? Thank you very much. Peter B. Function: eval.delta <- function(delta){ cat("VALUES\n") vlr <- NULL k <- 0 for(j in 1:length(delta)) if(delta[j] <= 2){

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

I have been trying to get the password sync working for Redhat 6.1 (samba 2.05a rpm). I got the PDC working this morning, so I'm happy. The FAQ section on the samba site was very helpful. The password sync was the only sticking point. The default passwd chat script didn't work well passwd chat = *New*UNIX*password* %n\n *ReType*new*UNIX*password* %n\n *passw

Hello, I'm trying to adjust a non linear model in which the biological response variable (ratio of germinated fungus spores) is dependent on 2 covariates (temperature and time). The response to temperature is modeled by a kind of beta function with 2 parameters (optimal and maximum temperatures) and the time function is a 2-parameter Weibull. Adjustments with nls or gnls work, but I need to

https://bugzilla.samba.org/show_bug.cgi?id=8913 Summary: Operation not supported (95) Product: rsync Version: 3.0.9 Platform: x64 OS/Version: Linux Status: NEW Severity: minor Priority: P5 Component: core AssignedTo: wayned at samba.org ReportedBy: martin.suc at googlemail.com

Hi all, I have some serious problems making Xen run. I have spent a couple of days looking in Google and trying a lot of different things, but none of them has worked. I hope you can help me... First of all: xen was running smoothly on my machine, but suddenly (I think I haven''t changed any configuration) it doesn''t work anymore. I am using a Debian 4.0 with Xen 3. When

Hi, I am rather new to R, but both myself and another much more experience user cannot figure this out. I have a collegue who has a 800+ nonlinear regressions to run for seed germination (different species, treatments, etc.) over time. I created a looping structure to extract the parameters from each regression; I will then use the parameters themselves for further analysis. I would like to

https://bugzilla.mindrot.org/show_bug.cgi?id=1280 Damien Miller <djm at mindrot.org> changed: What |Removed |Added ---------------------------------------------------------------------------- CC| |djm at mindrot.org Blocks| |1452 --- Comment #1 from Damien Miller

Hi I have been trying to replace my win 2003 DC with samba4 ADC , but i am getting the following error when i run samba-tool fsmo seize --role=all for all the 5 roles. I found that it is listed in Samba-Bugzilla ? Bug 9461 as on 28/1/2013.Is there any patch been released after that or is there any patch which is planned for release.Any workaround to overcome this [root at server1# samba-tool

Hi R-people, I performed controlled experiments to evaluated the seeds germination of two palms under four levels of water treatments. I conducted a generalized linear model (GLM) with a Poisson distribution to verify whether there were significant differences in the number of seed germination (NS-count variable) between treatments and species (explanatory variables). Thus, my model and output

Hi, Two days ago, I tried to turned my samba server, which worked perfectly into a PDC but despite of all my efforts, skimming thru the docs, mailing lists, hanging on #samba, I've never been able to logon to the new domain from the WinXP clients. Authentication is OK, but then access to the share is denied. Here's part of the log: [2006/09/28 19:03:36, 2]

Hi there,? I've been battling with an extension of this in my own data: getting appropriate error bars once data is stacked in a bar graph.? (original question: http://r.789695.n4.nabble.com/ggplot2-se-variable-in-geom-errorbar-s-limits- td3311176.html). It wouldn't let me reply to that thread. A modification of the earlier answer:? data(diamonds)? ?diamonds_df <- ddply(diamonds,