shoots: Height Growth Data of Bamboo Shoots

In biogeom: Biological Geometries

Height Growth Data of Bamboo Shoots

Description

The height data of four species of bamboo at Nanjing Forestry University campus in 2016.

Usage

data(shoots)

Details

In the data set, there are four columns of variables: Code, LatinName, x, and y. Code saves the number codes of different bamboo species; LatinName saves the Latin names of different bamboo species; x saves the investigation times (days from a specific starting time of growth, and where every bamboo has a different starting time of growth); and y saves the measured aboveground height values (cm).

Code = 1 represents Phyllostachys iridescens, and the starting time (namely time = 0) was defined as 12:00, 3rd April, 2016;

Code = 2 represents Phyllostachys mannii, and the starting time (namely time = 0) was defined as 12:00, 4th April, 2016;

Code = 3 represents Pleioblastus maculatus, and the starting time (namely time = 0) was defined as 12:00, 29th April, 2016;

Code = 4 represents Sinobambusa tootsik, and the starting time (namely time = 0) was defined as 12:00, 18th April, 2016.

References

Shi, P., Fan, M., Ratkowsky, D.A., Huang, J., Wu, H., Chen, L., Fang, S., Zhang, C. (2017) Comparison of two ontogenetic growth equations for animals and plants. Ecological Modelling 349, 1-10. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1016/j.ecolmodel.2017.01.012")}

Examples

data(shoots)
attach(shoots)
# Choose a species
# 1: Phyllostachys iridescens; 2: Phyllostachys mannii; 
# 3: Pleioblastus maculatus; 4: Sinobambusa tootsik
ind <- 4
x1  <- x[Code == ind]
y1  <- y[Code == ind] 

dev.new()
plot(x1, y1, cex=1.5, cex.lab=1.5, cex.axis=1.5, xlab="Time (days)", ylab="Height (cm)")


  delta0  <- c(0.5, 1, 2, 5, 10, 20)
  ini.val <- list(600, 25, 0, 40, delta0)
  resu1   <- fitsigmoid(MLRFE, x=x1, y=y1, ini.val=ini.val, simpver=NULL, 
               fig.opt=TRUE, control=list(reltol=1e-20, maxit=20000), 
               subdivisions = 100L, rel.tol=.Machine$double.eps^0.25,                
               abs.tol=.Machine$double.eps^0.25, stop.on.error=TRUE, 
               keep.xy=FALSE, aux=NULL)


graphics.off()

biogeom documentation built on May 29, 2024, 8:52 a.m.