ci_ht_regression_ch: Modelisation de Circonference en fonction de la Hauteur

In jprenaud-02/coopR: Fonctions R pour la Coop

Description

Cette fonction realise la modele Circonference - Hauteur a partir des tiges echantillons, en priorite selon le modele Coop de donnees, sinon a partir d'une regression lineaire.

Usage

ci_ht_regression_ch(Circ, Haut)

Arguments

Circ	: vecteur des circonferences (cm)
Haut	: vecteur des hauteurs totales (cm). Circ et Haut representent les memes individus, classes dans le meme ordre.

Value

liste de 7 elements : - C130.min : la circonference minimale des tiges utilisees pour la modelisation - C130.max : la circonference maximale des tiges utilisees pour la modelisation, - effectif : l'effectif de tiges utilise pour la modelisation, - modL : le modele retenu pour la modelisation (coop ou lineaire), - par1 : la valeur du parametre 1 du modele, - par2 : la valeur du parametre 2 du modele, - Observ. : les observations concernant le choix du modele.

Author(s)

Quentin Girard

References

Protocole Coop chene... document d Ingrid et Claudine

See Also

ci_ht_fonction_ch

Other coopR: accroissement_peuplement_ch, bd_2_xls_ch, calc_D0, calc_Dg, calc_Gha, calc_Ke, calc_RDI, calc_veg, ci_ht_fonction_ch, ci_ht_fonction, ci_ht_regression, graph_accroissement_ch, graph_age_ch, graph_ci_ht_ch, graph_evolution_ch, graph_mortalite, intervention_ch, liste_dispo_ch, load_bdd_to_R_ch, load_bdd_to_xls_ch, load_onf_ch, load_xls_to_R_ch, peuplement_ch

Examples

# la regression Coop fonctionne
Circ.ex1 = c(59,62,35,47,33,48,33,30,37,37,22,37,42,41,43,36,42,26,27,35,30,30,60,49,33,52,49,52,44)
Haut.ex1 = c(1570,1620,1400,1460,1410,1520,1270,1270,1370,1360,1230,1390,1470,1410,1520,1500,1510,1300,1260,1410,1300,NA,1590,1550,1430,1610,1500,1480,1570)
regr1 <- ci_ht_regression_ch(Circ.ex1, Haut.ex1)
regr1
plot(Circ.ex1, Haut.ex1) ; lines(sort(Circ.ex1), ci_ht_fonction_ch(sort(Circ.ex1), regr1$par1, regr1$par2, modL = regr1$modL))
# la regression Coop ne fonctionne pas
Circ.ex2 = c(1:20)
Haut.ex2 = 130 + c(1:20)
regr2 <- ci_ht_regression_ch(Circ.ex2, Haut.ex2)
regr2
plot(Circ.ex2, Haut.ex2) ; lines(sort(Circ.ex2), ci_ht_fonction_ch(sort(Circ.ex2), regr2$par1, regr2$par2, modL = regr2$modL))

jprenaud-02/coopR documentation built on May 3, 2019, 7:06 p.m.