similar to: Details related to Levene's test

Dear All, I sent already before a message concerning Levene's test in nested model, but I didn't get any answer. Optimistically I hope to get an answer this time. I also point a new question related to the whole model, because I haven't find any sure answer if I am analysing it in a suitable way or not. I really have tried to do my homework. I have response variable (y) and four

Hi, All! To calculate Levene's Test for Homogeneity of Variance I use R Commander, and this is the output: > levene.test(Dataset$age, Dataset$sex) Levene's Test for Homogeneity of Variance Df F value Pr(>F) group 1 0.8739 0.3567 33 I am not sure what means "Pr(>F)"? Can anyone explain/translate this? Regards, Iurie Malai Department of Psychology and

> From: Peter Dalgaard BSA <p.dalgaard at biostat.ku.dk> > Date: 01 Sep 2000 09:54:59 +0200 > > Prof Brian D Ripley <ripley at stats.ox.ac.uk> writes: Important omission: specification from Murray Jorgensen The test that I was thinking of basically does an anova on a modified response variable that is the absolute value of the difference between an observation

Hello all, I am low down on the learning curve of R but so far I have had little trouble using most of the packages. However, recently I have run into a wall when it comes to bootstrapping a Levene's test (from the car package) and thought you might be able to help. I have not been able to find R examples for the "boot" package where the test statistic specifically uses a

Dear All 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 4 factors plus response. Could someone please give me an example about how to use the command levene.test(y, group) ? Thanks in advance, Paul

Hello, I notice that when I do Levene's test to test equality of variances across levels of a factor, I get different answers in R and SPSS 16. e.g.: For the chickwts data, in R, levene.test(weight, feed) gives F=0.7493, p=0.5896. SPSS 16 gives F=0.987, p=0.432 Why this difference? Which one should I believe? (I would like to believe R :) Ravi -- View this message in context:

hello, I try to run simulation of levene's test to find the p-value but the error of replacement has length zero occur, could anyone help me to fix this problem? asim <- 1000 pv<-rep(NA,asim) for(i in 1:asim) {print(i) set.seed(i) g1 <- rnorm(20,0,2) g2 <- rnorm(20,0,2) g3 <- rnorm(20,0,2) x <- c(g1,g2,g3) group<-as.factor(c(rep(1,20),rep(2,20),rep(3,20))) library(Rcmdr)

Dear Contributors, I am asking help on the way how to solve a problem related to loops for that I always get confused with. I would like to perform the following procedure in a compact way. Consider that p is a matrix composed of 100 rows and three columns. I need to calculate the sum over some rows of each column separately, as follows: fa1<-(colSums(p[1:25,])) fa2<-(colSums(p[26:50,]))

Dear contributors I have tried this experiment: x<-c() for (i in 1:12){ x[i]<-list(cbind(x1[i],x2[i])) #this is a list of 12 couples of time series I am using to perform a test } # that compares them 2 by 2 # ################# #trace statistic test<-data.frame() cval<-array( , dim=c(2,3,12)) for (i in 2:12){ for (k in 1:2){ for (j in 1:3){ result[k,j,i]<-

Dear Contributors, I am back asking for help concerning the same type of dataset I was asking before, in a previous help request. I needed to sum data over subsample of three time series each of them made of 100 observations. The solution proposed were various, among which: db<-p dim( db ) <- c(25,4,3) db2 <- apply(db, c(2,3), sum) db3 <- t(apply(db2, 1, function(poff)

En innebygd og tegnsett-uspesifisert tekst ble skilt ut... Navn: ikke tilgjengelig Nettadresse: https://stat.ethz.ch/pipermail/r-help/attachments/20080214/5a98c799/attachment.pl

Please can you help me I have the following problem: I have selected an lm model through the step procedure which visualize for each step the AIC value; then I have calculated for the initial model and the selected one the AIC using the funnction AIC. The results are different.What's happened? Emilia Rocco Dipartimento di Statistica "G. Parenti" Università di Firenze e-mail:

Dear all, I need to compute tensor product of B-spline defined over equi-spaced break-points. I wrote my own program (it works in a 2-dimensional setting) library(splines) # set the break-points Knots = seq(-1,1,length=10) # number of splines M = (length(Knots)-4)^2 # short cut to splineDesign function bspline = function(x) splineDesign(Knots,x,outer.ok = T) # bivariate tensor product of

Dear R-users, I learned today that there exists an interesting topic in numerical analysis names "best polynomial approximation" (BSA). Given a function f the BSA of degree k, say pk, is the polynomial such that pk=arginf sup(|f-pk|) Although given some regularity condition of f, pk is unique, pk IS NOT calculated with least square. A quick google tour show a rich field of research

Hello, I am trying to plot a few correlograms on the same figure, with the function corrgram() from the package corrgram. However, the function does not seem to use the base graphic system, as setting out the multiple figure layout with, e.g., par(mfrow=c(2, 2,)) does not work. Does anybody know a workaround for this? Many thanks in advance for any advice best giuseppe -- Giuseppe Pagnoni,

Dear all, I am currently facing a problem related to the spatial autocorrelation of a sample of stations; these stations supply weekly data for a fixed time-window during the year (namely, 4-6 months per year). For this reason I'm trying to use the R package 'spdep' (specifically Moran's I) in order to get rid of it. Does anyone know how is it possible (if it is...) to

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 command levene.test(y, group) Thanks in advance, marta -- View this message in

Hello, I have two error messages on user.log as follow : 1. [2001/11/26 14:41:07, 0] smbd/oplock.c:oplock_break(813) oplock_break: no break received from client within 30 seconds. oplock_break failed for file aa1/EMI/SIGVM00 (dev = 302, inode = 194716). [2001/11/26 14:41:07, 0] smbd/oplock.c:oplock_break(859) oplock_break: client failure in oplock break in file aa1/EMI/SIGVM00 2.

Dear all I would like to show my audience that some variables are homogenous inside groups but different outside. I can use by with summary for all variables by(iris[,1:4], iris$Species, summary) what can be quite messy in case of more than few variables and about 8 groups or densityplot for one variable densityplot(~Petal.Length | Species, iris) I have two questions: 1. Is there

Dear R-helpers, I'm having a problem in using plot.design in Design Library. Tho following example code produce the error: > n <- 1000 # define sample size > set.seed(17) # so can reproduce the results > age <- rnorm(n, 50, 10) > blood.pressure <- rnorm(n, 120, 15) > cholesterol <- rnorm(n, 200, 25) > sex <-