psins.r

r软件 另一款可以计算核估计的软件包 需安装r软件

psi <-
-3*(2*a15^2*a35 + 2*a13*a15*a55 + a11*a35*a55)
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==0) 
& (r[5]==2) & (r[6]==1))
psi <-
-4*a15*a16*a35 - 2*a15^2*a36 - 2*a13*a16*a55 - a11*a36*a55 - 4*a13*a15*a56 - 
 2*a11*a35*a56
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==0) 
& (r[5]==1) & (r[6]==2))
psi <-
-2*a16^2*a35 - 4*a15*a16*a36 - 4*a13*a16*a56 - 2*a11*a36*a56 - 
 2*a13*a15*a66 - a11*a35*a66
if ((r[1]==2) & (r[2]==0) & (r[3]==1) & (r[4]==0) 
& (r[5]==0) & (r[6]==3))
psi <-
-3*(2*a16^2*a36 + 2*a13*a16*a66 + a11*a36*a66)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==4) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*a44*(4*a14^2 + a11*a44)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==3) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*(2*a14*a15*a44 + 2*a14^2*a45 + a11*a44*a45)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==3) 
& (r[5]==0) & (r[6]==1))
psi <-
-3*(2*a14*a16*a44 + 2*a14^2*a46 + a11*a44*a46)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==2) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15^2*a44 - 8*a14*a15*a45 - 2*a11*a45^2 - 2*a14^2*a55 - a11*a44*a55
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==2) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a15*a16*a44 - 4*a14*a16*a45 - 4*a14*a15*a46 - 2*a11*a45*a46 - 
 2*a14^2*a56 - a11*a44*a56
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==2) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16^2*a44 - 8*a14*a16*a46 - 2*a11*a46^2 - 2*a14^2*a66 - a11*a44*a66
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==1) 
& (r[5]==3) & (r[6]==0))
psi <-
-3*(2*a15^2*a45 + 2*a14*a15*a55 + a11*a45*a55)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==1) 
& (r[5]==2) & (r[6]==1))
psi <-
-4*a15*a16*a45 - 2*a15^2*a46 - 2*a14*a16*a55 - a11*a46*a55 - 4*a14*a15*a56 - 
 2*a11*a45*a56
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==1) 
& (r[5]==1) & (r[6]==2))
psi <-
-2*a16^2*a45 - 4*a15*a16*a46 - 4*a14*a16*a56 - 2*a11*a46*a56 - 
 2*a14*a15*a66 - a11*a45*a66
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==1) 
& (r[5]==0) & (r[6]==3))
psi <-
-3*(2*a16^2*a46 + 2*a14*a16*a66 + a11*a46*a66)
  if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==0) 
& (r[5]==4) & (r[6]==0))
psi <-
-3*a55*(4*a15^2 + a11*a55)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==0) 
& (r[5]==3) & (r[6]==1))
psi <-
-3*(2*a15*a16*a55 + 2*a15^2*a56 + a11*a55*a56)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==0) 
& (r[5]==2) & (r[6]==2))
psi <-
-2*a16^2*a55 - 8*a15*a16*a56 - 2*a11*a56^2 - 2*a15^2*a66 - a11*a55*a66
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==0) 
& (r[5]==1) & (r[6]==3))
psi <-
-3*(2*a16^2*a56 + 2*a15*a16*a66 + a11*a56*a66)
if ((r[1]==2) & (r[2]==0) & (r[3]==0) & (r[4]==0) 
& (r[5]==0) & (r[6]==4))
psi <-
-3*a66*(4*a16^2 + a11*a66)
if ((r[1]==1) & (r[2]==5) & (r[3]==0) & (r[4]==0) 
& (r[5]==0) & (r[6]==0))
psi <-
-15*a12*a22^2
if ((r[1]==1) & (r[2]==4) & (r[3]==1) & (r[4]==0) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*a22*(a13*a22 + 4*a12*a23)
if ((r[1]==1) & (r[2]==4) & (r[3]==0) & (r[4]==1) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*a22*(a14*a22 + 4*a12*a24)
if ((r[1]==1) & (r[2]==4) & (r[3]==0) & (r[4]==0) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*a22*(a15*a22 + 4*a12*a25)
if ((r[1]==1) & (r[2]==4) & (r[3]==0) & (r[4]==0) 
& (r[5]==0) & (r[6]==1))
psi <-
-3*a22*(a16*a22 + 4*a12*a26)
if ((r[1]==1) & (r[2]==3) & (r[3]==2) & (r[4]==0) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a13*a22*a23 + 2*a12*a23^2 + a12*a22*a33)
if ((r[1]==1) & (r[2]==3) & (r[3]==1) & (r[4]==1) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(a14*a22*a23 + a13*a22*a24 + 2*a12*a23*a24 + a12*a22*a34)
if ((r[1]==1) & (r[2]==3) & (r[3]==1) & (r[4]==0) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*(a15*a22*a23 + a13*a22*a25 + 2*a12*a23*a25 + a12*a22*a35)
if ((r[1]==1) & (r[2]==3) & (r[3]==1) & (r[4]==0) 
& (r[5]==0) & (r[6]==1))
psi <-
-3*(a16*a22*a23 + a13*a22*a26 + 2*a12*a23*a26 + a12*a22*a36)
if ((r[1]==1) & (r[2]==3) & (r[3]==0) & (r[4]==2) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a14*a22*a24 + 2*a12*a24^2 + a12*a22*a44)
if ((r[1]==1) & (r[2]==3) & (r[3]==0) & (r[4]==1) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*(a15*a22*a24 + a14*a22*a25 + 2*a12*a24*a25 + a12*a22*a45)
if ((r[1]==1) & (r[2]==3) & (r[3]==0) & (r[4]==1) 
& (r[5]==0) & (r[6]==1))
psi <-
-3*(a16*a22*a24 + a14*a22*a26 + 2*a12*a24*a26 + a12*a22*a46)
if ((r[1]==1) & (r[2]==3) & (r[3]==0) & (r[4]==0) 
& (r[5]==2) & (r[6]==0))
psi <-
-3*(2*a15*a22*a25 + 2*a12*a25^2 + a12*a22*a55)
if ((r[1]==1) & (r[2]==3) & (r[3]==0) & (r[4]==0) 
& (r[5]==1) & (r[6]==1))
psi <-
-3*(a16*a22*a25 + a15*a22*a26 + 2*a12*a25*a26 + a12*a22*a56)
if ((r[1]==1) & (r[2]==3) & (r[3]==0) & (r[4]==0) 
& (r[5]==0) & (r[6]==2))
psi <-
-3*(2*a16*a22*a26 + 2*a12*a26^2 + a12*a22*a66)
if ((r[1]==1) & (r[2]==2) & (r[3]==3) & (r[4]==0) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a13*a23^2 + a13*a22*a33 + 2*a12*a23*a33)
if ((r[1]==1) & (r[2]==2) & (r[3]==2) & (r[4]==1) 
& (r[5]==0) & (r[6]==0))
psi <-
-2*a14*a23^2 - 4*a13*a23*a24 - a14*a22*a33 - 2*a12*a24*a33 - 2*a13*a22*a34 - 
 4*a12*a23*a34
if ((r[1]==1) & (r[2]==2) & (r[3]==2) & (r[4]==0) 
& (r[5]==1) & (r[6]==0))
psi <-
-2*a15*a23^2 - 4*a13*a23*a25 - a15*a22*a33 - 2*a12*a25*a33 - 2*a13*a22*a35 - 
 4*a12*a23*a35
if ((r[1]==1) & (r[2]==2) & (r[3]==2) & (r[4]==0) 
& (r[5]==0) & (r[6]==1))
psi <-
-2*a16*a23^2 - 4*a13*a23*a26 - a16*a22*a33 - 2*a12*a26*a33 - 2*a13*a22*a36 - 
 4*a12*a23*a36
if ((r[1]==1) & (r[2]==2) & (r[3]==1) & (r[4]==2) 
& (r[5]==0) & (r[6]==0))
psi <-
-4*a14*a23*a24 - 2*a13*a24^2 - 2*a14*a22*a34 - 4*a12*a24*a34 - a13*a22*a44 - 
 2*a12*a23*a44
if ((r[1]==1) & (r[2]==2) & (r[3]==1) & (r[4]==1) 
& (r[5]==1) & (r[6]==0))
psi <-
-2*a15*a23*a24 - 2*a14*a23*a25 - 2*a13*a24*a25 - a15*a22*a34 - 
 2*a12*a25*a34 - a14*a22*a35 - 2*a12*a24*a35 - a13*a22*a45 - 2*a12*a23*a45
if ((r[1]==1) & (r[2]==2) & (r[3]==1) & (r[4]==1) 
& (r[5]==0) & (r[6]==1))
psi <-
-2*a16*a23*a24 - 2*a14*a23*a26 - 2*a13*a24*a26 - a16*a22*a34 - 
 2*a12*a26*a34 - a14*a22*a36 - 2*a12*a24*a36 - a13*a22*a46 - 2*a12*a23*a46
if ((r[1]==1) & (r[2]==2) & (r[3]==1) & (r[4]==0) 
& (r[5]==2) & (r[6]==0))
psi <-
-4*a15*a23*a25 - 2*a13*a25^2 - 2*a15*a22*a35 - 4*a12*a25*a35 - a13*a22*a55 - 
 2*a12*a23*a55
if ((r[1]==1) & (r[2]==2) & (r[3]==1) & (r[4]==0) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a16*a23*a25 - 2*a15*a23*a26 - 2*a13*a25*a26 - a16*a22*a35 - 
 2*a12*a26*a35 - a15*a22*a36 - 2*a12*a25*a36 - a13*a22*a56 - 2*a12*a23*a56
if ((r[1]==1) & (r[2]==2) & (r[3]==1) & (r[4]==0) 
& (r[5]==0) & (r[6]==2))
psi <-
-4*a16*a23*a26 - 2*a13*a26^2 - 2*a16*a22*a36 - 4*a12*a26*a36 - a13*a22*a66 - 
 2*a12*a23*a66
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==3) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a14*a24^2 + a14*a22*a44 + 2*a12*a24*a44)
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==2) 
& (r[5]==1) & (r[6]==0))
psi <-
-2*a15*a24^2 - 4*a14*a24*a25 - a15*a22*a44 - 2*a12*a25*a44 - 2*a14*a22*a45 - 
 4*a12*a24*a45
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==2) 
& (r[5]==0) & (r[6]==1))
psi <-
-2*a16*a24^2 - 4*a14*a24*a26 - a16*a22*a44 - 2*a12*a26*a44 - 2*a14*a22*a46 - 
 4*a12*a24*a46
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==1) 
& (r[5]==2) & (r[6]==0))
psi <-
-4*a15*a24*a25 - 2*a14*a25^2 - 2*a15*a22*a45 - 4*a12*a25*a45 - a14*a22*a55 - 
 2*a12*a24*a55
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==1) 
& (r[5]==1) & (r[6]==1))
psi <-
-2*a16*a24*a25 - 2*a15*a24*a26 - 2*a14*a25*a26 - a16*a22*a45 - 
 2*a12*a26*a45 - a15*a22*a46 - 2*a12*a25*a46 - a14*a22*a56 - 2*a12*a24*a56
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==1) 
& (r[5]==0) & (r[6]==2))
psi <-
-4*a16*a24*a26 - 2*a14*a26^2 - 2*a16*a22*a46 - 4*a12*a26*a46 - a14*a22*a66 - 
 2*a12*a24*a66
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==0) 
& (r[5]==3) & (r[6]==0))
psi <-
-3*(2*a15*a25^2 + a15*a22*a55 + 2*a12*a25*a55)
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==0) 
& (r[5]==2) & (r[6]==1))
psi <-
-2*a16*a25^2 - 4*a15*a25*a26 - a16*a22*a55 - 2*a12*a26*a55 - 2*a15*a22*a56 - 
 4*a12*a25*a56
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==0) 
& (r[5]==1) & (r[6]==2))
psi <-
-4*a16*a25*a26 - 2*a15*a26^2 - 2*a16*a22*a56 - 4*a12*a26*a56 - a15*a22*a66 - 
 2*a12*a25*a66
if ((r[1]==1) & (r[2]==2) & (r[3]==0) & (r[4]==0) 
& (r[5]==0) & (r[6]==3))
psi <-
-3*(2*a16*a26^2 + a16*a22*a66 + 2*a12*a26*a66)
if ((r[1]==1) & (r[2]==1) & (r[3]==4) & (r[4]==0) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*a33*(4*a13*a23 + a12*a33)
if ((r[1]==1) & (r[2]==1) & (r[3]==3) & (r[4]==1) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(a14*a23*a33 + a13*a24*a33 + 2*a13*a23*a34 + a12*a33*a34)
if ((r[1]==1) & (r[2]==1) & (r[3]==3) & (r[4]==0) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*(a15*a23*a33 + a13*a25*a33 + 2*a13*a23*a35 + a12*a33*a35)
if ((r[1]==1) & (r[2]==1) & (r[3]==3) & (r[4]==0) 
& (r[5]==0) & (r[6]==1))
psi <-
-3*(a16*a23*a33 + a13*a26*a33 + 2*a13*a23*a36 + a12*a33*a36)
if ((r[1]==1) & (r[2]==1) & (r[3]==2) & (r[4]==2) 
& (r[5]==0) & (r[6]==0))
psi <-
-2*a14*a24*a33 - 4*a14*a23*a34 - 4*a13*a24*a34 - 2*a12*a34^2 - 
 2*a13*a23*a44 - a12*a33*a44
if ((r[1]==1) & (r[2]==1) & (r[3]==2) & (r[4]==1) 
& (r[5]==1) & (r[6]==0))
psi <-
-(a15*a24*a33) - a14*a25*a33 - 2*a15*a23*a34 - 2*a13*a25*a34 - 
 2*a14*a23*a35 - 2*a13*a24*a35 - 2*a12*a34*a35 - 2*a13*a23*a45 - a12*a33*a45
if ((r[1]==1) & (r[2]==1) & (r[3]==2) & (r[4]==1) 
& (r[5]==0) & (r[6]==1))
psi <-
-(a16*a24*a33) - a14*a26*a33 - 2*a16*a23*a34 - 2*a13*a26*a34 - 
 2*a14*a23*a36 - 2*a13*a24*a36 - 2*a12*a34*a36 - 2*a13*a23*a46 - a12*a33*a46
if ((r[1]==1) & (r[2]==1) & (r[3]==2) & (r[4]==0) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15*a25*a33 - 4*a15*a23*a35 - 4*a13*a25*a35 - 2*a12*a35^2 - 
 2*a13*a23*a55 - a12*a33*a55
if ((r[1]==1) & (r[2]==1) & (r[3]==2) & (r[4]==0) 
& (r[5]==1) & (r[6]==1))
psi <-
-(a16*a25*a33) - a15*a26*a33 - 2*a16*a23*a35 - 2*a13*a26*a35 - 
 2*a15*a23*a36 - 2*a13*a25*a36 - 2*a12*a35*a36 - 2*a13*a23*a56 - a12*a33*a56
if ((r[1]==1) & (r[2]==1) & (r[3]==2) & (r[4]==0) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16*a26*a33 - 4*a16*a23*a36 - 4*a13*a26*a36 - 2*a12*a36^2 - 
 2*a13*a23*a66 - a12*a33*a66
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==3) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*(2*a14*a24*a34 + a14*a23*a44 + a13*a24*a44 + a12*a34*a44)
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==2) 
& (r[5]==1) & (r[6]==0))
psi <-
-2*a15*a24*a34 - 2*a14*a25*a34 - 2*a14*a24*a35 - a15*a23*a44 - a13*a25*a44 - 
 a12*a35*a44 - 2*a14*a23*a45 - 2*a13*a24*a45 - 2*a12*a34*a45
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==2) 
& (r[5]==0) & (r[6]==1))
psi <-
-2*a16*a24*a34 - 2*a14*a26*a34 - 2*a14*a24*a36 - a16*a23*a44 - a13*a26*a44 - 
 a12*a36*a44 - 2*a14*a23*a46 - 2*a13*a24*a46 - 2*a12*a34*a46
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==1) 
& (r[5]==2) & (r[6]==0))
psi <-
-2*a15*a25*a34 - 2*a15*a24*a35 - 2*a14*a25*a35 - 2*a15*a23*a45 - 
 2*a13*a25*a45 - 2*a12*a35*a45 - a14*a23*a55 - a13*a24*a55 - a12*a34*a55
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==1) 
& (r[5]==1) & (r[6]==1))
psi <-
-(a16*a25*a34) - a15*a26*a34 - a16*a24*a35 - a14*a26*a35 - a15*a24*a36 - 
 a14*a25*a36 - a16*a23*a45 - a13*a26*a45 - a12*a36*a45 - a15*a23*a46 - 
 a13*a25*a46 - a12*a35*a46 - a14*a23*a56 - a13*a24*a56 - a12*a34*a56
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==1) 
& (r[5]==0) & (r[6]==2))
psi <-
-2*a16*a26*a34 - 2*a16*a24*a36 - 2*a14*a26*a36 - 2*a16*a23*a46 - 
 2*a13*a26*a46 - 2*a12*a36*a46 - a14*a23*a66 - a13*a24*a66 - a12*a34*a66
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==3) & (r[6]==0))
psi <-
-3*(2*a15*a25*a35 + a15*a23*a55 + a13*a25*a55 + a12*a35*a55)
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==2) & (r[6]==1))
psi <-
-2*a16*a25*a35 - 2*a15*a26*a35 - 2*a15*a25*a36 - a16*a23*a55 - a13*a26*a55 - 
 a12*a36*a55 - 2*a15*a23*a56 - 2*a13*a25*a56 - 2*a12*a35*a56
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==1) & (r[6]==2))
psi <-
-2*a16*a26*a35 - 2*a16*a25*a36 - 2*a15*a26*a36 - 2*a16*a23*a56 - 
 2*a13*a26*a56 - 2*a12*a36*a56 - a15*a23*a66 - a13*a25*a66 - a12*a35*a66
if ((r[1]==1) & (r[2]==1) & (r[3]==1) & (r[4]==0) 
& (r[5]==0) & (r[6]==3))
psi <-
-3*(2*a16*a26*a36 + a16*a23*a66 + a13*a26*a66 + a12*a36*a66)
if ((r[1]==1) & (r[2]==1) & (r[3]==0) & (r[4]==4) 
& (r[5]==0) & (r[6]==0))
psi <-
-3*a44*(4*a14*a24 + a12*a44)
if ((r[1]==1) & (r[2]==1) & (r[3]==0) & (r[4]==3) 
& (r[5]==1) & (r[6]==0))
psi <-
-3*(a15*a24*a44 + a14*a25*a44 + 2*a14*a24*a45 + a12*a44*a45)
if ((r[1]==1) & (r