infl.growth(skewtools)
infl.growth()所属R语言包:skewtools
Local Influence Analysis
局部影响分析
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Diagnostic the influential observations of Heteroscedastic Nonlinear Regression Growth models using Local Influence Analysis of Cook with curvature correction
诊断有影响力的意见,异方差的非线性回归增长模型的局部影响分析库克与曲率校正
用法----------Usage----------
infl.growth(y, x, model, plot.it = TRUE)
参数----------Arguments----------
参数:y
a numeric vector of lengths
一个数值向量的长度
参数:x
a numeric vector of ages
的数字矢量年龄的
参数:model
a object related to fitted model of type nls
一个对象类型nls拟合模型
参数:plot.it
a logical value to plot. By default, is TRUE
一个逻辑值的图。默认情况下,是TRUE
值----------Value----------
参数:residuals
a numeric vector of model residuals
一个数值向量模型的残差
参数:outliers
a numeric vector with positions of outliers
一个数值向量与位置异常值
(作者)----------Author(s)----------
Javier E. Contreras-Reyes
参考文献----------References----------
参见----------See Also----------
conf.im, HNL.skew
conf.im,HNL.skew
实例----------Examples----------
data(merluzaChile)
x <- merluzaChile$edad
y <- merluzaChile$long
plot(x, y, main="Age-Length", ylab="Length (cm)", xlab="Age (years)")
beta <- c(80, 0.08, -0.187)
rho <- -0.1
sigma2 <- 3.2726
shape <- 0.1698
nu <- 11
modelVB1 <- HNL.skew(y, x, beta, rho, sigma2, shape, nu, loglik = TRUE,
model = "VB", type = "T", m.type = "power", error = 0.00001)
modelVB2 <- HNL.skew(y, x, beta, rho, sigma2, shape, nu, loglik = TRUE,
model = "VB", type = "N", m.type = "power", error = 0.00001)
I1=infl.growth(y, x, model=modelVB1, plot.it = TRUE)
I2=infl.growth(y, x, model=modelVB2, plot.it = TRUE)
plot(x, y, main="Age-Length", ylab="Length (cm)", xlab="Age (years)")
points(x[I1$outliers],y[I1$outliers],col="red")
legend(10,30,c("T"),col=c("red"),lty=1)
points(x[I2$outliers],y[I2$outliers],col="yellow")
legend(10,20,c("N: tau = 3"),col=c("yellow"),lty=1)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|