similar to: difference in results from R vs SPSS

Thanks Mike. My data has longitude and latitude coords and I used distAB {clim.pact} then mst {ape} to calculate my minimum spanning tree. The nodes are telecoms sites from all over Australia. My goal is to determine the minimum cost of linking them via cabling, and I'm starting by calculating the distance "as the crow flies", but will probably eventually need to calculate the

I am new to R and have been doing my utmost to learn it in order to teach my students how to use it as an alternative to programs you have to purchase ­ and I am very grateful for the program and what it can doŠ.but I am having a slight problem as I also am not very much into programming ­ though willing to learn. I have been trying to apply a mixed model ANOVA to some data I have. In short ­ one

R-Help Please can you tell me if there is a function within R that will return the day of any given date? Using 19 March 2006 as an example: given 20060319 (or any other date format) will return the day Sunday. Many thanks Fiona.

Dear All At institutes where there are PC labs used by Masters and undergraduate students and software is maintained /managed by a central IT services, do you always update R every time a new version is released - around beginning of April and October? Do you also ask them to install all patch releases as they turn up? Currently our IT department update once a year before the start of the Autumn

I've tried to download rw1090.exe 3 times and rw1081.exe once onto a Windows XP from the UK mirrors (I've tried both of them). When the download is complete and I double-click on the icon I get an error message telling me the set-up files are corrupted and to obtain a new copy of the program. I have spent quite a lot of time checking the message archives but I can't find a solution. If

Hi I have just upgraded from R2.6.0 to R2.7.1 (running on Windows) and a part of my code that previously ran ok now gives an error. The following is a simple example to demonstrate my problem. > a <- array(c(1,2,3,4,5,6,rep(NA,6)),dim=c(6,2)) > apply(a,2,sd,na.rm=T) In R2.6.0 this gives (which is what I would like) [1] 1.870829 NA In R2.7.1 it gives the following error

Hi, I've met a problem when I am trying llvm. I've added a small pass in it, and when I define or use an object of APInt class, it will get the following output while running: $ opt -load hello.so -hello test.bc -o test.hello.bc opt: CommandLine Error: Argument 'disable-opt' defined more than once! opt: CommandLine Error: Argument 'disable-inlining' defined more than

Is it possible to create a contour plot with the isolines labeled. I know you can do this with Matlab. Argh! I tried creating a filled contour plot, then using par(new=T), followed by overlaying the contour plot on top. However, the placement of the filled contour plot and the contour plot do not align correctly. Any suggestions would be appreciated. Thanks, Jon -- View this message in

Dear collegues, I have used R statistical program, package 'lmer', several times already. I never encountered major differences in the outcome between SPSS and R. ...untill my last analyses. Would some know were the huge differences come from. Thanks in advance, Ronald In SPSS the Pearson correlation between variable 1 and variable 2 is 31% p<0.001. In SPSS binary logistic

Dear colleagues, I know SPSS can not compute linear mixed models. I used 'R' before for computing multivariate analyses. But, I never encountered such a major difference in outcome between SPSS and 'R': In SPSS the Pearson correlation between variable 1 and variable 2 is 31% p<0.001. In SPSS binary logistic regression gives us an Odds Ratio (OR)=4.9 (95% CI 2.7-9.0),

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:

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

Hi I was wondering if someone could help me with an rpart problem. I can see that cross-validation is the default for tree selection in rpart -- has a bootstrap method been implemented anywhere? I think this is a different thing to 'bagging' or 'boosting' -- I still want 'one' tree at the end, I just would like it chosen using a bootstrap method. Any ideas??? Thanks

Hi, I am puzzeled with a differing result of princomp in R and FACTOR in SPSS. Regarding the amount of explained Variance, the two results are the same. However, the loadings differ substantially, in the unrotated as well as in the rotated form. In both cases correlation matrices are analyzed. The sums of the squared components is one in both programs. Maybe there is an obvious reason, but I

Dear list, I am calculating the 95th percentile of a set of values with R and with SPSS In R: > normal200<-rnorm(200,0,1) > qnorm(0.95,mean=mean(normal200),sd=sd(normal200),lower.tail =TRUE) [1] 1.84191 In SPSS, if I use the same 200 values and select Analyze -> Descriptive Statistics -> Frequencies and under "Statistics", I type in '95' under Percentiles,

First off, Marc Schwartz posted this link earlier today, read it. http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-are-p_002dvalues-not-di splayed-when-using-lmer_0028_0029_003f Second, your email is not really descriptive enough. I have no idea what OR is, so I have no reaction. Third, you're comparing estimates from different methods of estimation. lmer will give standard errors that

I am trying to stream various 320kbps MP3s across my 10Mbps LAN. Either locally or across the network, freeamp and xmms prefetch a chunk, then play it, then fetch another chunk etc. No continuous streaming. shout gets a lot of !s between its "."s when sending the file to icecast. Any ideas? -- Always hardwire the explosives -- Fiona Dexter quoting Monkey, J. Gregory Keyes, Dark

I am using the latest version of R 2.13.1 and need to load the following packages; library(maps) library(mapdata) library(mapproj) library(lattice) library(tgp) library(spatstat) library(akima) I get the follwing error when loading the packages, "maps", "mapdata" and "mapproj" Error: package 'maps' is not installed for 'arch=i386' Error: package

Hi all, Our section at work is looking at buying a high-powered computer (32 Gb RAM, 64-bit) to be able run models on large datasets and have more processing power for the software we commonly use. We mostly use R to fit GLMs and GAMs. Our department uses Windows as the standard OS (we are upgrading to Vista in the next few months) and so most people would be in favour of having the same OS on

I've just had polarity reversal provisioned by our telco to test hangup detect with a TDM400P I've set hanguponpolarityswitch=yes in zapata.conf When I start Asterisk I get "ignoring hanguponpolarityswitch" in /var/log/asterisk/messages I assume that the option is either not valid or conflicts with another setting somewhere. Any ideas?