Re: [R] loading igraph package on Solaris

2008-01-10 Thread Gabor Csardi
Peter, the changes you made should be ok. (No changes are needed for gcc/g++ anyway.) The problem is not the C compiler, which is (most likely) gcc, but the c++ compiler, which is not the GNU c++ compiler, but the Sun version. I guess that the GNU c++ compiler is called g++ or gxx. If you found

Re: [R] 4 dimensional graphics

2008-01-10 Thread Petr PIKAL
Thanks to all who replied. I will check all suggestions to see what is the best one in respect of readability. Suggestion of using LSD monitor from Alberto is quite attracting however I am not sure if legal :-) Thank you again and best regards. Petr [EMAIL PROTECTED] jiho <[EMAIL PROTECTED]>

[R] RODBC, postgresql on Linux

2008-01-10 Thread Poizot Emmanuel
Dear all, I facing pbs using RODBC library. I'm working on a Linux (Ubuntu Gutsy) os, with R version 2.5.1. I've got Postgresql 8.2.5 installed with the odbc-postgresql package. I try to connect in a R session to a postgresql database using: odbcConnect(dsn="thedatabasename",uid="theuser",pwf="

[R] hi

2008-01-10 Thread RAJA SEKHAR REDDY MAYALURU
Hi All , i want to use this R project to implement project mix model . i have seen some tutorials about R tool , i did nt fidn any thing related to implement product mix model . my requirement is : i want to give some data and some constraints , the tool need to give the output optimum slout

Re: [R] installing and using a package

2008-01-10 Thread jpolo
> I'm not sure what you have done but using a package requires two steps: > 1) installing or downloading the package, > 2) loading an installed package. > > Since you are using windows you can click on Packages, Install packages, > and > after choosing a CRAN mirror, selecting the package you want.

Re: [R] error message translation

2008-01-10 Thread jpolo
> Do you mean expansions on the terse messages, or translations from > another language? There aren't really any expansions available: the > messages are supposed to be clear in themselves. If you really want to > know why you got one you might be able to look it up in the source code, > but gen

Re: [R] changing the values in a list

2008-01-10 Thread dxc13
I will give it a try, thanks for your input! Peter McMahan wrote: > > > On Jan 10, 2008, at 9:48 PM, Peter McMahan wrote: > >> y <- matrix(rep(y,ncol(x)),nrow=nrow(y)) >> x[is.na(x)] <- y[is.na(x)] > > oops, i meant: > y <- matrix(rep(y,ncol(res)),nrow=nrow(y)) > res[is.na(res)] <- y[is.na(r

Re: [R] changing the values in a list

2008-01-10 Thread Peter McMahan
On Jan 10, 2008, at 9:48 PM, Peter McMahan wrote: > y <- matrix(rep(y,ncol(x)),nrow=nrow(y)) > x[is.na(x)] <- y[is.na(x)] oops, i meant: y <- matrix(rep(y,ncol(res)),nrow=nrow(y)) res[is.na(res)] <- y[is.na(res)] __ R-help@r-project.org mailing list h

Re: [R] changing the values in a list

2008-01-10 Thread Peter McMahan
try: res <- apply(res,2,function(x){x[is.na(x)] <- y[is.na(x)];x}) should work, though not the most efficient/elegant alternately: y <- matrix(rep(y,ncol(x)),nrow=nrow(y)) x[is.na(x)] <- y[is.na(x)] (also not the most efficient/elegant) On Jan 10, 2008, at 9:38 PM, dxc13 wrote: > > useR's, > >

[R] changing the values in a list

2008-01-10 Thread dxc13
useR's, Does anyone know of an efficient way to change the values in a list to different numbers stored in an object? For example, consider the following: y <- matrix(c(20, 13, 25, 10, 17), ncol = 1) > res [[1]] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [

Re: [R] installing and using a package

2008-01-10 Thread jpolo
> R distinguishes between "installing" and "loading" a package. It looks > like your installation of yags went fine. Now type > library(yags) > to load the package. > now help(yags) should bring up something, as should example(yags) . that helped. thank you john __

Re: [R] installing and using a package

2008-01-10 Thread jpolo
>> -Original Message- > You need to read the "An Introduction to R" manual which was probably > installed when you installed R on your un-named OS (Windows, I take it). yes, i did find what i needed there. thanks you. john __ R-help@r-project.

Re: [R] ggplot2, coord_equal and aspect ratio

2008-01-10 Thread hadley wickham
On Jan 10, 2008 6:16 PM, jiho <[EMAIL PROTECTED]> wrote: > Hi everyone, Hi Hadley, > > I am a heavy user of coord_equal() in ggplot2 since most of my data is > spatial, on x,y coordinates. Everything works. However by enforcing an > aspect ratio of 1 for the plotting region, coord_equal() usually >

Re: [R] cairo and rapache [was: Weird errors on MAC OS X Leopard when using Cairo/GDD through Rserve]

2008-01-10 Thread Simon Urbanek
Charles, Greg, the problem is that Apple's FreeType on Leopard uses CF calls (albeit just to get files), so it breaks on fork. (The message is moot, because the whole point of using fork is to *not* use exec). In any case the solution is to compile your own version of freetype, fontconfig

Re: [R] error message translation

2008-01-10 Thread Duncan Murdoch
On 10/01/2008 7:12 PM, [EMAIL PROTECTED] wrote: > sorry if my first msg was posted twice. this email system at school sucks. > > i've looked in the Admin and Main manual and the Intro manual for > something that tells me what error messages mean, but haven't found > anything. where can i find the

Re: [R] Extracting last time value

2008-01-10 Thread Gabor Grothendieck
Here is a slight variation of ones you have already received. This also uses zoo. In the DF<- line we read it in as a data frame. In the z<- line we convert it to zoo using chron date/times (see R News 4/1 for more). We calculate the positions of duplicate dates and in the last line we extract t

Re: [R] Cycle Regression Analysis in R?

2008-01-10 Thread Charles C. Berry
On Thu, 10 Jan 2008, Carson Farmer wrote: > Hello R community, > > Does anyone know of a package that will perform cycle regression > analysis? I have searched the R-help archives etc. but have come up with > nothing so far. > If I am unable to find an existing R package to do so, is there anyone

Re: [R] loading igraph package on Solaris

2008-01-10 Thread Peter McMahan
With a fresh copy of igraph 0.4.5 I made all of the changes at http://code.google.com/p/igraph/issues/detail?id=1 *except* for those to rinterface.c (I know very little about C and it wasn't clear what changes should be made). This still left some errors with pottsmodel_2.cpp ('Overloading ambi

Re: [R] pairwise cross tabulation tables

2008-01-10 Thread Charles C. Berry
On Thu, 10 Jan 2008, AndyZon wrote: > > Thank you so much, Chuck! > > This is brilliant, I just tried some dichotomous variables, it was really > fast. Yes, and if you are on a multicore system with multithreaded linear algebra, crossprod() will distribute the job across the cores making the e

Re: [R] installing and using a package

2008-01-10 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] > Sent: Thursday, January 10, 2008 4:06 PM > To: r-help@r-project.org > Subject: [R] installing and using a package > > i'm trying to learn R. i have no programing experience. i > tr

Re: [R] installing and using a package

2008-01-10 Thread Peter McMahan
R distinguishes between "installing" and "loading" a package. It looks like your installation of yags went fine. Now type library(yags) to load the package. now help(yags) should bring up something, as should example(yags) . __ R-help@r-project.org mai

Re: [R] installing and using a package

2008-01-10 Thread Charles Annis, P.E.
I'm not sure what you have done but using a package requires two steps: 1) installing or downloading the package, 2) loading an installed package. Since you are using windows you can click on Packages, Install packages, and after choosing a CRAN mirror, selecting the package you want. R will do

[R] ggplot2, coord_equal and aspect ratio

2008-01-10 Thread jiho
Hi everyone, Hi Hadley, I am a heavy user of coord_equal() in ggplot2 since most of my data is spatial, on x,y coordinates. Everything works. However by enforcing an aspect ratio of 1 for the plotting region, coord_equal() usually wastes a lot of space if the region of interest is not a perf

[R] error message translation

2008-01-10 Thread jpolo
sorry if my first msg was posted twice. this email system at school sucks. i've looked in the Admin and Main manual and the Intro manual for something that tells me what error messages mean, but haven't found anything. where can i find the meanings of the error messages? _

[R] installing and using a package

2008-01-10 Thread jpolo
i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/wi

[R] installing and using a package

2008-01-10 Thread jpolo
i'm trying to learn R. i have no programing experience. i tried to add on a package and just picked "yags". i can't get it to work or get any help. --- > install.packages("yags") Warning in install.packages("yags") : argument 'lib' is missing: using 'C:\Documents and Settings\j\My Documents/R/wi

Re: [R] Omnibus main effects in summary.lme?

2008-01-10 Thread Andrew Beckerman
Adam - Without resorting to the rather rich lmer/lme4 discussion realm, you need to base anova() comparisons of lme models with different fixed effects on maximum liklihood estimates rather tham REML. anova(update(l2,method="ML"),update(l2,~.-useful:nusience,method="ML")) should avoid the er

Re: [R] Cycle Regression Analysis in R?

2008-01-10 Thread roger koenker
Peter Bloomfield's Fourier Analysis of Time Series is a good reference for this sort of thing... On Jan 10, 2008, at 5:25 PM, Gabor Grothendieck wrote: > To do it from first principles using nonlinear optimization see: > http://finzi.psych.upenn.edu/R/Rhelp02a/archive/20100.html > > On Jan 10, 2

Re: [R] Cycle Regression Analysis in R?

2008-01-10 Thread Gabor Grothendieck
To do it from first principles using nonlinear optimization see: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/20100.html On Jan 10, 2008 5:27 PM, Carson Farmer <[EMAIL PROTECTED]> wrote: > Hello R community, > > Does anyone know of a package that will perform cycle regression > analysis? I have

Re: [R] Extracting last time value

2008-01-10 Thread Achim Zeileis
On Thu, 10 Jan 2008, Kondamani, Arjun (GMI - NY Corporate Bonds) wrote: > I have a dataframe as follows: > > Date timevalue > 20110620 11:18:007 > 20110620 11:39:009 > 20110621 11:41:008 > 20110621 11:40:006 > 20110622 1

Re: [R] Extracting last time value

2008-01-10 Thread jim holtman
This should do it: > x <- read.table(textConnection("Datetimevalue + 2011062011:18:007 + 2011062011:39:009 + 2011062111:41:008 + 2011062111:40:006 + 2011062214:05:008 + 2011062214:06:00

[R] Residual deviance in glm

2008-01-10 Thread Gang Chen
I'm running a categorical data analysis with a two-way design of nominal by ordinal structure like the Political Ideology Example (Table 9.5) in Agresti's book Categorical Data Analysis. The nominal variable is Method while the ordinal variable is Quality (Bad, Moderate, Good, Excellent). I

[R] Extracting last time value

2008-01-10 Thread Kondamani, Arjun (GMI - NY Corporate Bonds)
I have a dataframe as follows: Datetimevalue 2011062011:18:007 2011062011:39:009 2011062111:41:008 2011062111:40:006 2011062214:05:008 2011062214:06:006 For every date, I want t

[R] Plotting a smooth curve through points

2008-01-10 Thread Anderson, Walter
I have a set of Monte Carlo simulation results for a Quick Response Freight Model that I am currently preparing summary graphs. I want to show three things on the graphs; the model forecast, an approximate exponential best fit line, and a smooth line through mean+sd and mean-sd data points. As yo

[R] Omnibus main effects in summary.lme?

2008-01-10 Thread Adam D. I. Kramer
Hello, I've been running some HLMs using the lme function quite happily; it does what I want and I'm pretty sure I understand it. The issue is that I'm currently trying to estimate a model with a 14-level "nusiance" factor as an independent variable...which makes the output quite

[R] Cycle Regression Analysis in R?

2008-01-10 Thread Carson Farmer
Hello R community, Does anyone know of a package that will perform cycle regression analysis? I have searched the R-help archives etc. but have come up with nothing so far. If I am unable to find an existing R package to do so, is there anyone familiar with fitting sine functions to data. My p

Re: [R] A Contrast Question

2008-01-10 Thread Peter Dalgaard
(Ted Harding) wrote: > Hi Folks, > This should have been simple to answer, > but despite much chasing I don't seem able > to catch this particular mouse! > > Basically (somewhat simplified): I have a > binary response variable Y (0/1) and a > 2-level factor A (0/1). > > I want to assign a contrast

[R] A Contrast Question

2008-01-10 Thread Ted Harding
Hi Folks, This should have been simple to answer, but despite much chasing I don't seem able to catch this particular mouse! Basically (somewhat simplified): I have a binary response variable Y (0/1) and a 2-level factor A (0/1). I want to assign a contrast to A such that, when I run summary(glm

Re: [R] pairwise cross tabulation tables

2008-01-10 Thread AndyZon
Thank you so much, Chuck! This is brilliant, I just tried some dichotomous variables, it was really fast. Most categorical variables I am interested in are 3 levels, they are actually SNPs, I want to look at their interactions. My question is: after generating 0-1 codings, like 00, 01, 10, how s

Re: [R] question regarding kpss tests from urca, uroot and tseries packages

2008-01-10 Thread markleeds
>From: eugen pircalabelu <[EMAIL PROTECTED]> >Date: 2008/01/10 Thu PM 02:48:32 CST >To: R-help <[EMAIL PROTECTED]> >Subject: [R] question regarding kpss tests from urca, uroot and tseries >packages Schwert has an algorithm for deciding on the number of lags to use in an adf test but i can't say

[R] question regarding kpss tests from urca, uroot and tseries packages

2008-01-10 Thread eugen pircalabelu
Hi R users! I've come across using kpss tests for time series analysis and i have a question that troubles me since i don't have much experience with time series and the mathematical part underlining it. x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279,

Re: [R] general linear hypothesis glht() to work with lme()

2008-01-10 Thread array chip
Hi I updated my R version to 2.6.1, now glht() worked with lme() from nlmn package. thanks for responding. Now I have a questions: can I use the contrast() from contrast package to do the comparison that I am interested in? The comparison that I am intersted is to test the linear relationship amo

Re: [R] `[.data.frame`(df3, , -2) and NA columns

2008-01-10 Thread Tony Plate
Assigning a name vector to a dataframe that is shorter than the number of columns results in some columns having NA values for their names. "[.data.frame" has the following code in it: cols <- names(x) ... if (any(is.na(cols))) stop("undefined columns sele

Re: [R] Writing custom functions for rpart

2008-01-10 Thread R Help
I'm having a similar problem, and after looking through the test directory I believe I have it down pat, but would someone mind telling me if I have it correctly, and if not, what I need to fix? eval: this seems to be the evaluation of a node, used for for labeling it with whatever decision method

Re: [R] OS X binary: 32 or 64-bit?

2008-01-10 Thread Charles C. Berry
Robert, If you downloaded the binary from CRAN, it uses a 32-bit address space. Here are some links that may help with getting a 64-bit build: http://article.gmane.org/gmane.comp.lang.r.mac/2226 http://r.research.att.com/ Also, there is a list for MACs https://stat.eth

[R] Plotting Long Time Series

2008-01-10 Thread Vishal Belsare
I have some time series data which I am interested in plotting. The standard function for plotting ts, or zoo objects work fine, however, I am curious about plotting long time series with a dynamic aspect ratio. My point is this, say I have a time series of 4000 observations, so the image width is

Re: [R] diff in a dataframe

2008-01-10 Thread Gabor Grothendieck
This is consistent with how matrices and ts series in R work: they all use x[,j] only. On Jan 10, 2008 1:08 PM, Vishal Belsare <[EMAIL PROTECTED]> wrote: > Thanks Henrique, Don and Gabor. I did come around to solving it by > using the zoo library. Very useful stuff that one for handling a bunch >

Re: [R] diff in a dataframe

2008-01-10 Thread Vishal Belsare
Thanks Henrique, Don and Gabor. I did come around to solving it by using the zoo library. Very useful stuff that one for handling a bunch of long irregular time series. Gabor, thanks for your present and previous responses. The quickref was indeed helpful. I do have another question regarding zoo

Re: [R] how to calculate time offset between timezones?

2008-01-10 Thread dankelley
Since designations such as "AST" are used in the data files that I'm parsing, I guess I'll make up a list of my own conversions, perhaps based on http://www.timeanddate.com/library/abbreviations/timezones/ (I cite this website not because I know it to be authoritative, but just in case it provi

Re: [R] leave one out cross validation

2008-01-10 Thread Prof Brian Ripley
Please look at ?lm.influence -- this does all the work for you. On Thu, 10 Jan 2008, Anu Swatantran wrote: > Hi > > I am trying to validate my regression results using the leave one out crosss > validation method. Is any script available in R to use this method for a > linear regression equation?

[R] OS X binary: 32 or 64-bit?

2008-01-10 Thread Robert Terwilliger
Dear R Experts, I am using R.app (the Mac OS X binary) for neuroimage analysis, so I am loading in some large image files. I get the following error in the middle of my script: > source("3dLME.R") Read 1 record Read 1 record Read 1 record Read 1 record Read 1 record Error: cannot allocate vector

[R] question regarding hypothesis testing in contingency tables

2008-01-10 Thread eugen pircalabelu
Hi R-users! I have the following example: a<-data.frame(cat=c(5,10,15), dog=c(5,10, 15), mouse=c(10,10,20)) b<-data.frame(cat=c(15,10,5), dog=c(15, 10, 5), mouse=c(20,10,10)) rownames(b)<-c("scared", "happy", "sad") rownames(a)<-c("scared", "happy", "sad") Let's say that a and b are 2 contingency

[R] leave one out cross validation

2008-01-10 Thread Anu Swatantran
Hi I am trying to validate my regression results using the leave one out crosss validation method. Is any script available in R to use this method for a linear regression equation? Both R or SPLUS would do. any clues on how to write the script will also help. Thanks a lot, A [[alternativ

Re: [R] 4 dimensional graphics

2008-01-10 Thread jiho
On 2008-January-10 , at 17:41 , Petr PIKAL wrote: > Thank you > > Basically I have a rectangular space (like an aquarium) in which I > made > some analysis. > I can make > > image(lat, long, value) for each height but what I dream about is to > make > something like scatterplot3d(lat, long, h

Re: [R] unsupervised random forest classification

2008-01-10 Thread Liaw, Andy
Random Forests do not do clustering. You need to take the proximity matrix and feed it to algorithms of your choice for that. Best, Andy From: [EMAIL PROTECTED] > > Friends, > I would like to use Random Forest in unsupervised mode to > classify data that are > from a mixture of categorical an

Re: [R] 4 dimensional graphics

2008-01-10 Thread Greg Snow
Petr, You may want to look at the rgl package for doing the 3d plot, then code the points using different sizes/colors for your value. However, when looking at a 3 dimensional plot that has been projected onto 2 dimensions (computer screen or paper), the ability to really understand what you are

[R] Error: Bad value

2008-01-10 Thread Sven Garbade
Hi list, from time to time I got an "Error: bad value" and must restart R. The mail archives suggests memory corruption, but I do not run "special" C code, only base R calls. I use R release 2.6.22 on a Debian GNU/Linux (version "testing") on a i686 machine with 2 GB RAM. R was installed with d

Re: [R] diff in a dataframe

2008-01-10 Thread Gabor Grothendieck
Represent this as a time series. Using the zoo package: > library(zoo) > z <- zoo(cbind(price_g = c(0.34, 0.36), price_s = c(0.56, 0.76)), > as.Date(c("2000-01-01", "2000-01-05"))) > diff(log(z)) price_g price_s 2000-01-05 0.05715841 0.3053816 > diff(log(z), na.pad = TRUE)

Re: [R] diff in a dataframe

2008-01-10 Thread anna freni sterrantino
Hello Vishal Maybe this is what you want? x<-c(0.34,0.36,3) y<-c(0.56,0.76,4) a<-data.frame(x,y) >a xy 1 0.34 0.56 2 0.36 0.76 3 3.00 4.00 a$diff<-(log(a$x)) a$diff2<-(log(a$y)) >a xy diff diff2 1 0.34 0.56 -1.078810 -0.5798185 2 0.36 0.76 -1.021651 -0.2744368 3 3.00

Re: [R] Error on distance matrix

2008-01-10 Thread Marc
Thanks for your help. It is working now. Marc. Gavin Simpson wrote: > On Thu, 2008-01-10 at 10:48 +, Marc Moragues wrote: > >> Hi, >> >> I am trying to calculate a distance matrix on a binary data frame using >> dist.binary() {ade4}. This is the code I run and the error I get: >> >> >>

Re: [R] 4 dimensional graphics

2008-01-10 Thread Petr PIKAL
Thank you Basically I have a rectangular space (like an aquarium) in which I made some analysis. I can make image(lat, long, value) for each height but what I dream about is to make something like scatterplot3d(lat, long, height) with points set according to a value. Up to now i can do scatt

Re: [R] 4 dimensional graphics

2008-01-10 Thread Alberto Monteiro
Petr PIKAL wrote: > > I want to display 4 dimensional space by some suitable way. > You probably need an LSD monitor :-) > I searched CRAN and found miscellaneous 3 dim graphics > packages which I maybe can modify but anyway I am open > to any hint how to > efficiently display data like: > > lo

Re: [R] Subsetting array (NA:s)

2008-01-10 Thread Marc Schwartz
Lauri Nikkinen wrote: > Hi R-users, > > I have an array similar to this: > > tmp <- array(1:6, c(2,3,3)) > n1 <- c("one", "two") > n2 <- c("three", "four", "five") > n3 <- c("six", "seven", "eight") > dimnames(tmp) <- list(n1, n2, n3) > tmp[1,,1] <- NA > tmp[1,3,2] <- NA > tmp[2,,3] <- NA > tmp >

[R] Identifying x, y coordinates from periodogram generated by spec.pgram

2008-01-10 Thread stephen sefick
I would like to identify the x,y coordinates in a plot of the periodogram spec.pgram(a, taper=0, log="no") I have tried the identify() function to no avail I have also tried to pass the periodogram over to rggobi e<-plot.spec(spec.pgram(a, taper=0, log="no")) ggobi(e) to no avail I need to id

Re: [R] question on how to extract values from atomic vector

2008-01-10 Thread Jack Luo
Thanks a lot! It works. Jack On Jan 9, 2008 10:35 PM, jim holtman <[EMAIL PROTECTED]> wrote: > There are a couple of ways of getting the values: > > > x <- cbind(a=1:3, pi=pi) # simple matrix w/ dimnames > > attributes(x) > $dim > [1] 3 2 > > $dimnames > $dimnames[[1]] > NULL > > $dimnames[[2]]

Re: [R] load tcltk

2008-01-10 Thread Lorenz Gabmann
Problem solved: X11 from Mac OS X Tiger CD has to be installed then tcltk runs properly. Thanks! Lorenz Am 10.01.2008 um 15:38 schrieb Alberto Monteiro: > > Lorenz Gabmann wrote: >> >> Can someone help me with the following problem? >> > Let me see if I can guess what are the error messages. >

Re: [R] as.date - german date format

2008-01-10 Thread Konga
Hi Terry, thanks for the information, I should know this in advance but I´m not used to case sensitive syntax - thanks Delphi... regards, Eric Will Terry Therneau wrote: > >>> Konga wrote: Hi, I have a data.frame with the following variable: $ xx :Factor w/ 66 levels "0

Re: [R] how to calculate time offset between timezones?

2008-01-10 Thread Prof Brian Ripley
On Thu, 10 Jan 2008, dankelley wrote: > > #QUESTION# > > Is there a way to calculate the offset between timezones, e.g. from "AST" to > "GMT"? > > #DETAILS# > > I am working with data files that use local time, and that indicate the > timezone by e.g. "AST" (Atlantic Standard Time, my local time)

Re: [R] how to calculate time offset between timezones?

2008-01-10 Thread Armstrong, Whit
well, strangely, ISOdatetime(2008,1,1,1,0,0,tz="AST") creates a UTC timezone date on my system: > ISOdatetime(2008,1,1,1,0,0,tz="AST") [1] "2008-01-01 01:00:00 UTC" > and if you compare the numeric values of the UTC datetime and the GMT datetime, the are definitely the same: > as.numeric(ISOda

Re: [R] diff in a dataframe

2008-01-10 Thread Don MacQueen
So, what's the easiest way to add a column to a dataframe? Just do it. Here is a really simple example to illustrate: > foo <- data.frame(x=1:4,y=2:5) > foo x y 1 1 2 2 2 3 3 3 4 4 4 5 > foo$z <- c(NA,diff(foo$x)) > foo x y z 1 1 2 NA 2 2 3 1 3 3 4 1 4 4 5 1 Solutions using apply

Re: [R] Changing language of messages (was load tcltk)

2008-01-10 Thread Prof Brian Ripley
On Thu, 10 Jan 2008, Alberto Monteiro wrote: [...] > PS: is there any way to temporarily change the error messages > from one language to another? When I get error messages > in Portuguese it's a pain to translate them precisely to > something that I can google or RSiteSearch to know what it mean

Re: [R] plotting help needed

2008-01-10 Thread Don MacQueen
Try plot(1:2,all,type="b",col="blue",xlim=c(1,4),ylim=c(20,150),xaxt='n',yaxt='n') lines(2:3,all2,type='o',col='yellow') I don't really know what you mean by "4 intersects, just named by characters". By "no scale please" I assume you mean tick marks and tick mark labels on the axes, and that's

Re: [R] question about subset and join ts object(s)

2008-01-10 Thread tom soyer
Thanks Achim. Data manipulation in zoo and coerce back to ts. Sounds good! On 1/10/08, Achim Zeileis <[EMAIL PROTECTED]> wrote: > > On Thu, 10 Jan 2008, tom soyer wrote: > > > Hi, > > > > I have two questions about ts. > > > > (1) How do I subset a ts object and still preserve the time index? for

Re: [R] An "R is slow"-article

2008-01-10 Thread Tom Backer Johnsen
Gustaf Rydevik wrote: > Hi all, > > Reading the wikipedia page on R, I stumbled across the following: > http://fluff.info/blog/arch/0172.htm > > It does seem interesting that the C execution is that much slower from > R than from a native C program. Could any of the more technically > knowled

[R] Subsetting array (NA:s)

2008-01-10 Thread Lauri Nikkinen
Hi R-users, I have an array similar to this: tmp <- array(1:6, c(2,3,3)) n1 <- c("one", "two") n2 <- c("three", "four", "five") n3 <- c("six", "seven", "eight") dimnames(tmp) <- list(n1, n2, n3) tmp[1,,1] <- NA tmp[1,3,2] <- NA tmp[2,,3] <- NA tmp How to subset !is.na(x) rows resulting , , six

[R] unsupervised random forest classification

2008-01-10 Thread helen . mills
Friends, I would like to use Random Forest in unsupervised mode to classify data that are from a mixture of categorical and continuous variables. The examples I find online use only continuous data, followed by an mds plot. Is it correct to use RF to obtain proximities for my mixed dataset and then

[R] how to calculate time offset between timezones?

2008-01-10 Thread dankelley
#QUESTION# Is there a way to calculate the offset between timezones, e.g. from "AST" to "GMT"? #DETAILS# I am working with data files that use local time, and that indicate the timezone by e.g. "AST" (Atlantic Standard Time, my local time) in a header. I was guessing that > ISOdatetime(2008,

Re: [R] 4 dimensional graphics

2008-01-10 Thread hadley wickham
On Jan 10, 2008 8:36 AM, Petr PIKAL <[EMAIL PROTECTED]> wrote: > Dear all > > I want to display 4 dimensional space by some suitable way. I searched > CRAN and found miscellaneous 3 dim graphics packages which I maybe can > modify but anyway I am open to any hint how to efficiently display data > l

Re: [R] help with R

2008-01-10 Thread John Kane
--- [EMAIL PROTECTED] wrote: > (2) Where do I get to see the complete library of > packages offered by R? > In particular, I am interested in quantile > regression and logistic > regression. R webpage > CRAN > Packages link. Also have a look at the Task Views on the CRAN page. __

[R] Fwd: multinomial regression for clustered data

2008-01-10 Thread Niccolò Bassani
Hello dear R-users, does any of you know a way to perform a multinomial regression with clustered data (i.e. repeated measurements)? I made the first analysis with Stata option vce cluster in the mlogit command but was looking for a similar functionality in R too... thanks all! niccolò [[a

[R] distance between set of coordinates

2008-01-10 Thread Mafalda Viana
Hi! I have a set of 800 coordinates pairs and need to find the distance between each pair. I am able to find the distance between one pair using: <-geodist(LatStart,LongStart,LatEnd,LongEnd,units="km")) but doing one by one is nor feasible. Can you help me? Thanks, Mafalda --

Re: [R] `[.data.frame`(df3, , -2) and NA columns

2008-01-10 Thread John Kane
I don't know why this is happening but it has nothing to do with a negative index df[,-2] has not changed df. --- Dieter Menne <[EMAIL PROTECTED]> wrote: > Dear baseRs, > > I recently made a mistake when renaming data frame > columns, accidentally > creating an NA column. I found the following

Re: [R] how can i add a package that i want to use permanent

2008-01-10 Thread Bos, Roger
Another option is to modify the Rprofile.site file in the etc folder of your R installation. Then you don't have to save the workspace. You add as much startup scripting as you want to the Rprofile.site file. HTH, Roger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECT

Re: [R] load tcltk

2008-01-10 Thread Alberto Monteiro
Lorenz Gabmann wrote: > > Can someone help me with the following problem? > Let me see if I can guess what are the error messages. > Lade nötiges Paket: tcltk > Load (???) package: tcltk > Loading Tcl/Tk interface ... Error in dyn.load(file, ...) : >kann shared library > '/Library/Framewor

[R] 4 dimensional graphics

2008-01-10 Thread Petr PIKAL
Dear all I want to display 4 dimensional space by some suitable way. I searched CRAN and found miscellaneous 3 dim graphics packages which I maybe can modify but anyway I am open to any hint how to efficiently display data like: longitude, latitude, height, value Thank you Petr Pikal [EMAIL

Re: [R] question about subset and join ts object(s)

2008-01-10 Thread Achim Zeileis
On Thu, 10 Jan 2008, tom soyer wrote: > Hi, > > I have two questions about ts. > > (1) How do I subset a ts object and still preserve the time index? for > example: > > > x=ts(1:10, frequency = 4, start = c(1959, 2)) # the ts object > > x > Qtr1 Qtr2 Qtr3 Qtr4 > 1959 123 > 196

[R] load tcltk

2008-01-10 Thread Lorenz Gabmann
Hi All! Can someone help me with the following problem? Lade nötiges Paket: tcltk Loading Tcl/Tk interface ... Error in dyn.load(file, ...) : kann shared library '/Library/Frameworks/R.framework/Versions/2.6/ Resources/library/tcltk/libs/ppc/tcltk.so' nicht laden: dlopen(/Library/Frameworks

Re: [R] Error on distance matrix

2008-01-10 Thread Gavin Simpson
On Thu, 2008-01-10 at 10:48 +, Marc Moragues wrote: > Hi, > > I am trying to calculate a distance matrix on a binary data frame using > dist.binary() {ade4}. This is the code I run and the error I get: > > > sjlc.dist <- dist.binary(as.data.frame(data), method=2) #D = (a+d) / > (a+b+c+d) > Er

Re: [R] how can i add a package that i want to use permanent

2008-01-10 Thread Charles Annis, P.E.
?.First This will explain things. You make a hidden program called .First like this .First <- function() { library(lattice) } and then save your workspace. The next time you open R, lattice will load automatically. Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3

Re: [R] Error on distance matrix

2008-01-10 Thread Jari Oksanen
Jari Oksanen oulu.fi> writes: > > Marc Moragues scri.ac.uk> writes: > > > > > Hi, > > > > I am trying to calculate a distance matrix on a binary data frame using > > dist.binary() {ade4}. This is the code I run and the error I get: > > > > > sjlc.dist <- dist.binary(as.data.frame(data), met

[R] question about subset and join ts object(s)

2008-01-10 Thread tom soyer
Hi, I have two questions about ts. (1) How do I subset a ts object and still preserve the time index? for example: > x=ts(1:10, frequency = 4, start = c(1959, 2)) # the ts object > x Qtr1 Qtr2 Qtr3 Qtr4 1959 123 19604567 196189 10 I don't want the 1

[R] [R-pkgs] PwrGSD

2008-01-10 Thread Grant Izmirlian
Hello List: Please find uploaded to CRAN a new package, PwrGSD The package is intended for the design and analysis of group sequential trials There are two main functions, (1) GrpSeqBnds: computes group sequential stopping boundaries for interim analysis of a sequential trial based upon a

[R] [R-pkgs] prob package: elementary probability on finite sample spaces

2008-01-10 Thread G. Jay Kerns
Dear R Community, I am pleased to announce the beta-release of the prob package. The source code is now on CRAN, and binaries should be generated there before long. In the meantime, you can get it with install.packages("prob", repos = "http://r-forge.r-project.org";) The prob package gives a fr

[R] formatting ftable

2008-01-10 Thread Monica Pisica
Hi all, I am using ftable and xtabs to get total counts of species per location per hight category (understory, midstory and overstory). i can save the data in almost the format i would like it to be - but not quite. Here it is an example of what i get: x - is a dataframe with the following col

Re: [R] Error on distance matrix

2008-01-10 Thread Jari Oksanen
Marc Moragues scri.ac.uk> writes: > > Hi, > > I am trying to calculate a distance matrix on a binary data frame using > dist.binary() {ade4}. This is the code I run and the error I get: > > > sjlc.dist <- dist.binary(as.data.frame(data), method=2) #D = (a+d) / > (a+b+c+d) > Error in if (any(df

[R] as.date - german date format

2008-01-10 Thread Terry Therneau
>> Konga wrote: >>> Hi, >>> >>> I have a data.frame with the following variable: >>> $ xx :Factor w/ 66 levels "01.02.2006","01.03.2006",..: 3 3 3 3 3 3 3 3 >>> 3 >>> >>> now I?d like use "as.date" on it - if I comprehend the instructions on >>> http://finzi.psych.upenn.edu/R/library/base/html/fo

Re: [R] data.frame manipulation: Unbinding strings in a row

2008-01-10 Thread jim holtman
here is a quick hack: > x <- read.table(textConnection("ID ShopItems + ID1 A1 item1,item2,item3 + ID2 A2 item4,item5 + ID3 A1 item1,item3,item4"), header=TRUE) > y <- lapply(1:nrow(x), function(.row){ + .items <- strsplit(as.character(x$Items[.row]), ',')[[1

Re: [R] labels to values

2008-01-10 Thread Matthias Wendel
Hello, Henrique, thank you for the tip, but it was not quite what has been desired: > d[, 'Y6'] [1] 6 3 8 11 8 9 6 8 3 5 10 15 NA 9 8 3 8 16 6 6 NA 10 5 2 7 7 6 16 7 15 7 10 12 8 7 12 12 16 7 6 8 8 15 6 NA 8 99 7 12 8 9 16 7 16 8 7 7 1 15 [60] 12

[R] data.frame manipulation: Unbinding strings in a row

2008-01-10 Thread francogrex
Hi all, I have a data.frame I received with data that look like this (comma separated strings in last row): ID ShopItems ID1 A1 item1, item2, item3 ID2 A2 item4, item5 ID3 A1 item1, item3, item4 But I would like to unbind the strings in col(2) items so that

[R] `[.data.frame`(df3, , -2) and NA col umns

2008-01-10 Thread Dieter Menne
Dear baseRs, I recently made a mistake when renaming data frame columns, accidentally creating an NA column. I found the following strange behavior when negative indexes are used. Can anyone explain what happens here. No "workarounds" required, just curious. Dieter Version: Windows, R version

  1   2   >