similar to: GLM fitting in R and Statistica

I am a new R user attempting to convert files from Statistica to R. I can export from Statistica to SPSS .por format, but not to the SPSS .sav format. Is there a procedure for easily accomplishing this, which will allow me to keep variable short and long labels (big surveys LOTS of time to replace all this work). Many thanks for sharing your time and knowledge. steve Dr Steven M Burgess

Greetings, I've just put out the latest version of "The Popularity of Data Analysis Software" at http://r4stats.com/popularity. This update includes complete data for 2010, the addition of number of blogs for each software, more coverage of Statistica, and, where possible, measures regarding the implementations of the SAS Language: Carolina and the World Programming System (WPS).

Dear All, I am in the somewhat unfortunate position of having to reproduce the results previously obtained from (non-metric?) MDS on a "kinship" matrix using Statistica. A kinship matrix measures affinity between groups, and has its maximum values on the diagonal. Apparently, starting with a nxn kinship matrix, all it was needed to do was to feed it to Statistica flagging that the

Hello. I apologize in advance for the VERY lengthy e-mail. I endeavor to include enough detail. I have a question about survival curves I have been battling off and on for a few months. No one local seems to be able to help, so I turn here. The issue seems to either be how R calculates Kaplan-Meier Plots, or something with the underlying statistic itself that I am misunderstanding. Basically,

Better would be to use interval censored data. Create your data set so that you have (time1, time2) pairs, each of which describes the interval of time over which the tag was lost. So an animal first captured at time 10 sans tag would be (0,10); with tag at 5 and without at 20 would be (5,20), and last seen with tag at 30 would be (30, NA). Then survit(Surv(time1, time2,

Hello, this is a typical "how do you do it in Rails" question :) I would like to incorporate fit testing into a RoR app and I''d like to hear your recommendations on how to do it. In the ideal world I would like all the table test definitions and implementations stored in svn together with the app, ability to run the tests with a rake task + a hook to FitNesse, so that

Hi, I have a strange problem with testing classes which have a HABTM relation. It seems that while running my tests with rake the fixture for the association table is not loaded - as a result I get a failure on a test which simply checks the number of associated elements - it says it''s 0 although there should be one. The weird thing is that when I run the tests with rcov

Hi, I'd like to perform Dufrene-Legendre Indicator Species Analysis for a multivariate regression tree. However I have problems with arguments of duleg(veg,class,numitr=1000)function. How to obtain a vector of numeric class memberships for samples, or a classification object returned from mvpart? thanks in advance -- Best regards, Agnieszka Strzelczak

Hi everybody, I would like to compute the mean for 1 variable between the rows with the same levels. For example, with the dataset below: Factor1 Factor2 Value A X 1 A X 2 A Y 3 A Y 4 B X 5 B X 6 B Y

Hi, I'm trying to find out what threshold of indicator value in labadsv should be used to accept a specie as an indicator one? So far I assumed that indval=0.5 is high enough to avoid any mistakes but it was based only in my intuition. I'd be greatful for any advise best regards Agnieszka

Hello, I''m a novice to ferret, so far only used it via acts_as_ferret. My question is about a recommended pattern for an ''advanced search'', which would be searching by all fields of a model and some fields from related models, with range search, expression search and wildcards. The kind of search in which a user is presented with a huge form that allows them to set

Hello, I am a new R user and trying to learn how to implement the mahalanobis function to measure the distance between to 2 population centroids. I have used STATISTICA to calculate these differences, but was hoping to learn to do the analysis in R. I have implemented the code as below, but my results are very different from that of STATISTICA, and I believe I may not have interpreted the help

Probably, it's an obvious info, but I have not found anything in R FAQ concerning this feature/bug. The results of ks.test and wilcoxon.test (in the Mann-Whitney version, paired = 'FALSE') don't coincide with the results from the other statistical packages, e.g. Statistica, Medcalc, and (as for MW test) from the numerous online MW tests. E.g. Statistica p-value=0.0435353 Medcalc

Given this example #start code a<-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940, 760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430) b<-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90, 3220,490,20790,290,740,5350,940,3910,0,640,850,260) wilcox.test(a, b, paired=FALSE) #sum of rank for first sample sum.rank.a <-

Dear R-users, I used to give commands such as: > source(file="~/path/to/file.R", chdir=TRUE) but with the latest v. 2.1.0 it does not seem to work anymore. I tried to figure out what it was going on and it seems that the string for which > class(file) [1] "character" is changed to > class(file) [1] "file" "connection" when the connection is

Dear R gurus I am analysing data from a study of behaviour and shade utilization of chimpanzees. I am using GLMs in R (version 2.13.0) to test whether shade/sun utilization is predicted by behaviour observed. I am thus interested in whether an interaction of behaviour (as a predictor) and presence in the sun/shade (also predictor) predicts the counts I have for the respective categories. I have

Hello world, I am actually trying to transfer a lecture from Statistica to R and I ran into problems with spectral analysis, I think I just don't get it 8-( (The posting from "FFT, frequs, magnitudes, phases" from 2005 did not enlighten me) As a starter for the students I have a 10year data set of air temperature with daily values and I try to get a periodogram where the annual

Hi, Sorry to bother you again. I would like to estimate the effect of several categorical factors (two between subjects and one within subjects) on two continuous dependent variables that probably covary, with subjects as a random effect. *I want to control for the covariance between those two DVs when estimating the effects of the categorical predictors** on those two DVs*. The thing is, i

Hi, I cant understand where I am going wrong.Below is my code.I would really appreciate your help. Thanks. > genfile<-read.table("c:/tina/phd/bs871/hw/genfile.txt",skip=1) > > #read in SNP data > snp.dat <- as.matrix(genfile) > snp.name <- scan("c:/tina/phd/bs871/hw/genfile.txt",nline=1,what="character") Read 100 items

Hi everybody, I am starting my friendship with R. Now I want to do a Correspondence Analysis using the vegan library or multiv library. Both give me the same results which differ from the ones I get running other programs like Statistica... I wounder if there is a problem with my matrix. I have cualitataive variable in the columns and sampling units in the rows... I import my matrix into R