Re: [R] Data frame reordering to time series

2010-08-08 Thread steven mosher
In the real data the months are all complete, but the years can be missing. So years can be missing up front, in the middle, at the end. but if a year is present than every month has a value or NA. To create regular R ts I had to plow through the data frame, collect a year caluculate an index to

Re: [R] Sweave for inclusion of p value in a sentence of a LaTeX document

2010-08-08 Thread Matthieu Dubois
Dear Julia, my way to do that is to attribute the t.test to an object, and then refer to its p.value with the function \Sexpr e.g. \documentclass{article} \usepackage{Sweave} \begin{document} echo=FALSE= x-cbind(1,2,3) y-cbind(3,4,5) t - t.test(x,y) @

Re: [R] Error in as.environment(pos): using 'as.environment(NULL)' isdefunct

2010-08-08 Thread Immanuel Seeger
Hello, unfortunately I haven't gotten an answer up to now. You can see my question/problem below in the first mail. Is the given extract too unspecific? I haven't solved the problem but I think that it's in the handler.Berechnung function. Maybe I'm using wrong data-types (there are vectors

Re: [R] Data frame reordering to time series

2010-08-08 Thread Gabor Grothendieck
On Sun, Aug 8, 2010 at 2:01 AM, steven mosher mosherste...@gmail.com wrote: In the real data the months are all complete, but the years can be missing. So years can be missing up front, in the middle, at the end. but if a year is present than every month has a value or NA. To create regular R

Re: [R] metafor and meta-analysis at arm-level

2010-08-08 Thread Angelo Franchini
Dear Wolfgang, Is there any way for rma to add random effects only to each treatment arm, but not to the control one? Many thanks, Angelo On Thu, August 5, 2010 6:21 pm, Viechtbauer Wolfgang (STAT) wrote: Dear Angelo, rma(yi=o, sei=se, mods=~s+t-1, method=REML) is *a* way to run the

Re: [R] Latex errors when checking package

2010-08-08 Thread Christian Hennig
Dear Ben, thanks a lot! This did the trick (after adding install behind apt-get). Regards, Christian On Fri, 6 Aug 2010, Ben Bolker wrote: Christian Hennig chrish at stats.ucl.ac.uk writes: LaTeX errors when creating PDF version. This typically indicates Rd problems. LaTeX errors found: !

Re: [R] Error on random forest variable importance estimates

2010-08-08 Thread Pierre Dubath
Hello Andy, Thank you for your quick and helpful reply. I will try to follow your suggestions. Also, thank you for the R implementation of random forest. It is very useful for our work. Best, Pierre Liaw, Andy wrote: From: Pierre Dubath Hello, I am using the R randomForest package to

[R] Importing arguments for use by functions in a script

2010-08-08 Thread EvansA
QUESTION: Is there a way of passing arguments from an external file to a script so that they can be used directly by functions within the script? I have a series of interdependent functions. I wish to test the time for processing various datasets. I was initially doing something along the lines

[R] subset of string by index

2010-08-08 Thread david h shanabrook
How can I get a substring based on the index into the string? strM - c(abcde, cdefg) ind - c(1,3,5) I want to use ind to index into the strings so the result is: strMind - c(ace, ceg) __ R-help@r-project.org mailing list

Re: [R] plot the dependent variable against one of the predictors with other predictors as constant

2010-08-08 Thread Bert Gunter
I am thinking about using 'predict' command to generate the prediction of z with the new data.frame but there should be a better way. I'm puzzled. Why would you think that? This is the canonical way to do predictions in R. Bert Gunter Genentech Nonclinical Statistics

[R] subset of string by index

2010-08-08 Thread david h shanabrook
How can I get a substring based on the index into the string? strM - c(abcde, cdefg) ind - c(1,3,5) I want to use ind to index into the strings so the result is: strMind - c(ace, ceg) __ R-help@r-project.org mailing list

Re: [R] subset of string by index

2010-08-08 Thread Marc Schwartz
On Aug 8, 2010, at 9:50 AM, david h shanabrook wrote: How can I get a substring based on the index into the string? strM - c(abcde, cdefg) ind - c(1,3,5) I want to use ind to index into the strings so the result is: strMind - c(ace, ceg) Here is one way: apply(sapply(ind,

Re: [R] Data frame reordering to time series

2010-08-08 Thread steven mosher
Ok, I'm a bit confused by what you mean by regularly spaced After I do the do.call I do get a data structure with all the times present and every time has a NA or a data value. Steve On Sun, Aug 8, 2010 at 2:46 AM, Gabor Grothendieck ggrothendi...@gmail.comwrote: On Sun, Aug 8, 2010 at 2:01

[R] Integration under Tps

2010-08-08 Thread sam.e
Hi there, After many attempts at trying different methods I have used the Tps package to interpolate between a number of measurements of field data, in this case thickness of a sedimentary deposit. My input values are Longitude, Latitude and Thickness and I have produced a valid contour map

Re: [R] Data frame reordering to time series

2010-08-08 Thread Gabor Grothendieck
On Sun, Aug 8, 2010 at 11:21 AM, steven mosher mosherste...@gmail.com wrote: Ok, I'm a bit confused by what you mean by regularly spaced After I do the  do.call I do get a data structure with all the times present and every time has a NA or a data value. Steve regularly spaced means that

Re: [R] Data frame reordering to time series

2010-08-08 Thread Gabor Grothendieck
On Sun, Aug 8, 2010 at 11:55 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Sun, Aug 8, 2010 at 11:21 AM, steven mosher mosherste...@gmail.com wrote: Ok, I'm a bit confused by what you mean by regularly spaced After I do the  do.call I do get a data structure with all the times

Re: [R] package for measurement error models

2010-08-08 Thread David Winsemius
There is a MethComp package from Bendix Carstensen that does not seem to be on CRAN yet. http://staff.pubhealth.ku.dk/~bxc/MethComp/ But it is on R-Forge: https://r-forge.r-project.org/R/?group_id=364 -- David. On Aug 7, 2010, at 4:50 PM, Christos Argyropoulos wrote: I believe there was

[R] Bootstrap

2010-08-08 Thread marcellina79
I would like to bootstrap the difference between two mean costs (I have a different number of observations in the two group costs under consideration). So, this is my script: x - datos$direct_costsUS fun - function(x){m=(mean(x[63:136])-mean(x[1:62])) } boot.x- boot(x,

[R] Problems with Treemaps in VCD PORTFOLIO

2010-08-08 Thread David Epstein
Hi, I am trying to create a treemap for a non-symmetric tree where some branches have more sub-branches than others. I have tried tools in two packages, but have encountered some problems that I hope more experienced users on this list can help address. The PORTFOLIO package offers MAP.MARKET,

Re: [R] Bootstrap

2010-08-08 Thread Wu Gong
Hi, I assume you used boot package. Try this: library(boot) x - rnorm(136) fun - function(x,ind){ x - x[ind] m - mean(x[63:136])-mean(x[1:62]) m } boot.x- boot(x, fun, R = 5000, sim = ordinary) - A R learner. -- View this message in context:

[R] pdf device (and sweave) with localized characters

2010-08-08 Thread Roman Luštrik
Dear list, I'm trying to Sweave with localized characters in graph titles and axis labels through a pdf device. I set pdf encoding through pdf.options(encoding = CP1250). When I run my script through the command line like so: R CMD Sweave report.Rnw, the localized characters turn into two dots.

[R] Does anybody know how to control the appearance of the end of the line in lattice?

2010-08-08 Thread George Chen
Hi All, I am plotting vertical lines using xyplot in lattice and type=h. It works well, but the problem is that the tops of the lines are convex and the bottoms are concave. Is there a way to flatten the tops and bottoms? Here's my code: Source-matrix(1:30,10,3) colnames(Source)-c(x,y1,y2)

[R] p.adjust( , fdr)

2010-08-08 Thread qing chen
Hello, I am not sure about the p.adjust( , fdr). How do these adjusted p-values get? I have read papers of BH method. For independent case, we compare the ordered p-values with the alfa*i/m, where m is the number of tests. But I have checked that result based on the adjusted p-values is different

Re: [R] Does anybody know how to control the appearance of the end of the line in lattice?

2010-08-08 Thread Paul Murrell
Hi On 9/08/2010 7:49 a.m., George Chen wrote: Hi All, I am plotting vertical lines using xyplot in lattice and type=h. It works well, but the problem is that the tops of the lines are convex and the bottoms are concave. Is there a way to flatten the tops and bottoms? You want to control the

[R] Problem with installing a package in R!

2010-08-08 Thread gagea
Dear friends, I am having problem with installing some packages in R in Ubuntu like ISwR. Bellow is the outcome. Anybody can help me to solve this problem. Thanks a lot and sorry if it is very simple question. install.packages(ISwr) Loading Tcl/Tk interface ... done trying URL

Re: [R] Problem with installing a package in R!

2010-08-08 Thread Liviu Andronic
On Sun, 8 Aug 2010 13:55:28 -0700 (PDT) gagea gagzar...@yahoo.com wrote: ERROR: failed to lock directory ‘/home/ubuntu/R/i486-pc-linux-gnu-library/2.9’ for modifying Try removing ‘/home/ubuntu/R/i486-pc-linux-gnu-library/2.9/00LOCK’ Are you installing from user? Perhaps try from root. Liviu

Re: [R] Data frame reordering to time series

2010-08-08 Thread steven mosher
Thanks again, They worked for me as well. I did a simpler example with fewer years just to show that it worked...( shorted here for display) f - function(x) { +dat - x[-(1:2)] +tim - as.yearmon(outer(x$Year, seq(0, length = ncol(dat))/12, +)) +zoo(c(as.matrix(dat)), tim)

Re: [R] Data frame reordering to time series

2010-08-08 Thread Gabor Grothendieck
On Sun, Aug 8, 2010 at 5:54 PM, steven mosher mosherste...@gmail.com wrote: z-as.zooreg(as.ts(g)) z          X12345 X34567 X56789 1989(1)      NA      3      6 1989(2)      NA      3      6 1989(3)      NA      3      6 1989(4)      NA      3      6 1989(5)      NA      3      6 1989(6)  

Re: [R] Does anybody know how to control the appearance of the end of the line in lattice?

2010-08-08 Thread p_connolly
On Sun, 08-Aug-2010 at 12:49PM -0700, George Chen wrote: | Hi All, | | I am plotting vertical lines using xyplot in lattice and type=h. | It works well, but the problem is that the tops of the lines are | convex and the bottoms are concave. Is there a way to flatten the | tops and bottoms? The

[R] CRAN (and crantastic) updates this week

2010-08-08 Thread Crantastic
CRAN (and crantastic) updates this week New packages * cumSeg (1.0) Vito M.R. Muggeo http://crantastic.org/packages/cumSeg Estimation of number and location of change points in mean-shift (piecewise constant) models. Particularly useful to model genomic sequences of

Re: [R] p.adjust( , fdr)

2010-08-08 Thread Tal Galili
Hi Qing. I believe what you are referring to is: p.adjust(..., BH) The fdr option uses the method of Benjamini, Hochberg, and Yekutieli. (Not the original Benjamini, Hochberg article) It might be the method described here:

[R] efficient matrix element comparison

2010-08-08 Thread david h shanabrook
It is a simple problem in that I simply want to convert the For loop to a more efficient method. It simply loops through a large vector checking if the numeric element is not equal to the index of that element. The following example demonstrates a simplified example: rows - 10 collusionM -

Re: [R] efficient matrix element comparison

2010-08-08 Thread Henrique Dallazuanna
Try this: Matrix(diag(duplicated(matchM)), sparse = TRUE) On Sun, Aug 8, 2010 at 9:43 PM, david h shanabrook dhsha...@acad.umass.eduwrote: It is a simple problem in that I simply want to convert the For loop to a more efficient method. It simply loops through a large vector checking if the

[R] Need help on heatmap, K-means and hhierarchical clustering methods

2010-08-08 Thread meetsiddu1
Hi folks, I am new to the R software. I have been going through different materials to know more about R. I have the R software installed on my windows machine.I would like to know the R source code for the following problems on iris flower data set. I need to do the cluster analysis project

[R] Regular Expression

2010-08-08 Thread Steven Kang
Hi all, From a list of strings, I desire to filter out the followings: 1. Digits at the beginning of the strings 2. Character SPE following the digits (if it exists) 3. Any characters followed by hyphen The following produces the desired result, but would like to know whether this can be done

[R] Invalid plot type '1'

2010-08-08 Thread Stephen Liu
Hi folks, x [1] 1 2 3 4 5 y [1] 6 7 8 9 10 plot(x,y, type = 1 ) Error in plot.xy(xy, type, ...) : invalid plot type '1' I looked on man plot but unable to resolve. Please help. TIA B.R. Stephen L __ R-help@r-project.org mailing list

Re: [R] Need help on heatmap, K-means and hhierarchical clustering methods

2010-08-08 Thread Erik Iverson
On 08/08/2010 08:36 PM, meetsiddu1 wrote: Hi folks, I am new to the R software. I have been going through different materials to know more about R. I have the R software installed on my windows machine.I would like to know the R source code for the following problems on iris flower data set.

Re: [R] Invalid plot type '1'

2010-08-08 Thread Erik Iverson
On 08/08/2010 08:47 PM, Stephen Liu wrote: Hi folks, x [1] 1 2 3 4 5 y [1] 6 7 8 9 10 plot(x,y, type = 1 ) Error in plot.xy(xy, type, ...) : invalid plot type '1' I looked on man plot but unable to resolve. Please help. TIA You never say what you expect type = 1 (that's the

Re: [R] Regular Expression

2010-08-08 Thread Wu Gong
gsub(pattern = ^[0-9]+ (SPE )*(\\w+) - .*$, \\2, dat) - A R learner. -- View this message in context: http://r.789695.n4.nabble.com/Regular-Expression-tp2318086p2318101.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Invalid plot type '1'

2010-08-08 Thread Stephen Liu
Hi Erik, I followed following video; Statistics with R (part 3: plot and history tutorial) http://www.youtube.com/watch?v=NfH5peM1RtIfeature=related It is type = 1. A graph was generated there. But I got an error msg B.R. Stephen - Original Message From: Erik Iverson

Re: [R] efficient matrix element comparison

2010-08-08 Thread Nikhil Kaza
How about a - which(row(matchM)!=matchM) b - matchM[a] diag(collusionM[a,b]) -1 Nikhil Kaza Asst. Professor, City and Regional Planning University of North Carolina nikhil.l...@gmail.com On Aug 8, 2010, at 8:43 PM, david h shanabrook wrote: It is a simple problem in that I simply want to

Re: [R] Invalid plot type '1'

2010-08-08 Thread Erik Iverson
On 08/08/2010 09:14 PM, Stephen Liu wrote: Hi Erik, I followed following video; Statistics with R (part 3: plot and history tutorial) http://www.youtube.com/watch?v=NfH5peM1RtIfeature=related It is type = 1. A graph was generated there. But I got an error msg No, it was type = l (the

Re: [R] Regular Expression

2010-08-08 Thread Michael Bedward
Was going to suggest gsub(^[0-9]+ (SPE )?([^ -])( -.*)?, \\3, s) but I see Wu Gong beat me to the punch with a nicer one :) On 9 August 2010 12:13, Wu Gong w...@mtmail.mtsu.edu wrote: gsub(pattern = ^[0-9]+ (SPE )*(\\w+) - .*$, \\2, dat) __

Re: [R] Regular Expression

2010-08-08 Thread Michael Bedward
And my \\3 should have been a \\2 anyway ! On 9 August 2010 12:23, Michael Bedward michael.bedw...@gmail.com wrote: Was going to suggest gsub(^[0-9]+ (SPE )?([^ -])( -.*)?, \\3, s) but I see Wu Gong beat me to the punch with a nicer one :) On 9 August 2010 12:13, Wu Gong w...@mtmail.mtsu.edu

[R] R Base Code

2010-08-08 Thread nancy_shackelford
Hello all, I need to imitate the 'q' function (qnorm, qweibull, etc) for Clark's 2Dt distribution model. I'm not skilled enough in R to code it myself, so I thought I could find the base code for one of the existing 'q' functions and just modify it. However, I'm having trouble navigating

Re: [R] Invalid plot type '1'

2010-08-08 Thread Stephen Liu
Hi Erik, I got it. Thanks On R-Project The R Manuals http://www.r-project.org/ There are several manuals. Whether I should follow:- An Introduction to R http://www.r-project.org/ B.R. Stephen L - Original Message From: Erik Iverson er...@ccbr.umn.edu To: Stephen Liu

[R] Good Book To Work Through This Summer

2010-08-08 Thread TGS
Dear R users, I'm hoping to get a few suggestions about which books are good to follow along and learn R. I'm hoping to spend the summer going through a good R book as it is applied in linear regression. Thanks! __ R-help@r-project.org mailing list