Re: [R] read sas data into R

2007-10-08 Thread Peter Dalgaard
Mary Royerr wrote: > Hi > > I am having trouble using read.ssd. Can someone help? > > The code that I have written is > > *sashome<-"C:/Mary/Datasets"* > *read.ssd(read.ssd(file.path(sashome, "core", "sashelp"), "surv_1v",* > *sascmd = file.path(sashome, "sas.exe"))* > > Here the path that I have g

[R] simulation using ranks

2007-10-08 Thread Robinson Tomniavia
Hi there, how can i do simulations(100,000) using ranks. I mean i subtract main effects from the observations to form residuals, then i ranked the residuals. Next i allocate the ranks in place of the observations and do a nonparametric analysis on the ranks. The original data (observation) is ge

Re: [R] Linear regression with constraints on the parameters.

2007-10-08 Thread Moshe Olshansky
Hi Gopi, Simple linear regression minimizes sum of squares of the residuals. So in your case you can use Quadratic Programming (see quadprog package) to introduce linear constraints. Regards, Moshe. --- Gopi Goswami <[EMAIL PROTECTED]> wrote: > Hi there, > > > Is there an existing package in

[R] read sas data into R

2007-10-08 Thread Mary Royerr
Hi I am having trouble using read.ssd. Can someone help? The code that I have written is *sashome<-"C:/Mary/Datasets"* *read.ssd(read.ssd(file.path(sashome, "core", "sashelp"), "surv_1v",* *sascmd = file.path(sashome, "sas.exe"))* Here the path that I have given is correct, where the dataset su

[R] Greek in vector

2007-10-08 Thread Todd A. Gibson
Hello, I understand how to use substitute() to integrate Greek characters into plot labels, but I haven't been successful getting Greek to appear in the vector entries passed to legend() labels. Any hints? Thanks, -TAG __ R-help@r-project.org mailing l

[R] How to create something between a script and a package

2007-10-08 Thread Sergio Correia
(Before starting: I'm a total R noob so please bear with me in case of any error or faux pas). Hi, For a small project, I'm writing a few simple R functions and calling them from python (using RPy). I'm sharing the code with a couple of friends, using a subversion server. Now, I want something l

Re: [R] graph or svn

2007-10-08 Thread Rolf Turner
On 9/10/2007, at 12:04 PM, H. Paul Benton wrote: > Hello all, > > I'm having a spot of trouble. I have made a very nice distance > matrix from some comparison calculations. However, the matrix is > 3096 x > 3096, so it's relatively large. Currently I'm using the svn package as > it seems t

Re: [R] igraph and plotting connected components

2007-10-08 Thread elw
>> Hello there, >> >> I am using the igraph package to build graphs from my data. If I plot >> a graph though, it's not easy for me to see what's going on. Does >> anybody know how to rearrange a graph to get a plot without too many >> crossing lines? Maybe other packages? Edge-crossing mi

[R] Linear regression with constraints on the parameters.

2007-10-08 Thread Gopi Goswami
Hi there, Is there an existing package in R that does simple linear regression with linear constraints on the parameters? Here is the set up: y_i = \sum_{k = 1}^K \beta_k x_k + \epsilon_i where \sum_{k = 1}^K c_k \beta_k = c_0, for some known co

[R] saving simulation as .dat

2007-10-08 Thread Nigel . Walker
Hi group I have simulated some data using the garchSim and would like to save the series as a .dat file. I am not sure how to do this unfortunately any help would be appreciated. Thanks, Nigel __ R-help@r-project.org mailing list https://stat.ethz.ch/

Re: [R] igraph and plotting connected components

2007-10-08 Thread Duncan Murdoch
On 08/10/2007 8:27 PM, Dieter Best wrote: > Hello there, > > I am using the igraph package to build graphs from my data. If I plot a > graph though, it's not easy for me to see what's going on. Does anybody know > how to rearrange a graph to get a plot without too many crossing lines? Maybe

Re: [R] igraph and plotting connected components

2007-10-08 Thread Robert Gentleman
graph, RBGL and Rgraphviz, all available at www.bioconductor.org Dieter Best wrote: > Hello there, > > I am using the igraph package to build graphs from my data. If I plot a > graph though, it's not easy for me to see what's going on. Does anybody know > how to rearrange a graph to get a

Re: [R] Fill a polygon with a pattern

2007-10-08 Thread Duncan Murdoch
On 08/10/2007 5:42 PM, Alberto Monteiro wrote: > What is the simplest (namely: minimum number of libraries or conflicts) > way to fill a polygon with a pattern? > > For example, suppose I have (in a graphic file, like a jpeg) the > drawing of an "X". Then I fill a polygon (like a triangle) with >

[R] igraph and plotting connected components

2007-10-08 Thread Dieter Best
Hello there, I am using the igraph package to build graphs from my data. If I plot a graph though, it's not easy for me to see what's going on. Does anybody know how to rearrange a graph to get a plot without too many crossing lines? Maybe other packages? Thanks a lot in advance for

[R] variance explained by each term in a GAM

2007-10-08 Thread Julian M Burgos
Hello fellow R's, I do apologize if this is a basic question. I'm doing some GAMs using the mgcv package, and I am wondering what is the most appropriate way to determine how much of the variability in the dependent variable is explained by each term in the model. The information provided by sum

[R] graph or svn

2007-10-08 Thread H. Paul Benton
Hello all, I'm having a spot of trouble. I have made a very nice distance matrix from some comparison calculations. However, the matrix is 3096 x 3096, so it's relatively large. Currently I'm using the svn package as it seems to work for my data. The graph package errors and says the "adjacen

Re: [R] embedFonts rotates figure

2007-10-08 Thread Paul Smith
On 10/8/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > >> Consider the following code: > >> > >> pdf(file="figure.pdf",family="URWPalladio") > >> > >> curve(dlnorm(x,0,1.5),0,10,xlim=c(0,10),ylim=c(0,0.85),axes=F,xlab="",ylab="f") > >> segments(exp(-1.5^2),0,exp(-1.5^2),dlnorm(exp(-1.5^2),0,1.5

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Thomas Lumley
A much easier way to do this, and what Kevin probably was looking for, is replicate(). replicate(100, sum(sample(1:6,10, replace=TRUE))) gives 100 sums of 10 fair dice. -thomas On Mon, 8 Oct 2007, Alberto Monteiro wrote: > Kevin Zembower wrote: >> >> I'm trying to get R to simulate

[R] shared object for Rmpi

2007-10-08 Thread Erin Hodgess
Dear R yet again! I finally got the Rmpi package to install! Yay!. However, when I go into R to use it, there is an error with the dyn.load file. This should be simple (allegedly). Please help. Thanks, Sincerely, Erin mailto: [EMAIL PROTECTED] /home/faculty/hodgess/R-2.6.0/bin [EMAIL PROTE

Re: [R] do not plot polygon boundaries with spplot {sp}

2007-10-08 Thread roger koenker
Try the argument border -- if it gets passed to polygon then border: the color to draw the border. The default, 'NULL', means to use 'par("fg")'. Use 'border = NA' to omit borders. If not, then you may need to work around a bit. url:www.econ.uiuc.edu/~roger

Re: [R] Fill a polygon with a pattern

2007-10-08 Thread Greg Snow
The simplest approach using just R is probably to create the background by plotting the grid of your images, then draw with your background color over everything outside the polygon that you want. A possibly simpler approach is to draw the polygon in a solid color, then use an outside program such

Re: [R] How put description under the graph (using different lines)

2007-10-08 Thread Greg Snow
If your descriptions are just 1 line each then you can give yourself extra room at the bottom with par(mar=c(... Then use the mtext function to add the descriptions. If the descriptions are more detailed then you can try something like: > library(gplots) > layout( rbind( 1,2,c(3,4) ) ) > barplot(

[R] Fill a polygon with a pattern

2007-10-08 Thread Alberto Monteiro
What is the simplest (namely: minimum number of libraries or conflicts) way to fill a polygon with a pattern? For example, suppose I have (in a graphic file, like a jpeg) the drawing of an "X". Then I fill a polygon (like a triangle) with this pattern, and get something like: X XXX X

[R] How put description under the graph (using different lines)

2007-10-08 Thread akki
Hello, I need to use barplot' graphs, but under the graph I need to put different things. The structure of this graph must be: MY GRAPH -- description 1 --- - de

[R] do not plot polygon boundaries with spplot {sp}

2007-10-08 Thread Dylan Beaudette
Hi, Is there a simple way to suppress the plotting of polygon boundaries with spplot() ? # simple list of 12 colors cols <- brewer.pal(12, "Paired") # plot pile of polygons, with 12 classes: spplot(x, zcol='class2', col.regions=cols, scales=list(draw=T), xlab="Easting (m)", ylab="Northing (m)"

Re: [R] embedFonts rotates figure

2007-10-08 Thread Prof Brian Ripley
On Mon, 8 Oct 2007, Peter Dalgaard wrote: > Paul Smith wrote: >> Dear All, >> >> Consider the following code: >> >> pdf(file="figure.pdf",family="URWPalladio") >> >> curve(dlnorm(x,0,1.5),0,10,xlim=c(0,10),ylim=c(0,0.85),axes=F,xlab="",ylab="f") >> segments(exp(-1.5^2),0,exp(-1.5^2),dlnorm(exp(-1.

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Julian M Burgos
Well, if the dice are fair, then all rolls are independent. If you want to roll each dice 100 times, you can do something like rolls=matrix(sample(1:6,1000,replace=T),ncol=10) apply(rolls,2,sum) I'm sure that there must be a more elegant way to do it, though. Julian On Mon, 8 Oct 2007, Zemb

Re: [R] Linux editor for R+LaTeX, but not Emacs

2007-10-08 Thread Gregory. R. Warnes
I, personally, see tremendous slowdowns in XEmacs responsiveness when the buffer gets large, probably due to font-lock pretty-printing stuff. Since ESS can't / doesn't submit the next command until the fontification is done, things can slow down by a lot. -G On Oct 7, 2007, at 4:24PM , D

[R] "xlab" and "ylab" in clusplot

2007-10-08 Thread Haiyong Xu
Hi there, How can I change the "xlab" and "ylab" in clusplot instead of "Component 1" and "Component 2"? It always gives me the following error message. Error in plot.default(x1[, 1], x1[, 2], xlim = xlim, ylim = ylim, xlab = "Component 1", : formal argument "ylab" matched by multi

Re: [R] Getting intervals for within-group standard errors for each group using nlme and varIdent

2007-10-08 Thread Greg Snow
Just to follow-up on this for the archive and anyone else who is interested in this problem. I ended up solving this by running lme seperately on subsets of the data (one for group=='A' and one for group=='B'), the intervals for the statistics computed both ways were within rounding of each other,

Re: [R] embedFonts rotates figure

2007-10-08 Thread Peter Dalgaard
Paul Smith wrote: > Dear All, > > Consider the following code: > > pdf(file="figure.pdf",family="URWPalladio") > > curve(dlnorm(x,0,1.5),0,10,xlim=c(0,10),ylim=c(0,0.85),axes=F,xlab="",ylab="f") > segments(exp(-1.5^2),0,exp(-1.5^2),dlnorm(exp(-1.5^2),0,1.5),lty="dashed") > segments(1,0,1,dlnorm(1,0

[R] lmList function in lme4

2007-10-08 Thread Nathan Leon Pace, MD, MStat
I don¹t understand the following error message in lme4. > lmList(formula = oaas.bin ~ oaas.prob | subject, data = inductionPropRemiAll.df, family = binomial) Error in FUN(X[[1L]], ...) : unused argument(s) (family = function (link = "logit") > lmList(formula = oaas.bin ~ oaas.prob | subject, dat

[R] embedFonts rotates figure

2007-10-08 Thread Paul Smith
Dear All, Consider the following code: pdf(file="figure.pdf",family="URWPalladio") curve(dlnorm(x,0,1.5),0,10,xlim=c(0,10),ylim=c(0,0.85),axes=F,xlab="",ylab="f") segments(exp(-1.5^2),0,exp(-1.5^2),dlnorm(exp(-1.5^2),0,1.5),lty="dashed") segments(1,0,1,dlnorm(1,0,1.5),lty="dashed") segments(exp(

Re: [R] Mouse Location on Plots

2007-10-08 Thread Alberto Monteiro
Greg Snow wrote: > > Another option (which should work on all platforms) is to use tcltk and > bind a function to the "Motion" event. The playSudoku function in > the sudoku package has an example of this (as well as using > getGraphicsEvent). > > (wow, could that package actually be useful?)

Re: [R] R and FDA trials

2007-10-08 Thread Gregory Warnes
Hi All, I'm excited to see that R-Core has released the document "R: Regulatory Compliance and Validation Issues A Guidance Document for the Use of R in Regulated Clinical Trial Environments" (http://www.r-project.org/doc/R-FDA.pdf). I know it represents a great deal of effort on the part of

[R] variance explained by each term in a GAM

2007-10-08 Thread Julian M Burgos
Hello fellow R's, I do apologize if this is a basic question. I'm doing some GAMs using the mgcv package, and I am wondering what is the most appropriate way to determine how much of the variability in the dependent variable is explained by each term in the model. The information provided by

[R] get more leads with these healthcare lists

2007-10-08 Thread Hines K Efrain
Certified MDs in the United States 788,437 in total – 17,880 emails Many popular specialties like Emergency Medicine, Plastic Surgery, OBGYN, Oncology, Pediatrics and more Many unique fields like 'medical school attended' and 'location of residency training' Special Price: $321 *** FR

Re: [R] Mouse Location on Plots

2007-10-08 Thread Greg Snow
Another option (which should work on all platforms) is to use tcltk and bind a function to the "Motion" event. The playSudoku function in the sudoku package has an example of this (as well as using getGraphicsEvent). (wow, could that package actually be useful?) -- Gregory (Greg) L. Snow Ph.D.

Re: [R] Tart charts

2007-10-08 Thread Antony Unwin
Michael, > Try this alternative: > > # from http://research.microsoft.com/users/lamport/pubs/hair.pdf > hairsex <- matrix( >c(46, 45, 13, 12, > 1, 101, 0, 20), 2, 4, byrow=TRUE) > dimnames(hairsex) <- list("Gender"=c("Female", "Male"), > "Hair color"=c("Blond", "Brown", "Re

Re: [R] as.ordered

2007-10-08 Thread Friedrich Schuster
Hello, (Warning. This might not be the most complete or elegant solution ...) If you want a sorted dataframe: look here for example http://tolstoy.newcastle.edu.au/R/help/05/02/12391.html To convert the factors from a data frame, you have to call as.ordered for each factor separately (not for t

Re: [R] Simulate data based on correlation coefficient

2007-10-08 Thread Greg Snow
Similar questions have been asked here before. See http://finzi.psych.upenn.edu/R/Rhelp02a/archive/110042.html for one possible answer (your case will be simpler than this). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 >

Re: [R] as.ordered

2007-10-08 Thread Friedrich Schuster
Hello, (Warning. This might not be the most complete or elegant solution ...) If you want a sorted dataframe: look here for example http://tolstoy.newcastle.edu.au/R/help/05/02/12391.html To convert the factors from a data frame, you have to call as.ordered for each factor separately (not for

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Alberto Monteiro
Kevin Zembower wrote: > > I'm trying to get R to simulate the sum of the values on 10 fair dice > (yes, it's related to a homework problem, but is not the problem > itself). I tried to do this: > > rep(sum(sample(1:6,100,replace=T)), times=10) > [1] 341 341 341 341 341 341 341 341 341 341 > rep(s

[R] coxph() command design and data setup

2007-10-08 Thread Justin Montemarano
Hello all: I'm attempting to run a Cox proportional hazards function on survival data from insects and I have a few questions. My current command that I'm using to call the model is as follows (using coxph() from the survival library): coxph(Surv(day, censor) ~ treatment + room + chamber %in% tr

Re: [R] print Text on device

2007-10-08 Thread Greg Snow
Look at the textplot function in the gplots package. Does that do what you want? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Richard M. Heiberger
rep( sum(sample(1:6,100,replace=T)), times=10) Read carefully what you wrote. You asked R to sum a single sample, then make ten copies of the sum. You need to take ten samples, then sum each. apply() would be helpful. -Original Message- I'm trying to get R to simulate the sum

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Zembower, Kevin
Thanks so much, Chuck and Mark. Here's my script to simulate 10,000 rolls of 100 fair dice to demonstrate their conformity to a normal curve: > x<-replicate(1, sum(sample(1:6,100,replace=T))) > sdx<-sd(x) > sdx [1] 17.13966 > meanx<-mean(x) > meanx [1] 350.0451 > hist(x, freq=FALSE) > curve(dno

Re: [R] estfun & df

2007-10-08 Thread Douglas Bates
On 10/8/07, Achim Zeileis <[EMAIL PROTECTED]> wrote: > On Mon, 8 Oct 2007, Abdus Sattar wrote: > > > Hello EVERYONE, > > > > I need an URGENT help from you please! > > This type of requests is not considered to be very polite, please have a > look at the posting guide. > > > How can I see the "estf

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Prof Brian Ripley
See ?replicate, e.g. replicate(sum(sample(1:6,100,replace=TRUE)), n=10) Function arguments are (in gneral) evaluated as if they were evaluated when passed to the function. On Mon, 8 Oct 2007, Zembower, Kevin wrote: > I'm trying to get R to simulate the sum of the values on 10 fair dice > (yes,

Re: [R] problem with times

2007-10-08 Thread Gavin Simpson
On Mon, 2007-10-08 at 07:38 -0700, Samuel Okoye wrote: > Hello, I have got the following problem: > > times <- c("02.07.2007", "03.07.2007","03.09.2007", > > "04.07.2007","05.07.2007") > > mode(times) > [1] "numeric" > > tim <- as.character(times) > > mode(tim) > [1] "character" > > sort(times) >

[R] another try with Rmpi on Linux

2007-10-08 Thread Erin Hodgess
Dear R People: Here are my latest attempts for Rmpi: This is the list of locations for mpi.h: /home/faculty/hodgess/mpich2-1.0.5p4/src/include/mpi.h /home/faculty/hodgess/include/mpi.h /home/faculty/hodgess/lam-7.1.4/romio/adio/sgi/mpi3.1/mpi.h /home/faculty/hodgess/lam-7.1.4/share/incl

Re: [R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Charles C. Berry
See ?replicate which I think is what you are after. Chuck On Mon, 8 Oct 2007, Zembower, Kevin wrote: > I'm trying to get R to simulate the sum of the values on 10 fair dice > (yes, it's related to a homework problem, but is not the problem > itself). I tried to do this: >> rep(sum(s

[R] Dice simulation: Getting rep to re-evaluate sample()?

2007-10-08 Thread Zembower, Kevin
I'm trying to get R to simulate the sum of the values on 10 fair dice (yes, it's related to a homework problem, but is not the problem itself). I tried to do this: > rep(sum(sample(1:6,100,replace=T)), times=10) [1] 341 341 341 341 341 341 341 341 341 341 and noticed that sum(sample()) seems to

Re: [R] estfun & d

2007-10-08 Thread Douglas Bates
On 10/8/07, Abdus Sattar <[EMAIL PROTECTED]> wrote: > Hello EVERYONE, > I need an URGENT help from you please! > How can I see the "estfun" (empirical estimating function) and "df" (degree > of freedom) from the following mixed-model please? > (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sle

Re: [R] Reading distance matrices

2007-10-08 Thread Doran, Harold
Muri If it is lower triangular, then the upper triangle is empty. What do you expect it to have? We don't have any information on how you're reading in your data. But, read.table() has an argument for na.strings that you should look at. > -Original Message- > From: [EMAIL PROTECTED] > [m

Re: [R] estfun & df

2007-10-08 Thread Achim Zeileis
On Mon, 8 Oct 2007, Abdus Sattar wrote: > Hello EVERYONE, > > I need an URGENT help from you please! This type of requests is not considered to be very polite, please have a look at the posting guide. > How can I see the "estfun" (empirical estimating function) I guess (because you are not tel

[R] CompetingRiskFrailty modeling

2007-10-08 Thread Daniel Malter
Hi all can anybody give info what I do wrong when modeling a survival function with the CompetingRiskFrailty package in the following way? "accumulate", "reduce" and "closed" are dummy vectors with 0s indicating no event and 1s indicating the respective event in the appropriate column. ac

[R] problem with times

2007-10-08 Thread Samuel Okoye
Hello, I have got the following problem: > times <- c("02.07.2007", "03.07.2007","03.09.2007", "04.07.2007","05.07.2007") > mode(times) [1] "numeric" > tim <- as.character(times) > mode(tim) [1] "character" > sort(times) [1] "02.07.2007" "03.07.2007" "03.09.2007" "04.07.2007" "05.07.2007" Is it p

Re: [R] data.frame's and [2, c(T, F)] in R 2.6.0

2007-10-08 Thread data-ploner.com
On Oct 8, 2007, at 4:58 PM, Peter Dalgaard wrote: > data-ploner.com wrote: >> In R prior to 2.6.0 with matrizes as well as with data.frames it is >> possible to do: >> >> a[2, c(TRUE, FALSE)] >> >> In R 2.6.0 the mixed indexing works only for matrizes. Is this the >> intention, a bug, or did I get

[R] estfun & df

2007-10-08 Thread Abdus Sattar
Hello EVERYONE, I need an URGENT help from you please! How can I see the "estfun" (empirical estimating function) and "df" (degree of freedom) from the following mixed-model please? (fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)) Many thanks in advance for your kind help.

[R] Reading distance matrices

2007-10-08 Thread Muri Soares
Hi all, I've been having trouble reading in distance matrices (values in lower triangle only) that have been created by other programs. They load with NA values in the upper triangle. Is there an option to read the matrices properly? Thanks, Muri _

Re: [R] Applying function to data.frame

2007-10-08 Thread Gabor Grothendieck
Your function needs to be able to work on vector input, e.g. transform( y, time=Vectorize(str_to_millis)( as.character(time) ) ) or just do this: library(chron) y$time <- 24 * 60 * 60 * 1000 * as.numeric(times(y$time)) On 10/8/07, Rees, David <[EMAIL PROTECTED]> wrote: > Hi, > > If I have the

Re: [R] RTable width sum colum

2007-10-08 Thread Gabor Grothendieck
If it really is of class "table" (or is an array) or you convert it to such you can use addmargins. You may need to read ?addmargins carefully. On 10/8/07, Alessandra Marmo <[EMAIL PROTECTED]> wrote: > Hello > I need to print table width "colum sum" like > XV1V2TOT > A23

Re: [R] comparing matched proportions using glm

2007-10-08 Thread Charles C. Berry
On Mon, 8 Oct 2007, Corry Gellatly wrote: > > Thanks very much for your reply Chuck, I have a quick follow up > question. You mention putting the data into a 2x2x3 for log-linear > model, however my lists have many more than 3 strata, actually > thousands. I am trying to work out whether the propo

[R] as.ordered

2007-10-08 Thread Birgit Lemcke
Hello Members, I try to convert variables in a data.frame (bract.awn) in the class ordered. str(bract.awn) 'data.frame': 348 obs. of 2 variables: $ bracts.length.relative.to.flower...Min: Factor w/ 4 levels "1","2","3","4": 2 3 3 3 3 2 1 4 3 2 ... $ bract.awn.relative.to.body..

Re: [R] A rebel boxplot question

2007-10-08 Thread p.valdes
Thanks to all for the answers --- Thierry : a very interesting answer... but I'm still trying to figure how to install ggplot2 :-) install.packages("ggplot") Warning in install.packages("ggplot2") : argument 'lib' is missing: using '...R/i486-pc-linux-gnu-library/2.5' ?? also installing

Re: [R] R and FDA trials

2007-10-08 Thread Cody Hamilton
Ricardo, Tony Rossini (Novartis) provided a very good post on point (3) below (which you may have already read): http://tolstoy.newcastle.edu.au/R/e2/help/07/08/23820.html. Regards, -Cody -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ricardo Pietr

Re: [R] Specify plot size

2007-10-08 Thread Prof Brian Ripley
On Mon, 8 Oct 2007, Eric Thompson wrote: > Have you tried using pdf(), postscript(), or jpeg()? All of these have > arguments to specify the height and width of the device. And so does every other graphics device. Note though that the height/width of the device is not that of a plot (nor necess

Re: [R] Incompatible methods ("-.POSIXt", "Ops.difftime") for "-"

2007-10-08 Thread Prof Brian Ripley
On Mon, 8 Oct 2007, Albrecht, Dr. Stefan (APEP) wrote: Dear all, according to the Help-page of DateTimeClasses {base} I should be able to do time - z with timedate-time objects z a numeric vector (in seconds) or an object of class "difftime". However, on R version 2.6.0 (Windows X

[R] Error message on script execution

2007-10-08 Thread Sébastien
Dear R_users, I have some troubles with a visual basic application I have recently created. This application automatically creates and executes an R script based on a dataset and settings defined by the user. The .r file (Create_Diagnostic_plots.r) is saved first in a given folder and then cal

[R] pruning of Additive Regression Trees

2007-10-08 Thread Vadim Ogranovich
Hi, Does R provide functions for pruning of Additive Regression Tree Models? Suppose I have grown a large set of trees, say 3000 trees, and now want to find a 'computationally more economic' approximation of the function defined by the trees. In the case of a single tree this is pruning. Is t

Re: [R] RTable width sum colum

2007-10-08 Thread Marc Schwartz
On Mon, 2007-10-08 at 16:18 +0200, Alessandra Marmo wrote: > Hello > I need to print table width "colum sum" like > XV1V2TOT > A235 > B167 > C46 10 > > how i can do it by table function? > > Thanks > > Alessandra Presuming th

Re: [R] sort time

2007-10-08 Thread Marc Schwartz
On Mon, 2007-10-08 at 07:26 -0700, Samuel Okoye wrote: > Hello, I have got the following problem: > > times <- c("02.07.2007", "03.07.2007","03.09.2007", > "04.07.2007","05.07.2007") > > mode(times) > [1] "numeric" > > tim <- as.character(times) > > mode(tim) > [1] "character" > > sort(times)

Re: [R] sort time

2007-10-08 Thread Ben Bolker
Samuel Okoye wrote: > > Hello, I have got the following problem: > [snip] > > Assuming that these are in month/day/year format (which I conclude from your sorted values): times <- c("02.07.2007", "03.07.2007","03.09.2007", "04.07.2007","05.07.2007") times2 <- as.Date(times,format="%m.%d.%Y

Re: [R] data.frame's and [2, c(T, F)] in R 2.6.0

2007-10-08 Thread Peter Dalgaard
data-ploner.com wrote: > On Oct 8, 2007, at 4:58 PM, Peter Dalgaard wrote: >> >> >> As far as I can tell, It happens only when drop=TRUE and exactly one >> column is selected using a logical index. > > Thank you. > Maybe it is better letting [] work as in previous R ... ? Well, presumably the code

Re: [R] Odp: How to change row value based on previous row value

2007-10-08 Thread Peter Dalgaard
Petr PIKAL wrote: > Peter Dalgaard <[EMAIL PROTECTED]> napsal dne 08.10.2007 16:31:28: > > >> Petr PIKAL wrote: >> >>> Hi >>> >>> [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: >>> >>> >>> Hi, I would like to fix the data in the following data.frame, and I am h

Re: [R] linewidth in Rgraphviz

2007-10-08 Thread Ben Bolker
Robert Gentleman wrote: > > Hi Ben, >Sorry for the slow reply, but something like this works, > > [...] > Sorry not to have responded sooner. I had actually worked it out for myself (after finding another e-mail exchange on the Bioconductor list that indicated the problem had been

Re: [R] Bug or not

2007-10-08 Thread Ben Bolker
Darius Kasiulevičius wrote: > > I have this sample from help. I try use command erase.screen() but in > split mode sreen dont erase. Why? > > Not a bug. From the documentation: 'erase.screen' will appear not to work if the background colour is transparent (as it is by default on m

Re: [R] Odp: How to change row value based on previous row value

2007-10-08 Thread Petr PIKAL
Peter Dalgaard <[EMAIL PROTECTED]> napsal dne 08.10.2007 16:31:28: > Petr PIKAL wrote: > > Hi > > > > [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: > > > > > >> Hi, > >> > >> I would like to fix the data in the following data.frame, and I am > >> having trouble coding this in R - help please!

Re: [R] data.frame's and [2, c(T, F)] in R 2.6.0

2007-10-08 Thread Peter Dalgaard
data-ploner.com wrote: > In R prior to 2.6.0 with matrizes as well as with data.frames it is > possible to do: > > a[2, c(TRUE, FALSE)] > > In R 2.6.0 the mixed indexing works only for matrizes. Is this the > intention, a bug, or did I get something wrong? > Best regards > > Meinhard > > It

Re: [R] Mouse Location on Plots

2007-10-08 Thread Duncan Murdoch
On 08/10/2007 8:42 AM, Samuel Kemp wrote: > Hi, > > Does anyone know of a function that dynamically updates a variable with the > (x,y) co-ordinates of the mouse pointer on a plot? > > locator and identify only return the values when you click 'stop' on the > graphics device. > > Any help would

Re: [R] semivariogram

2007-10-08 Thread ONKELINX, Thierry
Have a look at the gstat package and it's documentation. Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodologie en kwaliteitszorg / Sec

[R] .C("myHadamaProduct", xx, yy, nr, nc, out=as.double(rep(0.0, n)))

2007-10-08 Thread bernardo lagos alvarez
Konw anybody as run the function psi on >> >> >> http://www.alglib.net/translator/dl/specialfunctions.psi.csharp.zip >> >> or >> >> http://www.alglib.net/translator/dl/specialfunctions.psi.cpp.zip >> >> using .C with R. >> ? >> >> Thank for your attention. >> >> Bernardo. >> >> >>> If anybody has a

[R] sort time

2007-10-08 Thread Samuel Okoye
Hello, I have got the following problem: > times <- c("02.07.2007", "03.07.2007","03.09.2007", "04.07.2007","05.07.2007") > mode(times) [1] "numeric" > tim <- as.character(times) > mode(tim) [1] "character" > sort(times) [1] "02.07.2007" "03.07.2007" "03.09.2007" "04.07.2007" "05.07.2007" I

Re: [R] (no subject)

2007-10-08 Thread Thomas Lumley
This is the same issue as FAQ 7.19 How do I produce PNG graphics in batch mode? The advice there may be helpful. -thomas On Fri, 5 Oct 2007, Zhang, Yongqing (NIH/NIA/IRP) [C] wrote: > Hi, All, > > I try to accessing the x11() device over cgi - namely, when trying to > create a graphi

[R] data.frame's and [2, c(T, F)] in R 2.6.0

2007-10-08 Thread data-ploner.com
In R prior to 2.6.0 with matrizes as well as with data.frames it is possible to do: a[2, c(TRUE, FALSE)] In R 2.6.0 the mixed indexing works only for matrizes. Is this the intention, a bug, or did I get something wrong? Best regards Meinhard Meinhard Ploner Althingstrasse 15 39031 Brun

[R] semivariogram

2007-10-08 Thread Duccio -
I need the most straightforward way to build semivariograms within R i am currently using s-gems software but i would like to even test R Thanks D [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

[R] "R graphical manuals" updated.

2007-10-08 Thread 小笠原理
Dear all I have updated "R graphical manuals" homepage. http://cged.genes.nig.ac.jp/RGM2/index.php Since obsolute homepage http://bg9.imslab.co.jp/Rhelp/ is no longer maintained, please change your bookmark. Osamu Ogasawara [EMAIL PROTECTED] [EMAIL PROTECTED] ___

[R] RTable width sum colum

2007-10-08 Thread Alessandra Marmo
Hello I need to print table width "colum sum" like XV1V2TOT A235 B167 C46 10 how i can do it by table function? Thanks Alessandra __ R-help@r-project.org mailing list https://stat.e

Re: [R] Odp: How to change row value based on previous row value

2007-10-08 Thread Peter Dalgaard
Petr PIKAL wrote: > Hi > > [EMAIL PROTECTED] napsal dne 08.10.2007 14:24:13: > > >> Hi, >> >> I would like to fix the data in the following data.frame, and I am >> having trouble coding this in R - help please! >> >> >>> x >>> >> AB xy z >> 1 1 10.0 100 1000 1 >>

Re: [R] Applying function to data.frame

2007-10-08 Thread jim holtman
try this: > x <- read.table(textConnection(" time val + 1 08:00:05.834 1 + 2 08:03:13.345 2 + 3 08:10:12.443 3"), header=TRUE) > z <- strsplit(as.character(x$time), ":") > newtime <- lapply(z, function(a) as.numeric(a) %*% c(360, 6, 1000)) > x$time <- unlist(newtime) > x

Re: [R] heatmap

2007-10-08 Thread Liaw, Andy
From: michael watson > > If you just want a coloured representation of your distance > matrix, use image() > > ?image() W/o the parens... > If you want to cluster your original data and show the > original data (not the distances) as a heatmap then use > heatmap(), but you should use so

Re: [R] .C("myHadamaProduct", xx, yy, nr, nc, out=as.double(rep(0.0, n)))

2007-10-08 Thread Bernardo Lagos Alvarez
> > Konw anybody as run the function psi on > > > http://www.alglib.net/translator/dl/specialfunctions.psi.csharp.zip > > or > > http://www.alglib.net/translator/dl/specialfunctions.psi.cpp.zip > > using .C with R. > ? > > Thank for your attention. > > Bernardo. > > >> If anybody has an idea what i

Re: [R] tcltk scrollbar

2007-10-08 Thread Peter Dalgaard
kapo coulibaly wrote: > This is approximately what I get: > -- --- > ! !!S! ! !S! > ! !!C! ! !C! > ! !

[R] Applying function to data.frame

2007-10-08 Thread Rees, David
Hi, If I have the following data.frame >y time val 1 08:00:05.834 1 2 08:03:13.345 2 3 08:10:12.443 3 > and the following function which converts the time string to the number of milliseconds since midnight > str_to_millis function( s ) { a <- as.numeric( unlist( strsplit(s,":

Re: [R] How to remove legend?

2007-10-08 Thread hadley wickham
On 10/8/07, ONKELINX, Thierry <[EMAIL PROTECTED]> wrote: > Hi Christophe, > > Although I'm a fan of ggplot2, it's lack of documentation is a drawback. > I was having a similar problem and Hadley suggested the "sc$legend <- > FALSE" solution. AFAIK it isn't documented (yet). Have a look at > Hadley'

Re: [R] Ambiguities in vector

2007-10-08 Thread Gavin Simpson
On Mon, 2007-10-08 at 15:35 +0200, Birgit Lemcke wrote: > Hello James, > > all of your suggestions work very well except of this: > > FemMal <- cbind(FemV1gezählt[2,], MalV1gezählt[2,]) > > colnames(FemMal) <- ("Females", "Males") > Fehler: syntax error The OP missed out c() above, hence the sy

Re: [R] Ambiguities in vector

2007-10-08 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 08.10.2007 15:35:01: > Hello James, > > all of your suggestions work very well except of this: > > FemMal <- cbind(FemV1gezählt[2,], MalV1gezählt[2,]) > > colnames(FemMal) <- ("Females", "Males") > Fehler: syntax error Syntax error means something is missing an

Re: [R] Ambiguities in vector

2007-10-08 Thread Birgit Lemcke
Hello James, all of your suggestions work very well except of this: FemMal <- cbind(FemV1gezählt[2,], MalV1gezählt[2,]) colnames(FemMal) <- ("Females", "Males") Fehler: syntax error FeMMal [,1] [ ,2] 1 133 79 2 203 237 3 51 76 But it works if I do that: Namen<-c("Female","Mal

[R] Odp: Specify plot size

2007-10-08 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 08.10.2007 15:15:58: > Hello, > > I have another (possibly easy) question: > > How to specify the size of a plot? When I draw a plot, I can freely change > the size of the window, which is nice for single plots to find the best > height/width ratio, but as I nee

  1   2   >