bandwidth(SNSequate)
bandwidth()所属R语言包：SNSequate
Automatic selection of the bandwidth parameter h
自动选择的带宽参数h
译者：生物统计家园网 机器人LoveR
描述----------Description----------
This functions implements the minimization of the combined penalty function described by Holland and Thayer (1989); Von Davier et al, (2004). It returns the optimal value of h for kernel continuization, according to the above mentioned criteria. Different types of kernels (others than the gaussian) are accepted.
这个功能实现了由荷兰和Thayer（1989）·冯·Davier等，（2004年）合并罚函数的最小化。它返回的最优值h内核continuization的，根据上述标准。不同类型的内核（别人比高斯）都可以接受。
用法----------Usage----------
bandwidth(scores, kert, degree, design, Kp = 1, scores2, degreeXA, degreeYA,
J, K, L, wx, wy, w, ...)
参数----------Arguments----------
参数：
Note that depending on the specified equating design, not all arguments are necessary as detailed below.
请注意，根据指定的等值设计，并非所有的参数是必要的，详情如下。
参数：scores
If the "EG" design is specified, a vector containing the raw sample frequencies coming from one group taking the test. If the "SG" design is specified, a matrix containing the (joint) bivariate sample frequencies for X (raws) and Y (columns). If the "CB" design is specified, a two column matrix containing the observed scores of the sample taking test X first, followed by test Y. The scores2 argument is then used for the scores of the sample taking test Y first followed by test X. If either the "NEAT_CB" or "NEAT_PSE" design is selected, a two column matrix containing the observed scores on test X (first column) and the observed scores on the anchor test A (second column). The scores2 argument is then used for the observed scores on test Y.
如果“EG”的设计，一个向量，包含原始采样频率从一组参加考试。如果指定的“SG”设计，矩阵中的（联合）二元采样频率为X（原糖）和Y（列）。如果“CB”的设计，一个两列的矩阵，包含的样本的观察分数测试X第一，其次是测试Y。 scores2参数，然后用于测试Xy首先其次是样品测试的分数。如果任一被选择时，“NEAT_CB”或“NEAT_PSE”设计一个两列的矩阵含有试验X（第一列上所观察到的分数）和观察分数上锚定测试A（第二列）。然后用所观察到的分数测试scores2Y参数。
参数：kert
A character string giving the type of kernel to be used for continuization. Current options include "gauss", "logis", and "uniform" for the gaussian, logistic and uniform kernels, respectively
的字符字符串，给出了用于continuization类型的内核。目前的选择包括“gauss”，“logis”，和“uniform”高斯，MF和统一的内核，分别
参数：degree
Either a number or vector indicating the number of power moments to be fitted to the marginal distributions, or the number or cross moments to be fitted to the joint distributions, respectively. For the "EG" design it will be a number (see Details).
一个数字或矢量的边缘分布，或数或交叉的时间来安装的联合分布，分别被安装到权力的时刻。对于“EG”的设计，这将是一个数（见详情）。
参数：design
A character string indicating the equating design (one of "EG", "SG", "CB", "NEAT_CE", "NEAT_PSE")
一个字符串，指示的等值设计（之一“EG”，“SG”，“CB”中，“NEAT_CE”，“NEAT_PSE”）
参数：Kp
A number which acts as a weight for the second term in the combined penalization function used to obtain h (see details).
作为第二项的权重的一个数字，它在合并惩罚函数用于获取h（见详情）。
参数：scores2
Only used for the "CB", "NEAT_CE" and "NEAT_PSE" designs. See the description of scores.
仅用于“CB”，“NEAT_CE”和“NEAT_PSE”设计。见的描述scores。
参数：degreeXA
A vector indicating the number of power moments to be fitted to the marginal distributions X and A, and the number or cross moments to be fitted to the joint distribution (X,A) (see details). Only used for the "NEAT_CE" and "NEAT_PSE" designs.
一个向量，表示电源瞬间被安装到的边缘分布的数量X和A，并数或交叉的瞬间被安装到联合分布(X,A)（见详情）。仅用于“NEAT_CE”和“NEAT_PSE”设计。
参数：degreeYA
Only used for the "NEAT_CE" and "NEAT_PSE" designs (see the description for degreeXA)
仅用于“NEAT_CE”和“NEAT_PSE的”设计（见说明degreeXA）
参数：J
The number of possible X scores. Only needed for "CB", "NEAT_CB" and "NEAT_PSE" designs
的数目可能X分数。只需要为“CB”，“NEAT_CB”和“NEAT_PSE”的设计
参数：K
The number of possible Y scores. Only needed for "CB", "NEAT_CB" and "NEAT_PSE" designs
的数目可能Y分数。只需要为“CB”，“NEAT_CB”和“NEAT_PSE”的设计
参数：L
The number of possible A scores. Needed for "NEAT_CB" and "NEAT_PSE" designs
的数目可能A分数。所需的“NEAT_CB”和“NEAT_PSE”设计
参数：wx
A number that satisfies 0<=w_x<=1 indicating the weight put on the data that is not subject to order effects. Only used for the "CB" design.
一个数字，满足0<=w_x<=1表示的重量把上的数据是没有顺序的影响。仅用于“CB”的设计。
参数：wy
A number that satisfies 0<=w_y<=1 indicating the weight put on the data that is not subject to order effects. Only used for the "CB" design.
一个数字，满足0<=w_y<=1表示的重量把上的数据是没有顺序的影响。仅用于“CB”的设计。
参数：w
A number that satisfies 0<=w<=1 indicating the weight given to population P. Only used for the "NEAT" design.
一个数字，满足0<=w<=1表示权重考虑到人口P。仅用于“干净”的设计。
参数：...
Further arguments currently not used.
目前没有进一步的论据。
Details
详细信息----------Details----------
To automatically select h, the function minimizes
要自动选择h，功能最大限度地减少
where PEN_1(h)=∑_j(\hat{r}_j-\hat{f}_h(x_j))^2, and PEN_2(h)=∑_jA_j(1-B_j). The terms A and B are such that PEN_2 acts as a smoothness penalty term that avoids rapid fluctuations in the approximated density (see Chapter 10 in Von Davier, 2011 for more details). The K term corresponds to the Kp argument of the bandwidth function. The \hat{r} values are assumed to be estimated by polynomial loglinear models of specific degree, which come from a call to loglin.smooth.
PEN_1(h)=∑_j(\hat{r}_j-\hat{f}_h(x_j))^2，PEN_2(h)=∑_jA_j(1-B_j)。的条款A和B是PEN_2作为一个平滑的惩罚，避免快速波动的近似的密度（参见第10章在冯·2011年Davier，为更多的细节）。 K项对应的Kpbandwidth函数的参数。 \hat{r}值假设要估计的具体degree，来自的呼叫loglin.smooth的多项式的对数线性模型。
值----------Value----------
A number which is the optimal value of h.
一个数，它的最优值h。
（作者）----------Author(s)----------
Jorge Gonzalez B.
参考文献----------References----------
New York, NY: Springer-Verlag.
参见----------See Also----------
loglin.smooth
loglin.smooth
实例----------Examples----------
#Example: The "Standard" column and firsts two rows of Table 10.1 in [例：“标准”列和第一行表10.1]
#Chapter 10 of Von Davier 2011[第10章冯Davier 2011]
data(Math20EG)
hx.logis<-bandwidth(scores=Math20EG[,1],kert="logis",degree=2,design="EG")$h
hx.unif<-bandwidth(scores=Math20EG[,1],kert="unif",degree=2,design="EG")$h
hx.gauss<-bandwidth(scores=Math20EG[,1],kert="gauss",degree=2,design="EG")$h
hy.logis<-bandwidth(scores=Math20EG[,2],kert="logis",degree=3,design="EG")$h
hy.unif<-bandwidth(scores=Math20EG[,2],kert="unif",degree=3,design="EG")$h
hy.gauss<-bandwidth(scores=Math20EG[,2],kert="gauss",degree=3,design="EG")$h
partialTable10.1<-rbind(c(hx.logis,hx.unif,hx.gauss),
c(hy.logis,hy.unif,hy.gauss))
dimnames(partialTable10.1)<-list(c("h.x","h.y"),c("Logistic","Uniform","Gaussian"))
partialTable10.1
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-30 11:16:37
