R/functions.R
boot.mcor.ic.Rd
Intervalles de confiance bootstrap pour une matrice de corrélation de Bravais-Pearson
boot.mcor.ic(mat, boot.mcor.res, conflevel = 0.95)
mat | Matrice des données |
---|---|
boot.mcor.res | Résultat du bootstrap de la matrice de corrélation de Bravais-Pearson |
conflevel | Niveau de confiance pour les intervalles |
Liste de quatre éléments : matrice des limites inférieures des intervalles de confiance bootstrap percentile, matrice des limites supérieures des intervalles de confiance bootstrap percentile, matrice des limites inférieures des intervalles de confiance bootstrap BCa, matrice des limites supérieures des intervalles de confiance bootstrap BCa
data(Mesures5,package="BioStatR") Mes5_red_lr = subset(Mesures5[,-5],subset=Mesures5$espece=="laurier rose") library(boot) boot.mcor <- boot(Mes5_red_lr[,c("masse","taille","masse_sec")], rho.mult, R=1000) boot.mcor#> #> ORDINARY NONPARAMETRIC BOOTSTRAP #> #> #> Call: #> boot(data = Mes5_red_lr[, c("masse", "taille", "masse_sec")], #> statistic = rho.mult, R = 1000) #> #> #> Bootstrap Statistics : #> original bias std. error #> t1* 0.8407223 0.0002719273 0.02902049 #> t2* 0.6705079 -0.0067273850 0.06751519 #> t3* 0.4150282 -0.0059770037 0.08885902boot.mcor.ic.res <- boot.mcor.ic(Mes5_red_lr[,c("masse", "taille","masse_sec")],boot.mcor) boot.mcor.ic.res#> $cor.ic.percentile.low #> [,1] [,2] [,3] #> [1,] 1.0000000 0.7748695 0.5073956 #> [2,] 0.7748695 1.0000000 0.2063711 #> [3,] 0.5073956 0.2063711 1.0000000 #> #> $cor.ic.percentile.up #> [,1] [,2] [,3] #> [1,] 1.0000000 0.8944618 0.7771115 #> [2,] 0.8944618 1.0000000 0.5726809 #> [3,] 0.7771115 0.5726809 1.0000000 #> #> $cor.ic.BCa.low #> [,1] [,2] [,3] #> [1,] 1.0000000 0.7721474 0.5072410 #> [2,] 0.7721474 1.0000000 0.2143602 #> [3,] 0.5072410 0.2143602 1.0000000 #> #> $cor.ic.BCa.up #> [,1] [,2] [,3] #> [1,] 1.0000000 0.8899512 0.7766753 #> [2,] 0.8899512 1.0000000 0.5797118 #> [3,] 0.7766753 0.5797118 1.0000000 #>