Re: [R] building packages on Windows

2007-10-01 Thread Gabor Grothendieck
This happened with older versions of R but its fixed in the more recent R 2.6.0 versions. On 10/1/07, Edna Bell <[EMAIL PROTECTED]> wrote: > Hi again. > > I'm sure that this is really simple. > > I'm trying to build a package on a Windows Vista machine. I use > Rcmd build --binary test > > but I

[R] building packages on Windows

2007-10-01 Thread Edna Bell
Hi again. I'm sure that this is really simple. I'm trying to build a package on a Windows Vista machine. I use Rcmd build --binary test but I get the "Please set TMPDIR to a valid temporary directory" I tried TMPDIR=c:\temp but to no avail. Please help. thanks __

Re: [R] plot graph with error bars trouble

2007-10-01 Thread Marcelo Laia
Dear Ben Bolker, Hadley Wickham, Jim Lemon and all I'm very thank to you! I tried all suggestion and I wiil comments it next. Ben Bolker suggestion is very nice and I did a little changes in col=c(1,1),pch=c(1,16) and I include a xlab, ylab and a legend(5, 4, c("02H02","306"), pch = c(1,16), co

[R] problems with glm

2007-10-01 Thread stephenc
I am having a couple of problems someone may be able to cast some light on. Question 1: I am making a logistic model but when i do this: glm.model = glm(as.factor(form$finished) ~ ., family=binomial, data=form[1:15,]) I get this: Error in model.frame(formula, rownames, variables, varname

[R] convert a tar.gz to a Windows zip

2007-10-01 Thread Edna Bell
Dear R Gurus; Is there a simple way to convert a Linux produced tar.gz file (a package) to a Windows binary zip package, please? Thanks in advance __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the post

Re: [R] Ordering of names on X- and Y-axis

2007-10-01 Thread Matthew Keller
Hi Michael, This type of thing is pretty simple to do in R (take a look at ?plot and then ?par - very useful... necessary actually... for graphing). One way to do it is to use the "axis" function: my.labels <- c("a","b","c","d","e","f","g","h","i","j") plot(1:10,rnorm(10),axes=F,ylab="whatev",xl

Re: [R] Ordering of names on X- and Y-axis

2007-10-01 Thread jim holtman
Try ordering the factors: x <- read.table(textConnection("SemType, Length GeoLocation, Sentence GeneralInfo, Paragraphs GeneralInfo, Paragraphs GeneralInfo, Sentence GeneralInfo, Paragraphs NatLang, Phrase Advice, Article GeneralInfo, Article Advice, Article Resource, Sentence"), header=TRUE) # o

Re: [R] RODBC and Oracle

2007-10-01 Thread Thomas Harte
jack, i had a similar problem with Oracle at work. i don't have access to my yahoo email account at work but i'll take a stab at directing you to the solution here and now FWIW. open Data Sources then click on the System DSN (i think?) tab. highlight Oracle 10g and hit configure. there's a defaul

[R] Ordering of names on X- and Y-axis

2007-10-01 Thread Michael Kaisser
Hi, I am new to R. I have a bit of data looking like this: SemType, Length GeoLocation, Sentence GeneralInfo, Paragraphs GeneralInfo, Paragraphs GeneralInfo, Sentence GeneralInfo, Paragraphs NatLang, Phrase Advice, Article GeneralInfo Advice, Article Resource, Sentence ... (roughly 40,000 lin

Re: [R] merging, interpolating two simulataneous time series

2007-10-01 Thread Steve Powers
No, it's a little more complex. I do not want to interpolate var2 as a function of time. I want interpolate var2 as a function of var1, but in an iterative chronological sequence that does require matching times first. So the two var1 values for the interpolation are chosen such that var1 times

[R] Dealing with multiple observation per time point in time series

2007-10-01 Thread adschai
Hi - I think this is a classic question. I'm using fCalendar time series. I have observations that have multiple instances at second level. Is there anyway to handle this in R? I can inject millisecond unit arbitrarily inside there but I don't feel it's a correct way to do. Any ideas would be re

Re: [R] R CMD build not excluding .svn

2007-10-01 Thread Michael Lawrence
Thanks for the help. I looked into it further and noticed that the src/RGtk2 directory was actually a symlink. The perl find() does not follow symlinks, as invoked, but tar, as invoked, does. Thus, when the tar/untar effectively copies the files, there is a mismatch with the exclude paths. Is ther

Re: [R] Diff. between time series and stochastic process

2007-10-01 Thread Rolf Turner
On 2/10/2007, at 12:52 AM, Duncan Murdoch wrote: > On 10/1/2007 8:39 AM, Megh Dal wrote: >> Hi, Can anyone give me a good explanation about what is the >> difference between Stochastic Process and Time Series? In my >> knowledge Time series process is one type of Stochastic process. >> Am I

Re: [R] tricky problem with "if" function -

2007-10-01 Thread Tim Calkins
This is confusing. try this. number.runs <- 1 if (number.runs == 1) run.number <- 1 if (number.runs > 1) { for (i in 1:number.runs) { print(i) } } perhaps the easiest way to print 1 and then 1 2 would be for (i in 1:number.runs) { print(i) } but i'm not sure that is what you want. i

[R] tricky problem with "if" function -

2007-10-01 Thread Matthew Keller
Hi all, This question involves using a "for" loop to make a "decision" in a script. I've written a rather intricate script, and near the start of it, I want it either to do a loop (if a variable called "number.runs" > 1) or not do a loop (if "number.runs" is 1). This is probably trivial but I can

Re: [R] sub or gsub function

2007-10-01 Thread Tim Calkins
see ?toupper On 10/2/07, Michael Sun <[EMAIL PROTECTED]> wrote: > Dear list, > > x <- c("year","date","time","starting price","end price","max","min") > > how can I change all these lower-case letter to capital letter? > > how can I use sub or gsub function??? > > sub("year","YEAR",x) > will only

[R] sub or gsub function

2007-10-01 Thread Michael Sun
Dear list, x <- c("year","date","time","starting price","end price","max","min") how can I change all these lower-case letter to capital letter? how can I use sub or gsub function??? sub("year","YEAR",x) will only change one each time... Appreciate for any helps Cheers mam [[alterna

Re: [R] Version control and R package development (was: plot graph with error bars trouble)

2007-10-01 Thread jiho
On 2007-October-01 , at 23:13 , jiho wrote: > [...] > The end result looks like this: > http://cbetm.univ-perp.fr/irisson/svn/ > to see what an inner view looks like: > http://cbetm.univ-perp.fr/irisson/svn/distribution_data/tetiaroa/ > trunk/data/?rev=0&sc=1 > You can play around wi

[R] Version control and R package development (was: plot graph with error bars trouble)

2007-10-01 Thread jiho
First a quick summary of the beginning of an interesting discussion (not all indented correctly but in the correct order at least): > From: "Gabor Grothendieck" <[EMAIL PROTECTED]> > On 9/30/07, hadley wickham <[EMAIL PROTECTED]> wrote: >> On 9/30/07, jiho <[EMAIL PROTECTED]> wrote: >>> [...] >>

Re: [R] lpSolve doesn't compile because of a malloc.h error

2007-10-01 Thread Prof Brian Ripley
Isn't the message rather clear? You need to edit the package sources and replace malloc.h by stdlib.h (or remove the former if the latter is already there, as in this case). And BTW, 'gcc-4.2.2' is unreleased, but imminent (as R 2.6.0 is). On Mon, 1 Oct 2007, vittorio wrote: > Under freebsd 6

[R] lpSolve doesn't compile because of a malloc.h error

2007-10-01 Thread vittorio
Under freebsd 6.2-p7 i386, R 2.5.1,gcc-4.2.2 I'm unable to compile package lpSolve because: hpbsd# R CMD INSTALL lpSolve_5.5.8.tar.gz * Installing to library '/usr/local/lib/R/library' * Installing *source* package 'lpSolve' ... ** libs cc -std=gnu99 -I/usr/local/lib/R/include -I/usr/local/lib/R/

Re: [R] xyplot

2007-10-01 Thread Thomas Lumley
It's a FAQ (7.22) -thomas On Mon, 1 Oct 2007, Samor Gandhi wrote: > Hello, > I am calling the following code with the loop! It makes 3 out graphs but > empty! COuld you help me plase on that? Thank you in advance > > z1 <- dbConnect(MyData, "something", "A1", "A2") > for (tt in c("x

Re: [R] orientlib

2007-10-01 Thread Duncan Murdoch
On 10/1/2007 1:27 PM, Abu Naser wrote: > Hi Duncan, > > Let me explain again about what I want to do. > > I have a three coordinates (x1,y1,z1),(x2,y2,z2), and (x3,y3,z3). I want to > use (x1,y1,z1) > and (x2,y2,z2) for defining axis of rotation Okay so far: the axis is (x1-x2, y1-y2, z1-z2)

Re: [R] how to plot a graph with different pch

2007-10-01 Thread hadley wickham
You might find it easier to use the ggplot2 (http://had.co.nz/ggplot2) package: install.packages("ggplot2") library(ggplot2) qplot(mpg, wt, data=mtcars, colour=factor(cyl)) # or qplot(mpg, wt, data=mtcars, shape=factor(cyl)) ggplot2 takes care of mapping the variable values to colours and sizes,

Re: [R] how to plot a graph with different pch

2007-10-01 Thread Charles C. Berry
On Mon, 1 Oct 2007, Letticia Ramlal wrote: > I am trying to plot a graph but the points on the graph should be > different symbols and colors. It should represent what is in the legend. > I tried using the points command but this does not work. WHAT "does not work" ??? points() puts points on a

Re: [R] how to plot a graph with different pch

2007-10-01 Thread Henrique Dallazuanna
Perhaps plot(mtcars$wt,mtcars$mpg,xlab= "Weight(lbs/1000)", ylab="Miles per Gallon", pch=c(3,17,19), col=c("blue", "red", "purple4")) On 01/10/2007, Letticia Ramlal <[EMAIL PROTECTED]> wrote: > > I am trying to plot a graph but the points on the graph should be > different symbols and colors. I

Re: [R] Shapiro-Welch W value interpretation

2007-10-01 Thread Omar Baqueiro
Thanks for your answer, I will try to make myself clearer and do not worry to point out what I've got wrong. > There is some confusion in your query. > First, how do you know that your data are indeed normally distributed? In this specific case I *know* that my data is normally distributed becaus

Re: [R] Idioms for empty dataframes

2007-10-01 Thread Thomas Lumley
On Mon, 1 Oct 2007, Ranjan Bagchi wrote: > I'm fairly new to R, coming from a programming background -- it's quite > nice to work with dataframes, though, as opposed to explicit iteration. > > One thing I've found, which is surprising is that zero-length dataframes > seem to cause errors: It's no

[R] how to plot a graph with different pch

2007-10-01 Thread Letticia Ramlal
I am trying to plot a graph but the points on the graph should be different symbols and colors. It should represent what is in the legend. I tried using the points command but this does not work. Is there another command in R that would allow me to use different symbols and colors for the points?

Re: [R] saving and loading complex objects

2007-10-01 Thread Gabor Grothendieck
formulas have environments and since most model structures include the formula the entire environment of the formula will be pulled in. On 10/1/07, hadley wickham <[EMAIL PROTECTED]> wrote: > Why are you using eval and not: > > path <- paste( "C://Program Files//R//R-2.5.1//arima//",nn, sep="") >

Re: [R] plot graph with error bars trouble

2007-10-01 Thread Gabor Grothendieck
On 10/1/07, hadley wickham <[EMAIL PROTECTED]> wrote: > > These seem nearly identical to what you can get with R-Forge or with > > TortoiseSVN (and likely other svn clients too). Since any developer > > is likely to have an svn client a web interface more sophisticated than > > what is already ava

Re: [R] saving and loading complex objects

2007-10-01 Thread hadley wickham
Why are you using eval and not: path <- paste( "C://Program Files//R//R-2.5.1//arima//",nn, sep="") save(x, file = path) ? It's possible that extra environments are getting saved. (Someone with more knowledge of eval might suggest a better solution, but that might be a good place to start) Had

Re: [R] plot graph with error bars trouble

2007-10-01 Thread hadley wickham
> These seem nearly identical to what you can get with R-Forge or with > TortoiseSVN (and likely other svn clients too). Since any developer > is likely to have an svn client a web interface more sophisticated than > what is already available via the net has less utility than if this info were > n

[R] plots

2007-10-01 Thread amna khan
DEAR SIR How can I use ts.plt() function for time series of unequal length? Regards -- AMINA SHAHZADI Department of Statistics GC University Lahore, Pakistan. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://

[R] saving and loading complex objects

2007-10-01 Thread Lloyd Lubet
Dear Sir: When I try to save large and very complex recursive objects with some components containing models (such as the output from arima or lm ), the resulting file sizes increase by 4 meg per save. example directory: amex8 meg argentina 12 meg australia 16 m

Re: [R] Clustering techniques using R

2007-10-01 Thread Prof Brian Ripley
On Mon, 1 Oct 2007, Maura E Monville wrote: > Now that I've loaded a file into an R data.frame and played with > linear regression until I got a good model, my next step is clustering > using the coefficients of the regression model (I have many files) > Thanks to some R experts' guidelines I co

Re: [R] orientlib

2007-10-01 Thread Abu Naser
Hi Duncan, Let me explain again about what I want to do. I have a three coordinates (x1,y1,z1),(x2,y2,z2), and (x3,y3,z3). I want to use (x1,y1,z1) and (x2,y2,z2) for defining axis of rotation and then calculate the angle of rotation of a line from (x1,y1,z1) and (x3,y3,z3) around that axis o

[R] Clustering techniques using R

2007-10-01 Thread Maura E Monville
Now that I've loaded a file into an R data.frame and played with linear regression until I got a good model, my next step is clustering using the coefficients of the regression model (I have many files) Thanks to some R experts' guidelines I could find plenty of documentation on regression analys

Re: [R] Plot for two Series

2007-10-01 Thread Gabor Grothendieck
See ?ts.plot Also ?plot.zoo and ?xyplot.zoo in the zoo package. On 10/1/07, amna khan <[EMAIL PROTECTED]> wrote: > Dear Sir > > I want to plot the two time series having unequal number of observations, > say, one series having values from year 1960 to year 2000 and other having > values from year

Re: [R] plot graph with error bars trouble

2007-10-01 Thread Gabor Grothendieck
On 10/1/07, hadley wickham <[EMAIL PROTECTED]> wrote: > On 10/1/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > > On 10/1/07, hadley wickham <[EMAIL PROTECTED]> wrote: > > > > > > The biggest drawback (to me) to both google code and R-forge, is their > > > failure to offer a nice interface to b

[R] Plot for two Series

2007-10-01 Thread amna khan
Dear Sir I want to plot the two time series having unequal number of observations, say, one series having values from year 1960 to year 2000 and other having values from year 1975 to year 2000. How to plot them on a single graph sheet. Regards -- AMINA SHAHZADI Department of Statistics GC Univ

Re: [R] mean of subset of rows

2007-10-01 Thread Jeffrey Robert Spies
You were on the right track with the for loop, but often you can do the same thing looplessly (I know, it's not really a word) in R: If your data is like this: data<-data.frame(ID=rep(letters[1:4], 5), size=runif(20)) then apply either tapply(data$size, data$ID, mean) or aggregate(data$size

Re: [R] Odp: Idioms for empty dataframes

2007-10-01 Thread Ranjan Bagchi
Thanks... that's very helpful. Sorry about the typo too > t1 <- t[t$bob < 50, ,drop=F] > t1 [1] bob <0 rows> (or 0-length row.names) > order(t$bob) [1] 1 > On Mon, 1 Oct 2007, Petr PIKAL wrote: > Hi > > [EMAIL PROTECTED] napsal dne 01.10.2007 18:01:13: > >> I'm fairly new to R, coming from a pr

Re: [R] mean of subset of rows

2007-10-01 Thread John Kane
--- [EMAIL PROTECTED] wrote: > Dear list, > this must be an easy one: > > I have a data.frame of two columns, "ID" with four > different levels (A > to D) and numerical "size", and each of the 4 > different IDs is > repeated a > different number of times. I would like to get the > mean size

Re: [R] orientlib

2007-10-01 Thread Duncan Murdoch
On 10/1/2007 11:09 AM, Abu Naser wrote: > Hi Duncan, > Sorry, I forgot to mention. I meant (x1,y1,z1) and (x3,y3,z3). I'm afraid I still don't understand what angle you are trying to calculate. Duncan Murdoch > > With regards, > Abu > > > Date: Mon, 1 Oc

[R] function for special matrix design

2007-10-01 Thread kevinchang
Hi All, I am trying to make a matrix with a particular value for all the elements above and including diagonal and another particular value for all the elements below diagonal. Obviously , Matrix function in Matrix package does not work since it only allow filling by either row or column. So I am

[R] Clustering literature was Re: nonlinear regression

2007-10-01 Thread Petr PIKAL
Hi It is preferable to echo your posts to r-help, you usually get more answers and some definitelly superb to mine. It is also better to start a new mail if your question has nothing to do with original subject "Maura E Monville" <[EMAIL PROTECTED]> napsal dne 01.10.2007 17:44:43: > Unluckil

Re: [R] mean of subset of rows

2007-10-01 Thread joris . dewolf
data <- data.frame(ID = rep(letters[1:4],5),size=rnorm(20,0,1)) aggregate(data$size, by = list(data$ID),mean) <[EMAIL PROTECTED] e.ehu.es>

[R] Odp: Idioms for empty dataframes

2007-10-01 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 01.10.2007 18:01:13: > I'm fairly new to R, coming from a programming background -- it's quite > nice to work with dataframes, though, as opposed to explicit iteration. > > One thing I've found, which is surprising is that zero-length dataframes > seem to cause

[R] Interpretation of residual variance components and scale parameters in GLMMs

2007-10-01 Thread Emily H. DuVal
Dear R-listers, I am working with generalized linear mixed models to quantify the variance due to two nested random factors, but have hit a snag in the interpretation of variance components. Despite my best efforts with Venables & Ripley 2002, Fahrmeir & Tutz 2001, R-help archives, Google, and ot

[R] Idioms for empty dataframes

2007-10-01 Thread Ranjan Bagchi
I'm fairly new to R, coming from a programming background -- it's quite nice to work with dataframes, though, as opposed to explicit iteration. One thing I've found, which is surprising is that zero-length dataframes seem to cause errors: > t <- data.frame(bob=c(100)) > order(t$bob) [1] 1 > t1

[R] mean of subset of rows

2007-10-01 Thread darteta001
Dear list, this must be an easy one: I have a data.frame of two columns, "ID" with four different levels (A to D) and numerical "size", and each of the 4 different IDs is repeated a different number of times. I would like to get the mean size for each ID as another data.frame. I have tried th

Re: [R] plot graph with error bars trouble

2007-10-01 Thread hadley wickham
On 10/1/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > On 10/1/07, hadley wickham <[EMAIL PROTECTED]> wrote: > > > > The biggest drawback (to me) to both google code and R-forge, is their > > failure to offer a nice interface to browser the svn repository and > > view the timeline of changes.

Re: [R] Graphics and LaTeX documents with the same font - multiple y-axes

2007-10-01 Thread hadley wickham
It may not be misleading, but does it really help you understand the problem more than a scatterplot of rugosity vs. elevation would? Hadley On 10/1/07, Monica Pisica <[EMAIL PROTECTED]> wrote: > > Well, > > I agree that plots with multiple y-axes can be very confusing but i am > not sure

[R] Multiple Regression

2007-10-01 Thread Ultan
Hi, I have only recently started to use the R package and I was wondering how to perform a multiple regression on a set of data were there is no intercept and a condition that none of the coefficients can be less than zero. I would appreciate if anybody could help! Thanks Ultan -- View this me

[R] Re : how to create this graph?

2007-10-01 Thread justin bem
Try this plot.new() plot.window(xlim=c(0,10), ylim=c(0,15)) rect(1,0,2,10) rect(2,0,4,4.1,col="grey") rect(2,4.1,4,10,col="grey") segments(0, 0, 5, 0) segments(0, 10, 5, 10) text(3,3.5,"41") text(3,6, "51") text(4.5,3..5,"G/V") text(4.5,6, "U") Ju

Re: [R] Disentagling formulas

2007-10-01 Thread Eik Vettorazzi
maybe I didn't understand the problem, but you can do sth like formula<-y~X1+X2+X3 lm(update.formula(formula,f(Y)~.)) maybe ?all.vars or ?terms will help you too. hth. Rebecca Sela schrieb: > I am writing a program in which I would like to take in a formula, change the > response (Y) variable

Re: [R] center option of basehaz in survfit

2007-10-01 Thread Thomas Lumley
On Thu, 27 Sep 2007, David Koons wrote: > I have a very general question about what the centering option in > basehaz does to factors. (basehaz computes the baseline cumulative > hazard for a coxph object using the Breslow estimator). > > Lets say I'm interested in a survival model with two (di

[R] Error loading DCluster

2007-10-01 Thread Birgit Lemcke
When I try to load the DCluster package I get the following error message (R 2.5.1; PowerBook G4; Mac OS X 10.4.10) (Load needed package) Lade nötiges Paket: spdep (error: package spdep could not been loaded) Fehler: Paket 'spdep' konnte nicht geladen werden (additional warning message) Zusätzli

Re: [R] plot graph with error bars trouble

2007-10-01 Thread Gabor Grothendieck
On 10/1/07, hadley wickham <[EMAIL PROTECTED]> wrote: > > The biggest drawback (to me) to both google code and R-forge, is their > failure to offer a nice interface to browser the svn repository and > view the timeline of changes. I particularly like trac (e.g. > http://src.ggobi.org/) despite it'

Re: [R] Adding circles or ellipses to graphs

2007-10-01 Thread Greg Snow
For circles look at the symbols function. For ellipses see the ellipse package and the my.symbols function in the TeachingDemos package. -Original Message- From: "Friedman, Steven" <[EMAIL PROTECTED]> To: "r-help@r-project.org" Sent: 10/1/07 7:10 AM Subject: [R] Adding circles or ellips

Re: [R] orientlib

2007-10-01 Thread Abu Naser
Hi Duncan, Sorry, I forgot to mention. I meant (x1,y1,z1) and (x3,y3,z3). With regards, Abu > Date: Mon, 1 Oct 2007 10:48:12 -0400> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]> CC: r-help@r-project.org> Subject: Re: [R] orientlib>> On 10/1/2007 9:44 A

Re: [R] Graphics and LaTeX documents with the same font - multiple y-axes

2007-10-01 Thread Monica Pisica
Well, I agree that plots with multiple y-axes can be very confusing but i am not sure that always they are also misleading One time i was asked to make a 2 y-axes plot, one y axes was for elevation (heights in meters) and one for rugosity (values between 1 and 2 - unitless). The ide

Re: [R] plot graph with error bars trouble

2007-10-01 Thread hadley wickham
> >> It's something I have thought a little bit about, but I haven't made > >> much progress. Ideally, if it's something that I do for ggplot2, I > >> should do it for all my other R packages too. > > I don't see why doing it for *some* package(s) should entail > doing it for all of the

Re: [R] export csv

2007-10-01 Thread Petr PIKAL
Hi > Did you try reading the help pages or looking at the examples in the > help pages? Did you even make an attempt at using the command? If so, > did you get an error or unexpected results? > > I don't know how you can expect people to help you when you appear to > have made no effort yourself.

Re: [R] Disentagling formulas

2007-10-01 Thread peter salzman
Rebecca, you can keep formula as a character vector that can be easily manipulated formula1 <- "Y~X1+X2+X3" formula2 <- gsub("Y","newY",formula1) lm(as.formula(formula1),data=mydata) lm(as.formula(formula2),data=mydata) best, peter On 10/1/07, Rebecca Sela <[EMAIL PROTECTED]> wrote: > > I am

Re: [R] export csv

2007-10-01 Thread Eric Thompson
Did you try reading the help pages or looking at the examples in the help pages? Did you even make an attempt at using the command? If so, did you get an error or unexpected results? I don't know how you can expect people to help you when you appear to have made no effort yourself. On 10/1/07,

Re: [R] Ambiguities in vector

2007-10-01 Thread Birgit Lemcke
Hello James, first I have to thank you for your help but there are some things I don´t understand now. I am not sur if I understand what this example gives me back: ratings <- data.frame(id = c(1,2,3,4), att1 = c(1,1,0,1), att2 = c (1,0,0,1), att3 = c(0,1,1,1)) ratings id att1 att2 att3

Re: [R] orientlib

2007-10-01 Thread Duncan Murdoch
On 10/1/2007 9:44 AM, Abu Naser wrote: > Hi Duncan, > > Thanks for your replay. > I am confused as I dont have the quaternion backgorund. > > Let me explain what I am trying to do: > I have 3 coordinates in space (x1,y1,z1),(x2,y2,z2), and (x3,y3,z3) . Now, I > use (x1,y1,z1) and (x2,y2,z2) to

Re: [R] optimize() stuck in local plateau ?

2007-10-01 Thread Ravi Varadhan
The local maximum (x* = 0.106) found by the golden-section search does lie in the interval [x_1, x_2]. So, it is consistent with the explanation on the help page. --- Ravi Varadhan, Ph.D. Assistant Professor, The C

Re: [R] optimize() stuck in local plateau ?

2007-10-01 Thread Mike Lawrence
I may be misunderstanding, but my example seems to violate the "local minimum inside [x_1,x_2] will be found as solution, even when f is constant in there" rule since the search in my example continues on towards +1. On 1-Oct-07, at 10:40 AM, Ravi Varadhan wrote: > Please read the help for

Re: [R] xyplot

2007-10-01 Thread Wayne.W.Jones
When you call xyplot in a for loop you have to use the print command. For instance modifying the xyplot example: The following wont work: Depth <- equal.count(quakes$depth, number=8, overlap=.1) for(i in 1:3){xyplot(lat ~ long | Depth, data = quakes)} But this will: for(i in 1:3){print(x

Re: [R] export csv

2007-10-01 Thread John Kane
?write.csv --- elyakhlifi mustapha <[EMAIL PROTECTED]> wrote: > Hello, > I wanna know how to export a data frame to Excel. > For example I wanna export this data frame > > > M[1:5,] > var1 var2 distance > 1 41 42 0.2 > 2 41 43 0.304347826086957 > 3 41 46 0.

[R] xyplot

2007-10-01 Thread Samor Gandhi
Hello, I am calling the following code with the loop! It makes 3 out graphs but empty! COuld you help me plase on that? Thank you in advance z1 <- dbConnect(MyData, "something", "A1", "A2") for (tt in c("xyz", "abc", "m1")) { message(paste("Here", tt, "!!!")) mydata <- dbReadTable(z1, t

Re: [R] plot graph with error bars trouble

2007-10-01 Thread Martin Maechler
> "GaGr" == Gabor Grothendieck <[EMAIL PROTECTED]> > on Sun, 30 Sep 2007 21:15:37 -0400 writes: GaGr> On 9/30/07, hadley wickham <[EMAIL PROTECTED]> wrote: >> On 9/30/07, jiho <[EMAIL PROTECTED]> wrote: >> > On 2007-September-30 , at 22:40 , hadley wickham wrote: >> >

Re: [R] Apparently Conflicting Results with coxph

2007-10-01 Thread Kevin E. Thorpe
Peter Dalgaard wrote: > Kevin E. Thorpe wrote: >> Dear List: >> >> I have a data frame prepared in the couting process style for including >> a binary time-dependent covariate. The first few rows look like this. >> >> PtNo StartEnd Status Imp >> 1 1 0 608.0 0 0 >> 2 2

Re: [R] Apparently Conflicting Results with coxph

2007-10-01 Thread Peter Dalgaard
Kevin E. Thorpe wrote: > Dear List: > > I have a data frame prepared in the couting process style for including > a binary time-dependent covariate. The first few rows look like this. > > PtNo StartEnd Status Imp > 1 1 0 608.0 0 0 > 2 2 0 513.0 0 0 > 3

Re: [R] R CMD build not excluding .svn

2007-10-01 Thread Prof Brian Ripley
It does this via an exclude list which it uses to remove files from a staged build. Printing out the exclude list (at ca line 232) should be informative. One possibility is simply that there is a permissions issue in deleting that directory. On Sun, 30 Sep 2007, Michael Lawrence wrote: > Hi,

Re: [R] Reenviar: variables selection in Aalen and Cox model

2007-10-01 Thread Ravi Varadhan
Check out the function coxpath() in the package "glmpath". This is a better approach than doing stepwise covariate selection. Best, Ravi. --- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Healt

Re: [R] export csv

2007-10-01 Thread Wayne.W.Jones
Following on from your example: write.table(M,file="Myfile.csv",sep=",",row.names=F) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Yves Moisan Sent: 01 October 2007 14:45 To: r-help@r-project.org Subject: Re: [R] export csv Hello, I wanna know how

Re: [R] orientlib

2007-10-01 Thread Abu Naser
Hi Duncan, Thanks for your replay. I am confused as I dont have the quaternion backgorund. Let me explain what I am trying to do: I have 3 coordinates in space (x1,y1,z1),(x2,y2,z2), and (x3,y3,z3) . Now, I use (x1,y1,z1) and (x2,y2,z2) to define axis of rotaion. Then, I wanted to use 3rd co

Re: [R] export csv

2007-10-01 Thread Yves Moisan
Hello, I wanna know how to export a data frame to Excel. For example I wanna export this data frame > M[1:5,] var1 var2 distance 1 41 42 0.2 2 41 43 0.304347826086957 3 41 46 0.19047619047619 4 41 47 0.156 5 41 48 0.209302325581395 I do

Re: [R] optimize() stuck in local plateau ?

2007-10-01 Thread Ravi Varadhan
Please read the help for optimize() carefully. The following excerpted from there should help explain your problem: "The first evaluation of f is always at x_1 = a + (1-phi)(b-a) where (a,b) = (lower, upper) and phi = (sqrt 5 - 1)/2 = 0.61803.. is the golden section ratio. Almost always, the seco

Re: [R] Disentagling formulas

2007-10-01 Thread Gabor Grothendieck
Try this: fo <- y ~ x fo[[2]] <- as.name("z") fo On 10/1/07, Rebecca Sela <[EMAIL PROTECTED]> wrote: > I am writing a program in which I would like to take in a formula, change the > response (Y) variable into something else, and then pass the formula, with > the new Y variable to another funct

[R] Apparently Conflicting Results with coxph

2007-10-01 Thread Kevin E. Thorpe
Dear List: I have a data frame prepared in the couting process style for including a binary time-dependent covariate. The first few rows look like this. PtNo StartEnd Status Imp 1 1 0 608.0 0 0 2 2 0 513.0 0 0 3 2 513 887.0 0 1 4 3

[R] Disentagling formulas

2007-10-01 Thread Rebecca Sela
I am writing a program in which I would like to take in a formula, change the response (Y) variable into something else, and then pass the formula, with the new Y variable to another function. That is, I am starting with formula <- Y~X1+X2+X3 and I'd like to do something like

[R] export csv

2007-10-01 Thread elyakhlifi mustapha
Hello, I wanna know how to export a data frame to Excel. For example I wanna export this data frame > M[1:5,] var1 var2 distance 1 41 42 0.2 2 41 43 0.304347826086957 3 41 46 0.19047619047619 4 41 47 0.156 5 41 48 0.209302325581395 I don'

Re: [R] plot graph with error bars trouble

2007-10-01 Thread Jim Lemon
Marcelo Laia wrote: > Hi, > > I have a data set like this: > > MutantRepTime OD > 02H02100.029 > 02H02200.029 > 02H02300.023 > 02H02180.655 > 02H02280.615 > 02H02380.557 > 02H021121.776 > 02H0221

[R] Adding circles or ellipses to graphs

2007-10-01 Thread Friedman, Steven
Hello, I'm developing an ordination using metaMDS (package vegan). The analysis identifies 3 distinct groups that I'd like to define by either adding circles or ellipses to help identify the groups. The data set is a spatial temporal data base depicting change in each of 4 areas over 3 time pe

Re: [R] 3-dimensional graph

2007-10-01 Thread bbolker
John Sorkin wrote: > > Windows XP > R 2.3.1 > > I have a funciton > fit1<-lm(y~x+z) > Is there a function that will produce a 3-dimensional plot of y,x,z? > > I looked at the help files, but did not find a clean answer to my questio > Check out ?scatter3d in the Rcmdr package. Ben Bolker

Re: [R] Errorlog

2007-10-01 Thread Duncan Murdoch
On 10/1/2007 8:47 AM, Bernd Stampfl wrote: > ERRORLOG > > Just wanted to import a bulk of data into R via BLOOMBERG. > > And it crashed. > > Is there any useful "something" like an ERRORLOG? > > Checked the web but did not get useful information. If by "crashed" you don't mean crashed, then t

Re: [R] Concatenating one character vector into one string

2007-10-01 Thread Duncan Murdoch
On 10/1/2007 8:42 AM, Bos, Roger wrote: > Rainer, > > This works: paste(x, sep="", collapse="") > > I would have thought this 'paste(x) would work, but it didn't. > > Perhaps someone else can tell us why we need collapse="". What would > break if the default was collapse="" instead of collapse

Re: [R] Concatenating one character vector into one string

2007-10-01 Thread michael watson (IAH-C)
We need both sep and collapse to deal with situations where we are pasting two character vectors and also where we are pasting just one. paste(c("T","i","o","s"), c("his ", "s ", "ne ", "entence."),sep="",collapse="") -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] O

[R] corMatrix crashes R 2.5.1 (windows XP) with corARMA structure

2007-10-01 Thread Benjamin Tyner
R-helpers, n <- 100 arcoefs <- c(0.8) macoefs <- c(-0.6) p <- length(arcoefs) q <- length(macoefs) require(nlme) tmp <- corARMA(value=c(arcoefs,macoefs), form=~1, p=p, q=q) Sigma <- corMatrix(tmp, covariate = 1:n) # results in segfault Have I used these commands in an improper way? Thanks Ben _

Re: [R] Concatenating one character vector into one string

2007-10-01 Thread Rainer M. Krug
Thanks a million Rainer [EMAIL PROTECTED] wrote: > x <- c("This ", "is ", "one ", "sentence.") > paste(x,collapse="") > > > > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] Behalf Of Rainer M. Krug > Sent: 01 October 2007 13:23 > To: r-help > Subject: [R] Con

Re: [R] Diff. between time series and stochastic process

2007-10-01 Thread Duncan Murdoch
On 10/1/2007 8:39 AM, Megh Dal wrote: > Hi, Can anyone give me a good explanation about what is the difference > between Stochastic Process and Time Series? In my knowledge Time series > process is one type of Stochastic process. Am I right? I need a further > explanation. You're partly right.

[R] Errorlog

2007-10-01 Thread Bernd Stampfl
ERRORLOG Just wanted to import a bulk of data into R via BLOOMBERG. And it crashed. Is there any useful "something" like an ERRORLOG? Checked the web but did not get useful information. Cheers, Bernd -- View this message in context: http://www.nabble.com/Errorlog-tf4548140.html#a12978580 S

Re: [R] simple matching with R

2007-10-01 Thread Birgit Lemcke
Hello Jeff, thanks a lot for your help. It seems to work well now. Greetings Birgit Am 28.09.2007 um 20:33 schrieb Jeffrey Robert Spies: > Hi Birgit, > > I've updated the recipe here, including a change to the > dissimilarity function (making it more efficient): > > http://www.r-cookbook.com

Re: [R] Concatenating one character vector into one string

2007-10-01 Thread Bos, Roger
Rainer, This works: paste(x, sep="", collapse="") I would have thought this 'paste(x) would work, but it didn't. Perhaps someone else can tell us why we need collapse="". What would break if the default was collapse="" instead of collapse=NULL? Thanks, Roger -Original Message- From:

[R] Diff. between time series and stochastic process

2007-10-01 Thread Megh Dal
Hi, Can anyone give me a good explanation about what is the difference between Stochastic Process and Time Series? In my knowledge Time series process is one type of Stochastic process. Am I right? I need a further explanation. Thanks and regards, -

  1   2   >