[R] levene.test

2010-07-14 Thread martanair
I am trying to use Levene's test (of package car), but I do not understand quite well how to use it. '?levene.test' does not unfortunately provide any example. My data are in a data frame and correspond to 1 factor plus response. Could someone please give me an example about how to use the com

[R] levene.test

2010-07-14 Thread martanair
I am trying to use Levene's test (of package car), but I do not understand quite well how to use it. '?levene.test' does not unfortunately provide any example. My data are in a data frame and correspond to 1 factor plus response. Could someone please give me an example about how to use the com

Re: [R] levene.test

2010-07-14 Thread martanair
the value of the F-test is uncorrect. If I calculate the Levene's Test for Homogeneity of Variance with SPSS I get F = 5.128 Pr (> F) = 0.033 Why? what is the exact formula which R calculate Levene's test for two groups? Thanks marta -- View this message in context: http://r.789695.n4.nabble.com

Re: [R] levene.test

2010-07-14 Thread martanair
how do I calculate the differences from the mean? levene.test(filedati$y , filedati$fertilizzante, location="mean") Is it right? regards marta -- View this message in context: http://r.789695.n4.nabble.com/levene-test-tp2288632p2288832.html Sent from the R help mailing list archive at Nabble.c

[R] Multilevel

2011-01-29 Thread martanair
Why this? I write: radon.data <- list ("n", "J", "x", "y", "county") radon.inits <- function (){ list (a=rnorm(J), b=rnorm(1), mu.a=rnorm(1), sigma.y=runif(1), sigma.a=runif(1)) } radon.parameters <- c ("a", "b", "mu.a", "sigma.y", "sigma.a") # with 10 iterations radon.1 <- bugs (rad

Re: [R] package R2WinBUGS question; was: Multilevel

2011-01-29 Thread martanair
I have R2WinBUGS and this is my radon1.bug # Bugs code for multilevel model for radon # with bsmt as an individual predictor # varying-intercept model model { for (i in 1:n){ y[i] ~ dnorm (y.hat[i], tau.y) y.hat[i] <- a[county[i]] + b*x[i] } b ~ dnorm (0, .0001) tau.y <- pow(si

[R] greedy.wilks

2012-06-20 Thread martanair
I have used the greedy.wilks to stepwise discriminant analysis, but it doesn't work with my dataset. I don't understand which is the problem Could you help me please? Thanks in advance Marta > str(data_indiciN2) 'data.frame': 200 obs. of 36 variables: $ gruppo: Factor w/ 2 levels "0","1": 2

[R] Stepwise Discriminant Analysis - greedy.wilks

2012-06-20 Thread martanair
I don't understand which is the problem Could you help me please? Thanks in advance Marta > str(data_indiciN2) 'data.frame': 200 obs. of 36 variables: $ gruppo: Factor w/ 2 levels "0","1": 2 2 2 2 2 2 2 2 2 2 ... $ I001N2: num 19.32 8.22 28.35 7.24 14.7 ... $ I002N2: num 2.92 2.5