Thank you very much, I'll look into that package! Kirsten _____________________________________________________________ Kirsten Martin Von Holle Lab Department of Biology, BL 405 University of Central Florida Orlando, FL 32816 ________________________________ From: Jagat.K.Sheth [via R] [ml-node+s789695n4651207...@n4.nabble.com] Sent: Wednesday, November 28, 2012 4:07 PM To: Kirsten Martin Subject: Re: Conditional model in R
Sounds like a finite mixture model. I haven't read your references but an overall model for such an approach could be f(Y=0; pi, kappa) = 1- pi + pi*f(Y=0|Z=1; kappa) where pi=Pr(Z=1) is the probability of an event, z, and y is the value observed when the event occurs and f is the probability density of Y with parameters kappa. You could try 'fmr' in Jim Lindsey's gnlm package (available at http://www.commanster.eu/rcode.html ) which fits generalized nonlinear regression models with two or three point mixtures using maximum likelihood. > -----Original Message----- > From: [hidden email]<UrlBlockedError.aspx> [mailto:r-help-bounces@r- > project.org] On Behalf Of Kirsten Martin > Sent: Wednesday, November 28, 2012 1:33 PM > To: [hidden email]<UrlBlockedError.aspx> > Subject: [R] Conditional model in R > > Hello all, > > I have a data set where the response variable is the percent cover of a > specific plant (represented in cover classes 0,1,2,3,4,5, or 6). This > data > set has a lot of zeros (plots where the plant was not present). > I am trying to model cover class of the plant as a function of both > total > nitrogen and shrub cover. > > After quite a bit of research I have come across a conditional approach > to > modeling data with a lot of zeros (Fletcher et al. 2005, Welsh et al. > 1996). > In this approach you model the presence/absence data using a logistic > regression and then model the presence only data using ordinary (least > squares) regression. > > I have successfully come up with both a logistic model and an ols model > with > good fits. I am running into trouble combining the two (as outlined in > the > third step of the Fletcher et al. 2005 paper). > > Does anyone have any experience or any advice on doing this? How does > one > come up with an overall model for the data using this approach? > > Thanks for your help! > Kirsten > > > > -- > View this message in context: > http://r.789695.n4.nabble.com/Conditional-model-in-R-tp4651188.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > [hidden email]<UrlBlockedError.aspx> mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting- > guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ [hidden email]<UrlBlockedError.aspx> mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ________________________________ If you reply to this email, your message will be added to the discussion below: http://r.789695.n4.nabble.com/Conditional-model-in-R-tp4651188p4651207.html To unsubscribe from Conditional model in R, click here<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=4651188&code=a21tYXJ0aW5Aa25pZ2h0cy51Y2YuZWR1fDQ2NTExODh8LTgwNzYyNDM3Nw==>. NAML<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml> -- View this message in context: http://r.789695.n4.nabble.com/Conditional-model-in-R-tp4651188p4651288.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
