Re: [R] Importing Sas files to R

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of Rekha > Sent: Monday, December 27, 2010 8:33 PM > To: r-help@r-project.org > Subject: [R] Importing Sas files to R > > Dear All, > > Am trying to import SAS files to R. > > For

[R] Importing Sas files to R

Dear All, Am trying to import SAS files to R. For that i used library "Hmisc" and try<-sas.get("C:\\ex.sas7bdat") . But it shows error that Error in sas.get("C:\\ex.sas7bdat") : SAS member name is required In addition: Warning message: In sas.get("C:\\ex.sas7bdat") : C:\ex.sas7bdat/formats.sc

Re: [R] Reading sas7bdat files into R

You might post to r-help (see www.r-project.org for how to get to the R mailing lists). There are perhaps 3 programs that will read sas7bdat files. One is the SAS dataset reader (I forgot its actual name) which produces csv files (watch out for bugs such as unmatched quote marks in character var

Re: [R] linear regression with dates

On Dec 27, 2010, at 8:56 PM, Entropi ntrp wrote: Thanks for the response. I proivded the necessary details below, and also have a general question for how to deal with dates in R. Is there a way to make R read dates as numbers? Here is the details of the R code: egfr <- read.csv(file.ch

Re: [R] levelplot blocks size

Thanks for your help. Might you be able to explain in a little more detail how to use those functions to solve this specific problem? I'm happy to put in the work myself and have looked up those functions but am new to R and still a little unsure about how I would go about using those functions

Re: [R] linear regression with dates

Thanks for the response. I proivded the necessary details below, and also have a general question for how to deal with dates in R. Is there a way to make R read dates as numbers? Here is the details of the R code: egfr <- read.csv(file.choose(), header=TRUE, sep=",") #egfr is a matrix read fr

Re: [R] Reading sas7bdat files into R

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of John Sorkin > Sent: Monday, December 27, 2010 5:18 PM > To: cha...@ohsu.edu; r-help@r-project.org > Subject: Re: [R] Reading sas7bdat files into R > > Try sas.get to read the libr

Re: [R] another superscript problem

On 2010-12-27 16:51, David Winsemius wrote: On Dec 27, 2010, at 6:40 PM, T.D. Rudolph wrote: I've exceeded the maximum time I am willing to accept for solving simple problems so I thank all in advance for your assistance. I am trying to plot text combined with an object value and a superscri

Re: [R] Reading sas7bdat files into R

Try sas.get to read the library. You will need the Hmisc package. Let me know if you have success. John John Sorkin Chief Biostatistics and Informatics Univ. of Maryland School of Medicine Division of Gerontology and Geriatric Medicine jsor...@grecc.umaryland.edu -Original Message- From: S

Re: [R] Reading sas7bdat files into R

Thanks for your reply John. I guess you are correct, sasxport.get function is for files in SAS transfer format. Do you know how I can convert the .sas7bdat file into SAS export file that I suppose will be in SAS transfer format? Thanks Sashi -Original Message- From: John Sorkin [mail

Re: [R] Reading sas7bdat files into R

I can't say for sure as I have never used sasxport.get, but I believe that a SAS export file is not the same as a sas library. I believe you are trying to read a SAS library file with an R program that may have been written to read a SAS export file. John John Sorkin Chief Biostatistics and Info

Re: [R] another superscript problem

On Dec 27, 2010, at 6:40 PM, T.D. Rudolph wrote: I've exceeded the maximum time I am willing to accept for solving simple problems so I thank all in advance for your assistance. I am trying to plot text combined with an object value and a superscript. obv = 5 text = "Population mean ="

[R] Reading sas7bdat files into R

Hi All, I am trying to import a .sas7bdat file into R. I tried using Hmisc package's sasxport.get() function. temp <- sasxport.get("path\abcd.sas7bdat") I get an error that says "Error in lookup.xport(file) : file not in SAS transfer format" I am not familiar with SAS transfer format. Could

Re: [R] Finding indexes of minum and maximum elements of an array

Here's one way: H = array(c(8:5,1:4),dim=c(2,2,2)) which(H==min(H),arr.ind=TRUE) dim1 dim2 dim3 [1,]112 which(H==max(H),arr.ind=TRUE) dim1 dim2 dim3 [1,]111 - Phil Spector Statis

[R] another superscript problem

I've exceeded the maximum time I am willing to accept for solving simple problems so I thank all in advance for your assistance. I am trying to plot text combined with an object value and a superscript. obv = 5 text = "Population mean =" ss = ^o # degrees Something like this (very naive so yo

Re: [R] linear regression with dates

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of Entropi ntrp > Sent: Monday, December 27, 2010 3:05 PM > To: r-help@r-project.org > Subject: [R] linear regression with dates > > Hi, > I am trying to do simple linear regression

[R] linear regression with dates

Hi, I am trying to do simple linear regression using dates in R but receiving error messages. With the data shown below, I would like to regress x on y. x y 11/12/1999 56.8 11/29/1999 17.9 01/04/2000 27.4 1/14/2000 96.8 1/31/2000 49.5 R gives the following erro

Re: [R] Fitting mixed effects Baseline category logit models

Hi: Try the mlogit() function in package mlogit. It appears to handle random effects terms, but I don't know offhand if that includes hierarchical random effects. You can always try it out... HTH, Dennis On Mon, Dec 27, 2010 at 2:43 PM, DHIMAN BHADRA wrote: > Hello everyone, > > I want to fit a

Re: [R] Finding indexes of minum and maximum elements of an array

On 10-12-27 5:26 PM, Eduardo de Oliveira Horta wrote: Hello there I wish to get the "coordinates" of the minimum element of an array. For example, if the array were H = array(c(8:5,1:4),dim=c(2,2,2)) H , , 1 [,1] [,2] [1,]86 [2,]75 , , 2 [,1] [,2] [1,]1

[R] Fitting mixed effects Baseline category logit models

Hello everyone, I want to fit a baseline category logit model (with 3-4 categories) with nested random effects. (For example, I have clusters(i) and households within clusters (j) resulting in the nested random effects structure : b_i +d_j(i)). Is there a R function/package that I can use ? Any h

[R] Finding indexes of minum and maximum elements of an array

Hello there I wish to get the "coordinates" of the minimum element of an array. For example, if the array were > H = array(c(8:5,1:4),dim=c(2,2,2)) > H , , 1 [,1] [,2] [1,]86 [2,]75 , , 2 [,1] [,2] [1,]13 [2,]24 then > min(H) [1] 1 and > max(H) [1]

Re: [R] Gamma & Lognormal Model

On 10-12-27 3:08 PM, Louisa wrote: Dear, I'm very new to R Gui and I have to make an assignment on Gamma Regressions. Surfing on the web doesn't help me very much so i hope this forum may be a step forward. This forum isn't for homework. You should ask your instructor for help. Duncan Murdo

[R] Gamma & Lognormal Model

Dear, I'm very new to R Gui and I have to make an assignment on Gamma Regressions. Surfing on the web doesn't help me very much so i hope this forum may be a step forward. The question sounds as follows: The data set is in the library MASS first install library(MASS) then type data(mammals)

[R] R-code to generate random rotation matrix for rotation testing

Dear list, I am looking for an implementation of random rotation matrix generation in R to do a rotation test: I want to use the matrices to create random multivariate normal matrices with common covariance structure and mean based on an observed data matrix. The rRotationMatrix-function in th

[R] Any functions to manipulate (merge, cut, remove) hclust objects? (maybe through phylo?)

Hello all, I'm now working with hclust objects and was hoping to perform some basic editing on them like: - Joining = the merging of two hclust objects (so they will share one root) - Splicing = So to cut/extract a branch out of an hclust object - that by itself will be an hclust obje

Re: [R] Problem using pkg "survival"

On Dec 27, 2010, at 19:48 , James Steinberg wrote: > Hello all. > > I've been attempting to utilize the "survival" pkg ( > http://cran.r-project.org/web/packages/survival/index.html), while reading > through this guide (http://www.ms.uky.edu/~mai/Rsurv.pdf). I figured working > through the guide

Re: [R] matrix looping accessing previous column

Using which() improves my code, but now I'm receiving some data.frame error and still does not convert the values to negative. #TAKE PART OF MATRIX CEM1_PARTIAL <- CEM1[1:10,1:10] #CEM1_PARTIAL=apply(CEM1_PARTIAL,2,as.character) #CEM1_PARTIAL <- as.numeric(as.character(CEM1_PARTIAL)) #CEM1_PART

[R] Problem using pkg "survival"

Hello all. I've been attempting to utilize the "survival" pkg ( http://cran.r-project.org/web/packages/survival/index.html), while reading through this guide (http://www.ms.uky.edu/~mai/Rsurv.pdf). I figured working through the guide would be the best way to go, before attempting my own data. I t

Re: [R] can't install R with *local* gcc

Hello, I infer from your answer that piece of information which is most important for us now: That my configuration instructions for R should work (should be correct). Given all these unspecified build-systems we have to fight (at least Gmp, Mpfr, Gcc with everything in it, R, and finally our own

Re: [R] matrix looping accessing previous column

That just converts my values in V3 to NA -- View this message in context: http://r.789695.n4.nabble.com/matrix-looping-accessing-previous-column-tp3165308p3165340.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.or

Re: [R] matrix looping accessing previous column

Convert your column to numeric: CEM1_PARTIAL$V3 <- as.numeric(as.character(CEM1_PARTIAL$V3)) On Mon, Dec 27, 2010 at 4:28 PM, R_novice wrote: > > Thank you for the quick response :-). I've applied your suggestion to my > code, but I still receive an error: > > > CEM1_PARTIAL$V3[which(CEM1_PART

Re: [R] matrix looping accessing previous column

Thank you for the quick response :-). I've applied your suggestion to my code, but I still receive an error: > CEM1_PARTIAL$V3[which(CEM1_PARTIAL$V4 == "DOWN") - 1] <- > CEM1_PARTIAL$V3[which(CEM1_PARTIAL$V4 == "DOWN") - 1] * -1 Warning message: In Ops.factor(CEM1_PARTIAL$V3[which(CEM1_PARTIAL$V4

Re: [R] need help with data management

On Dec 25, 1:36 pm, Gabor Grothendieck wrote: > On Sat, Dec 25, 2010 at 8:08 AM, analys...@hotmail.com > > > > > > wrote: > > I have a data frame that reads > > > client ID date transcations > > > 323232   11/1/2010 22 > > 323232   11/2/2010 0 > > 323232   11/3/2010 missing > > 121212   11/10/2

Re: [R] Calculation of BIC done by leaps-package

On 2010-12-26 06:47, Jan Henckens wrote: Hi Folks, I've got a question concerning the calculation of the Schwarz-Criterion (BIC) done by summary.regsubsets() of the leaps-package: Using regsubsets() to perform subset-selection I receive an regsubsets object that can be summarized by summary.reg

Re: [R] Help on lme output. what are the estimates in $modelStruct

Lei wrote: > > > I have been playing with lme fitting different kinds of model. I am > puzzled by the estimates associated with $modelStruct and I have no idea > what they are? > > e.g if I fit a random intercept model using the following > > fit.0<-lme(y~time,random=list(subject=pdIdent(~1))

Re: [R] matrix looping accessing previous column

Try this: > x V1 V2 V3 V4 1 PROBE 1 2.5 UP 2 PROBE 2 1.0 UP 3 PROBE 3 1.4 DOWN 4 PROBE 4 2.0 UP 5 PROBE 5 1.3 DOWN x$V3[which(x$V4 == "DOWN") - 1] <- x$V3[which(x$V4 == "DOWN") - 1] * -1 On Mon, Dec 27, 2010 at 4:11 PM, R_novice wrote: > > Hi, > > I have a matrix with number

[R] matrix looping accessing previous column

Hi, I have a matrix with numbers and character. I want to evaluate each cell and change the value of the cell before it depending on the evaluation. My evaluation: if a cell had the word "down" change the cell preceding it to a negative number by multiplying that value by a -1. I am have trouble

Re: [R] Where to find "rcompgen" package?

On 25.12.2010 15:44, David Winsemius wrote: On Dec 25, 2010, at 8:08 AM, Megh Dal wrote: Dear all, can somebody point me from where to download "rcompgen" package? CRAN does not seem to hold that. http://lmgtfy.com/?q=rcompgen Or more specific (since the first Google hits do not tell t

Re: [R] read.data? without separator

Would the read.fwf function work for you? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Fror > Sent: Sunday

[R] [R-pkgs] First release of package "cg" for comparison of groups

Version 0.9.0, a first release of cg, is published on CRAN and its mirrors. The "cg" name stands for "compare groups." Its genesis and evolution are driven by common "in-practice" needs to compare samples, treatments, administrations, conditions, etc. in medicine research & development. The current

Re: [R] Can't merge on datetime?

Convert to POSIXct: x$date <- as.POSIXct(x$date, tz = "GMT") y$date <- as.POSIXct(y$date, tz = "GMT") merge(x, y, all = TRUE) On Mon, Dec 27, 2010 at 2:14 PM, Dimitri Shvorob wrote: > > x = structure(list(date = structure(list(sec = c(0, 0, 0, 0, 0, 0, > 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,

[R] [R-pkgs] mcga 1.1 (machine coded genetic algorithms) package released

mcga 1.1 (machine coded genetic algorithms) package implements a genetic algorithm optimisation tool with machine coded chromosomes. The machine coded chromosomes stand for chromosomes that are not decoded and encoded. The byte representation of 'double' type variables are crossed-over and muta

[R] Can't merge on datetime?

x = structure(list(date = structure(list(sec = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), min = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), hour = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,

Re: [R] Nagelkerke R square for Prediction data

Hello I found some small postings dated to 22 Oct 2008 on the message subject. Recently, I have been working with binary logistic regressions. I didn't use the design package. Yet, I needed the "fit" indices. Therefore, I wrote a small function to output the Nagelkerke's R, and the Cox-&-Snell

Re: [R] Drop column from a data frame

On Dec 26, 2010, at 8:22 PM, John Sorkin wrote: I am trying to drop a column of a data frame. The code below attempts to drop a numeric column (which does not work but gives no error or warning) and a factor column (which does not work but gives an error). I would appreciate someone tellin

Re: [R] levelplot blocks size

Look at the functions cut, findInterval, tapply, and aggregate. Sent from my iPod On Dec 26, 2010, at 4:34 PM, "jonathan" wrote: > > Thanks for your advice, but my data is not decimals, so I don't need to round > the values. Instead, what I need to really do is "group" the values into > large

Re: [R] filled.contour colors

Great, I will have a look at it! Thanks! On Mon, Dec 27, 2010 at 12:52 PM, Jim Lemon [via R] wrote: > On 12/27/2010 02:21 PM, randhindi wrote: >> >> Hi, >> >> I am trying to set the color scale in filled.contour based on a specific >> value instead of a relative position. >> Specifically, I wan

Re: [R] Parsing a Simple Chemical Formula

> From: han...@depauw.edu > To: dwinsem...@comcast.net > Date: Sun, 26 Dec 2010 22:36:49 -0500 > CC: r-h...@stat.math.ethz.ch > Subject: Re: [R] Parsing a Simple Chemical Formula > > Hi David & others... > > I did find the function you recommended, plus, it's even easier (but a > little hidden i

[R] Heteroskedasticity and autocorrelation of residuals

Hello everyone, I'm working on a current linear model Y = a0 + a1* X1 + ... + a7*X7 + residuals. And I know that this model presents both heteroskedasticity (tried Breusch-Pagan test and White test) and residuals autocorrelation (using Durbin Watson test). Ultimately, this model being meant to be

Re: [R] how to replace my double for loop which is little efficient!

Thank Berend, It seems like that it is better to attach a PDF file for avoiding messy code. Yes, I want to obtain is Tanimoto coefficient and your web site "wikipedia" is about this coefficient. I also search R site about tanimoto coefficient and learn it more. About your code, I has saved and

Re: [R] Parsing a Simple Chemical Formula

> Date: Sun, 26 Dec 2010 20:24:23 -0800 > From: spencer.gra...@structuremonitoring.com > To: han...@depauw.edu > CC: r-h...@stat.math.ethz.ch > Subject: Re: [R] Parsing a Simple Chemical Formula > > Mike Marchywka's post mentioned a CRAN package, "rpu

Re: [R] modifying user agent strings in http requests

> From: soumen...@gmail.com > Date: Mon, 27 Dec 2010 11:17:49 +0530 > To: r-h...@stat.math.ethz.ch > Subject: [R] modifying user agent strings in http requests > > Hi all. > > How does one change user agent strings in http requests made in R? And > how do I figure out what my current user agent st

Re: [R] Doing a mixed-ANOVA after accounting for a covariate

Thank you David and Rich for your reply. I will try to get more concrete with an artificial example that replicates the structure of the data sets I'm analyzing. It will include analysis suggestions based on what I was able to figure out from the demo(MMC.apple). At the end there is my current puzz

Re: [R] filled.contour colors

On 12/27/2010 02:21 PM, randhindi wrote: Hi, I am trying to set the color scale in filled.contour based on a specific value instead of a relative position. Specifically, I want the values below 0 to be in a gradient of green, and those above 0 to be red. 0 would be white. I tried: posZero = a

Re: [R] how to replace my double for loop which is little efficient!

Found this: http://en.wikipedia.org/wiki/Jaccard_index#Tanimoto_coefficient_.28extended_Jaccard_coefficient.29 and an R site search with "tanimoto" yielded some more interesting stuff. The rest is up to you. Berend -- View this message in context: http://r.789695.n4.nabble.com/how-to-replac

Re: [R] Question about mars() -function

On Dec 26, 2010, at 17:54 , Tiina Hakanen wrote: > Hi! > > I have some questions about MARS model's coefficient of determination. I use > the MARS method in my master's thesis and I have noticed some problems with > the MARS model's R^2. > > You can see the following example that the MARS mode

Re: [R] how to replace my double for loop which is little efficient!

bbslover wrote: > > thanks for your help. I am sorry I do not full understand your code, so i > can not correct using your code to my data. here is the attachment of my > data, and what I want to compute is the equation in the word document of > the attachment: > > the code form Berend can get