找回密码
 注册
查看: 359|回复: 0

R语言 skewtools包 infl.growth()函数中文帮助文档(中英文对照)

[复制链接]
发表于 2012-9-30 09:52:37 | 显示全部楼层 |阅读模式
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:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2025-5-26 04:53 , Processed in 0.027418 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表