Displaying 1 result from an estimated 1 matches for "stderrr1".
Did you mean:
stderrr2
2005 Nov 07
1
Newbie on functions
...ction (m, cov, Q, R) {
theta <- (acos(R/sqrt(Q^3)))
beta <- (-2)*sqrt(Q)*cos(theta/3)+m[1]/3
rho1 <- (-2)*sqrt(Q)*cos((theta+2*pi)/3)+m[1]/3
rho2 <- (-2)*sqrt(Q)*cos((theta-2*pi)/3)+m[1]/3
stderrb <- deltamethod( ~(-2)*sqrt(Q)*cos(theta/3)+x1/3,m,cov)
stderrr1 <- deltamethod( ~(-2)*sqrt(Q)*cos((theta+2*pi)/3)+x1/3, m,
cov) stderrr2 <- deltamethod( ~(-2)*sqrt(Q)*cos((theta-2*pi)/3)+x1/3,
m, cov) stderr <- c(stderrb,stderrr1,stderrr2)
results <- c(beta,rho1,rho2,stderr)
results2 <- t(results)
results2
}
When I call the funct...