search for: b_hat

...important if I want to do things like u' * A * u, where u' is a row vector and u is a column vector, regardless of what originall format the "u" is... I want to recast it to column vector or row vector... How can I do that? ---------------------------------------------------- b_hat=solve(t(X) %*% X) %*% t(X) %*% Y; attr(b_hat, "dim") [1] 4 1 rbind(b_hat) [,1] -4.814763 V -5.804245 -5.122668 -4.308326 [[alternative HTML version deleted]]

Hi, I try to execute the seven lines of code below to plot a graph. But I am failing as the messages below show. Where am I going wrong? install.packages("rgl") library(rgl) y_hat = X%*%B_hat open3d(windowRect = c(100,100,900,900),family = "serif") color = rainbow(length(y_hat))[rank(y_hat)] plot3d(educ,exper,wage,col = color,type = "s",size = 0.5,xlim = c(0,20),ylim = c(0,60),zlim = c(-10,70),box = FALSE,axes = TRUE) planes3d(B_hat[2],B_hat[3],-1,B_hat[1],alpha = 0....

In a multiple linear regression with two independent variables is there any function in R to test for the coefficients being different than some given values? Example: x1<-rnorm(100) x2<-rnorm(100) y<-3+0.6*x1+0.3*x2 lm(y~x1+x2) Obtain a test for the coefficients for x1 being different than 0.6 and for x2 different than 0.3 Thanks Manuel

...arameters. a_guess<-2 b_guess<-3 q_guess<-7 q_guess-1/a_guess p0<-c(log(a_guess), log(b_guess), log(q_guess)) #Optimize to find minimum of the negative likelihood function est_p_A<-optim(par=p0, fn=neg.LL_sinmad, y=y) est_p_A$par a_hat1<-exp(est_p_A$par[1]) b_hat1<-exp(est_p_A$par[2]) q_hat1<-exp(est_p_A$par[3]) q_hat1-1/a_hat1 est_p_B<-nlm(f=neg.LL_sinmad, p=c(est_p_A$par[1], est_p_A$par[2], est_p_A$par[3]), print.level=1, y=y) a_hat<-exp(est_p_A$par[1]) b_hat<-exp(est_p_A$par[2]) q_hat<-exp(est_p_A$par[3]) q_h...

...l guess of the parameters. a_guess<-2 b_guess<-3 q_guess<-7 q_guess-1/a_guess p0<-c(log(a_guess), log(b_guess), log(q_guess)) #Optimize to find minimum of the negative likelihood function est_p_A<-optim(par=p0, fn=neg.LL_sinmad, y=y) est_p_A$par a_hat1<-exp(est_p_A$par[1]) b_hat1<-exp(est_p_A$par[2]) q_hat1<-exp(est_p_A$par[3]) q_hat1-1/a_hat1 est_p_B<-nlm(f=neg.LL_sinmad, p=c(est_p_A$par[1], est_p_A$par[2], est_p_A$par[3]), print.level=1, y=y) a_hat<-exp(est_p_A$par[1]) b_hat<-exp(est_p_A$par[2]) q_hat<-exp(est_p_A$par[3]) q_hat-1/a_hat CONFIDENTIA...