similar to: sas.get problem

Hello again I have a file from SPSS in .sav format. when I run library(foreign) cvar<-as.data.frame(read.spss("c:\\NDRI\\cvar\\data\\cvar2rev3.sav")) I get a warning Warning message: c:\NDRI\cvar\data\cvar2rev3.sav: Compression bias (0) is not the usual value of 100. The data appear to be OK, but I am concerned. (I tried searching the archives and the documenation for data

I have a dataframe called cvar, with two variables (among many others) called MSA and ACTUP. Both are numeric. This was working fine. Then I found out that for two MSAs, ACTUP should be 1, not 0. so I tried cvar$ACTUP[cvar$MSA == 6840] <- 1 cvar$ACTUP[cvar$MSA == 5360] <- 1 but when I try table(cvar$MSA, cvar$ACTUP) the level of ACTUP for those two MSAs has not changed, and is still

Hello Using R 2.2.1 on a Windows machine. Has anyone programmed the Welch test for equality of variances? I tried RSiteSearch, but this gave references to t test and oneway.test, which are not quite what I need.....I need the Welch test itself, for use in a meta-analysis (to determine if variances are equal). TIA Peter Peter L. Flom, PhD Assistant Director, Statistics and Data Analysis

Hello I am sure this is elementary, but I can't figure it out..... I am using SAS v8.2 and R 1.4.1 on a Windows platform. I have a large (6 megabyte) file in .sd2 format. I want to import it into R. If necessary, I can create a smaller file, as I don't need all the variables that are in this file I downloaded the package foreign. I then tried to use PROC CPORT in SAS to create a

Hi List, Merlin Magix hardware v02 I'm trying to get asterisk to act as a voicemail server for a lucent merlin magix PBX that we purchased used. We have 4 FXO channels between the two PBXs on a Sangoma A200 card. The 770 dialgroup is working properly, in that calls to 770 are answered by Asterisk. The magix is sending mode codes in the format #XX#XXX#, where the 2nd block of digits

# Hi all, # consider the following code (please, run it: # it's fully working and requires just few minutes # to finish): require(CreditMetrics) require(clusterGeneration) install.packages("Rdonlp2", repos= c("http://R-Forge.R-project.org", getOption("repos"))) install.packages("Rsolnp2", repos= c("http://R-Forge.R-project.org",

Hello I have a dataframe that looks like this MSA CITY HIVEST YEAR YR CAT 1 0200 Albuquerque 0.50 1996 1996 5 2 0520 Atlanta 13.00 1997 1997 5 3 0720 Baltimore 29.10 1994 1994 1 4 0720 Baltimore 13.00 1995 1995 5 5 0720 Baltimore 3.68

Dear R-users, I am seeking advises on the sas.get function from the Hmisc package. I have tried to import some of our SAS files using the syntax presented in the function help example but the importation always failed. The function does not seem to recognize our sas files and complains about the lack of format library files (I am not SAS proficient, but I guess that is what the

Good afternoon Running R2.10.0 on Windows I have a data frame that includes (among much else) a factor (In_2006) and a continuous variable (math_3_4). I would like to find the 2 cases for In_2006 = 0 that are closest to each case where In_2006 = 1. My data looks like In_2006 math_3_4 0 55.1 1 51.6 1 18.1 1 26.6 1 14.1

Hello I am running R 1.8.1 on a Windows platform I am attempting to fit an ordinal logistic regression model, using the polr function, as described in Venables and Ripley. But when I try model4 <- polr(ypsxcat~committed + as.factor(sex) + as.factor(drugusey) + anycsw + as.factor(sex)*committed + as.factor(sex)*as.factor(drugusey)+as.factor(sex)*anycsw, data = duhray) I get a message

Dear R Users, I have just installed R 2.8 today and also updated the Hmisc package to the latest 3.4-4 version. However, the sas.get function no longer works in the new Hmsic 3.4-4. Below is a sample code and output from the old and new, but the old one worked. Can someone tell me how sas.get should be used in Hmisc 3.4-4. Thanks, Larry Package Hmisc version 3.0-12 > a1 <-

Pursing my earlier question, when I tried loading Lindsey's gnlm, I got a message Loading required package: rmutil Warning message: There is no package called 'rmutil' in: library(package, character.only = TRUE, logical = TRUE, warn.conflicts = warn.conflicts, According to the R documentation http://finzi.psych.upenn.edu/R/doc/html/packages.html rmutil is in the standard

Running R 1.8.1 on a Windows machine In dotchart, I would like to shrink the labels on the tick marks (that is, the numbers) without shrinking anything else. I could not find this in either the Rhelp archives or in ?dotchart, which recmmends cex to avoid 'label overlap', but cex shrinks all the characters in the plot. Is there a way to do this? Thanks Peter Peter L. Flom, PhD

Hello, Is there any way we convert SAS dataset into R dataset? Kindest Regards, Peppy Adi-Purnomo ------ Peppy Adi-Purnomo Energy Market Analyst Energy Link Ltd Dunedin - New Zealand Ph.: +64 3 479 2475 Fax: +64 3 477 8424 Email: s.adi.purnomo at energylink.co.nz www.EnergyLink.co.nz -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Dear R-helpers, I am helping a SAS user run some analyses in R that she cannot do in SAS and she is complaining about R's peculiar (to her!) way of recoding variables. In particular, she is wondering if there is an R package that allows this kind of SAS recoding: IF TYPE='TRUCK' and count=12 THEN VEHICLES=TRUCK+((CAR+BIKE)/2.2); Thanks for any help or suggestions you might be able

Hello I am writing an article demonstrating how quantile plots work. In one figure, I generate 3 sets of normal data, and add different sorts of outliers to each. I then make 9 qqnorm plots, in a 3x3 array: Outlier1 Outler2 Outlier3 Dataset1 Dataset2 Dataset3 which I did using mfrow = c(3,3). I'd like to use the space on the page as efficiently as possible,so I

I used garchFit function to fit 1600 observations of EURO/USD 2-day returns in GARCH(1,1) model. As part of the summary I got warning message: NaNs produced in: sqrt(diag(fit$cvar)) And didn't get any estimates for 3 params' std.error, t value or probability: Error Analysis: Estimate Std. Error t value Pr(>|t|) mu -0.004827 0.020141 -0.240 0.811 ar1 0.010311

Hello, I am trying to import SAS files into R using sas.get() in Hmisc. I am running both SAS and R on a UNIX server. I've pasted my command and ensuing errors below: sas.get(lib="EOG", mem="eog3pub00_1") Error in if (status != 0) { : argument is of length zero In addition: Warning message: In sas.get(lib = "EOG", mem = "eog3pub00_1") :

Thanks for all the help on my earlier questions..... How do you plot a simple time series with unequal intervals? I have the following dateofpoll <- as.ts(c("6/1/02", "7/1/02", "10/1/02", "1/4/03", "1/25/03", "6/7/03", "7/16/03", "8/17/03", "9/4/03")) reelect <- c(51, 47, 49, 51, 49, 49, 46, 45, 40)

Hello all I have a feeling this is very simple......but I am not sure how to do it My boss has two variables, one is an average of 4 numbers, the other is an average of 3 of those numbers i.e var1 = (X1 + X2 + X3 + X4)/4 var2 = (X1 + X2 + X3)/3 all of the X variables are supposed to be measuring similar constructs not surprisingly, these are highly correlated (r = .98), the question is how