Re: [R] lattice: overlap histogram

2008-10-28 Thread Ferry
Thank you Deepayan. Still on lattice, can the same overlap function be applied to categorical variable plot such as barchart? I mean, can I create an overlap barchart plot? On Tue, Oct 28, 2008 at 5:41 PM, Deepayan Sarkar <[EMAIL PROTECTED]>wrote: > On Tue, Oct 28, 2008 at 4:31 PM, Ferry <[EMAIL

Re: [R] Plotting weather data

2008-10-28 Thread Daniel Malter
This certainly involves arranging multiple plots and plotting multiple variables within the same plot (in the upper graph). The lower plots looks rather standard and just have to be arranged right. Further, you have to put the labels (if they are desired). There is a book on R-graphics which may he

Re: [R] definition of a function

2008-10-28 Thread Berwin A Turlach
G'day George, On Wed, 29 Oct 2008 06:21:37 +0200 KARAVASILIS GEORGE <[EMAIL PROTECTED]> wrote: > Hello, R subscribers. > I would like to define a function like this one: > Fun <- function(x)sin(x)/x > with value 1 at x=0. How can I do that? Fun <- function(x) ifelse(x==0, 1, sin(x)/x) > Is

Re: [R] definition of a function

2008-10-28 Thread Jorge Ivan Velez
Dear George, Perhaps redefining Fn as Fun2 <- function(x) ifelse(x==0,1,sin(x)/x) could help. Now, for the symbolic part (I'm not sure if understood correctly), try ?curve as follows: curve(Fun2(x),-10,10) HTH, Jorge On Wed, Oct 29, 2008 at 12:21 AM, KARAVASILIS GEORGE <[EMAIL PROTECTED]>wr

[R] definition of a function

2008-10-28 Thread KARAVASILIS GEORGE
Hello, R subscribers. I would like to define a function like this one: Fun <- function(x)sin(x)/x with value 1 at x=0. How can I do that? Is there a way to plot it "symbolically", without using x as a vector? e.g. x <- seq(from=-10, to=10, length=100) plot(x, Fun(x)) Thank you in advance. ___

[R] Plotting weather data

2008-10-28 Thread Rachael Howard
Hello. I am relatively new to R and trying to figure out how go about graphing something like this: http://www.math.yorku.ca/SCS/Gallery/images/NYweather.jpg Specifically, I am having a hard time graphing the mean, min, and max temperatures as they are depicted on that graph. When I plot th

Re: [R] behavior of "by"

2008-10-28 Thread Jeff Laake
Thanks to those that replied. It was reproducible on my system but it was nested within other code which would have been obvious had I left the browse> portions in it. The reason for the behavior was obvious when I used str(samples). Even though Effort looked like a vector in samples, it was

Re: [R] help with doing a manipulation on a column of a data frame based on another column

2008-10-28 Thread Rolf Turner
On 29/10/2008, at 2:41 PM, stephen sefick wrote: #this is my stab at - I am sure that I am missing something. If this doesn't make sense then please ask for more details. #This may show my low level of programing knowledge hester. <- c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4) value <- rnorm(16) x <-

[R] help with doing a manipulation on a column of a data frame based on another column

2008-10-28 Thread stephen sefick
#this is my stab at - I am sure that I am missing something. If this doesn't make sense then please ask for more details. #This may show my low level of programing knowledge hester. <- c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4) value <- rnorm(16) x <- data.frame(value, hester.) z <- (if(x[,"hester."]==

Re: [R] behavior of "by"

2008-10-28 Thread Sebastian P. Luque
On Tue, 28 Oct 2008 18:04:57 -0700, Jeff Laake <[EMAIL PROTECTED]> wrote: > Any insight into the behavior of "by" in the following case would be > appreciated. There is a note in the help details for "by" about > documenting behavior since v2.7 but I don't entirely understand what > it is saying.

Re: [R] behavior of "by"

2008-10-28 Thread Rolf Turner
On 29/10/2008, at 2:04 PM, Jeff Laake wrote: Any insight into the behavior of "by" in the following case would be appreciated. There is a note in the help details for "by" about documenting behavior since v2.7 but I don't entirely understand what it is saying. I'm using R2.7.2 Windows. I'

[R] behavior of "by"

2008-10-28 Thread Jeff Laake
Any insight into the behavior of "by" in the following case would be appreciated. There is a note in the help details for "by" about documenting behavior since v2.7 but I don't entirely understand what it is saying. I'm using R2.7.2 Windows. I'm interested if the following behavior was a cha

Re: [R] how to get the value of aov summary into another variable

2008-10-28 Thread Jorge Ivan Velez
Dear Waverley, Try # Data x=rnorm(10) y=rep(1:2,5) aov.ex = aov(x~y) tail(unlist(summary(aov.ex)),2) # or anova(lm(x~y))$"Pr(>F)" HTH, Jorge On Tue, Oct 28, 2008 at 8:48 PM, Waverley <[EMAIL PROTECTED]>wrote: > Hi, > > I have a question of aov. e.g. > > aov.ex = aov(x~y) > summary(aov.ex

Re: [R] how to get the value of aov summary into another variable

2008-10-28 Thread Peter Alspach
Waverley In general, examine the structure of the object and go from there. Thus: str(summary(aov.ex)) will reveal that summary(aov.ex) is a list of 1. Thus you probably want to examine that single element in more detail: str(summary(aov.ex)[[1]]) # the [[1]] to get the first element of the l

[R] how to get the value of aov summary into another variable

2008-10-28 Thread Waverley
Hi, I have a question of aov. e.g. aov.ex = aov(x~y) summary(aov.ex) The aov summary will print to the screen. How can I extract the aov result, in particular the values of Pr(>F) and F value into a vector so that I can use them for other use? Thanks. -- Waverley @ Palo Alto __

Re: [R] repeated measure one way anova

2008-10-28 Thread Chris Evans
JohnLi sent the following at 28/10/2008 19:56: > Thanks a lot ! I learn a lot about R while I solve my problem according to > those provided informatin. Me too: thanks all. This particular question has suddenly become relevant to me as I'm trying to do a bit of simulating to model power for a r

Re: [R] lattice: overlap histogram

2008-10-28 Thread Deepayan Sarkar
On Tue, Oct 28, 2008 at 4:31 PM, Ferry <[EMAIL PROTECTED]> wrote: > Dear R users, > > Is it possible to have an overlap histogram plot? For example: > > stuff <- data.frame(Mode = c("Land", "Air"), AgeGroup = c("Young", "Adult", > "Old"), Value = sample(1:300)) > histogram( ~Value | AgeGroup * Mode

Re: [R] outputting (writing) output into a dataframe

2008-10-28 Thread Henrique Dallazuanna
Try this: do.call(rbind.data.frame, lapply(p1, function(i)power.prop.test(p1=i,p2=0.5*i,power=0.8, sig.level=0.05))) On Tue, Oct 28, 2008 at 7:16 PM, Farrel Buchinsky <[EMAIL PROTECTED]> wrote: > I have solved this problem once before but don't recall exactly how. > Is there a url t

Re: [R] acf(): meaning of the blue horizontal lines

2008-10-28 Thread Oliver Bandel
Zitat von David Scott <[EMAIL PROTECTED]>: > On Tue, 28 Oct 2008, Oliver Bandel wrote: > > > Zitat von David Scott <[EMAIL PROTECTED]>: > > > >> On Tue, 28 Oct 2008, Oliver Bandel wrote: > >> > >>> Hello, > >>> > >>> what are they meaning? It could be something that > >>> would show a threshhold a

Re: [R] acf(): meaning of the blue horizontal lines

2008-10-28 Thread David Scott
On Tue, 28 Oct 2008, Oliver Bandel wrote: Zitat von David Scott <[EMAIL PROTECTED]>: On Tue, 28 Oct 2008, Oliver Bandel wrote: Hello, what are they meaning? It could be something that would show a threshhold above which the result is indicating different meanings then just random noise. But

[R] lattice: overlap histogram

2008-10-28 Thread Ferry
Dear R users, Is it possible to have an overlap histogram plot? For example: stuff <- data.frame(Mode = c("Land", "Air"), AgeGroup = c("Young", "Adult", "Old"), Value = sample(1:300)) histogram( ~Value | AgeGroup * Mode, data = stuff, auto.key = TRUE) Instead of having 2 * 3 panel, I want to hav

Re: [R] acf(): meaning of the blue horizontal lines

2008-10-28 Thread Oliver Bandel
Zitat von David Scott <[EMAIL PROTECTED]>: > On Tue, 28 Oct 2008, Oliver Bandel wrote: > > > Hello, > > > > what are they meaning? It could be something that > > would show a threshhold above which the result is > > indicating different meanings then just random noise. > > But there is no descript

Re: [R] Distributions Comparison

2008-10-28 Thread Rubén Roa-Ureta
Igor Telezhinsky wrote: Dear all, I have recently found out about the R project. Could anyone tell me if it is possible to make the comparison of two distributions using R packages? The problem is TRICKY because I have the distributions of measurements and each measurement in the given distribut

Re: [R] Saving a 3d array into a matlab file

2008-10-28 Thread Henrik Bengtsson
Hi, On Tue, Oct 28, 2008 at 2:25 PM, Minho Chae <[EMAIL PROTECTED]> wrote: > Dear R users, > > I am tryting to save an 3d array to a matlab file like the following. > > A <- array(1:24, c(2,3,4)) > writeMat(filename, A=A) > > But if I load the mat file from Matlab, it is not 3d matrix anymore. > D

Re: [R] Repository missing hmisc

2008-10-28 Thread stephen sefick
I couldn't download this either with install.packages from two or three mirrors- I haven't tried since thinking that the binaries were being built for windows and mac ... If this isn't right I would be interested. On Tue, Oct 28, 2008 at 6:10 PM, David Huffer <[EMAIL PROTECTED]> wrote: > I'm try

Re: [R] Distributions Comparison

2008-10-28 Thread Rolf Turner
On 29/10/2008, at 10:53 AM, Igor Telezhinsky wrote: Dear all, I have recently found out about the R project. Could anyone tell me if it is possible to make the comparison of two distributions using R packages? The problem is TRICKY because I have the distributions of measurements and eac

Re: [R] Sweave, Bibtex, package references

2008-10-28 Thread ONKELINX, Thierry
Paolo, Try using \nocite{*}. That will cite every entry in your bibliography. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest Cel biometrie, methodo

[R] Repository missing hmisc

2008-10-28 Thread David Huffer
I'm trying to install the Hmisc package, however it doesn't appear in the list after "utils:::menuInstallPkgs()". Any help? -- David __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide htt

Re: [R] How to add a data line (series) to a plot using add=TRUE

2008-10-28 Thread stephen sefick
add is not a graphical parameter On Tue, Oct 28, 2008 at 3:49 PM, Mark Na <[EMAIL PROTECTED]> wrote: > Hello, > > I'd like to use the add=TRUE parameter to add a second data line (series) to > an existing plot, but R is giving me an error (see below). > try (untested): plot(aspen_sort,ylim=c(1,10

[R] Saving a 3d array into a matlab file

2008-10-28 Thread Minho Chae
Dear R users, I am tryting to save an 3d array to a matlab file like the following. A <- array(1:24, c(2,3,4)) writeMat(filename, A=A) But if I load the mat file from Matlab, it is not 3d matrix anymore. Does anyone know how I can preserve the 3d structure? Any help will be greatly appreciated.

[R] Distributions Comparison

2008-10-28 Thread Igor Telezhinsky
Dear all, I have recently found out about the R project. Could anyone tell me if it is possible to make the comparison of two distributions using R packages? The problem is TRICKY because I have the distributions of measurements and each measurement in the given distribution has its own error, whi

[R] outputting (writing) output into a dataframe

2008-10-28 Thread Farrel Buchinsky
I have solved this problem once before but don't recall exactly how. Is there a url that shows how? What I want to do now is quite specific but my query is actually very general There are many functions in which one specifies several parameters and an output is generated. Well what happens if one

Re: [R] Fitting weibull and exponential distributions to left censoring data

2008-10-28 Thread Fernando Marmolejo Ramos
Dear Barja Have you looked at "gamlss" y "gamlss.dist" libaries? There some functions called WEI, WEI2, and WEI3 Cheers, Fer __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www

Re: [R] does qplot works with Sweave?

2008-10-28 Thread hadley wickham
You need to explicitly print ggplot objects: print(qplot(carat,price,data=dsmall,colour = color)) Hadley On Tue, Oct 28, 2008 at 3:43 PM, Felipe Carrillo <[EMAIL PROTECTED]> wrote: > Hi Hadley: > I'm practicing writing a document using MikTex but can't make qplot to > execute. It works when usin

[R] get coordinates to multi file

2008-10-28 Thread Alessandro
Hi All, I need to understand How I could to get coordinate to multi-files in a some code: I import X,Y,Z txt files (=vegetation_2.txt, vegetation_3.txt, vegetation_4.txt, etc etc.) with for(i in 2:4) assign(paste('vegetation',i,sep='_'),read.delim(paste('vegetation_',i,'.txt' ,sep='')))

[R] does qplot works with Sweave?

2008-10-28 Thread Felipe Carrillo
Hi Hadley: I'm practicing writing a document using MikTex but can't make qplot to execute. It works when using 'plot' though. Is this a normal behaviour? \documentclass[9pt]{article} \title{ggplot2 example} \begin{document} \maketitle \section*{Examples of using ggplot2} The goal is to be ab

Re: [R] acf() plots of a data.frame: what meaning?

2008-10-28 Thread David Scott
On Tue, 28 Oct 2008, Oliver Bandel wrote: Hello, what is the meaning of the plots of an acf()-plot, when using a data frame as argument? The result is NOT obvious for me. There are combinations of the columns of a dataframe in the reulting plot. But an acf() is just defined for onyl one time-

Re: [R] acf(): meaning of the blue horizontal lines

2008-10-28 Thread David Scott
On Tue, 28 Oct 2008, Oliver Bandel wrote: Hello, what are they meaning? It could be something that would show a threshhold above which the result is indicating different meanings then just random noise. But there is no description on the definition of those lines, so it means nothing, if it is

Re: [R] slightly OT: (un)supervised clustering?

2008-10-28 Thread Dylan Beaudette
On Tuesday 28 October 2008, viktoras didziulis wrote: > Hi, > > my question is not exactly about R... What I am looking for are hints > and directions on suitable methods (available in R or elsewhere) to > solve a grouping (or pattern recognition) problem of environmental > features in an environm

Re: [R] repeated measure one way anova

2008-10-28 Thread JohnLi
Thanks a lot ! I learn a lot about R while I solve my problem according to those provided informatin. John > Dieter Menne wrote: >> JohnLi site.uottawa.ca> writes: >>> I am a new comer for Statistics R. I am using R for one way repeated >>> measure anova, for example, on the following data cons

[R] How to add a data line (series) to a plot using add=TRUE

2008-10-28 Thread Mark Na
Hello, I'd like to use the add=TRUE parameter to add a second data line (series) to an existing plot, but R is giving me an error (see below). This code: > rap<-plot(aspen_sort,ylim=c(1,1),log="y") ...produces the plot to which I'd like to add the second line. But this code: > rap<-plot

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Peter Dalgaard
David M Smith wrote: [] > "David" is the sixth most common first name in the US (with a frequency of > 2.363%) and "Smith" is the most common last name (with a frequency of > 1.006%). Now, presumably first and last names aren't chosen independently > (choice of a first name that in the same

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Ben Bolker
David M Smith revolution-computing.com> writes: > > David Andrew Smith from XL Solutions emailed me privately to clear up the > confusion, too (thanks, David A). To avoid any future confusion, I'll use > my middle initial (M) from now on. > > To get this thread back to R and statistics, there

Re: [R] color individual bar of histogram?

2008-10-28 Thread Peter Dalgaard
Henrique Dallazuanna wrote: > Other option without rect is: > > h <- hist(x, plot = FALSE) > bin <- as.numeric(cut(latest.ob, h$breaks)) > plot(h, col = replace(rep("white", length(h$breaks) - 1), bin, "red")) Nice. Some small improvements seem possible: h <- hist(x, plot = FALSE) bin <- cut(lat

Re: [R] Need help

2008-10-28 Thread JohnLi
Thank you very much! I solve my problem according these information. Guichong > Hello John.After a simple search in google, I came up with this: > http://www.ats.ucla.edu/stat/R/seminars/Repeated_Measures/repeated_measures.htm > > Enjoy, > Tal. > > > > On Mon, Oct 27, 2008 at 6:15 PM, JohnLi <[EM

Re: [R] How to export text into separate text files

2008-10-28 Thread Lauri Nikkinen
I think I found the solution: for (i in 1:ncol(r)) { cat("Some text,", "\n", file=paste("temp", i, ".txt"), append=T) cat("\n", file=paste("temp", i, ".txt"), append=T) cat("More text, More text, More text", file=paste("temp", i, ".txt"), append=T) cat("\n", file=paste("temp", i, "

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread David M Smith
David Andrew Smith from XL Solutions emailed me privately to clear up the confusion, too (thanks, David A). To avoid any future confusion, I'll use my middle initial (M) from now on. To get this thread back to R and statistics, there might be an interesting counting problem here somewhere. :) Wh

[R] slightly OT: (un)supervised clustering?

2008-10-28 Thread viktoras didziulis
Hi, my question is not exactly about R... What I am looking for are hints and directions on suitable methods (available in R or elsewhere) to solve a grouping (or pattern recognition) problem of environmental features in an environmental gradient as described below. Given environmental samp

Re: [R] color individual bar of histogram?

2008-10-28 Thread Henrique Dallazuanna
Other option without rect is: h <- hist(x, plot = FALSE) bin <- as.numeric(cut(latest.ob, h$breaks)) plot(h, col = replace(rep("white", length(h$breaks) - 1), bin, "red")) On Tue, Oct 28, 2008 at 4:48 PM, Whit Armstrong <[EMAIL PROTECTED]>wrote: > Anyone know a quick way to color one bar of a hi

Re: [R] color individual bar of histogram?

2008-10-28 Thread Whit Armstrong
That's great, Peter. Thanks very much. -Whit On Tue, Oct 28, 2008 at 3:13 PM, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > Whit Armstrong wrote: >> Anyone know a quick way to color one bar of a histogram? >> >> I want to mark the bar in which the most recent observation falls. >> >> So, for inst

Re: [R] color individual bar of histogram?

2008-10-28 Thread Peter Dalgaard
Whit Armstrong wrote: > Anyone know a quick way to color one bar of a histogram? > > I want to mark the bar in which the most recent observation falls. > > So, for instance: > > x <- rnorm(100) > latest.ob <- x[100] > hist(x) > ## how do I mark the bucket that latest.ob falls into? > Something

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread eugene dalt
I attended a course at XLSolutions and glad they can publish a list of recommended R books. IMHO they shouldn't post it here though. This name issue is hilarious given thousands of David Smith out there but to cut the story short, my class at XLSolutions was taught by a David Andrew Smith w

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread eugene dalt
I attended a course at XLSolutions and glad they can publish a list of recommended R books. IMHO they shouldn't post it here though. This name issue is hilarious given thousands of David Smith out there but to cut the story short, my class at XLSolutions was taught by a David Andrew Smith w

Re: [R] Computational problems in R

2008-10-28 Thread Steven McKinney
I misinterpreted the suggestion from Duncan Murdoch, and I think Amy did too, judging from her comment about B possibly being negative. Duncan Murdoch kindly explained the idea off-list, many thanks to him for taking the time. To avoid further confusion, I've set up a basic example explaining the

[R] change data into array

2008-10-28 Thread Suyan Tian
Sorry to bother everyone again with one question about array. I have one data set saved in this way patienttime x1 x2 response 1 1 18 1 1 2 27 1 2 1 36 1 2

[R] color individual bar of histogram?

2008-10-28 Thread Whit Armstrong
Anyone know a quick way to color one bar of a histogram? I want to mark the bar in which the most recent observation falls. So, for instance: x <- rnorm(100) latest.ob <- x[100] hist(x) ## how do I mark the bucket that latest.ob falls into? Thanks, Whit ___

Re: [R] License Question

2008-10-28 Thread James MacDonald
>>> "Barry Rowlingson" <[EMAIL PROTECTED]> wrote: > Oh, btw, I'd fix your blog links: go to: > http://www.experience-rplus.com/blog.asp?b=act and click a > screenshot... Now that is hilarious. A fine example of defensive coding for the web circa 2008 ;-D Best, Jim > > Barry > > __

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Bert Gunter
David et. al: My sentiments are the same, but my apology to you (and I'm glad to learn that it's not you!). -- Bert _ From: David Smith [mailto:[EMAIL PROTECTED] Sent: Tuesday, October 28, 2008 11:09 AM To: Bert Gunter Subject: Re: [R] Recommended R books by XLSolutions Corporation

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Steve_Friedman
I agree with Bert, The R user group lists should not be a vehicle for private commercial vendors. Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 Office (305) 224 - 4282 [EMAIL PROTECTED]

Re: [R] Sweave Error

2008-10-28 Thread Martin Morgan
[EMAIL PROTECTED] wrote: dear R users, I am using sweave to generate report for my data analysis. I recently updated R ro 2.8.0, and now I have the following results when compile the the tex file generated from R. Hi Claudio How do you invoke latex? On many systems, this works % R CMD pdfla

Re: [R] abind

2008-10-28 Thread Tony Plate
It looks like you are trying to construct a ragged array, where the extent of the dimensions varies. However, in R, ordinary arrays have a regular structure, e.g., the rows of a matrix always have the same number of columns. This is the kind of object abind() constructs. So, to bind your two

Re: [R] Dose Profile

2008-10-28 Thread Bert Gunter
Why don't you ask your "someone at work" to help, as you appear to have local resources handy? -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Laura Bonnett Sent: Tuesday, October 28, 2008 9:33 AM To: r-help@r-project.org Subject: [R] Dose

Re: [R] Pkg rgl: installation fails because of x11

2008-10-28 Thread megha patnaik
This works. After installing the dependencies and installing RGL from inside R, library(rgl) rgl.open() should open a display window for you. 2008/10/27 Ben Bolker <[EMAIL PROTECTED]> > Duncan Murdoch stats.uwo.ca> writes: > > > > > Matthieu Stigler wrote: > > > Hello > > > > > > Im trying to

Re: [R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread Bert Gunter
David: IMHO, this is beginning to be commercial abuse of this list. Please no more advertising. (Others may have different opinions). -- Bert -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Tuesday, October 28, 2008 9:33 AM To:

Re: [R] abind

2008-10-28 Thread Henrique Dallazuanna
You can use lists for this: list(a, b) On Tue, Oct 28, 2008 at 3:39 PM, Suyan Tian <[EMAIL PROTECTED]>wrote: > I am trying to combine two arrays with different dimensions into one. For > example > > The first one is > 1 2 3 > 4 5 6 > > The second one is > 7 8 > 9 10 > > The resulted one woul

[R] abind

2008-10-28 Thread Suyan Tian
I am trying to combine two arrays with different dimensions into one. For example The first one is 1 2 3 4 5 6 The second one is 7 8 9 10 The resulted one would be like , , 1 1 2 3 4 5 6 , , 2 7 8 9 10 I used abind to do this, but failed. Could somebody please let me know how to do thi

[R] Sweave Error

2008-10-28 Thread [EMAIL PROTECTED]
dear R users, I am using sweave to generate report for my data analysis. I recently updated R ro 2.8.0, and now I have the following results when compile the the tex file generated from R. This is pdfTeXk, Version 3.141592-1.40.3 (Web2C 7.5.6) %&-line parsing enabled. entering extended mode (

Re: [R] Marginal effects in negative binomial

2008-10-28 Thread John Fox
Dear Roberto, It depends on what you mean by "marginal effects," but the effects package will calculate what I would prefer to call "partial effects" for terms in a generalized linear model; it should work fine with objects produced by glm.nb. I hope this helps, John ---

[R] Random Forest Bug

2008-10-28 Thread Roeder Jens (CR/AEM5)
Dear help list, I think I found a bug a the R Random Forest. Hopefully, you are able to reproduce it. I use R version 2.7.2 and RF version 4.5-27. This is a minimal code to describe the problem: library(randomForest) tries <- 20 dimension <- 20 n <- 200 outlyingness <- rep(NaN,tries) for (o_numbe

[R] Recommended R books by XLSolutions Corporation

2008-10-28 Thread dsmith
We've listed books we've recommended to our introductory R/R-PLUS courses attendees [1]http://www.experience-rplus.com/books.asp If you'd like to recommend a book, email [EMAIL PROTECTED] R-PLUS Rocks! R-PLUS Team 24 Hours Tech Support Think Fast, Think B

Re: [R] Fixing an only one coefficient in an ARIMA model

2008-10-28 Thread Kellie Wills
Hello Yohann, Try arima(lh, order=c(3,0,0), fixed=c(0, NA, NA, NA)) The NA entries in "fixed" (for AR2, AR3, and intercept) will be allowed to vary. Kellie Wills Engineering Service Manager REvolution Computing [EMAIL PROTECTED] On Tue, Oct 28, 2008 at 7:10 AM, Yohann MOREAU <[EMAIL PROTECTED]

Re: [R] Sweave, Bibtex, package references

2008-10-28 Thread Paolo Sonego
Dear Thierry, Thank you very much for you very useful advice! I managed to solve my specific problem with this code: <>= package <- grep("^package:", search(), value = TRUE) keep <- sapply(package, function(x) x == "package:base" || !is.null(attr(as.environment(x), "path"))) package <- sub("^p

[R] acf(): meaning of the blue horizontal lines

2008-10-28 Thread Oliver Bandel
Hello, what are they meaning? It could be something that would show a threshhold above which the result is indicating different meanings then just random noise. But there is no description on the definition of those lines, so it means nothing, if it is not clearly defined. Where can I find a deta

Re: [R] License Question

2008-10-28 Thread rplus
Thank you Simon, We are working on the linux version as well and will roll out more products during our conference in Las Vegas next year. Regards R-PLUS Team 24 Hours Tech Support Think Fast, Think Big, Without The High Cost! [1]www.Experience-Rplus.com/Conference.asp

[R] acf() plots of a data.frame: what meaning?

2008-10-28 Thread Oliver Bandel
Hello, what is the meaning of the plots of an acf()-plot, when using a data frame as argument? The result is NOT obvious for me. There are combinations of the columns of a dataframe in the reulting plot. But an acf() is just defined for onyl one time-sreies. So what do the plots mean? I did not

Re: [R] odbcConnectExcel2007

2008-10-28 Thread Chrischizinski
Agnolucci, Paolo wrote: > > Hi > > > > sorry if this is really basic but I am just starting on R. > > > > Can anyone point me at how to write R objects into Excel 2007 files. I > have seen how to set up a connection to a file through > odbcConnectExcel2007(xls.file, readOnly = FALSE, ...

[R] Dose Profile

2008-10-28 Thread Laura Bonnett
Hi Everyone, I have data in a long format e.g. there is one row per patient but each follow-up appointment is included in the row. So, a snippet of the data looks like this: TrialNo Drug Sex Rand Adate1 Date1 Dose1 Time1 Adate2 Date2 Dose2 Time2 B1001029 LTG M 15719 30/04/2003 15825 150 10

[R] R training Courses

2008-10-28 Thread Sharon Lazenby
All, There are still places available on the following courses. For more detailed information and a registration form please contact [EMAIL PROTECTED], or visit our website at www.mango-solutions.com . The R Language

Re: [R] How to export text into separate text files

2008-10-28 Thread Lauri Nikkinen
Thanks for the response. I have used the method Jim Holtman suggested when exporting data.frames. In the current situation, I'm combining text and data from data.frame, and I don't know how to separate the results into the separate text files. I can output everything from the loop into a single fil

Re: [R] Source code for ppr (Projection Pursuit Regression)

2008-10-28 Thread Erin Hodgess
To get the references, check ?ppr For the code: > methods(ppr) [1] ppr.default* ppr.formula* Non-visible functions are asterisked > getAnywhere(ppr.default) A single object matching 'ppr.default' was found It was found in the following places registered S3 method for ppr from namespace stat

Re: [R] License Question

2008-10-28 Thread David Smith
[replying to r-help as well] Hi Tom, > Thanks. I am looking forward to more information. We're putting together the academic program right now, but our website will be updated with more details in the next few weeks. As I mentioned though our general attitude is that RPro should be available t

Re: [R] Best way of figuring out whether graphical elements overlap?

2008-10-28 Thread Greg Snow
There is also the spread.labs function in the TeachingDemos package that uses a different method from the plotrix function and should not move any labels that are not overlapping. There are also the dynIdentify and TkIdentify functions in the same package that allow you to interactively move la

[R] Source code for ppr (Projection Pursuit Regression)

2008-10-28 Thread Arvind Iyer
Dear R users, I am looking for the source code of the implementation of ppr (Projection Pursuit Regression) in R. It will be great if citations of the source papers on which the implementation is based, are also provided. Thank you, Arvind Iyer, Grad student, Deptt. of Biomedical Engineering Vit

[R] group sequential analysis - stopping for futility

2008-10-28 Thread Alice . Pressman
Hello - I am wondering if anyone has written some R code to calculate futility bounds for a group sequential analysis of clinical trials data. The library ldbounds has a function 'bounds' which calculates the effectiveness stopping bounds for various spending functions, but it does not app

Re: [R] Transferring results from R to MS Word

2008-10-28 Thread Marc Schwartz
on 10/28/2008 10:23 AM [EMAIL PROTECTED] wrote: > John Kane wrote: > >> --- On Mon, 10/27/08, Tom Backer Johnsen <[EMAIL PROTECTED]> >> wrote: >> >>> From: Tom Backer Johnsen <[EMAIL PROTECTED]> Subject: Re: [R] >>> Transferring results from R to MS Word To: "Greg Snow" >>> <[EMAIL PROTECTED]> Cc

[R] Fitting weibull and exponential distributions to left censoring data

2008-10-28 Thread Borja Soto Varela
Dear R-users I have some datasets, all left-censoring, and I would like to fit distributions to (weibull,exponential, etc..). I read one solution using the function survreg in the survival package. i.e survreg(Surv(...)~1, dist="weibull") but it returns only the scale parameter. Does anyone know h

Re: [R] Transferring results from R to MS Word

2008-10-28 Thread heimdal
John Kane wrote: > --- On Mon, 10/27/08, Tom Backer Johnsen <[EMAIL PROTECTED]> wrote: > > > From: Tom Backer Johnsen <[EMAIL PROTECTED]> > > Subject: Re: [R] Transferring results from R to MS Word > > To: "Greg Snow" <[EMAIL PROTECTED]> > > Cc: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]> > > Receive

Re: [R] ggplot2: how to combine position=stack and position=dodge in a single graph?

2008-10-28 Thread hadley wickham
On Fri, Oct 24, 2008 at 5:18 PM, Elena Schulz <[EMAIL PROTECTED]> wrote: > Hi Hadley, > > thanks a lot for your quick answer. >> You should be able to replicate any dodging layout with facetting > You mean instead of facetting by years, facetting by months? I will try this > an see how the plot loo

[R] how to easily access S-PLUS data sets in R?

2008-10-28 Thread Mark Daniel Ward
Dear R-help colleagues, Greetings! I'm interested in using some of the data sets from S-PLUS in a class that I'm teaching on R. Our university has a license for S-PLUS, and I recognize that data sets from S-PLUS can be individually exported from S-PLUS and then imported into R. Is there a

[R] odbcConnectExcel2007

2008-10-28 Thread Agnolucci, Paolo
Hi sorry if this is really basic but I am just starting on R. Can anyone point me at how to write R objects into Excel 2007 files. I have seen how to set up a connection to a file through odbcConnectExcel2007(xls.file, readOnly = FALSE, ...) but it doesn't say anything on how to write data

[R] Best way of figuring out whether graphical elements overlap?

2008-10-28 Thread Johannes Graumann
Hi all, I'm plotting impulses, where some of them should have labels hovering above them. I know of plotrix' spread.labels function, but would like to save that for instances where there truely is to little space for the label. Does anybody have any hints what' the most efficient way might be to

Re: [R] Stuck with FFT

2008-10-28 Thread stephen sefick
#look at dput this is an easy way to share #data or another great way is to make a #fake data series x <- (c(105L, 115L, 140L, 180L, 230L, 300L, 360L, 418L, 450L, 480L, 495L, 500L, 495L, 480L, 450L, 410L, 360L, 300L, 230L, 180L, 140L, 115L, 105L, 100L, 105L, 115L, 140L, 180L, 230L, 300L, 360L, 410L

Re: [R] Transferring results from R to MS Word

2008-10-28 Thread John Kane
--- On Mon, 10/27/08, Tom Backer Johnsen <[EMAIL PROTECTED]> wrote: > From: Tom Backer Johnsen <[EMAIL PROTECTED]> > Subject: Re: [R] Transferring results from R to MS Word > To: "Greg Snow" <[EMAIL PROTECTED]> > Cc: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]> > Received: Monday, October 27, 2008, 5:4

Re: [R] How to export text into separate text files

2008-10-28 Thread John Kane
"jim holtman" <[EMAIL PROTECTED]> just answered a similar question for someone who wanted to input txt files. His suggestion is below. It should work for output as well Jim Holtman's suggestion. - for (i in 2:4){ input <- r

[R] How to export text into separate text files

2008-10-28 Thread Lauri Nikkinen
Hello, I'm producing text from my data.frame using cat function. I would like to use for loop to export each column in my data.frame into separate text files. Here is the example code r <- t(Indometh) for (i in 1:ncol(r)) { cat("Some text,", "\n") cat("\n") cat("More text, More tex

[R] Fixing an only one coefficient in an ARIMA model

2008-10-28 Thread Yohann MOREAU
Good afternoon, I would like fitting an ARIMA model without the first coefficient. For example, I want to fit an AR(3) like this : y[t]=a[1]*y[t-1]+a[2]*y[t-2]+a[3]*y[t-3], where a[1]=0. How can I specify it in the function "arima", if it is possible ? Thank you in advance. Yohann Moreau

Re: [R] Transferring results from R to MS Word

2008-10-28 Thread Eric Lecoutre
Tom, Sorry to react only at tis point ; it makes months I didn"t read R-help, having unfortunalety to work for SAS. Though I am not anymore the maintainer of the package, I am the creator of it and may provide some insights. - There are a lot of HTML methods, effectively hidden (namespace) ; you

[R] bug calling setRepositories in R 2.8.0

2008-10-28 Thread Kjetil Halvorsen
The following happened when calling setRepositories on R 2.8.0 compiled as a shared library on ubuntu intrepid RC: > setRepositories(gra=TRUE) > The program 'R' received an X Window System error. This probably reflects a bug in the program. The error was 'BadWindow (invalid Window parameter)'. (

Re: [R] gsubfn, strapply, REGEX Problem

2008-10-28 Thread Johannes Graumann
Thank you! Joh Gabor Grothendieck wrote: > Have just made an improvement to the development > version to ignore escaped left parens in the regexp > in setting the backref default. This improvement > should address your problem so that this now > works without errors: > > library(gsubfn) > # o

Re: [R] gsubfn, strapply, REGEX Problem

2008-10-28 Thread Gabor Grothendieck
Have just made an improvement to the development version to ignore escaped left parens in the regexp in setting the backref default. This improvement should address your problem so that this now works without errors: library(gsubfn) # overwrite relevant function with devel version of it source("h

  1   2   >