Learn R Programming

embryogrowth (version 5.0)

resultNest_4p_weight: Fit using the nest database with weight

Description

Fit using the nest database with weight

Usage

resultNest_4p_weight

Arguments

format

A list with fitted information about data(nest)

Details

Result of the fit using the nest database with weight

References

Girondot, M. & Kaska, Y. Submitted. A model to predict temperature dependency on embryo growth rate and incubation duration from field data.

Examples

Run this code
library(embryogrowth)
data(nest)
formated <- FormatNests(nest)
w <- weightmaxentropy(formated, control_plot=list(xlim=c(20,36)))
x <- structure(c(118.768297442004, 475.750095909406, 306.243694918151,
116.055824800264), .Names = c("DHA", "DHH", "T12H", "Rho25"))
# pfixed <- c(K=82.33) or rK=82.33/39.33
pfixed <- c(rK=2.093313)
# K or rK are not used for dydt.linear or dydt.exponential
resultNest_4p_weight <- searchR(parameters=x,
	fixed.parameters=pfixed, temperatures=formated,
	derivate=dydt.Gompertz, M0=1.7, test=c(Mean=39.33, SD=1.92),
	method = "BFGS", maxiter = 200, weight=w)
data(resultNest_4p_weight)
plotR(resultNest_4p_weight, ylim=c(0,0.50), xlim=c(15, 35))

Run the code above in your browser using DataLab