elev.com(smatr)
elev.com()所属R语言包:smatr
Test for equal elevation among several lines fitted with (standardised) major axes of common slope.
配有共同的斜率(标准)的主轴线的几行之间的平等海拔试验。
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Test if several major axis or standardised major axis lines share a common elevation. This can now be done via sma(y~x+groups), see help on the sma function.
测试,如果几大轴或标准化的主要轴线都有一个共同的高程。现在可以通过sma(y~x+groups)信息,请参阅上的sma功能的帮助。
用法----------Usage----------
elev.com(y, x, groups, data = NULL, method = 'SMA', alpha = 0.05,
V = array(0, c(2, 2, length(unique(groups)))),
group.names = sort(unique(groups)))
参数----------Arguments----------
参数:y
The Y-variable for all observations (as a vector).
Y-的所有观测的变量(矢量)。
参数:x
The X-variable for all observations (as a vector).
X-变量的所有观测值(矢量)。
参数:groups
Coding variable identifying which group each observation belongs to (as a factor or vector).
编码变量确定每个观测值属于哪个组的一个因素或矢量。
参数:data
(Optional) data frame containing the data.
(可选)数据框包含的数据。
参数:method
The line fitting method:
该生产线拟合的方法:
'SMA' or 1 standardised major axis (this is the default)
“SMA”或标准化的长轴(这是默认设置)
'MA' or 2 major axis
MA或2轴
参数:alpha
The desired confidence level for the 100(1-alpha)% confidence interval for the common elevation. (Default value is 0.05, which returns a 95% confidence interval.)
所需的置信水平为100(1-α)%的置信区间的共同高程。 (默认值是0.05,它返回一个95%的置信区间)。
参数:V
The estimated variance matrices of measurement error, for each group. This is a 3-dimensional array with measurement error in Y in the first row and column, error in X in the second row and column,and groups running along the third dimension. Default is that there is no measurement error.
测量误差的估计方差矩阵,每个组。与Y中的测量误差的第一行和列中,X中的第二行和列中的错误,并沿第三维度上运行的组中,这是一个3维阵列。默认是不存在测量误差。
参数:group.names
(optional: rarely required). A vector containing the labels for "groups". (Only actually useful for reducing computation time in simulation work).
(可选:很少需要)。一个向量“组”的标签。 (实际上只用在模拟工作,减少计算时间)。
Details
详细信息----------Details----------
Calculates a Wald statistic to test for equal elevation of several MA's or SMA's with a common slope. This is done by testing for equal mean residual scores across groups.
计算Wald统计量进行测试的几个MA或SMA的一个共同的斜率等于抬高。这样做等于平均剩余分数,各组的测试。
Note that this test is only valid if it is reasonable to assume that the axes for the different groups all have the same slope.
注意,这个测试是有效的,如果它是合理的假设,为不同的组所有轴具有相同的斜率。
The test assumes the following:
测试假设如下:
each group of observations was independently sampled
观察各组的独立采样
the axes fitted to all groups have a common slope
轴安装到所有群体都有一个共同的斜率
y and x are linearly related within each group
y和x是线性相关的每个组内的
residual scores independently follow a normal distribution with equal variance at all points along the line, within each group
剩余的分数,独立遵循沿线的所有点等于方差的正态分布,各组内
Note that we do not need to assume equal variance across groups, unlike in tests comparing several linear regression lines.
请注意,我们并不需要承担平等的差异,各组不同,在测试中,比较了几种线性回归线。
The assumptions can be visually checked by plotting residual scores against fitted axis scores, and by constructing a Q-Q plot of residuals against a normal distribution, available using the plot.sma function. On a residual plot, if there is a distinct increasing or decreasing trend within any of the groups, this suggests that all groups do not share a common slope.
假设可以肉眼检查剩余的分数对装轴的分数,并通过构建的残留物对正常的分布,可使用plot.sma功能的QQ图的绘制。剩余的图,如果有一个明显的增加或减少的趋势在任何一组,这表明,所有群体都有一个共同的斜率。
The common slope (b) is estimated from a maximum of 100 iterations, convergence is reached when the change in b < 10^-6.
共同的斜率(b),估计从最高的100次迭代时,达到收敛的变化b < 10^-6。
值----------Value----------
参数:stat
The Wald statistic testing for no shift along the common axis
Wald统计量的共同轴线沿线无移位测试
参数:p
The P-value of the test. This is calculated assuming that stat has a chi-square distribution with (g-1) df, if there are g groups
P-值的测试。这是计算假设统计,卡方分布(G-1),DF,如果有G功能组
参数:a
The estimated common elevation
估计常见的海拔
参数:ci
A 100(1-alpha)% confidence interval for the true common elevation
A 100(1-α)%置信区间为真正的共同抬高
参数:as
Separate elevation estimates for each group
各组独立的高程估计
(作者)----------Author(s)----------
Warton, D. <a href="mailto avid.Warton@unsw.edu.au">David.Warton@unsw.edu.au</a>, translated to R by Ormerod, J. 2005-12-08
参考文献----------References----------
<h3>See Also</h3>
实例----------Examples----------
# Load leaf longevity data[加载叶寿命数据]
data(leaflife)
# Test for common SMA slope amongst species at low soil nutrient sites[测试常见的物种SMA斜率之间的在低土壤养分网站]
# with different rainfall:[不同的降雨:]
leaf.low.soilp <- subset(leaflife, soilp == 'low')
slope.com(log10(longev), log10(lma), rain, data=leaf.low.soilp)
# Now test for common elevation of the groups fitted with an axis[立即测试装有一个轴的共同的海拔的基团]
# of common slope, at low soil nutrient sites:[共同坡度,在低土壤养分网站:]
elev.com(log10(longev), log10(lma), rain, data = leaf.low.soilp)
# Or test for common elevation amongst the MA's of common slope,[或测试之间的共同海拔MA的共同的斜率,]
# for low soil nutrient sites, and construct 99% a confidence interval[为低土壤养分网站,并构建99%的置信区间]
# for the common elevation:[的共同海拔:]
elev.com(log10(longev), log10(lma), rain, method='MA',
data = leaf.low.soilp, alpha=0.01)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|