Displaying 2 results from an estimated 2 matches for "obeta".
Did you mean:
beta
2011 Jul 20
2
Bootstrap
...17 4.74
14 5.49
11 4.12
17 5.92"), header=TRUE)
data <- as.matrix(xx)
closeAllconnections()
Nt <- NULL
for (Ncount in 1:100)
{
y <- data[,1]
x <- data[,2]
n <- length(x)
X <- cbind(rep(1,n),x) #covariate/design matrix
obeta<- c(1,1) #previous/starting values of beta
nbeta <- c(0,0) #new beta
iter=0
while(crossprod(obeta-nbeta)>10^(-12))
{
nbeta <- obeta
eta <- X%*%nbeta
mu <- eta
mu1 <- 1/eta
W <- diag(...
2003 May 21
1
Code Help
...(clover$Line),var))
om <- order(mn.tab)
plot(mn.tab,sd.tab,las=1,xlab="mean",ylab="s.d.")
lines(lowess(mn.tab,sd.tab))
model.i <- glm(Y ~ C(clover$Line,Cline),data=clover,family=gaussian),
Beta <- summary(model.i)$coefficients
new.order <- order(Beta[-1,1]) + 1
Obeta <- Beta[c(1,new.order),]
detect.diff <- quantile(Obeta,probs=seq(0.2,0.8,0.2))
std.dev <- (min(Obeta[-1,2]) + max(Obeta[-1,2]))/2*sqrt(5)
sink("ClovAttr.txt")
print(paste("_____________",vnames[(i+2)],"_____________"))
cat("__________ mean table _______...