similar to: R ANOVA gives diferent results than SPSS

An embedded and charset-unspecified text was scrubbed... Name: nem el?rhet? URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110420/4913654b/attachment.pl>

Dear R community, I am trying to resolve a discrepancy between the way SPSS and R handle 2-way, repeated measures ANOVA. An experiment was performed in which samples were drawn before and after treatment of four groups of subjects (control and disease states 1, 2 and 3). Each group contained five subjects. An experimental measurement was performed on each sample to yield a

An embedded and charset-unspecified text was scrubbed... Name: nem el?rhet? URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110421/2b2b582b/attachment.pl>

Thanks, I was getting to try this, but got side tracked by actual work... Your analysis reproduces the SPSS unscaled estimates. It still remains to figure out how Nick got > coefficients(lm(ZDEPRESSION ~ ZMEAN_PA * ZDIVERSITY_PA, data=s1)) (Intercept) ZMEAN_PA ZDIVERSITY_PA ZMEAN_PA:ZDIVERSITY_PA 0.07342198 -0.39650356

Hi John, Thanks for the comment... but that appears to mean that SPSS has a big problem. I have always been told that to include an interaction term in a regression, the only way is to do the multiplication by hand. But then it seems to be impossible to stop SPSS from re-standardizing the variable that corresponds to the interaction term. Am I missing something? Is there a way to perform the

Hi Simon, Yes, if I uses coefficients() I get the same results for lm() and lm.ridge(). So that's consistent, at least. Interestingly, the "wrong" number I get from lm.ridge()$coef agrees with the value from SPSS to 5dp, which is an interesting coincidence if these numbers have no particular external meaning in lm.ridge(). Kind regards, Nick ----- Original Message -----

Hi, I have some problems with my repeated measures analysis. When I compute it with SPSS I get different results than with R. Probably I am doing something wrong in R. I have two groups (1,2) both having to solve a task under two conditions (1,2). That is one between subject factor (group) and one within subject factor (task). I tried the following: aov(Score

Hallo, I hope I am posting to the right place. I was advised to try this list by Ben Bolker (https://twitter.com/bolkerb/status/859909918446497795). I also posted this question to StackOverflow (http://stackoverflow.com/questions/43771269/lm-gives-different-results-from-lm-ridgelambda-0). I am a relative newcomer to R, but I wrote my first program in 1975 and have been paid to program in about

Hi all, 1. I am doing some data analysis using both R and SPSS, to get used to both software packages. I had performed linear regression in R and SPSS for a number of times before this last one and the resulting coefficient values always matched. However, this last data set I was analyzing using simple linear regression and using the command lm(y~x), gave me different readings from R and

I have two questions I would appreciate assistance with: (1) I believe "by" is the command used to split a file. In the following example, "mydat" is the dataframe , "group" the variable I want to split the analysis by & "WK1FREQ,WK2FREQ" the variables attach(mydat) by (GROUP) cor.test (WK1FREQ,WK2FREQ) I have also tried: (by (GROUP) cor.test

Dear List, I apologize in advance if this is silly. I tried to replicate an analysis I did previously in SPSS using R, and was surprised to find different results. So my question is: shouldn't the following SPSS syntax REGRESSION DEPENDENT INC89 /METHOD=ENTER hiedyrs experien SE93rec. Yeld the same results of the following R command modelB<-lm(INC89~HIEDYRS+EXPERIEN+SE93REC) I

dear all, i'm trying to reproduce an spss-anova in R. It is an 2x3x3 repeated measures desingn with repeated contrasts. In R i've coded a contrast matrix for all factors and made a split in the aov summary - but I can't get the repeated interaction contrasts. The output from SPSS looks like this: TaskSw * CongNow * CongBefore: SS df Mean Square F Sig. 1 vs. 2 1 vs. 2 1 vs. 2

I asked you before, but in case you missed it: Are you looking at the right place in SPSS output? The UNstandardized coefficients should be comparable to R, i.e. the "B" column, not "Beta". -pd > On 5 May 2017, at 01:58 , Nick Brown <nick.brown at free.fr> wrote: > > Hi Simon, > > Yes, if I uses coefficients() I get the same results for lm() and

Dear all, I have found some difference of the results between multinom() function in R and multinomial logistic regression in SPSS software. The input data, model and parameters are below: choles <- c(94, 158, 133, 164, 162, 182, 140, 157, 146, 182); sbp <- c(105, 121, 128, 149, 132, 103, 97, 128, 114, 129); case <- c(1, 3, 3, 2, 1, 2, 3, 1, 2, 2); result <- multinom(case ~ choles

I had no problems running regression models in SPSS and R that yielded the same results for these data. The difference you are observing is from fitting different models. In R, you fitted: res <- lm(DEPRESSION ~ ZMEAN_PA * ZDIVERSITY_PA, data=dat) summary(res) The interaction term is the product of ZMEAN_PA and ZDIVERSITY_PA. This is not a standardized variable itself and not the same as

Hi Nick, I think that the problem here is your use of $coef to extract the coefficients of the ridge regression. The help for lm.ridge states that coef is a "matrix of coefficients, one row for each value of lambda. Note that these are not on the original scale and are for use by the coef method." I ran a small test with simulated data, code is copied below, and indeed the output from

Dear Nick, On 2017-05-05, 9:40 AM, "R-devel on behalf of Nick Brown" <r-devel-bounces at r-project.org on behalf of nick.brown at free.fr> wrote: >>I conjecture that something in the vicinity of >> res <- lm(DEPRESSION ~ scale(ZMEAN_PA) + scale(ZDIVERSITY_PA) + >>scale(ZMEAN_PA * ZDIVERSITY_PA), data=dat) >>summary(res) >> would reproduce the

SPSS seems to have changed its default datafile format, resulting in issues for read.spss(). In Windows this results in a warning, in Debian the import completely fails: Debian (R version 2.8.0 (2008-10-20) i486-pc-linux-gnu, foreign_0.8-29) > read.spss("/home/jeroen/samples/Tomato.sav") Error in iconv(names(rval), cp, "") : unsupported conversion from 'CP65001'

Hello everyone! My data is composed of 277 individuals measured on 8 binary variables (1=yes, 2=no). I did two similar cluster analyses, one on SPSS 18.0 and one on R 2.9.2. The objective is to have the means for each variable per retained cluster. 1) the R analysis ran as followed: > call data > dist=dist(data,method="euclidean") >

Hi, Here is (I hope) all the relevant output from R. > mean(s1$ZDEPRESSION, na.rm=T) [1] -1.041546e-16 > mean(s1$ZDIVERSITY_PA, na.rm=T) [1] -9.660583e-16 > mean(s1$ZMEAN_PA, na.rm=T) [1] -5.430282e-15 > lm.ridge(ZDEPRESSION ~ ZMEAN_PA * ZDIVERSITY_PA, data=s1)$coef ZMEAN_PA ZDIVERSITY_PA ZMEAN_PA:ZDIVERSITY_PA -0.3962254 -0.3636026