2。方法让<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M135">
F米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i mathvariant="bold-italic">
η米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
,米米l:米o><米米l:米i mathvariant="bold-italic">
x米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
是一个1参数统计模型(
4)gydF4y2Ba参数化<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M136">
θ米米l:米我><米米l:mo>
∈米米l:米o><米米l:米i>
Θ米米l:米我><米米l:mo>
⊂米米l:米o><米米l:米i mathvariant="double-struck">
R米米l:米我>米米l:math>
对于一些参数空间<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M137">
Θ米米l:米我>米米l:mrow>
居住在<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M138">
(米米l:米o><米米l:米o>
- - - - - -米米l:米o><米米l:米i>
∞米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
∞米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
。也让<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M139">
问米米l:米我>米米l:mrow>
y米米l:米我>米米l:mrow>
τ米米l:米我><米米l:mrow>
∣米米l:米o>米米l:mrow>
x米米l:米我>米米l:mrow>
对应的分位数函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M140">
τ米米l:米我><米米l:mo>
∈米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
0 1米米l:米n><米米l:米o stretchy="false">
)米米l:米o>米米l:math>
。至关重要的是,注意,(
5)gydF4y2Ba处理的合适的转换<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M141">
θ米米l:米我>米米l:mrow>
在线性预测参数链接功能。相比之下我们的提议小说直接建模<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M142">
问米米l:米我>米米l:mrow>
y米米l:米我>米米l:mrow>
τ米米l:米我><米米l:mrow>
∣米米l:米o>米米l:mrow>
x米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
通过平滑和一对一的函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M143">
G米米l:米我>米米l:mrow>
的形式<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M144">
(13)米米l:米text>米ml:mtd>
η米米l:米我>米米l:mrow>
τ米米l:米我>米米l:mrow>
=米米l:米o><米米l:米i mathvariant="script">
G米米l:米我><米米l:mfenced separators="|">
问米米l:米我>米米l:mrow>
y米米l:米我>米米l:mrow>
τ米米l:米我><米米l:mrow>
∣米米l:米o>米米l:mrow>
x米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
=米米l:米o><米米l:米年代up>
G米米l:米我>米米l:mrow>
∗米米l:米我>米米l:mrow>
θ米米l:米我><米米l:mrow>
∣米米l:米o>米米l:mrow>
x米米l:米我><米米l:mo>
,米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我>米米l:mrow>
,米米l:米o>米米l:mtd>
也被纳入到VGLM / VGAM对数似,即(
9)gydF4y2Ba和(
12)gydF4y2Ba。在这里,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M145">
τ米米l:米我><米米l:mo>
=米米l:米o><米米l:米年代up>
τ米米l:米我>米米l:mrow>
1米米l:米n>米米l:mrow>
,米米l:米o><米米l:米o>
…米米l:米o><米米l:米o>
,米米l:米o><米米l:米年代ub>
τ米米l:米我>米米l:mrow>
l米米l:米我>米米l:mrow>
T米米l:米我>米米l:mrow>
是一个指定向量分位数的利息。的例子(
13)<米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M146">
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
ξ米米l:米我>米米l:mrow>
τ米米l:米我>米米l:mrow>
,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M147">
l米米l:米我><米米l:mi mathvariant="normal">
o米米l:米我><米米l:mi mathvariant="normal">
g米米l:米我><米米l:mi mathvariant="normal">
我米米l:米我><米米l:mi mathvariant="normal">
t米米l:米我><米米l:mi>
ξ米米l:米我>米米l:mrow>
τ米米l:米我>米米l:mrow>
,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M148">
ξ米米l:米我>米米l:mrow>
τ米米l:米我>米米l:mrow>
。
gydF4y2B一个方程(
13)gydF4y2Ba是这项工作的核心。它允许通过造型的选择<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M149">
G米米l:米我>米米l:mrow>
分位数函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M150">
问米米l:米我>米米l:mrow>
y米米l:米我>米米l:mrow>
,它代表了一个新的修改VGLM / VGAM框架。请注意,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M151">
G米米l:米我>米米l:mrow>
就像一个<我t一个l我c>链接功能我t一个l我c>VGLM / VGAM框架内表
1gydF4y2B一个。两个音符:首先,没有普遍性,(
13)gydF4y2Ba可以看到(严格)的函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M152">
θ米米l:米我>米米l:mrow>
自从分位数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M153">
τ米米l:米我>米米l:mrow>
协变量和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M154">
x米米l:米我>米米l:mrow>
是已知的。其次,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M155">
G米米l:米我>米米l:mrow>
∗米米l:米我>米米l:mrow>
是单调和一对一的,作为一个复合的结果吗<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M156">
G米米l:米我>米米l:mrow>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M157">
问米米l:米我>米米l:mrow>
y米米l:米我>米米l:mrow>
也具有这样的特性。然而,在拟合过程中,内部irl算法需要的倒数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M158">
G米米l:米我>米米l:mrow>
∗米米l:米我>米米l:mrow>
。处理1参数分布在这个阶段通过费舍尔得分,因为逆简化实现<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M159">
G米米l:米我>米米l:mrow>
∗米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:mrow>
可以手动,然后纳入irl算法。在一些情况下,逆的<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M160">
G米米l:米我>米米l:mrow>
∗米米l:米我>米米l:mrow>
没有一个封闭的形式,如1参数γ分布,和另一个迭代法是用来近似<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M161">
G米米l:米我>米米l:mrow>
∗米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:mrow>
。为了实现这一有效,两个选择是可用的。这些都是(gydF4y2B一个一个)new来nR一个ph年代on.basic ()米ono年代p一个ce>从<年代一个n年代- - - - - -年代er我f>VGAMextra和(b)VGAM:b我年代ect我on.basic ()米ono年代p一个ce>,两个<我t一个l我c>矢量化我t一个l我c>著名的牛顿迭代的实现和二分算法,来解决一个实值函数的根在给定的时间间隔<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M162">
(米米l:米o><米米l:米i>
一个米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
b米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
。进一步给出了部分细节
2.2gydF4y2B一个。
gydF4y2B一个这项工作的一个优点是,VGLM / VGAM框架可以规避分位数穿越的问题(例如,
2,
11),gydF4y2Ba教派。2.5)通过选择<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M163">
H米米l:米我>米米l:mrow>
1米米l:米n>米米l:mrow>
=米米l:米o><米米l:米年代ub>
我米米l:米我>米米l:mrow>
米米米l:米我>米米l:mrow>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M164">
H米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
=米米l:米o><米米l:米年代ub>
H米米l:米我>米米l:mrow>
3米米l:米n>米米l:mrow>
=米米l:米o><米米l:米o>
⋯米米l:米o><米米l:米o>
=米米l:米o><米米l:米年代ub>
H米米l:米我>米米l:mrow>
d米米l:米我>米米l:mrow>
=米米l:米o><米米l:米年代ub>
1米米l:米n>米米l:mrow>
米米米l:米我>米米l:mrow>
(一个<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M165">
米米米l:米我>米米l:mrow>
向量的)。在这种并行性假设的方法借力量在整个数据集,以便添加剂模型<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M166">
η米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
关于<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M167">
x米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
,米米l:米o><米米l:米o>
…米米l:米o><米米l:米o>
,米米l:米o><米米l:米年代ub>
x米米l:米我>米米l:mrow>
d米米l:米我>米米l:mrow>
是平行的。每个家庭功能表
2gydF4y2B一个和
3gydF4y2B一个有一个<米ono年代p一个ce>平行米ono年代p一个ce>论点是<米ono年代p一个ce>假米ono年代p一个ce>默认情况下。使用的语法<年代一个n年代- - - - - -年代er我f>VGAM(基于钱伯斯和黑斯蒂
12),gydF4y2Ba设置<米ono年代p一个ce>平行= TRUE米ono年代p一个ce>(或<米ono年代p一个ce>平行= FALSE ~ 1米ono年代p一个ce>)的结果<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M168">
H米米l:米我>米米l:mrow>
1米米l:米n>米米l:mrow>
=米米l:米o><米米l:米年代ub>
我米米l:米我>米米l:mrow>
米米米l:米我>米米l:mrow>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M169">
H米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
=米米l:米o><米米l:米年代ub>
H米米l:米我>米米l:mrow>
3米米l:米n>米米l:mrow>
=米米l:米o><米米l:米o>
⋯米米l:米o><米米l:米o>
=米米l:米o><米米l:米年代ub>
1米米l:米n>米米l:mrow>
米米米l:米我>米米l:mrow>
;即。,我t我年代f一个l年代e为only the intercept.
理想情况下,链接转换的支持<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M195">
Y米米l:米我>米米l:mrow>
来<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M196">
R米米l:米我>米米l:mrow>
因为<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M197">
η米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
应该是无限的。三种最常见的情况如下。为<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M198">
Y米米l:米我><米米l:mo>
∈米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
0米米l:米n><米米l:米o>
,米米l:米o><米米l:米i>
∞米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
推荐一个日志链接,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M199">
Y米米l:米我><米米l:mo>
∈米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
0 1米米l:米n><米米l:米o stretchy="false">
)米米l:米o>米米l:math>
logit链接是一个不错的选择<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M200">
Y米米l:米我><米米l:mo>
∈米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米o>
- - - - - -米米l:米o><米米l:米i>
∞米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
∞米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
意味着一个身份是自然的链接。这些病例实施七1参数分布。函数的选择<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M201">
G米米l:米我>米米l:mrow>
为每一个<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M202">
问米米l:米我>米米l:mrow>
第五列所示的表吗
2,gydF4y2Ba而产生的分位数链接的函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M203">
θ米米l:米我>米米l:mrow>
最后一列所示。
gydF4y2B一个我们现在描述的分位数链接指数和Topp-Leone分布为例。首先,对<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M204">
Y米米l:米我><米米l:mo>
~米米l:米o><米米l:米i mathvariant="normal">
e米米l:米我><米米l:mi mathvariant="normal">
x米米l:米我><米米l:mi mathvariant="normal">
p米米l:米我><米米l:mi mathvariant="normal">
o米米l:米我><米米l:mi mathvariant="normal">
n米米l:米我><米米l:mi mathvariant="normal">
e米米l:米我><米米l:mi mathvariant="normal">
n米米l:米我><米米l:mi mathvariant="normal">
t米米l:米我><米米l:mi mathvariant="normal">
我米米l:米我><米米l:mi mathvariant="normal">
一个米米l:米我><米米l:mi mathvariant="normal">
l米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
,速度参数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M205">
θ米米l:米我><米米l:mo>
>米米l:米o><米米l:米n mathvariant="normal">
0米米l:米n>米米l:math>
,密度和CDF给出<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M206">
f米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米i>
θ米米l:米我><米米l:msup>
e米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米i>
θ米米l:米我><米米l:mi>
y米米l:米我>米米l:mrow>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M207">
F米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米年代up>
e米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米i>
θ米米l:米我><米米l:mi>
y米米l:米我>米米l:mrow>
。有轻微改变符号,分位数函数是由<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M208">
F米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:mrow>
,也就是说,<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M209">
(14)米米l:米text>米ml:mtd>
问米米l:米我>米米l:mrow>
Y米米l:米我>米米l:mrow>
τ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
=米米l:米o><米米l:米o>
- - - - - -米米l:米o><米米l:米frac>
1米米l:米n>米米l:mrow>
θ米米l:米我>米米l:mrow>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我>米米l:mrow>
,米米l:米o>米米l:mtd>
坐落在<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M210">
(米米l:米o><米米l:米n mathvariant="normal">
0米米l:米n><米米l:米o>
,米米l:米o><米米l:米i>
∞米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
不管的值<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M211">
τ米米l:米我>米米l:mrow>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M212">
θ米米l:米我>米米l:mrow>
。考虑到的值<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M213">
τ米米l:米我>米米l:mrow>
(由用户指定),(
14)gydF4y2Ba的函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M214">
θ米米l:米我>米米l:mrow>
。因此,新分位数链接指数分布如表所示
2gydF4y2B一个只是获得了吗<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M215">
G米米l:米我>米米l:mrow>
对数变换,如下:<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M216">
(15)米米l:米text>米ml:mtd>
η米米l:米我><米米l:mfenced separators="" open="(" close=")">
θ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我>米米l:mrow>
=米米l:米o><米米l:米row>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
- - - - - -米米l:米o><米米l:米frac>
1米米l:米n>米米l:mrow>
θ米米l:米我>米米l:mrow>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我>米米l:mrow>
=米米l:米o><米米l:米row>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
。米米l:米o>米米l:mtd>
这个分位数链接已经实现<年代一个n年代- - - - - -年代er我f>VGAMextra通过函数<米ono年代p一个ce>经验值问l我nk(),如表所示
3gydF4y2B一个。它的逆矩阵(表示<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M217">
θ米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
η米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
从逆(获得)可以手动
15)gydF4y2Ba。注意,相应的家庭功能(<米ono年代p一个ce>指数()米ono年代p一个ce>)中实现<年代一个n年代- - - - - -年代er我f>VGAM包括一个(已知)位置参数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M218">
一个米米l:米我>米米l:mrow>
,提供密度<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M219">
f米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米i>
θ米米l:米我><米米l:msup>
e米米l:米我>米米l:mrow>
- - - - - -米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
- - - - - -米米l:米o><米米l:米i>
一个米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:mrow>
。默认情况下<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M220">
一个米米l:米我><米米l:mo>
=米米l:米o><米米l:米n mathvariant="normal">
0米米l:米n>米米l:math>
,它是由参数<米ono年代p一个ce>位置米ono年代p一个ce>。
gydF4y2B一个其次,考虑到Topp-Leone分布<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M221">
Y米米l:米我><米米l:mo>
~米米l:米o><米米l:米i mathvariant="normal">
T米米l:米我><米米l:mi mathvariant="normal">
o米米l:米我><米米l:mi mathvariant="normal">
p米米l:米我><米米l:mi mathvariant="normal">
p米米l:米我><米米l:mo>
- - - - - -米米l:米o><米米l:米i mathvariant="normal">
l米米l:米我><米米l:mi mathvariant="normal">
e米米l:米我><米米l:mi mathvariant="normal">
o米米l:米我><米米l:mi mathvariant="normal">
n米米l:米我><米米l:mi mathvariant="normal">
e米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
是谁的支持<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M222">
(米米l:米o><米米l:米n mathvariant="normal">
0 1米米l:米n><米米l:米o stretchy="false">
)米米l:米o>米米l:math>
和<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M223">
(16)米米l:米text>米ml:mtd>
问米米l:米我>米米l:mrow>
Y米米l:米我>米米l:mrow>
年代米米l:米我><米米l:mo>
;米米l:米o><米米l:米年代ub>
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
=米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米年代qrt>
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米年代ubsup>
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
年代米米l:米我>米米l:mrow>
,米米l:米o>米米l:mtd>
与<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M224">
0米米l:米n><米米l:米o>
<米米l:米o><米米l:米年代ub>
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:math>
。在这里,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M225">
θ米米l:米我><米米l:mo>
=米米l:米o><米米l:米i>
年代米米l:米我>米米l:math>
。来验证这个限制注意<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M226">
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo>
>米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:math>
,对于任何形状参数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M227">
年代米米l:米我><米米l:mo>
∈米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
0 1米米l:米n><米米l:米o stretchy="false">
)米米l:米o>米米l:math>
,因此对于任何<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M228">
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
∈米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
0 1米米l:米n><米米l:米o stretchy="false">
)米米l:米o>米米l:math>
,<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M229">
(17)米米l:米text>米ml:mtd>
0米米l:米n><米米l:米o>
<米米l:米o><米米l:米年代ubsup>
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
年代米米l:米我>米米l:mrow>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
⟺米米l:米o>米米l:mtd>
0米米l:米n><米米l:米o>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米年代ubsup>
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
年代米米l:米我>米米l:mrow>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
⟺米米l:米o>米米l:mtd>
0米米l:米n><米米l:米o>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米年代qrt>
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米年代ubsup>
τ米米l:米我>米米l:mrow>
j米米l:米我>米米l:mrow>
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
年代米米l:米我>米米l:mrow>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
。米米l:米o>米米l:mtd>
因此,允许分位数函数来模拟,我们以分对数转换为<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M230">
G米米l:米我>米米l:mrow>
。这个分布是产生的分位数链接<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M231">
η米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
τ米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米row>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
问米米l:米我>米米l:mrow>
Y米米l:米我>米米l:mrow>
(米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
τ米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:mrow>
,如表所示
2gydF4y2B一个。分配CDF<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M232">
F米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米年代up>
y米米l:米我><米米l:mo>
·米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
2米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i>
y米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:mrow>
年代米米l:米我>米米l:mrow>
为<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M233">
0米米l:米n><米米l:米o>
<米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:math>
,密度<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M234">
f米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米n mathvariant="normal">
2米米l:米n><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i>
y米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
·米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
2米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i>
y米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米年代up>
]米米l:米o>米米l:mrow>
年代米米l:米我><米米l:mo>
- - - - - -米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:mrow>
。分位数函数来源于解方程<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M235">
τ米米l:米我>米米l:mrow>
0米米l:米n>米米l:mrow>
=米米l:米o><米米l:米i>
F米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
y米米l:米我><米米l:mo>
;米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米fenced open="[" close="]" separators="|">
y米米l:米我><米米l:mo>
·米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米n mathvariant="normal">
2米米l:米n><米米l:米o>
- - - - - -米米l:米o><米米l:米i>
y米米l:米我><米米l:msup>
)米米l:米o>米米l:mrow>
年代米米l:米我>米米l:mrow>
,因为<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M236">
0米米l:米n><米米l:米o>
<米米l:米o><米米l:米年代ub>
τ米米l:米我>米米l:mrow>
0米米l:米n>米米l:mrow>
<米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:math>
,从而导致二次方程<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M237">
y米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
- - - - - -米米l:米o><米米l:米n mathvariant="normal">
2米米l:米n><米米l:米i>
y米米l:米我><米米l:mo>
+米米l:米o><米米l:米年代ubsup>
τ米米l:米我>米米l:mrow>
0米米l:米n>米米l:mrow>
1米米l:米n><米米l:米o>
/米米l:米o><米米l:米i>
年代米米l:米我>米米l:mrow>
=米米l:米o><米米l:米n mathvariant="normal">
0米米l:米n>米米l:math>
。解决方案必须在撒谎<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M238">
(米米l:米o><米米l:米n mathvariant="normal">
0 1米米l:米n><米米l:米o stretchy="false">
)米米l:米o>米米l:math>
和实际上是(
16),gydF4y2Ba的函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M239">
年代米米l:米我>米米l:mrow>
。家庭的功能<米ono年代p一个ce>推翻()米ono年代p一个ce>从<年代一个n年代- - - - - -年代er我f>VGAM估计<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M240">
年代米米l:米我>米米l:mrow>
,默认的链接<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M241">
η米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
)米米l:米o><米米l:米o>
=米米l:米o><米米l:米i mathvariant="monospace">
l米米l:米我><米米l:mi mathvariant="monospace">
o米米l:米我><米米l:mi mathvariant="monospace">
g米米l:米我><米米l:mi mathvariant="monospace">
我米米l:米我><米米l:mi mathvariant="monospace">
t米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
年代米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
。
年代ec><年代ec id="sec2.2">
2.2。软件实现实际使用的他人,我们实现了七VGLM-quantile链接,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M242">
η米米l:米我>米米l:mrow>
τ米米l:米我>米米l:mrow>
在<年代一个n年代- - - - - -年代er我f>R包<年代一个n年代- - - - - -年代er我f>VGAMextra。总结在表
2gydF4y2B一个。这个包<年代一个n年代- - - - - -年代er我f>VGAM是一个要求<年代一个n年代- - - - - -年代er我f>VGAMextra因为造型功能<米ono年代p一个ce>vglm ()米ono年代p一个ce>和<米ono年代p一个ce>vgam ()米ono年代p一个ce>,但最后家庭表的函数
2,gydF4y2Ba住在那里。本文对<年代一个n年代- - - - - -年代er我f>VGAMextra 0。<年代一个n年代- - - - - -年代er我f>0- - - - - -- - - - - - - - - -<年代一个n年代- - - - - -年代er我f>2和<年代一个n年代- - - - - -年代er我f>VGAM 1。<年代一个n年代- - - - - -年代er我f>1- - - - - -- - - - - - - - - -<年代一个n年代- - - - - -年代er我f>0或更高要求;他们是可用的<米ono年代p一个ce>www。年代t一个t.auckland.ac.nz ~ vmir178米ono年代p一个ce>和<米ono年代p一个ce>www。年代t一个t.auckland.ac.nz ~仪/ VGAM /预映/米ono年代p一个ce>在旧版本上可用的都是凹口(<米ono年代p一个ce>http://CRAN.R-project.org)。
gydF4y2B一个一种特殊情况是<米ono年代p一个ce>g一个米米一个1问link ()米ono年代p一个ce>,1参数(形状)伽马分布,定义为<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M243">
(18)米米l:米text>米ml:mtd>
η米米l:米我><米米l:mfenced separators="" open="(" close=")">
θ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我>米米l:mrow>
=米米l:米o><米米l:米row>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
g米米l:米我><米米l:mi mathvariant="monospace">
一个米米l:米我><米米l:mi mathvariant="monospace">
米米米l:米我><米米l:mi mathvariant="monospace">
米米米l:米我><米米l:mi mathvariant="monospace">
一个米米l:米我><米米l:mfenced separators="" open="(" close=")">
τ米米l:米我><米米l:mo>
,米米l:米o><米米l:米i mathvariant="monospace">
年代米米l:米我><米米l:mi mathvariant="monospace">
h米米l:米我><米米l:mi mathvariant="monospace">
一个米米l:米我><米米l:mi mathvariant="monospace">
p米米l:米我><米米l:mi mathvariant="monospace">
e米米l:米我><米米l:mo>
=米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
,米米l:米o>米米l:mtd>
的主要参数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M244">
τ米米l:米我>米米l:mrow>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M245">
θ米米l:米我>米米l:mrow>
。它的逆矩阵(表
3)gydF4y2Ba不承认一个封闭的形式,它是近似的函数<米ono年代p一个ce>VGAM:new来nR一个phson.basic ()米ono年代p一个ce>,一个<我t一个l我c>矢量化我t一个l我c>牛顿迭代算法的实现。几乎所有其他实现的是一个标量参数,但是我们使用向量的长度<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M246">
n米米l:米我>米米l:mrow>
。它的工作原理如下。我们的数据是有效的<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M247">
y米米l:米我>米米l:mrow>
我米米l:米我>米米l:mrow>
=米米l:米o><米米l:米o stretchy="false">
{米米l:米o><米米l:米年代ub>
y米米l:米我>米米l:mrow>
我米米l:米我>米米l:mrow>
,米米l:米o><米米l:米年代ub>
x米米l:米我>米米l:mrow>
我米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
d米米l:米我>米米l:mrow>
}米米l:米o><米米l:米o>
,米米l:米o><米米l:米o>
我米米l:米我><米米l:mo>
=米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n><米米l:米o>
,米米l:米o><米米l:米o>
…米米l:米o><米米l:米o>
,米米l:米o><米米l:米i>
n米米l:米我>米米l:math>
感兴趣的,而分位数,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M248">
τ米米l:米我>米米l:mrow>
或<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M249">
p米米l:米我>米米l:mrow>
,必须由用户输入。形状参数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M250">
θ米米l:米我>米米l:mrow>
由irl估计,因此它可以在每个迭代。因此,对于每一个<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M251">
η米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
p米米l:米我><米米l:msub>
)米米l:米o>米米l:mrow>
0米米l:米n>米米l:mrow>
,“逆”是由根,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M252">
θ米米l:米我>米米l:mrow>
的函数,<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M253">
(19)米米l:米text>米ml:mtd>
f米米l:米我><米米l:mfenced separators="" open="(" close=")">
θ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
p米米l:米我><米米l:mo>
,米米l:米o><米米l:米i>
η米米l:米我>米米l:mrow>
=米米l:米o><米米l:米i>
η米米l:米我><米米l:msub>
θ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
p米米l:米我>米米l:mrow>
0米米l:米n>米米l:mrow>
- - - - - -米米l:米o><米米l:米row>
日志米米l:米我>米米l:mrow>
米米l:米o><米米l:米row>
g米米l:米我><米米l:mi mathvariant="monospace">
一个米米l:米我><米米l:mi mathvariant="monospace">
米米米l:米我><米米l:mi mathvariant="monospace">
米米米l:米我><米米l:mi mathvariant="monospace">
一个米米l:米我><米米l:mfenced separators="|">
τ米米l:米我><米米l:mo>
,米米l:米o><米米l:米i mathvariant="monospace">
年代米米l:米我><米米l:mi mathvariant="monospace">
h米米l:米我><米米l:mi mathvariant="monospace">
一个米米l:米我><米米l:mi mathvariant="monospace">
p米米l:米我><米米l:mi mathvariant="monospace">
e米米l:米我><米米l:mo>
=米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
。米米l:米o>米米l:mtd>
最后,所有的逆VGLM-quantile链接表所示
3gydF4y2B一个的名称,以及相应的实现<年代一个n年代- - - - - -年代er我f>VGAMextra。需要inverse-links在不同阶段的irl费舍尔得分,而内部交换机之间<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M254">
η米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
θ米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
(即表
2),<米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M255">
θ米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i>
η米米l:米我><米米l:mo>
;米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
(即表
3)gydF4y2Ba。具体来说,该算法需要得分向量和各irl eim迭代,给出的链式法则公式如下:<米米l:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M256">
(20)米米l:米text>米ml:mtd>
∂米米l:米o><米米l:米i mathvariant="script">
l米米l:米我>米米l:mrow>
∂米米l:米o><米米l:米i>
η米米l:米我>米米l:mrow>
=米米l:米o><米米l:米frac>
∂米米l:米o><米米l:米i mathvariant="script">
l米米l:米我>米米l:mrow>
∂米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
·米米l:米o><米米l:米frac>
∂米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
∂米米l:米o><米米l:米i>
η米米l:米我>米米l:mrow>
,米米l:米o>米米l:mtd>
- - - - - -米米l:米o><米米l:米i mathvariant="double-struck">
E米米l:米我><米米l:mfenced open="[" close="]" separators="|">
∂米米l:米o>米米l:mrow>
2米米l:米n>米米l:mrow>
l米米l:米我>米米l:mrow>
∂米米l:米o><米米l:米年代up>
η米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
=米米l:米o><米米l:米o>
- - - - - -米米l:米o><米米l:米i mathvariant="double-struck">
E米米l:米我><米米l:mfenced open="[" close="]" separators="|">
∂米米l:米o>米米l:mrow>
2米米l:米n>米米l:mrow>
l米米l:米我>米米l:mrow>
∂米米l:米o><米米l:米年代up>
θ米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
∂米米l:米o><米米l:米i>
θ米米l:米我>米米l:mrow>
∂米米l:米o><米米l:米i>
η米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
。米米l:米o>米米l:mtd>
在内部,利用计算逆函数<米ono年代p一个ce>VGAM:埃塔2th埃塔()米ono年代p一个ce>或<米ono年代p一个ce>VGAM:th埃塔2et一个 ()米ono年代p一个ce>。的<年代一个n年代- - - - - -年代er我f>VGAMextra手册和Miranda-Soberanis [
13)gydF4y2Ba提供进一步的细节的推导分位数链接,而绮(
6)gydF4y2Ba描述irl和费舍尔评分算法估算VGLMs和VGAMs。补充第二作者的主页给额外的细节链接功能。
年代ec><年代ec id="sec2.3">
2.3。软件使用像往常一样为用户,这种方法通过调用造型功能<米ono年代p一个ce>VGAM:vglm ()米ono年代p一个ce>和<米ono年代p一个ce>VGAM:VGAM()下面描述的,除了两个修改。
gydF4y2B一个首先,我们给出以下输出显示中央参数处理<米ono年代p一个ce>VGAM:vglm ()米ono年代p一个ce>:
<我nl我ne- - - - - -graphic xlink:href="//www.newsama.com/downloads/journals/jps/2019/3493628.text-inline.001">
第一次调整发生争论<米ono年代p一个ce>公式米ono年代p一个ce>,一个象征性的描述模型。通常,一个表达式<米ono年代p一个ce>y~x2 + x3响应应该足够了吗<米ono年代p一个ce>y协变量和<米ono年代p一个ce>x2和<米ono年代p一个ce>x3。这有效地适用于单变量甚至多个回答说<米ono年代p一个ce>y1, y2,和<米ono年代p一个ce>y3,唯一的变化是设置<米ono年代p一个ce>cb我nd(y1、y2、y3)米ono年代p一个ce>~x2 + x3。这里,右边(RHS)公式应用于每个线性预测。
gydF4y2B一个分位数模型使用VGLMs VGAMs,<米ono年代p一个ce>Q。reg()必须纳入公式的参数见表
4gydF4y2B一个。对于一个给定的分位数的兴趣,进入<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M257">
τ米米l:米我><米米l:mo>
=米米l:米o><米米l:米o stretchy="false">
(米米l:米o><米米l:米年代ub>
τ米米l:米我>米米l:mrow>
1米米l:米n>米米l:mrow>
,米米l:米o><米米l:米o>
…米米l:米o><米米l:米o>
,米米l:米o><米米l:米年代ub>
τ米米l:米我>米米l:mrow>
l米米l:米我>米米l:mrow>
)米米l:米o>米米l:mrow>
T米米l:米我>米米l:mrow>
,<米ono年代p一个ce>Q。reg()复制的响应矩阵YgydF4y2B一个成<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M258">
N米米l:米我><米米l:mi>
O米米l:米我><米米l:mi>
年代米米l:米我><米米l:mo>
·米米l:米o><米米l:米i mathvariant="normal">
d米米l:米我><米米l:mi mathvariant="normal">
我米米l:米我><米米l:mi mathvariant="normal">
米米米l:米我><米米l:mo stretchy="false">
(米米l:米o><米米l:米i mathvariant="bold-italic">
τ米米l:米我><米米l:mo stretchy="false">
)米米l:米o>米米l:math>
列,<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M259">
N米米l:米我><米米l:mi>
O米米l:米我><米米l:mi>
年代米米l:米我>米米l:math>
表示的列数YgydF4y2B一个。然后,RHS列的公式适用于每组根据感兴趣的分位点的数量。通常是一个矢量,这样的反应<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M260">
N米米l:米我><米米l:mi>
O米米l:米我><米米l:mi>
年代米米l:米我><米米l:mo>
=米米l:米o><米米l:米n mathvariant="normal">
1米米l:米n>米米l:math>
和<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M261">
l米米l:米我><米米l:mo>
=米米l:米o><米米l:米i>
米米米l:米我>米米l:math>
。
我们使用仿真生成<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M268">
n米米l:米我><米米l:mo>
=米米l:米o><米米l:米n mathvariant="normal">
200年米米l:米n>米米l:math>
从麦克斯韦分布随机变量的速度参数是一个协变量的函数<我nl我ne- - - - - -for米ul一个><米ml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M269">
x米米l:米我>米米l:mrow>
2米米l:米n>米米l:mrow>
。考虑到非线性趋势数据集,添加剂与立方模型在线性平滑样条函数似乎是一个更好的选择方案,诸如VGLMs。在这个例子中,我们执行以下步骤确认方法的性能。(1)
生成随机偏离麦克斯韦分布。
(2)运行条件VGAM-quantile建模使用<米ono年代p一个ce>米一个xwellQlink ()米ono年代p一个ce>基于VGAM家庭功能<米ono年代p一个ce>VGAM:麦克斯韦()米ono年代p一个ce>费舍尔得分,估计麦克斯韦分布。