Re: [R] Path to R script

Johannes W. Dietrich wrote: > > At 19:29 Uhr -0400 12.05.2010, Gabor Grothendieck wrote: >>If you are running the file via source("myfile.R") then just put this >>in the file: >> >>this.dir <- dirname(sys.frame(1)$ofile) >> >>This is a bit fragile since changes to the internals of source could >

Re: [R] read table for Fisher Exact

Thanks, David! - Original Message > From: David Winsemius > To: Frank E Harrell Jr > Cc: "Shi, Tao" ; r-help@r-project.org > Sent: Thu, May 13, 2010 8:29:31 PM > Subject: Re: [R] read table for Fisher Exact > > On May 13, 2010, at 9:24 PM, Frank E Harrell Jr wrote: > On > 05/13/2

Re: [R] Extract model matrix from plm

There is actually a model.matrix method for plm objects. Unfortunately, it is not currently exported and therefore can't be used by the end user of plm. I'll fix that in the next update. However, you won't find in this matrix the dummies because demeaned values are used. If you really want the matr

Re: [R] How to draw a graphic using data with coordinates and production rate?

Hi, assuming your production figures are in a square matrix, where points with no production take zero and points with production take the production figure, you could use image() This is the most basic approach I would think. Daniel -- View this message in context: http://r.789695.n4.nabbl

Re: [R] More complex historgram

Are you asking us to do your homework? This is not a homework list. For the t-test look in any introductory R manual. For the histograms, look in the lattice library. HTH Daniel -- View this message in context: http://r.789695.n4.nabble.com/More-complex-historgram-tp2197823p2216076.html Sent f

Re: [R] Split data frame by conditional and column at the same time

Does anything speak against selecting only x where x$B<20 and then splitting it by x$C? split(x[!x$B<20,],x$C) HTH, Daniel -- View this message in context: http://r.789695.n4.nabble.com/Split-data-frame-by-conditional-and-column-at-the-same-time-tp2197908p2216072.html Sent from the R help ma

Re: [R] Add 1 week to date with hours included for subset

On May 13, 2010, at 11:25 PM, emorway wrote: Forum, I've got hourly data for roughly 5 months of time that I want to subset into successive weeks. The first two elements of the dataset are well.80.2$date[1] [1] 6/12/2000 2:00 9256 Levels: 10/1/2001... well.80.2$date[2] [1] 6/12/2000

Re: [R] read table for Fisher Exact

On May 13, 2010, at 9:24 PM, Frank E Harrell Jr wrote: On 05/13/2010 08:16 PM, Shi, Tao wrote: Hi Prof. Harrell, Could you please elaborate on why chi-square test is more appropriate in this case? Thank you very much! ...Tao Exact tests tend to not be very accurate. Typically their P-

[R] Add 1 week to date with hours included for subset

Forum, I've got hourly data for roughly 5 months of time that I want to subset into successive weeks. The first two elements of the dataset are well.80.2$date[1] >[1] 6/12/2000 2:00 >9256 Levels: 10/1/2001... well.80.2$date[2] >[1] 6/12/2000 3:00 >9256 Levels: 10/1/2001 and so on until m

[R] Text Editing - UltraEdit and R

Hi, I use the text editor UltraEdit. Does anyone have any idea how to (a) run bloc lines of R code (rather than the entire file) from within UE and (b) execute graphics that will open up inside of R. I can run entire files of R code from within UltraEdit just fine via a power tool. But I would

Re: [R] Multiple plots; single x(y) labels

Try this: x1 <- rnorm(7) x2 <- rnorm(27) x3 <- rnorm(47) x4 <- rnorm(67) xyplot(value ~ i | id, do.call(rbind, lapply(ls(patt = "x"), function(x)data.frame(id = x, value = get(x), i = seq_along(get(x), type='l') On Thu, May 13, 2010 at 1:04 PM, Xin Ge wrote

Re: [R] R cannot access the web?

Check if your network has a proxy. see ?setInternet2 On Thu, May 13, 2010 at 7:29 PM, Trojan wrote: > > Hi, > > I had originally posted regarding an error when trying to install package - > GenABEL - it has now become clear that R is not able to connect to the net. > Below are a couple of things

Re: [R] read table for Fisher Exact

On 05/13/2010 08:16 PM, Shi, Tao wrote: Hi Prof. Harrell, Could you please elaborate on why chi-square test is more appropriate in this case? Thank you very much! ...Tao Exact tests tend to not be very accurate. Typically their P-values are too large. See @Article{cra08how, author =

Re: [R] read table for Fisher Exact

Hi Prof. Harrell, Could you please elaborate on why chi-square test is more appropriate in this case? Thank you very much! ...Tao - Original Message > From: Frank E Harrell Jr > To: r-help@r-project.org > Sent: Thu, May 13, 2010 5:35:07 AM > Subject: Re: [R] read table for Fisher

Re: [R] R cannot access the web?

I was also able to download the NYTimes page, using R2.10.1 on Windows XP. -- View this message in context: http://r.789695.n4.nabble.com/R-cannot-access-the-web-tp2215881p2216005.html Sent from the R help mailing list archive at Nabble.com. __ R-he

Re: [R] Multiple plots; single x(y) labels

Not sure why you want to do that b/c this can be done with one function call of "xyplot" in lattice. Anyway, you can always set axes=F and then replot each axis using 'axis' for each plot. Make sure before you do that, set your xlim and ylim, so all four plots have the same xlim and ylim. Lik

Re: [R] How to rank matrix data by deciles?

Vincent - I'm afraid there's no solution other than artificially modifying the zeroes: vec [1] 26.58950617 5.73074074 5.9622 5.6478 20.95728395 0. 0.0700 12.8689 [9] 3.64543210 0.05049383 25.6089 3.53246914 0. 31.39049383 3.77641975 13.196172

Re: [R] How to rank matrix data by deciles?

Dear Phil, You helped me with a request to rand matrix columns by deciles two weeks ago. This really un-blocked me on this project but I found a little bug. As in before, my data is in a matrix: > madebt[1:16,1:2] X4.19.2010 X4.16.2010 [1,] 26.61197531 26.58950617 [2,] 5.72765432

Re: [R] How to rank matrix data by deciles?

Dear Phil, You helped me with a request to rand matrix columns by deciles two weeks ago. This really unblocked me on this project but I found a little bug. As in before, my data is in a matrix: > madebt[1:16,1:2] X4.19.2010 X4.16.2010 [1,] 26.61197531 26.58950617 [2,] 5.72765432 5.7

Re: [R] R cannot access the web?

On May 13, 2010, at 6:29 PM, Trojan wrote: Hi, I had originally posted regarding an error when trying to install package - GenABEL - it has now become clear that R is not able to connect to the net. Below are a couple of things I've tried with the resulting errors I am running R 2.6.2 o

[R] R cannot access the web?

Hi, I had originally posted regarding an error when trying to install package - GenABEL - it has now become clear that R is not able to connect to the net. Below are a couple of things I've tried with the resulting errors I am running R 2.6.2 on Linux e15 X86 How can I fix this? Thank you! 1

[R] S-plus Finmetrics "cpredict" command

Hi, I am currently transferring a portion of S+ code into R. Do you guys know if there is an alternative to the "cpredict" command, which does hard/soft conditional forecasting on VARs. I've looked at the vars package, but it does not inclue this type of function. Thank you.

[R] Where has the stats-rosuda-devel mailing list gone?

ESET NOD32 Antivirus, version of virus signature database 5113 (20100513) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Simple question on binning data

I second it! Made me lol :) Dennis On Thu, May 13, 2010 at 3:03 PM, Rolf Turner wrote: > > On 14/05/2010, at 9:54 AM, Frank E Harrell Jr wrote: > > > > > Binning is seldom needed and usually distorts. It is > > the statistical equivalent of a former governor from Alaska. > > >

[R] ARMA(1,1)-GARCH(1,1) rolling estimation question

Hi all, I got the daily stock return data from 2005 - 2008, calculated from HF minute data. (Thanks to Jeff and Josh). Now, I set 05 - 07yr as the carlibration period for estimating the parameters of ARMA(1,1)-GARCH(1,1) model, aqnd leave 08 for backtesting. So I use the return data observations

Re: [R] Simple question on binning data

On 14/05/2010, at 9:54 AM, Frank E Harrell Jr wrote: > Binning is seldom needed and usually distorts. It is > the statistical equivalent of a former governor from Alaska. I nominate this for a fortune. cheers, Rolf Turner

Re: [R] Simple question on binning data

On 05/13/2010 04:05 PM, Bert Gunter wrote: ?hist However, I'm going to go way out on a limb here, as I know absolutely nothing of what you're up to, and suggest that you do NOT need to bin your data. -- Bert Bert is correct. Binning is seldom needed and usually distorts. It is the statisti

Re: [R] a question about "latex" in Hmisc

Tao Shi, Here is the repair for your query. Use this version of show.dvi temporarily (until Hmisc includes it in its official version of show.dvi). This function and option setting will let the command latex(x) display the dvi file in a YAP window on Windows. This version is not to be used

Re: [R] Simple question on binning data

It's very simple to write a "binit()" function. If all you want to do is e.g., bin 107 values into sums of 10 at a time, then write a loop that sums x[10*i:11*i-1] (not tested and not syntactically correct). The one I wrote for myself discards any partial bin (101-107 in my example) and leav

Re: [R] aggregate.zoo

Here is a second solution where dly and year are as in the original post. The function returns the minimum index and applying that to the sequence numbers of the days in each year we get a vector of the indexes of the yearly minimums and then can use that to subscript dly: f <- function(ix) ix[wh

Re: [R] tune svm

Hi, On Thu, May 13, 2010 at 3:54 PM, jparga wrote: > > Hello, I hope you can help me! > I`m trying to tune svm parameters: cost and gamma for a landsat image > classification, but I get an error and I can't understand it. > I write this: >> tune(svm, Class~., data = mdt01bis, ranges = list(gamma

Re: [R] Simple question on binning data

?hist However, I'm going to go way out on a limb here, as I know absolutely nothing of what you're up to, and suggest that you do NOT need to bin your data. -- Bert Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-b

[R] GAM, GAMM and numerical integration, help please

I am trying to apply methods used by Chaloupka & Limpus (1997) ( http://www.int-res.com/articles/meps/146/m146p001.pdf) to my own turtle growth data. I am having trouble with two things... 1) After the GAM is fit, the residuals are skewed. >m1 <- gam(growth~s(mean.size, bs="cr")+s(year,bs="cr",k

Re: [R] frequencies of time series data

Hello, Adam Sean wrote: I have an excel file in following shape: date Category 01/21/2010turf 01/21/2010jack 01/21/2010jack 01/22/2010 turf 01/22/2010 jack 01/22/2010 psi 01/22/2010 psi 01/22/2010 psi 01/23/2010 turf 01/23/2010 turf 01/23/2010 jack 01/23/

Re: [R] Adding name to nth row

We need a self-contained, reproducible example of what you have, and what you want. Is this close? tmp <- data.frame(a = rnorm(278)) tmp$newcol <- "" tmp[seq(97,278, by = 12), "newcol"] <- "A Name" ecvet...@uwaterloo.ca wrote: Hello, I have a data frame with many rows, and I want to create

Re: [R] Simple question on binning data

There would be several people who could help if you gave us a minimal, reproducible example like the posting guide asks for. If you have a vector of continuous data, and need to create a categorical variable (in R, a factor) from that continuous variable, then ?cut can help you. george5000

[R] Adding name to nth row

Hello, I have a data frame with many rows, and I want to create a column with a name only at every 12th row, starting from 97 to 278. Thanks in advance! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

[R] Simple question on binning data

Hello everyone, I have a data set, and I need to bin my data using a bin width of say g(n). Would anyone be willing to tell me how to do this in R? Thanks -- View this message in context: http://r.789695.n4.nabble.com/Simple-question-on-binning-data-tp2202644p2202644.html Sent from the R help

[R] frequencies of time series data

I have an excel file in following shape: date Category 01/21/2010turf 01/21/2010jack 01/21/2010jack 01/22/2010 turf 01/22/2010 jack 01/22/2010 psi 01/22/2010 psi 01/22/2010 psi 01/23/2010 turf 01/23/2010 turf 01/23/2010 jack 01/23/2010 psi I need to outp

[R] tune svm

Hello, I hope you can help me! I`m trying to tune svm parameters: cost and gamma for a landsat image classification, but I get an error and I can't understand it. I write this: > tune(svm, Class~., data = mdt01bis, ranges = list(gamma = 2^(-15:3), cost > = 2^(-5:15))) and R gives: Error en predic

Re: [R] aggregate.zoo

Achim, Thanks for pointing me in the right direction! I broke your user-defined function down into two separate aggregate calls on which.min and length. library(zoo) n <- 1000 dly <- zoo( x = rnorm(n), order.by = seq.Date(as.Date('1970-01-01'), by = 'day', length.out = n) ) year <- as.numer

Re: [R] a question about "latex" in Hmisc and .dvi file

I didn't see that in the Adobe help file. There was some reference in the Adobe online forum http://forums.adobe.com/ but the first few hits were queries, not responses to the query. If you pursue it, please report success back to the list. Rich On Thu, May 13, 2010 at 4:08 PM, Felipe Carrill

Re: [R] a question about "latex" in Hmisc and .dvi file

I was used to see the dvi file with menus and options to rotate  the document and was expecting the pdf file  to landscape and be rotated automatically. Is there anywhere in your function that I can have the pdf to be rotated clockwise automatically? > >From: RICHARD M. HEIBERGER >To: Felipe

Re: [R] Path to R script

In this case, deploy them as a R package is the way to go. - Original Message > From: Johannes W. Dietrich > To: "Shi, Tao" ; johannes.dietr...@rub.de > Cc: r-help@r-project.org > Sent: Thu, May 13, 2010 9:14:24 AM > Subject: Re: [R] Path to R script > > At 13:19 Uhr -0700 12.05.2010

Re: [R] need help in igraph package of R

Hi, See comments in line. On Thu, May 13, 2010 at 1:07 PM, bababasi wrote: > > hi > > I am struck with a problem in igraph package of R. My problem is as follows > > I want to plot a power law fit for my data (in .net format --- pajek format) > > syntax for that in R is > > g <- read.graph("file

Re: [R] a question about "latex" in Hmisc and .dvi file

Your example showed up in landscape for me on the pdf file. I needed to rotate the image in the Adobe Reader. I like the appearance better with show.dvi(dvi(x.tex <- latex(report,landscape=TRUE), width=8.5, height=11)) where I changed the paper size in the dvi() command. Rich [[altern

[R] ESP Ghostscrip ...

Hello All: Could anybody tell me please what is the meaning of this error message: ESP Ghostscript 815.02: Unrecoverable error, exit code 1 Error in plot.xy(xy, type, ...) : ignoring SIGPIPE signal Calls: plot -> plot.default -> plot.xy Execution halted My code has to go to 18 loops for three ti

Re: [R] a question about "latex" in Hmisc and .dvi file

That function works nicely generating the pdf on the fly instead of the dvi file, however, consider the example below. I can't get the pdf to landscape. report <- structure(list(Date = c("2/26/2010", "2/27/2010", "2/28/2010", "3/1/2010", "3/2/2010", "3/3/2010", "3/4/2010", "3/5/2010", "3/6/2010"

Re: [R] aggregate.zoo

On Thu, 13 May 2010, Hutchinson,David [PYR] wrote: I have been using aggregate function to determine the annual minimum in a daily time series. However, I was also hoping to determine the date of occurrence without having to loop through each year to identify the date of each minimum returned fr

Re: [R] Comparing histograms?

On 13-May-10 16:00:51, Jonathan Greenberg wrote: > Rhelpers: > I'm curious what the appropriate analysis to use for testing the > hypothesis that two histograms are statistically different from one > another? Thanks! > > --j That's potentially several questions in one, and each question might ha

Re: [R] a question about "latex" in Hmisc

Hi Prof. Harrell, Originally, I was thinking that by putting tables (especially those long tables) into separate .tex files would make my 'big' .tex file more organized (i.e. rather than having an one giant .tex file). But now I'm having a second thought, b/c I don't really have to worry how m

Re: [R] a question about "latex" in Hmisc

Hi Duncan, I have no problems running latex or pdflatex. I re-ran my code and it seems the problem is that R (or MiKTex ?) put the .tex file (the one wrapping everything in a document format for displaying purpose) in a R tmp folder, but the corresponding .log, .aux, and .dvi in my working dir

Re: [R] How to display data values for points in a plot?

Nish wrote: Hello, I would like to know how to display values for points in a plot funtion. For example, plot( y=dat$a, x=dat$b, main="plot1", ylab="a", xlab="b", ylim=c(-10, 10), xlim=c(-10, 10), type = "p", pch=17, col=, c

Re: [R] Access to files downloaded

Joey Zhou wrote: Dear All: I downloaded "gee" source code and now I am having problem accessing the file. download.packages("gee", type="source", destdir="~") --- Please select a CRAN mirror for use in this session --- trying URL ' http://www.ibiblio.org/pub/languages/R/CRAN/src/contrib/g

[R] install.packages Error

Hi all! I am trying to install package - GenABEL in a Linux (RedHat) environment using R 2.6.2. The CRAN Mirror starts up, I select the site, but I end up with the following error: ** Warning: unable to access index for repository

[R] need help in igraph package of R

hi I am struck with a problem in igraph package of R. My problem is as follows I want to plot a power law fit for my data (in .net format --- pajek format) syntax for that in R is g <- read.graph("filename.net", "pajek") d <- degree (g, mode="in") power.law.fit (d+1, 2) it gives me desired o

[R] need help in igraph package of R

hi I am struck with a problem in igraph package of R. My problem is as follows I want to plot a power law fit for my data (in .net format --- pajek format) syntax for that in R is g <- read.graph("filename.net", "pajek") d <- degree (g, mode="in") power.law.fit (d+1, 2) it gives me desired o

[R] How to display data values for points in a plot?

Hello, I would like to know how to display values for points in a plot funtion. For example, plot( y=dat$a, x=dat$b, main="plot1", ylab="a", xlab="b", ylim=c(-10, 10), xlim=c(-10, 10), type = "p", pch=17, col=, cex=1.5

Re: [R] Creating data frames in a loop

Thank yoy. this works. -- View this message in context: http://r.789695.n4.nabble.com/Creating-data-frames-in-a-loop-tp2195527p2197780.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz

[R] Access to files downloaded

Dear All: I downloaded "gee" source code and now I am having problem accessing the file. > download.packages("gee", type="source", destdir="~") --- Please select a CRAN mirror for use in this session --- trying URL ' http://www.ibiblio.org/pub/languages/R/CRAN/src/contrib/gee_4.13-14.tar.gz' Co

[R] Comparing histograms?

Rhelpers: I'm curious what the appropriate analysis to use for testing the hypothesis that two histograms are statistically different from one another? Thanks! --j __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

[R] More complex historgram

CID DISTRICTSCHOOL YRofSTUDY DATE.ENTRY SEX AGE EGGS.M EGGS.MIIEGGS.H POBADN0130 PO Badongo 0 18.1.08 2 11 0 -1 73 POBADN0117 PO Badongo 0 18.1.08 1 10 0 -1 59 POBADN0160 PO

Re: [R] Whiskers on the default boxplot {graphics}

Hi Robert, Your points are well taken. However, I reserve mine, b/c I think without this detailed discussion, an average R user would simply confused the "interquartile range" said in boxplot help file with the results of "IQR". Changing it to "length of box" makes it more exact and consisten

[R] aggregate.zoo

I have been using aggregate function to determine the annual minimum in a daily time series. However, I was also hoping to determine the date of occurrence without having to loop through each year to identify the date of each minimum returned from aggregate(). Is there a more convenient way of doin

Re: [R] Multiple plots; single x(y) labels

On May 13, 2010, at 1:24 PM, Xin Ge wrote: > Hi David (and Others), > > I think I messed it up. Lets start afresh, I do not want to use > lattice for this. I'm using multiple "plot" commands and then > eventually would like to get a *combine* x-label and y-label for > this plot. > > Like fo

Re: [R] Multiple plots; single x(y) labels

Hello Xin, If you need to recreate plots from different systems you need full control over the graphics output. This is not always easy with traditional graphics, better have a look at the plotting functions from lattice (based on grid graphics system in R) or maybe ggplot2. http://www.googl

Re: [R] Whiskers on the default boxplot {graphics}

I agree. I was convinced by Ehlers' example that type =2 was a better match to fivenum's result -- David.. On May 13, 2010, at 1:36 PM, Joshua Wiley wrote: On Thu, May 13, 2010 at 7:55 AM, David Winsemius > wrote: Yes, and experimentation leads me to the conclusion that the only possible c

Re: [R] Whiskers on the default boxplot {graphics}

On Thu, May 13, 2010 at 7:55 AM, David Winsemius wrote: > Yes, and experimentation leads me to the conclusion that the only possible > candidate for matching up the results of fivenum[c(2,4] with quantile(y, > c(1,3)/4, type=i) is for type=5. I'm not able to prove that to myself from > mathematic

Re: [R] Multiple plots; single x(y) labels

Hi David (and Others), I think I messed it up. Lets start afresh, I do not want to use lattice for this. I'm using multiple "plot" commands and then eventually would like to get a *combine* x-label and y-label for this plot. Like for example, the following plot has a combine x-label ("Height") an

Re: [R] a question about "latex" in Hmisc and .dvi file

> > At least I thought it did, since I believe show.dvi is an exported > function, can't you just overwrite it in your global environment and > anything that calls it will use your copy? I'll have to look into how I > got your code working on my system, because I know I didn't use that trick! >

Re: [R] Help with reading information of "summary"-Object

Tino - When you call summary for an lm object, the function that actually gets called is "summary.lm". This function has a help page, and, if you type ?summary.lm and look at the "Value" section, you'll see the following: coefficients: a p x 4 matrix with columns for the estimated

Re: [R] a question about "latex" in Hmisc and .dvi file

Notice that I didn't quite get my function into the Hmisc environment and that is why you need the awkward phrase show.dvi(dvi(x.tex <- latex(x))) instead of the smooth phrase latex(x) Somehow with my own data and not using the awkward bit, it worked on my Linux machine? At least I

Re: [R] Multiple plots; single x(y) labels

On May 13, 2010, at 12:59 PM, Xin Ge wrote: Hi David, Thanks for your reply. By single x and y-labels I meant something like this: http://zoonek.free.fr/blosxom//R/2006-08-10_lattice_xyplot_quakes.png which lattice gives by default. The code you sent doesn't seem to solve the problem, I

Re: [R] Help with reading information of "summary"-Object

Hello, Tino Schöllhorn wrote: Hi, I am quite new to R - but quite expierience in programming. Welcome to R! It's easiest if you give reproducible examples so we can help you. You can include code to create objects, or use the ?dput function. Nonetheles I have some problemes in access

[R] pnmath on windows

Hello All, I'd like to enable threaded computing on my system (windows 7 pro x64) if possible. pnmath seems to be the most accepted solution out there for the time being, but it requires compiling. Are there binaries for windows available that don't require manual compiliation? If not, does any

Re: [R] a question about "latex" in Hmisc and .dvi file

On Thu, May 13, 2010 at 12:46 PM, Erik Iverson wrote: > Looks good. Thank you. > I came up with something on my own which was basically to re-write > print.latex. I call pdflatex a few times because with the longtable > package, sometimes things don't line up right until you run pdflatex >

Re: [R] Multiple plots; single x(y) labels

Hi David, Thanks for your reply. By single x and y-labels I meant something like this: http://zoonek.free.fr/blosxom//R/2006-08-10_lattice_xyplot_quakes.png which lattice gives by default. The code you sent doesn't seem to solve the problem, I'm sorry if I havent' explained it clearly before. A

[R] Help with reading information of "summary"-Object

Hi, I am quite new to R - but quite expierience in programming. Nonetheless I have some problemes in accessing information of the "summary" object. Here is what I do: model <- lm ( y ~ myVariable ) "summary(model)" gives me an object which has a lot of information about the regression. Now

Re: [R] Multiple plots; single x(y) labels

On May 13, 2010, at 12:35 PM, David Winsemius wrote: ?plot # ylim and you need to have the data in a form (before plotting) where you can determine the shared max and min for the y limits On May 13, 2010, at 12:04 PM, Xin Ge wrote: Hi All, Can anyone please help me with getting a sing

Re: [R] a question about "latex" in Hmisc and .dvi file

Looks good. I came up with something on my own which was basically to re-write print.latex. I call pdflatex a few times because with the longtable package, sometimes things don't line up right until you run pdflatex multiple times. I think mine is less flexible though, and I like your soluti

Re: [R] Whiskers on the default boxplot {graphics}

On May 13, 2010, at 12:18 PM, Robert Baer wrote: And try this (which seems to leave us with type=2) and is listed in ? quantile as "Discontinuous sample quantile types 1, 2, and 3" quantile(1:101, c(1,3)/4, type=2) 25% 75% 26 76 I think Peter may be right,. If I do it with the rnorm functi

Re: [R] Multiple plots; single x(y) labels

?plot # ylim and you need to have the data in a form (before plotting) where you can determine the shared max and min for the y limits On May 13, 2010, at 12:04 PM, Xin Ge wrote: Hi All, Can anyone please help me with getting a single x and y-axis label while plotting muliple plots. H

Re: [R] Path to R script

At 19:29 Uhr -0400 12.05.2010, Gabor Grothendieck wrote: If you are running the file via source("myfile.R") then just put this in the file: this.dir <- dirname(sys.frame(1)$ofile) This is a bit fragile since changes to the internals of source could break it but it does currently work. Great,

Re: [R] Whiskers on the default boxplot {graphics}

And try this (which seems to leave us with type=2) and is listed in ?quantile as "Discontinuous sample quantile types 1, 2, and 3" quantile(1:101, c(1,3)/4, type=2) 25% 75% 26 76 David, try this: fivenum(1:101) quantile(1:101, c(1,3)/4, type=5) -Peter On 2010-05-13 8:55, David Winsemius

Re: [R] Help with interpolation

... and what package is interpNA in, pray tell? (Such information is requested in the posting guide). Base R does it with a bit of fiddling: ?approx But please read the Help carefully. Here is a simple example of how it is used: > x <- 1:5 > y <- c(2,rep(NA,3),7) > approx(x,y, xout =1:5)## Note

Re: [R] Path to R script

At 13:19 Uhr -0700 12.05.2010, Shi, Tao wrote: Don't quite understand your question, but it looks like a more IT issue to me. I guess you store your R scripts in a central location (e.g. a server) and everybody call them from their own workstation, right? Thank you, but this would be an easy

Re: [R] Whiskers on the default boxplot {graphics}

David, try this: fivenum(1:101) quantile(1:101, c(1,3)/4, type=5) -Peter On 2010-05-13 8:55, David Winsemius wrote: On May 13, 2010, at 10:25 AM, Robert Baer wrote: Hi Peter, You're absolutely correct! The description for 'range' in 'boxplot' help file is a little bit confusing by using

Re: [R] Help with interpolation

On Wed, May 12, 2010 at 5:14 PM, Jorge Nieves wrote: > I have a two Colum matrix ( I named it holder) as shown bellow. It seems > to be three columns, but the first column contains the row names. You can't have a matrix whose columns are different classes. > > I am trying to fill all  NAs by lin

[R] Multiple plots; single x(y) labels

Hi All, Can anyone please help me with getting a single x and y-axis label while plotting muliple plots. Here is the code: par(mfcol=c(2,2)) plot(x <- sort(rnorm(7)), type = "s", main = "", ylab="", xlab="") plot(x <- sort(rnorm(27)), type = "s", main = "", ylab="", xlab="") plot(x <- sort(rnorm(

Re: [R] a question about "latex" in Hmisc and .dvi file

Here is the full repair for the latex functions in Hmisc to make pdflatex work in Windows. This version is still slightly awkward. I hope that Charles and Frank will smooth it out and put it in their next release. I added two new options() and revised show.dvi so it will use them. Rich librar

[R] Help with interpolation

> Hi, > > I have a two Colum matrix ( I named it holder) as shown bellow. It > seems to be three columns, but the first column contains the row > names. > > I am trying to fill all NAs by linear interpolation the known values > in the matrix. I tried using the interpNA function, but I get an e

Re: [R] read table for Fisher Exact

How much would that differ? And how should I build my table if I want to do a Pearson Chi Sqaure? - Original Message - From: "Frank E Harrell Jr [via R]" Date: Thursday, May 13, 2010 2:36 pm Subject: Re: read table for Fisher Exact To: visser > On 05/12/2010 03:31 PM, visser wrote:

Re: [R] how to make monthly time series out of daily

Reto Baumgartner gmail.com> writes: > I am using the Rmetrics package and would like to convert a daily > price time serie into a monthly one. In SPlus I could use: > aggregateSeries(timeSerie, by="months",FUN=first). I wrote now my own function. But is there an easier way? #---

[R] Questions about ggplot2

Hi I have two questions about using ggplot2. First, I have multiple columns of data that I would like to combine into one histogram where each column of data would correspond to one bar in the histogram. Each column has 0 or 1s and I want my bars in the histogram to correspond to the sum of th

Re: [R] Compiling R with --enable-R-shlib for rpy2 error

On 13.05.2010 11:45, Padma TAN wrote: Hi, I am trying to compile R with the command below in order to install rpy2. ./configure --prefix=/usr/local/R/R-2.9.2 --with-gnu-ld --with-cairo --with-x --enable-R-shlib However, error the I have gotten was: /usr/bin/ld: CConverters.o: relocation R_X

Re: [R] Source code for gee

On 13.05.2010 16:29, Joey Zhou wrote: Dear all: I am looking for the source code for gee(generalized estimation equations by Liang and Zeger) package. Can you help? download.packages("gee", type="source", destdir="~") Uwe Ligges Thanks much, Joey [[alternative HTML version del

Re: [R] R version of SAS/IntrNet

Hi, more links in the FAQ: http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-Web-Interfaces On Thursday 13 May 2010 01:46:51 am Scott Czepiel wrote: > I recall hearing about -- or was I just dreaming? -- a way to compile > R as an apache module so it could stay in memory while the webserver > was

Re: [R] Help with Loops

On May 13, 2010, at 10:49 AM, Amit Patel wrote: Hi I have tried many attempts but cant get the loop right, as I am not a strong programmer. What I am basically trying to do is compare 2 spreadsheets. The problem is that one of them only contain a portion of the overall data (TESTSAMP), w

Re: [R] A primitive OO in R -- where next?

This is not really OO at all, in my opinion. It's an example of the amazing flexibility of the language. I'd like to add on to what Erik said, with an example: sum(1:10) [1] 55 foo <- sum foo(1:10) [1] 55 junk <- list(a=sum) junk$a(1:10) [1] 55 sum is an R object; it happens to be

Re: [R] How ls() only functions or anything else but functions?

On 5/13/2010 11:04 AM, Chuck Cleland wrote: > On 5/13/2010 10:02 AM, John Edwards wrote: >> Hello, >> >> How ls() only functions or only data objects (basically anything other than >> functions) such as data.frame, numeric ...? > > c(ls.str(mode = "function")) > > ls()[!(ls() %in% c(ls.str(mode="

  1   2   >