Nothing
TMLjac21.G <-
function(d.Beta,d.sigma,rs0,delta,X,cu) {
# Jacobian of sigma-equation wrt Beta.hat
n <- length(rs0); cl <- -cu; zero <- 1.e-6
D1 <- D2 <- D3 <- D <- rep(0,n); p <- ncol(X)
rsd <- (rs0-X%*%d.Beta)/d.sigma
Fo <- pnorm(rsd)
den <- 1-Fo; ok <- den > zero
fo <- dnorm(rsd)
ai <- (pmax(rs0,cl)-X%*%d.Beta)/d.sigma
bi <- (cu - X%*%d.Beta)/d.sigma
foai <- dnorm(ai)
fobi <- dnorm(bi)
Foai <- pnorm(ai)
Fobi <- pnorm(bi)
fopai <- foai*(-ai)
fopbi <- fobi*(-bi)
D1 <- - delta*ww(rs0,cl,cu)*2*rsd
D2[ok]<- - ( (1-delta)*fo/(1-Fo)^2*(foai*ai-fobi*bi + Fobi - Foai) )[ok]
D3[ok]<- - ( (1-delta)/(1-Fo)*( fopai*ai - fopbi*bi ) )[ok]
D <- D1 + D2 + D3
Jac <- t(X)%*%(as.vector(D))/d.sigma/(n-p)
Jac}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.