search for: nfinal

Displaying 6 results from an estimated 6 matches for "nfinal".

Did you mean: final
2006 Jan 06
1
lmer p-vales are sometimes too small
...: ", sum(p.aov<0.001)/n.sims, "\n") cat("lmer.1: ", sum(p.lmer.1<0.001)/n.sims, "\n") cat("lmer.2: ", sum(p.lmer.2<0.001)/n.sims, "\n") cat("lmer.3: ", sum(p.lmer.3<0.001)/n.sims, "\n") cat("\nFinal aov analysis: \n") print(summary(fm.aov)$"Error: Block:Treat") cat("\nFinal lmer analysis 1: \n") print(anova(fm.lmer.1)) cat("\nFinal lmer analysis 2: \n") print(anova(fm.lmer.2)) cat("\nFinal lmer analysis 3: \n") print(anova(fm.lm...
2002 Jul 18
1
sem: incorrect parameter estimates
...;, 1e6, 'D->B', 'p4', 0.2, 'C<>C', 'p5', 1e6 ), ncol=3, byrow=T) p.sem <- sem(p.ram, p.cov, 293, fixed.x=c('A','D'), debug=T, analytic.gradient=F) print(summary(p.sem)) cat('\n\nFinal parameter values:\n') print(p.sem$coeff,digits=2) cat('\nExpected final values (from AMOS software):\n') cat('5.0e-01 8.4e-01 7.9e+05 4.3e-02 1.3e+06\n\n') cat('\n\nThe chisquare value from the AMOS fit is 69.7.\n') ====== end R code ======= -- Baxter Rogers (rogers...
2010 Nov 03
0
bad optimization with nnet?
...matrix(rnorm(10),10,1) rx <- max(x)-min(x) rz <- max(z)-min(z) x <- x/rx z <- z/rz erreur <- 10^9 for(i in 1:100){ temp.mod <- nnet(x=x,y=z,size=10,rang=1,maxit=1000) if(temp.mod$value<erreur){ res.mod <- temp.mod erreur <- res.mod$value } } cat("\nFinal error : ",res.mod$value,"\n") Normaly it is easy task for an MLP with 10 hidden units to reduce the final error to almost 0 (althougt there is nothing to predict). But the smallest error that I get is : 0.753895 (very poor result) Maybe this problem is already known?? Maybe th...
2011 Jul 20
0
Analysis of unbalanced data in nlme or car
I am analyzing a dataset on the effects of six pesticides on population growth rate of a predatory mite. The response variable is the population growth rate of the mite expressed as ln(Nfinal/Nstarting) of the mite, where N final the population of the mite at the end of the experiment and N starting the population of the mite at the beginning of the experiment. Each of the six treatments was ran in three time blocks with 2 replicates per block. In the third block I lost 1 replicate for...
2009 May 05
0
stepAICc function (based on MASS:::stepAIC.default)
...") ddf <- c(NA, abs(diff(rdf))) AICc <- sapply(models, "[[", "AICc") heading <- c("Stepwise Model Path \nAnalysis of Deviance Table", "\nInitial Model:", deparse(as.vector(formula(object))), "\nFinal Model:", deparse(as.vector(formula(fit))), "\n") aod <- if (usingCp) data.frame(Step = change, Df = ddf, Deviance = dd, `Resid. Df` = rdf, `Resid. Dev` = rd, Cp = AICc, check.names = FALSE) else data.frame(Ste...
2002 Aug 07
2
Constructing titles from list of expressions
...9 7 heading 1040 8 1033 47 Stepwise\040Model\040Path\040\nAnalysis\040of\040Deviance\040Table 1033 15 \nInitial\040Model: 1033 61 gwRZ\040~\040I(f^2)\040+\040I(vS^2)\040+\040I(o^2)\040+\040f\040+\040vS\040+\040 o\040+\040f:vS\040+\040f:o\040+\040 1033 17 \040\040\040\040vS:o\040+\040f:vS:o 1033 13 \nFinal\040Model: 1033 61 gwRZ\040~\040I(f^2)\040+\040I(vS^2)\040+\040I(o^2)\040+\040f\040+\040vS\040+\040 o\040+\040f:vS\040+\040f:o\040+\040 1033 8 \040\040\040\040vS:o 1033 1 \n 254 1026 767 1040 14 1033 12 coefficients 1033 9 residuals 1033 7 effects 1033 4 rank 1033 13 fitted.values 1033 6 assign 1033...