It need to apply nonlinear modellingto physiological data collected on three NZ tree fern species (cyadea =Cyatheadealbata, cyasmi = Cyathea smithii, dicsqu= Dicksoniasquarrosa) . The data was collected with a fluorometerinstrument which was used to measure the electron transport rate in photosystemII (etr) (a proxy for photosynthesis) as a function of light (par).Please find the data in theprovided
etr_lc.csvdata set. There were three individuals of each species (ind) and those are true statistical replicates, so no need to aggregate the data in this case. You will see that the etrdeclines at high light levels indicating photoinhibition. Use the attached self-starter function of theEilers& Peeters(1988) photosynthesis model, which accounts for this photoinhibitory effect. Run a varying coefficient model allowing individual parameter estimates for each tree fern species. Check the model assumptions and apply variance modelling techniques if required. Then formulate a model with common parameter estimates for all species using the same variance structure as in the varying coefficient model and compare the two models using a likelihood ratio test. If the varying coefficient model is supported by the likelihood ratio test, follow up with a
post-hoccomparison using the attached custom pwcfunction (you also find the code for this function in the nonlinear regression course PDF).
"spec","ind","par","etr" "cyadea","1",0,0 "cyadea","1",33,10.20096 "cyadea","1",47,14.2128 "cyadea","1",67,19.52916 "cyadea","1",113,28.7133 "cyadea","1",171,36.34092 "cyadea","1",274,44.07564 "cyadea","1",392,45.60528 "cyadea","1",587,44.13066 "cyadea","1",926,37.33632 "cyadea","1",1410,23.0958 "cyadea","2",0,0 "cyadea","2",33,10.56132 "cyadea","2",47,14.94318 "cyadea","2",67,20.9643 "cyadea","2",113,32.5101 "cyadea","2",171,44.38476 "cyadea","2",274,57.77016 "cyadea","2",392,63.88032 "cyadea","2",587,61.635 "cyadea","2",926,53.28204 "cyadea","2",1410,NA "cyadea","3",0,0 "cyadea","3",33,10.96326 "cyadea","3",47,15.1998 "cyadea","3",67,20.45778 "cyadea","3",113,29.4252 "cyadea","3",171,36.48456 "cyadea","3",274,41.31372 "cyadea","3",392,40.00752 "cyadea","3",587,36.981 "cyadea","3",926,33.0582 "cyadea","3",1410,NA "cyadea","4",0,0 "cyadea","4",33,9.9792 "cyadea","4",47,13.83774 "cyadea","4",67,18.82566 "cyadea","4",113,27.33696 "cyadea","4",171,35.1918 "cyadea","4",274,41.77404 "cyadea","4",392,41.16 "cyadea","4",587,39.93948 "cyadea","4",926,38.11416 "cyadea","4",1410,23.688 "cyadea","5",0,0 "cyadea","5",33,9.95148 "cyadea","5",47,13.34424 "cyadea","5",67,16.85586 "cyadea","5",113,21.68922 "cyadea","5",171,25.42428 "cyadea","5",274,28.07952 "cyadea","5",392,27.65952 "cyadea","5",587,26.13324 "cyadea","5",926,28.00224 "cyadea","5",1410,21.3192 "cyadea","6",0,0 "cyadea","6",33,9.68814 "cyadea","6",47,13.3245 "cyadea","6",67,17.5875 "cyadea","6",113,23.91984 "cyadea","6",171,29.15892 "cyadea","6",274,34.06368 "cyadea","6",392,34.08048 "cyadea","6",587,34.02252 "cyadea","6",926,32.28036 "cyadea","6",1410,NA "cyasmi","1",0,0 "cyasmi","1",33,10.69992 "cyasmi","1",47,15.01227 "cyasmi","1",67,20.80953 "cyasmi","1",113,32.05923 "cyasmi","1",171,42.80472 "cyasmi","1",274,54.26022 "cyasmi","1",392,55.89528 "cyasmi","1",587,53.12937 "cyasmi","1",926,43.55904 "cyasmi","1",1410,30.2022 "cyasmi","2",0,0 "cyasmi","2",33,10.51974 "cyasmi","2",47,14.78526 "cyasmi","2",67,20.5422 "cyasmi","2",113,31.65582 "cyasmi","2",171,43.52292 "cyasmi","2",274,56.73444 "cyasmi","2",392,63.71568 "cyasmi","2",587,59.41614 "cyasmi","2",926,50.5596 "cyasmi","2",1410,33.1632 "cyasmi","3",0,0 "cyasmi","3",33,10.49202 "cyasmi","3",47,14.72604 "cyasmi","3",67,20.34522 "cyasmi","3",113,31.3236 "cyasmi","3",171,41.44014 "cyasmi","3",274,49.94472 "cyasmi","3",392,51.20304 "cyasmi","3",587,47.82876 "cyasmi","3",926,38.892 "cyasmi","3",1410,37.9008 "cyasmi","4",0,0 "cyasmi","4",33,10.68606 "cyasmi","4",47,14.86422 "cyasmi","4",67,20.2608 "cyasmi","4",113,29.56758 "cyasmi","4",171,37.92096 "cyasmi","4",274,45.4566 "cyasmi","4",392,45.76992 "cyasmi","4",587,44.3772 "cyasmi","4",926,38.11416 "cyasmi","4",1410,27.2412 "cyasmi","5",0,0 "cyasmi","5",33,10.86624 "cyasmi","5",47,15.08136 "cyasmi","5",67,20.48592 "cyasmi","5",113,30.1371 "cyasmi","5",171,39.1419 "cyasmi","5",274,48.21852 "cyasmi","5",392,50.87376 "cyasmi","5",587,48.32184 "cyasmi","5",926,42.39228 "cyasmi","5",1410,NA "dicsqu","1",0,0 "dicsqu","1",33,14.60592 "dicsqu","1",47,20.87316 "dicsqu","1",67,28.12992 "dicsqu","1",113,41.83368 "dicsqu","1",171,54.06912 "dicsqu","1",274,65.142 "dicsqu","1",392,70.32102 "dicsqu","1",587,61.0134 "dicsqu","1",926,57.33756 "dicsqu","1",1410,36.00702 "dicsqu","2",0,0 "dicsqu","2",33,14.72184 "dicsqu","2",47,20.90088 "dicsqu","2",67,28.12992 "dicsqu","2",113,42.02016 "dicsqu","2",171,53.16192 "dicsqu","2",274,64.1718 "dicsqu","2",392,66.40284 "dicsqu","2",587,58.8672 "dicsqu","2",926,49.81788 "dicsqu","2",1410,NA "dicsqu","3",0,0 "dicsqu","3",33,14.64456 "dicsqu","3",47,20.76228 "dicsqu","3",67,27.55662 "dicsqu","3",113,39.59592 "dicsqu","3",171,48.9888 "dicsqu","3",274,54.8856 "dicsqu","3",392,53.20476 "dicsqu","3",587,49.3626 "dicsqu","3",926,39.00834 "dicsqu","3",1410,NA "dicsqu","4",0,0 "dicsqu","4",33,14.81844 "dicsqu","4",47,20.81772 "dicsqu","4",67,27.67128 "dicsqu","4",113,39.84456 "dicsqu","4",171,49.98672 "dicsqu","4",274,56.133 "dicsqu","4",392,60.62868 "dicsqu","4",587,55.8012 "dicsqu","4",926,45.11808 "dicsqu","4",1410,39.53712 "dicsqu","5",0,0 "dicsqu","5",33,14.95368 "dicsqu","5",47,21.09492 "dicsqu","5",67,28.24458 "dicsqu","5",113,41.0256 "dicsqu","5",171,50.98464 "dicsqu","5",274,58.4892 "dicsqu","5",392,64.75308 "dicsqu","5",587,59.1738 "dicsqu","5",926,46.05804 "dicsqu","5",1410,38.8311 "dicsqu","6",0,0 "dicsqu","6",33,14.76048 "dicsqu","6",47,20.42964 "dicsqu","6",67,26.60112 "dicsqu","6",113,36.86088 "dicsqu","6",171,43.81776 "dicsqu","6",274,49.3416 "dicsqu","6",392,48.66792 "dicsqu","6",587,45.6834 "dicsqu","6",926,37.12842 "dicsqu","6",1410,33.18294 ##=========================================## #### Eilers-Peeters light response model #### ##=========================================## ## Self-starter Funtion for the Eilers-Peeters light response model accounting for photoinhibition (Eilers & Peeters, 1988) ## Model parameters: ## ETRmax = maxmimum electron transport rate derived from chlorophyll fluorescence ## measurements (to be replaced by PSmax, the maximal rate of photosynthesis when ## dealing with gas exchange measurements) ## PPFDsat = light saturation point, corresponds to the x-value at ETRmax (PPFD = Photosynthetic photon flux density) ## a = light-use-efficiency, initial slope of the curve ep <- function(ppfd,="" etrmax,="" ppfdsat,="" alpha)="" {="" etrmax="" *="" 2="" *="" (1="" +="" alpha)="" *="" (ppfd/ppfdsat)/((ppfd/ppfdsat)^2="" +="" 2="" *="" alpha="" *="" ppfd/ppfdsat="" +="" 1)="" }="" #="" ep=""><- function(ppfd,="" etrmax,="" ppfdsat,="" alpha)="" {="" ppfd="" (="" ppfd^2/(alpha="" *="" ppfdsat^2)="" +="" ppfd/etrmax="" -="" 2="" *="" ppfd/(alpha="" *="" ppfdsat)="" +="" (1/alpha))="" }="" ep_init=""><- function(mcall,="" lhs,="" data)="" #="" initial="" value="" routine="" for="" parameter="" estimation="" {="" xy=""><- sortedxydata(x="mCall[["ppfd"]]," y="LHS," data="data)" #="" xy=""><- sortedxydata("ppfd",="" "etr",="" data="treefern)" #="" etrmax=""><- max(xy$x)="" #="" max="" etr="" rate="" #="" ppfdsat=""><- xy$x[which.max(xy$y)]="" #="" light="" saturation="" #="" alpha=""><- coef(lm(y="" ~="" x,="" data="xy[1:3," ]))[[2]]="" #="" initial="" slope="quantum" yield="" (light="" use="" efficiency)="" sp=""><- smooth.spline(x="xy$x," y="xy$y)" #="" plot(y="" ~="" x,="" data="xy)" xv=""><- seq(min(xy$x),="" max(xy$x),="" length.out="200)" #="" lines(xv,="" predict(object="sp," x="xv)$y)" preds=""><- as.data.frame(predict(sp,="" x="xv," deriv="1))" #="" first="" derivative="" of="" the="" smoothing="" spline="" alpha=""><- mean(preds[preds$x="">< 40,="" "y"])="" #="" initial="" slope="quantum" yield="" etrmax=""><- nlsstrtasymptote(xy)="" #="" max="" etr="" rate="" ppfdsat=""><- nlsstclosestx(xy,="" yval="NLSstRtAsymptote(xy))" #="" light="" saturation="" init.pars=""><- c(etrmax,="" ppfdsat,="" alpha)="" names(init.pars)=""><- mcall[c("etrmax",="" "ppfdsat",="" "alpha")]="" init.pars="" }="" ssep=""><- selfstart(ep, ep_init, c("etrmax", "ppfdsat", "alpha")) selfstart(ep,="" ep_init,="" c("etrmax",="" "ppfdsat",="">