search for: zvalu

Suppose I create an R program called myTest.R with only one line like the following: type <- as.integer(readline("input type (1: type1; 2: type2)? ")) Then I'd like to run myTest.R in batch mode by constructing an input file called answers.R with the following: source("myTest.R") 1 When I ran the following at the terminal: R CMD BATCH answer.R output.Rout it failed

...<- summary(model1)$tTable mod2.out <- summary(model2)$tTable mod3.out <- summary(model3)$tTable indir <- mod3.out[2, 1] * mod2.out[3, 1] effvar <- (mod3.out[2, 1])^2 * (mod2.out[3, 2])^2 + (mod2.out[3, 1])^2 * (mod3.out[2, 2])^2 serr <- sqrt(effvar) zvalue = indir/serr out <- list(Model.1 = mod1.out, Model.2 = mod2.out, Model.3 = mod3.out, Indirect.Effect = indir, SE = serr, z.value = zvalue, N = nrow(NEWDAT),n.grps=length(unique(NEWDAT$grpid))) return(out) } Any help much appreciated! Best, Gavril -- Gavril Bilev, P...

...create a grid of circles with the right number of dimensions: Read_data=read.table("C:/files/circle_data.txt", head = T) number_rows = nrow(Read_data) number_cols = ncol(Read_data) xaxis = rep(seq(1,ncol(Read_data),1),nrow(Read_data)) yaxis = rep(seq(1,nrow(Read_data),1),ncol(Read_data)) zvalues = rep(0.1,(number_rows*number_cols)) #This is just a dummy vector symbols(xaxis, yaxis, circles = zvalues, inches = FALSE, xlab="", ylab="") axis(2, at=c(1:number_rows), rownames(Read_data), las = 1) axis(1, at=c(1:number_cols), colnames(Read_data)) BUT... this only uses dumm...

Dear R users, Recent changes to the MuMIn package now means that the model averaging command (model.avg) no longer returns confidence intervals, but instead returns zvalues and corresponding pvalues for fixed effects included in models. Previously I have used this package for model selection/averaging following Greuber et al (2011) where it suggests that one should use confidence intervals from model averaging to assess whether your fixed effects have an affect or...

Dear all, how can I use vcovHC() to get robust/corrected standard errors from an arima() output? I ran an arima model with AR(1) and got the estimate, se, zvalue and p-value using coeftest(arima.output). However, I cannot use vcovHC(arima.output) to get corrected standard errors. It seems vcovHC works only with lm and plm objects? Is there another way I can get robust/corrected standard errors, or am I missing something? Thank you! Nicole -...