Nothing
fit.saery.indep <-
function(X, ydi, D, md, sigma2edi, conf.level){
fit <- REML.saery.indep(X, ydi, D, md, sigma2edi)
sigmau1.hat <- fit[[1]][1]
sigmau2.hat <- fit[[1]][2]
stderr <- stderr.saery.indep(fit)
alpha <- 1-conf.level
k <- 1-alpha/2
z <- qnorm(k)
lower.s <- c(sigmau1.hat-z*stderr[[2]][1], sigmau2.hat-z*stderr[[2]][2])
upper.s <- c(sigmau1.hat+z*stderr[[2]][1], sigmau2.hat+z*stderr[[2]][2])
sigmas <- data.frame(Estimate=c(sigmau1.hat, sigmau2.hat), lower=lower.s, upper=upper.s, row.names= c("sigmau1", "sigmau2"))
beta.u.hat <- BETA.U.saery.indep(X, ydi, D, md, sigma2edi, sigmau1.hat, sigmau2.hat)
pv <- pvalue(beta.u.hat[[1]], fit)
lower.b <- beta.u.hat[[1]]-z*stderr[[1]]
upper.b <- beta.u.hat[[1]]+z*stderr[[1]]
betas <- data.frame(Estimate=beta.u.hat[[1]], std.err=stderr[[1]], t=beta.u.hat[[1]]/stderr[[1]], 'p-value'=pv, lower=lower.b, upper=upper.b)
return(list("Regression"=betas, "SIGMA"=sigmas, Fsig=fit[[2]], iteration=fit[[3]]))
}
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.