Fernando, For background reading you may check
Stevens. M.H.H. (2009) A Primer of Ecology with R. Use R! Series. Springer. and associated 'primer' package in CRAN. The book discusses functional responses. Cheers, Jari Oksanen ________________________________________ From: r-sig-ecology-boun...@r-project.org [r-sig-ecology-boun...@r-project.org] on behalf of Pierre THIRIET [pierre.d.thir...@gmail.com] Sent: 30 May 2013 09:58 To: r-sig-ecology@r-project.org Subject: Re: [R-sig-eco] Functional response type II question, curve fitting and Dear Fernando, I would suggest the function nls() that fit non-linear model. You can specify by hand your equation and starting values of your parameters (your guessed value) R include some Self-Starting functions for the model the most used, this could be easier for you to work with it. check ?selfStart() According to wikipedia, Type II functional response is similar to the Michaelis--Menten equation, which is implemented in R as self starting function:http://stat.ethz.ch/R-manual/R-patched/library/stats/html/SSmicmen.html By adpting examples, you could easily fit your model, if the choosen equation is suitable. About statistical comparisons of estimated parameters among curves, I have no idea but I found, among other this post: http://stats.stackexchange.com/questions/26611/how-to-test-the-effect-of-a-grouping-variable-with-a-non-linear-model Good luck, Pierre Le 30/05/2013 07:49, Luis Fernando García Hernández a écrit : > Dear all, > > I am new on the list and on the more complex applications of R, so I ask > you to excuse me if my mail is too long. > > This time I have to do several tasks related tu functional reponse. I want > to evaluate the relationship between prey density (x) vs prey consumed(y) > the idea is to fit the data set to the functional reponse equation: > Na/TP=a/(1+aThN). I have idea about the value of most of the equation > parameters, but I am completely lost about how to fit the data points to > this equation. > > On another hand I have to graph and compare if there are significative > differences on estimated parameters of the curves for lepidopterans and > ants, the final data should look like the attacher picture.If any of you > have some idea about how to do this, would be really appreciated. > > Below you can find the data and the pictures. > > Thanks in advance > > Prey Density Eaten Lepidopteran 1 2 Lepidopteran 1 3 Lepidopteran 1 3 > Lepidopteran 3 4 Lepidopteran 3 5 Lepidopteran 3 4 Lepidopteran 5 7 > Lepidopteran 5 5 Lepidopteran 5 6 Lepidopteran 10 8 Lepidopteran 10 10 > Lepidopteran 10 7 Ant 1 3 Ant 1 1 Ant 3 4 Ant 3 2 Ant 5 4 Ant 5 6 Ant > 10 5 Ant 10 6 > > > _______________________________________________ > R-sig-ecology mailing list > R-sig-ecology@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-ecology [[alternative HTML version deleted]] _______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology