Re: [R] including images in a plot

2007-09-28 Thread Frede Aakmann Tøgersen
Perhaps http://finzi.psych.upenn.edu/R/Rhelp02a/archive/94069.html is a good starting point? Med venlig hilsen Frede Aakmann Tøgersen > -Oprindelig meddelelse- > Fra: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] På vegne af Steve C. Wang > Sendt: 28. september 2007 01:15 > Til:

Re: [R] TukeyHSD doubts

2007-09-28 Thread Prof Brian Ripley
This is a repost of https://stat.ethz.ch/pipermail/r-help/2007-September/141727.html Please do study the posting guide to see why you did not get an answer and what to do when you do not. It is nothing to do with TukeyHSD, as the differences are there in the means. It is clear that your cova

Re: [R] ifelse and dates do not work together: What workaround?

2007-09-28 Thread Gabor Grothendieck
See ?replace to do it in one line. On 9/27/07, Farrel Buchinsky <[EMAIL PROTECTED]> wrote: > I encountered the above problem. I went to the help files and > discovered the reason why. My insight as to why it was happening did > not immediately provide me with a solution by which I could accomplish

Re: [R] converting numbers in "YYYYMM" format to last calendar day and last exchange trading day of the month

2007-09-28 Thread Gabor Grothendieck
Using as.yearmon and as.Date.yearmon with frac=1 (both from the zoo package), the last calendar day of the month is: library(zoo) d <- c("200701", "200702") lastday <- as.Date(as.yearmon(d, "%Y%m"), frac = 1) and the zoo quickref vignette("zoo-quickref")

Re: [R] Legend

2007-09-28 Thread S Ellison
>>> P Ehlers <[EMAIL PROTECTED]> 26/09/2007 11:06:54 >>> >a <- .33 >b <- .55 >legend("bottom", fill=c("red","blue"), > legend=c(bquote(p == .(a)), bquote(p == .(b))), bty="n") paste should do it: legend("bottom", fill=c("red","blue"), legend=paste("p=", c(a,b), sep=""), bty="n") and would be

[R] RMySQL installing problems

2007-09-28 Thread andreiabb
Dear forum, I'm trying to install RMySQL but I'm having problems in loading it, here is the message that I'm getting in R install.packages("RMySQL", dependencies=TRUE) trying URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.5/RMySQL_0.6-0.zip' Content type 'application/zip' length 39

[R] BATCH run library problem

2007-09-28 Thread Flash
I have installed R-2.5.1 on my linux machine (Ubuntu LTS 6.01) by compiling the source. I have no problem to install packages on my local folder, /home/flash/R_LIBS And I have no problem to load the package inside of R. for example, >library(cluster) However, when I'm trying to run batch comma

[R] Converting to mode numeric within a matrix

2007-09-28 Thread Sergey Goriatchev
Hello, I create a matrix: best <- matrix(0, ncol=2, nrow=num.selected, dimnames=list(the.best$.Name, c("Probability(%)", "Inside"))) best[,1] <- as.numeric(the.best$Total*100) best[,2] <- ifelse(the.best$Weight==0, "No", "Yes") What I want is the second column of mode numeric, but it is of mod

Re: [R] ifelse and dates do not work together: What workaround?

2007-09-28 Thread Peter Dalgaard
Gabor Grothendieck wrote: > See ?replace to do it in one line. > Also, don't miss the bleeding obvious (or so one might think...): > z <- ifelse(TRUE, as.Date("2007-01-01"), as.Date("2007-01-02")) > class(z) <- "Date" > z [1] "2007-01-01" In the special case of censoring, you can also get away

[R] Writing successful self-starting models for nlme

2007-09-28 Thread I R Cleasby
Hi fellow R users, My problem is the following: I wish to use the nlme package to analyse the growth curves of house sparrow chicks, especially in relation to factors such as sex, position in the weight hierarchy etc. as fixed effects and chick ID nested within BroodID as random effects.

[R] Printing the contents of a list

2007-09-28 Thread Michael Stegh
Dear List, is there a way to send the contents of a variable (e.g. a list) directly to a laser printer. I am aware that I could save the variable to a ascii-file and then print the file, but I would prefer to print from R directly. I have been looking for a way to do so for the better part o

[R] Creating nice looking lists: how?

2007-09-28 Thread Sergey Goriatchev
Hello, For my functions I want to create output similar in appearance to that of what you get when you print a summary of lm model: Residuals: Min1Q Median3Q Max -0.209209 -0.043133 0.001793 0.044105 0.234750 Coefficients: Estimate Std

Re: [R] New R website: R-Cookbook.com

2007-09-28 Thread hadley wickham
Hi Jeff, That looks like a nice initiative. However, if you are interested in getting contributions from the community, it might be good to spell out how others might use the content of the site. Currently you have copyright r-cookbook.com, but maybe you could consider a creative commons (http:/

Re: [R] kurtosis

2007-09-28 Thread bbolker
S Bina wrote: > > Hi, > > I cannot find the function kurtosis. Is it sth additional I am meant to > download? I use the MacOS X version of R. > > Many thanks > Samira > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listin

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread jiho
On 2007-September-28 , at 16:57 , Frank E Harrell Jr wrote: > jiho wrote: >> On 2007-September-28 , at 15:18 , Paul Smith wrote: >>> On 9/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > I know how to export graphics as pdf files and then how to include > them in LaTeX documents. Ho

Re: [R] New R website: R-Cookbook.com

2007-09-28 Thread Jeff Spies
Hi Philippe, I don't want to be competing in any way with any of the fantastic official R resources that exist; I only want to supplement them. Although this list is probably not the best place to discuss the proper use of wiki's for documentation/learning, I'll make a few comments on th

[R] orientlib

2007-09-28 Thread Abu Naser
Hi All user, I have been using R-2.5.1. dose orientlib support this version? I would like to try. it. I have been wondering how to install the library. With regards, abu _ Celeb spotting – Play CelebMashup and win cool prizes

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Frank E Harrell Jr
Paul Smith wrote: > On 9/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: >>> I know how to export graphics as pdf files and then how to include >>> them in LaTeX documents. However, I do not know how to do in order to >>> have the text of the graphics written with the font selected for the >>>

Re: [R] Creating nice looking lists: how?

2007-09-28 Thread Ptit_Bleu
Hi, I don't know if it will help you but to retrieve the slope of lm(y~x), I use : coefficients(lm(y~x))[2] ([1] for the intercept) I should not tell this but It took me a long time to find this obvious thing. Have a nice week-end, Ptit Bleu. ---

[R] off topic: Job advert - computational biologist

2007-09-28 Thread Crispin Miller
Hi, I hope you don't mind me posting this to the list: We are seeking a postdoctoral scientist to join the Applied Computational Biology and Bioinformatics group at the Paterson Institute for Cancer Research. The role will be to develop novel approaches for the analysis of high throughput genomic

[R] simple matching with R

2007-09-28 Thread Birgit Lemcke
Hello! I am R beginner and I have a question obout a simple matching. I have to datasets that i read in with: MalVar29_37<-read.table("MalVar29_37.csv", sep = ";") FemVar29_37<-read.table("FemVar29_37.csv", sep = ";") They look like this and show binary variables: V1 V2 V3 V4 V5 V6 V7 V8

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Frank E Harrell Jr
jiho wrote: > On 2007-September-28 , at 15:18 , Paul Smith wrote: >> On 9/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: I know how to export graphics as pdf files and then how to include them in LaTeX documents. However, I do not know how to do in order to have the text

[R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Paul Smith
Dear All, I know how to export graphics as pdf files and then how to include them in LaTeX documents. However, I do not know how to do in order to have the text of the graphics written with the font selected for the LaTeX document. Is that possible? Thanks in advance, Paul _

Re: [R] Help with "old" function in "base"

2007-09-28 Thread Prof Brian Ripley
Kernel density estimation is in package stats, and has been for over 3 years. That package has a C entry point registered as R_band_den_bin which may be what you are looking for. On Fri, 28 Sep 2007, Van Wyk, Jaap wrote: > Hi > > I think there was a built-in function (perhaps C-code) called > "

[R] Help with "old" function in "base"

2007-09-28 Thread Van Wyk, Jaap
Hi I think there was a built-in function (perhaps C-code) called "band_den_bin" in the base library (used for kernel density estimators). Is there a newer version of this function. I have an old routine using this but it won't run in the newer version of R. Any help is much appreciated. Jacob

Re: [R] Plots with discontinuity balls

2007-09-28 Thread Duncan Murdoch
Paul Smith wrote: > On 9/28/07, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > >>> Can R plot graphs like the one at >>> >>> http://www.mathwords.com/f/f_assets/floor_graph.gif >>> >>> with the balls at the discontinuity points? >>> >> You can use segments() to draw the segments and symbols(

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread jiho
On 2007-September-28 , at 15:18 , Paul Smith wrote: > On 9/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: >>> I know how to export graphics as pdf files and then how to include >>> them in LaTeX documents. However, I do not know how to do in >>> order to >>> have the text of the graphics w

Re: [R] Is there a model like that in R?

2007-09-28 Thread Gustaf Rydevik
On 9/28/07, Manuel Ramon <[EMAIL PROTECTED]> wrote: > Hi to everyone, > > > > I am starting to work with a model that is not familiar to me. > > The model would be like that: > > > >y = . - a*max(Q(i) - Q(0), 0) + . > > > > where Q(i) is the accumulated effect of a variable at time i and Q(0) a

[R] lmer giving negative, or no, estimated standard errors

2007-09-28 Thread Chaloner, Kathryn
R Users, Emine Bayman sent this out earlier and we do not think it went through. Appologies if it did. We want to fit GLMM with lmer with binomial data and a one-way random effects model (overall mean is a fixed effect and there are random effects for each binomial). We are using the La

[R] Drawing functions on Cartesian coordinate systems

2007-09-28 Thread Paul Smith
Dear All, Can R draw plots of functions on a Cartesian coordinate system with axes like the ones shown at http://en.wikipedia.org/wiki/Image:Cartesian-coordinate-system-with-circle.svg ? I have already searched the R web-site, but found nothing. Thanks in advance, Paul __

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread hadley wickham
> Yes there is harm. But to make bold lines, easy to read titles is fine. > See the spar function in > http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see > the setps, ps.slide, and setpdf functions in the Hmisc package. I was interested to see that you have code for drawin

Re: [R] simple matching with R

2007-09-28 Thread Jeffrey Robert Spies
> The next problem for me is now to deal with the NAs. I thought > perhaps it is possible to exclude the variable from the row > comparison if in one of the rows is an NA? If you exclude the NAs in one dataset, you'll need to exclude the exact same row in the other dataset. The question to

[R] error message from dbConnect, using ROracle

2007-09-28 Thread Muller, John H
I am trying to use ROracle (v 0.5-8) with R on Windows (v 2.5.1) I get an error message when I run the following code > library(DBI) > library(ROracle) > m <- dbDriver("Oracle") > con <- dbConnect(m, dsn='', uid='', pwd='') Note masked real values with ###

[R] Elasticity

2007-09-28 Thread Giusy
Hello, I'm starting with data analysis with R, and I'd like to know if there is in any package a function to have price elasticity of an item. For example I have information about sale quantity and price, and I want to know if the series is sensitive to the price. Sale Quantity Price 34

Re: [R] RMySQL NA/NULL value storage error

2007-09-28 Thread Chris Stubben
Adam Wilson-4 wrote: > > I am running R 2.5.1, RMySQL 0.6 , and DBI 0.2-3 on Windows XP > > Like others, I am having trouble with NA/Null value conversions between R > and a MySQL database via DBI, but I could not find my exact problem in the > archives. Most of the time NA values in R get tra

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Paul Smith
On 9/28/07, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote: > >> I know how to export graphics as pdf files and then how to include > >> them in LaTeX documents. However, I do not know how to do in > >> order to > >> have the text of the graphics written with the font selected for the

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Paul Smith
On 9/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > > I know how to export graphics as pdf files and then how to include > > them in LaTeX documents. However, I do not know how to do in order to > > have the text of the graphics written with the font selected for the > > LaTeX document. Is t

Re: [R] kurtosis

2007-09-28 Thread Erich Neuwirth
According to the docs package e1071 has kurtosis. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-containe

Re: [R] plot or boxplot!

2007-09-28 Thread John Kane
--- Samuel Okoye <[EMAIL PROTECTED]> wrote: > Hello, > > if we suppose that > > times <- c("2006-05-14", "2006-06-12", > "2006-06-12", "2006-05-14", "2006-05-14", > "2006-06-12") > value <- c(2,3,1,4,3,1) > > then with > > plot(times, value) Have you tried this? I thin

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Prof Brian Ripley
On Fri, 28 Sep 2007, Paul Smith wrote: > Dear All, > > I know how to export graphics as pdf files and then how to include > them in LaTeX documents. However, I do not know how to do in order to > have the text of the graphics written with the font selected for the > LaTeX document. Is that possibl

Re: [R] different colors for two wireframes in same plot

2007-09-28 Thread deepayan . sarkar
On 9/27/07, marcg <[EMAIL PROTECTED]> wrote: > Thanks a lot > > This already looks nice and I already checked the ?wireframe, but with no > examples and as newcommer its hard to find out a correct code. > > If we set drape=F in the example: > > g <- expand.grid(x = 1:10, y = 5:15, gr = 1:2) > g$z <

Re: [R] Creating nice looking lists: how?

2007-09-28 Thread Charles C. Berry
On Fri, 28 Sep 2007, Sergey Goriatchev wrote: > Hello, > > For my functions I want to create output similar in appearance to that > of what you get when you print a summary of lm model: > > Residuals: > Min1Q Median3Q Max > -0.209209 -0.043133 0.001793 0.04

Re: [R] simple matching with R

2007-09-28 Thread Jeffrey Robert Spies
Not sure how you want to handle the NAs, but you could try the following: #start MalVar29_37 <- read.table(textConnection("V1 V2 V3 V4 V5 V6 V7 V8 V9 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 NA NA NA NA NA NA NA NA NA 0 1 0 0 0 1 0 0 0"), header=TRUE)

[R] errorbar il lattice plot

2007-09-28 Thread Federico Calboli
Hi Everyone, I would like to add errorbars to a lattice plot. I already have a function that adds error bars to a plot (originally written by George Gilchrist) which takes as arguments the coordinates of the point and then the size of the error bar. I would like to integrate it in my lattice co

[R] kurtosis

2007-09-28 Thread S Bina
Hi, I cannot find the function kurtosis. Is it sth additional I am meant to download? I use the MacOS X version of R. Many thanks Samira __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Converting to mode numeric within a matrix

2007-09-28 Thread R Help
Your problem is that a matrix can't be a mix of data types, but a data frame can. Try the following code. newBest = as.data.frame(best) Now you can convert the column of numbers in newBest to numbers from the factors they start as. Check the help file for data.frame(...) for more details of how

Re: [R] Plots with discontinuity balls

2007-09-28 Thread Paul Smith
On 9/28/07, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > > Can R plot graphs like the one at > > > > http://www.mathwords.com/f/f_assets/floor_graph.gif > > > > with the balls at the discontinuity points? > > You can use segments() to draw the segments and symbols() to draw the > balls. For example

Re: [R] New R website: R-Cookbook.com

2007-09-28 Thread Philippe Grosjean
Hello Jeff, Good initiative,... but why not to put this in the official R Wiki (http://wiki.r-project.org)? There is a section named 'tips' dedicated to such little recipes (http://wiki.r-project.org/rwiki/doku.php?id=tips:tips). It should be better to centralize all these little tips, don't y

[R] plot graph with error bars trouble

2007-09-28 Thread Marcelo Laia
Hi, I have a data set like this: MutantRepTime OD 02H02100.029 02H02200.029 02H02300.023 02H02180.655 02H02280.615 02H02380.557 02H021121.776 02H02212 1.859 02H023121.668 02H02

[R] Is there a model like that in R?

2007-09-28 Thread Manuel Ramon
Hi to everyone, I am starting to work with a model that is not familiar to me. The model would be like that: y = . - a*max(Q(i) - Q(0), 0) + . where Q(i) is the accumulated effect of a variable at time i and Q(0) a threshold above it there is effect on y. The value of Q(i) could b

[R] File selection by condition

2007-09-28 Thread Svempa
I have a large number of textfiles, and one matix in every texfile. I now want to find an easy way to select those textfiles that contain at least one value over a certain limit, and collect those matrices in a vector or something alike. How do I do that? -- View this message in context: http:/

Re: [R] Cairo on windows

2007-09-28 Thread Michael Lawrence
On 9/27/07, Moisan Yves <[EMAIL PROTECTED]> wrote: > > > My experience is that cairoDevice is a lot slower than Cairo, > especially > on Windows, and about equally flaky. If you see how many layers are > involved with Cairo on Windows you will not be surprised. > > Actually, I ended up trying to l

Re: [R] Drawing functions on Cartesian coordinate systems

2007-09-28 Thread Charles Annis, P.E.
Yes, R can do that. Well, actually YOU can do that using R. But it is hard to believe that you looked very hard before writing. Did you look at these R functions? ?plot ?line ?points ?arrows Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.Statistical

Re: [R] Expressing number in percentage

2007-09-28 Thread S Ellison
Sorry! paste(x*100,"%") >>> Marc Schwartz <[EMAIL PROTECTED]> 27/09/2007 15:56:05 >>> On Thu, 2007-09-27 at 14:36 +0100, stat stat wrote: > I am wondering if there is any procedure to write a particular value > in Percentage format, still maintaining it's numeric character. for > example I want

Re: [R] simple matching with R

2007-09-28 Thread Birgit Lemcke
Thanks a lot for both solutions of my problem. I tried it immediately and I understood how they are working. The next problem for me is now to deal with the NAs. I thought perhaps it is possible to exclude the variable from the row comparison if in one of the rows is an NA? Furthermore it wou

Re: [R] Expressing number in percentage

2007-09-28 Thread S Ellison
Or paste(x,"%") >>> Marc Schwartz <[EMAIL PROTECTED]> 27/09/2007 15:56:05 >>> On Thu, 2007-09-27 at 14:36 +0100, stat stat wrote: > I am wondering if there is any procedure to write a particular value > in Percentage format, still maintaining it's numeric character. for > example I want to write

Re: [R] plot or boxplot!

2007-09-28 Thread Uwe Ligges
Samuel Okoye wrote: > Hello, > > if we suppose that > > times <- c("2006-05-14", "2006-06-12", "2006-06-12", "2006-05-14", > "2006-05-14", "2006-06-12") > value <- c(2,3,1,4,3,1) > > then with > > plot(times, value) > > we have two boxplots in one graph for 2006-0

Re: [R] ELF file OS ABI invalid yes?????

2007-09-28 Thread Peter Dalgaard
vittorio wrote: > Compilation of MCMCpack under freebsd 6.2 i386 fails because of the following > cryptic error: > > * Installing *source* package 'MCMCpack' ... > checking for C++ compiler default output file name... a.out > checking whether the C++ compiler works... yes > checking whether we ar

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Frank E Harrell Jr
jiho wrote: > > On 2007-September-28 , at 16:57 , Frank E Harrell Jr wrote: >> jiho wrote: >>> On 2007-September-28 , at 15:18 , Paul Smith wrote: On 9/28/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: >> I know how to export graphics as pdf files and then how to include >> them i

[R] transparency of one layer in multiple wireframe plot

2007-09-28 Thread mdgi
Sorry sending this again - is anyone familiar with multiple transparent wireframe plots? I already checked the ?wireframe, but with no examples and as newcommer its hard to find out a correct code. If we set drape=F in the example: g <- expand.grid(x = 1:10, y = 5:15, gr = 1:2) g$z <- log((g$x

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Dietrich Trenkler
Paul Smith schrieb: > Dear All, > > I know how to export graphics as pdf files and then how to include > them in LaTeX documents. However, I do not know how to do in order to > have the text of the graphics written with the font selected for the > LaTeX document. Is that possible? > > Thanks in adv

Re: [R] Password-protect script files

2007-09-28 Thread S Ellison
> Is there any way to password-protect script files (either within R or > otherwise)? The question seemed to me to be more about password protection against modification, rather than encryption. I'd have thought this was something a decent OS could take care of. It seems as daft to try to get a

[R] Barnard's exact test

2007-09-28 Thread João Fadista
Dear all, I would like to know if there is any R package that has the Barnard´s exact test. I also like to have any opinions about which test is more powerful: Fisher´s exact test or Barnard´s exact test. Best regards João Fadista Ph.d. student UNIVERSITY OF AARHUS F

Re: [R] ELF file OS ABI invalid yes?????

2007-09-28 Thread Rainer Hurling
Hello Vittorio, I am running FreeBSD 7.0-CURRENT (i386) from yesterday and I am able to compile MCMCpack_0.9-1.tar.gz without any problem on R-2.6.0 RC (2007-09-27 r43002). Native FreeBSD devel/pcre-7.3 and lang/gcc42 packages are installed. Please look if your library versions are ok. Regard

Re: [R] errorbar il lattice plot

2007-09-28 Thread deepayan . sarkar
On 9/28/07, Federico Calboli <[EMAIL PROTECTED]> wrote: > Hi Everyone, > > I would like to add errorbars to a lattice plot. I already have a function > that > adds error bars to a plot (originally written by George Gilchrist) which > takes > as arguments the coordinates of the point and then the si

Re: [R] RMySQL installing problems

2007-09-28 Thread Prof Brian Ripley
You do need the RMySQL client DLLs in your PATH. Also make sure that your MySQL version is not one listed as problematic at http://www.stats.ox.ac.uk/pub/RWin/ReadMe On Fri, 28 Sep 2007, andreiabb wrote: > > Dear forum, > > I'm trying to install RMySQL but I'm having problems in loading it, he

Re: [R] Plots with discontinuity balls

2007-09-28 Thread Paul Smith
On 9/28/07, Duncan Murdoch <[EMAIL PROTECTED]> wrote: > >>> Can R plot graphs like the one at > >>> > >>> http://www.mathwords.com/f/f_assets/floor_graph.gif > >>> > >>> with the balls at the discontinuity points? > >>> > >> You can use segments() to draw the segments and symbols() to draw the > >>

Re: [R] orientlib

2007-09-28 Thread Duncan Murdoch
On 9/28/2007 11:20 AM, Abu Naser wrote: > Hi All user, > > I have been using R-2.5.1. dose orientlib support this version? I would like > to try. it. > I have been wondering how to install the library. As far as I know it is fine. You install it as you would install any other package, but that

Re: [R] Drawing functions on Cartesian coordinate systems

2007-09-28 Thread Greg Snow
I think he wants the axes crossing at 0,0 not on the outer edges like the default. You can put the axes in the plot (though it tends to distract rather than help in many cases) by: > axis(1, pos=0) > axis(2, pos=0) You will need to draw the arrowheads yourself. There are options (under ?par) fo

Re: [R] Creating nice looking lists: how?

2007-09-28 Thread Katharine Mullen
The source code for print.summary.lm and summary.lm is in the file /src/library/stats/R/lm.R of the R source code, which you can download from CRAN if you don't have it already. The file lm.R is also at https://svn.r-project.org/R/trunk/src/library/stats/R/lm.R On Fri, 28 Sep 2007, Sergey Goriatc

[R] fitted values in LMER for the fixed-effects only

2007-09-28 Thread Anouk Simard
Hi, I would like to extract the fitted values from a model using LMER but only for the fix portion of the model and not for the fix and random portion (e.g it is the procedure outpm or outp in SAS). I am aware of the procedure fitted() but I not sure it give the fitted values both for the fixe

Re: [R] Drawing functions on Cartesian coordinate systems

2007-09-28 Thread Paul Smith
On 9/28/07, Greg Snow <[EMAIL PROTECTED]> wrote: > I think he wants the axes crossing at 0,0 not on the outer edges like > the default. > > You can put the axes in the plot (though it tends to distract rather > than help in many cases) by: > > > axis(1, pos=0) > > axis(2, pos=0) > > You will need t

Re: [R] File selection by condition

2007-09-28 Thread jim holtman
Read the files in, select those that meet the criteria and store them in a list match.list <- list() for (i in list.files(pattern="some pattern")){ x.in <- read.table(i) if (your match criteria) match.list[[i]] <- x.in } On 9/28/07, Svempa <[EMAIL PROTECTED]> wrote: > > I have a large num

[R] P-P plot

2007-09-28 Thread Maura E Monville
Sorry for my silly questions. I'm a beginner with R and most statistics concepts. I carried out a simple linear regression where the dependent variable is explained through a combination of powers of cos(independent variable). I can see R returned a good R^2 factor (> 0.99) but I have a hard time a

[R] Découvrez Bordeaux inscrit au patrimoine de l'UNESCO !

2007-09-28 Thread redaction
[1]voyages et vacances Destination vacances voyages et vacances [2]Merveilles de l'Egypte voyages et vacances [3]La Sicile mystérieuse voyages et vacances [4]Athènes et ses joyaux vacances et séjour chateau de bonagui

[R] boot.ci and NA values

2007-09-28 Thread Lucia Costanzo
Hi, I have been trying to run the code below. In the event of non-convergence, the statistic in the boot function returns NA as a value. To obtain a BCA confidence interval, I use boot.ci but, if NA recorded as a value for the statistic for one the replicates the following error appears: Error

Re: [R] Creating nice looking lists: how?

2007-09-28 Thread Greg Snow
Others have shown you where to find the actual printing code for print.summary.lm, but the short answer is: Use the cat function for the general text (the word "coefficients:" and the signif codes at the bottom) And put the actual coefficients in a matrix, use colnames and rownames to add column a

Re: [R] Drawing functions on Cartesian coordinate systems

2007-09-28 Thread Greg Snow
Look at the las argument in ?par for the easiest solution (can be passed to axis). -- 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

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Henrik Bengtsson
On 9/28/07, hadley wickham <[EMAIL PROTECTED]> wrote: > > Yes there is harm. But to make bold lines, easy to read titles is fine. > > See the spar function in > > http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see > > the setps, ps.slide, and setpdf functions in the Hmisc p

Re: [R] Drawing functions on Cartesian coordinate systems

2007-09-28 Thread Paul Smith
On 9/28/07, Greg Snow <[EMAIL PROTECTED]> wrote: > Look at the las argument in ?par for the easiest solution (can be passed > to axis). Thanks a lot, Greg! That is it! Paul > > -Original Message- > > From: [EMAIL PROTECTED] > > [mailto:[EMAIL PROTECTED] On Behalf Of Paul Smith > > Sent:

Re: [R] P-P plot

2007-09-28 Thread Greg Snow
Does this give you what you want? fit <- lm( Petal.Width ~ Petal.Length, data=iris) tmp1 <- resid(fit) tmp2 <- pnorm( tmp1, 0, summary(fit)$sigma ) par(mfrow=c(2,1)) qqnorm(tmp1) qqline(tmp1) plot( ppoints(length(tmp1)), sort(tmp2), xlab='Theoretical Percentiles', ylab='Sample Percentiles') ab

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Austin, Matt
I find them helpful for displaying the differences in coupled biomarkers between therapies with different therapeutic targets. You can quickly see/show which marker is affected quicker and then the other compensating in response and how they can be reversed depending on target. Having each measure

Re: [R] fitted values in LMER for the fixed-effects only

2007-09-28 Thread Douglas Bates
On 9/28/07, Anouk Simard <[EMAIL PROTECTED]> wrote: > Hi, > > I would like to extract the fitted values from a model using LMER but > only for the fix portion of the model and not for the fix and random > portion (e.g it is the procedure outpm or outp in SAS). I am aware of > the procedure fitted()

[R] plot x-axis at 0

2007-09-28 Thread James
Hi, I am trying to create a plot of a simple P&L graph for an option I'm pricing using Rmetrics. I'm not trying to do anything fancy. I just want to figure out if it's possible to put the x-axis at 0 rather than at the bottom the the entire plot. It seems like something that would be re

Re: [R] plot x-axis at 0

2007-09-28 Thread Greg Snow
Use xaxt='n' in the original plot to suppress the default axis (or axes=FALSE to supress both), then axis(1, pos=0) to draw the axis at 0. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Messa

Re: [R] plot x-axis at 0

2007-09-28 Thread James
On Sep 28, 2007, at 3:22 PM, Greg Snow wrote: > Use xaxt='n' in the original plot to suppress the default axis (or > axes=FALSE to supress both), then axis(1, pos=0) to draw the axis > at 0. That was what I was looking for. Thank you. I just need to do a bit of refining now. James ___

[R] Selecting values

2007-09-28 Thread Brian Perron
Hello all, An elementary question that I am sure can be easily cracked by an R enthusiast. Let's say I have multiple scores (y) on subjects (x.sample). Some subjects have a few more scores than others. Can somebody suggest some code that will select the first score for each subject? For example

Re: [R] Selecting values

2007-09-28 Thread Marc Schwartz
On Fri, 2007-09-28 at 17:48 -0400, Brian Perron wrote: > Hello all, > > An elementary question that I am sure can be easily cracked by an R > enthusiast. Let's say I have multiple scores (y) on subjects (x.sample). > Some subjects have a few more scores than others. Can somebody suggest some > c

Re: [R] Selecting values

2007-09-28 Thread Matthew Keller
Is this easier? x.index <- duplicated(x.sample)==FALSE cbind(x.sample[x.index],y[x.index]) - Matt On 9/28/07, Marc Schwartz <[EMAIL PROTECTED]> wrote: > On Fri, 2007-09-28 at 17:48 -0400, Brian Perron wrote: > > Hello all, > > > > An elementary question that I am sure can be easily cracked by a

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Frank E Harrell Jr
hadley wickham wrote: >> Yes there is harm. But to make bold lines, easy to read titles is fine. >> See the spar function in >> http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see >> the setps, ps.slide, and setpdf functions in the Hmisc package. > > I was interested to see

Re: [R] Selecting values

2007-09-28 Thread Marc Schwartz
Here is yet another approach using aggregate(), which internally, basically does what my first solution did: > aggregate(z[, 2], list(z[, 1]), "[", 1) Group.1 x 1 1 -1.2006469 2 2 -0.1614918 3 3 -0.5717729 4 4 -0.2398887 5 5 1.1690564 See ?aggregate Note

Re: [R] Graphics and LaTeX documents with the same font

2007-09-28 Thread Daniel Lakeland
On Fri, Sep 28, 2007 at 12:13:46PM -0500, hadley wickham wrote: > > Yes there is harm. But to make bold lines, easy to read titles is fine. > > See the spar function in > > http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see > > the setps, ps.slide, and setpdf functions in t

[R] Importing only one function from a package

2007-09-28 Thread Søren Højsgaard
Dear List In a package I want to import the mApply function from the Hmisc package, and I would like to import only that function. 1) If I write "Depends: Hmisc" in the DESCRIPTION file I get the whole Hmisc package, so that is not the way to go ahead. 2) According to "Writing R extensions"

Re: [R] Importing only one function from a package

2007-09-28 Thread Matthew Keller
Hi Soren, What I do in cases like this is just copy the function and place it into my script at the top (or write it into its own source file and call it from script). Best, Matt On 9/28/07, Søren Højsgaard <[EMAIL PROTECTED]> wrote: > Dear List > > In a package I want to import the mApply funct

Re: [R] Importing only one function from a package

2007-09-28 Thread Deepayan Sarkar
On 9/28/07, Søren Højsgaard <[EMAIL PROTECTED]> wrote: > Dear List > > In a package I want to import the mApply function from the Hmisc package, and > I would like to import only that function. > > 1) If I write "Depends: Hmisc" in the DESCRIPTION file I get the whole Hmisc > package, so that is

[R] RODBC and Oracle

2007-09-28 Thread John McHenry
Hi WizaRds, I'm experiencing a problem connecting to an Oracle 10g database via RODBC (I'm getting this on Microsoft XP). The same SQL queries via PL/SQL Developer work just fine, but when I pump the query through sqlQuery in RODBC then I get a data frame back with 0 rows. I cut the query down

[R] Problem compiling R - "use of NULL environment is defunct"

2007-09-28 Thread Saptarshi Guha
Hello, I downloaded R-2.5.1 and configured it on a AMD x86-74 running Redhat. During compilation i get this error configure --prefix=$HOME/mine make mkdir -p -- ../../../../library/methods/libs make[5]: Leaving directory `/home/sguha/tmp/R-2.5.1/src/library/ methods/src' make[4]: Leaving direc

Re: [R] New R website: R-Cookbook.com

2007-09-28 Thread Jeff Spies
Hadley, This was something I forgot to change before making the site public. Thanks for pointing out the creative commons license--it's definitely what I was thinking: share-alike content. I'd be glad to hear any other suggestions you might have, Jeff. On Sep 28, 2007, at 9:15 AM, hadley w

[R] seq() question

2007-09-28 Thread Derek Ogle
Hello all, I am attempting to use the seq() function to define breaks for the cut() function but am finding performance that I do not understand. The root of my problem appears to be contained in the following simple example. > brks <- seq(0.8,1.2,0.1) > brks [1] 0.8 0.9 1.0 1.1 1.2 > brks==1.

Re: [R] seq() question

2007-09-28 Thread Marc Schwartz
On Fri, 2007-09-28 at 21:58 -0500, Derek Ogle wrote: > Hello all, > > I am attempting to use the seq() function to define breaks for the > cut() function but am finding performance that I do not understand. > The root of my problem appears to be contained in the following simple > example. > > >

  1   2   >