Re: [R] PGF Device

Dear Cameron, thanks for that message and your effort. I'm really glad to hear this. Lasse Am Donnerstag, den 07.05.2009, 15:09 -0700 schrieb cameron.bracken: > > baptiste auguie-2 wrote: > > > > I think the pgfSweave project on R-forge is working on this (as far as > > i know it currently re

[R] dyn.load a package for dev purposes?

Hi, I don't seem to be able to get a package to detach and reload in my debugging R session once I have already loaded it. R CMD REMOVE -l /usr/lib64/R/library quantmod R CMD build /home/user1/eclipse/workspace/quantmod/; R CMD INSTALL -l /usr/lib64/R/library quantmod_0.3-7.tar.gz in R... de

Re: [R] I'm offering $300 for someone who know R-programming to do the assignments for me.

Homework?... Regards, Yihui -- Yihui Xie Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin University of China, Beijing, 100872, China On Fri, May 8, 2009 at 12:19 PM, ayaku1..

[R] I'm offering $300 for someone who know R-programming to do the assignments for me.

There are six assignments in total. It won't take you long if you were familiar with R. For those who are interested, please send me an email with your profile (your experience with R, how long and how often have you been using it.) I will be paying through paypal. Thanks!

[R] display data ( text) on leveplot / Superpose panel

I am trying to display data ion top of colored rectangle on the same levelplot. While using superpose the shape and symbol for levelplot changes to small circle. I am unable to figure what exactly I am doing wrong... It may be straightforward solution... .I will appreciate help to fix this issu

[R] if ((x >.2 || x<(-.2)) && (col(x)!=row(x))) {x=x[,-col(x)]}

Hi. I have a correlation matrix 'x' which is of size 923x923 I need to remove variables that are highly correlated. I don't have a sophisticated way of selecting which of the two in a highly correlated pair to remove. I thought I would just go through each entry of the correlation matrix and if i

[R] mgarch

The mgarchBEKK package can be downloaded from the author's website at http://www.vsthost.com/vstDocs/mgarchBEKK/release/ Cheers, Joye [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

Re: [R] Do you use R for data manipulation?

>> Farrel Buchinsky wrote: >>> Is R an appropriate tool for data manipulation and data reshaping and data >>> organizing? I think so but someone who recently joined our group thinks >>> not. >>> The new recruit believes that python or another language is a far better >>> tool for developing data ma

[R] why doesn't t.test value include standard error

Hello, everybody! I'm back to ask the obvious on behalf of the silent majority :) Today a student asked me "what standard error was used in this t.test output". I looked into it and was a little surprised that a t.test output object does not have a "slot" for the standard error. Of course, we c

[R] DPpackage help: DPdensity

Hi all, in the function DPdensity( ) example -- ?DPdensity --, we can extract the density estimate by: # Extracting the density estimate cbind(fit1.1$x1,fit1.1$dens) cbind(fit1.2$x1,fit1.2$dens) How is density calculated? I can imagine it is a mixture of normal densities, but where is the in

Re: [R] paste with apply, spaces and NA

Hi, On Thu, May 7, 2009 at 6:36 PM, William Dunlap wrote: > Do you get similar results rif you use 10 instead of NA > in your examples, with more spaces if you use 1? No. Only that particular combination of NA values and strings with spaces gives the odd result. > I think this has to do wi

[R] problem with normalize.quantiles

Hi All, I am using normalize.quantiles (package preprocessCore) to process a 8300 row x 96 col matrix, all numerical data. When I limit the input matrix to 8300 x 67, the function works fine and my data is normalized. But when I increase the input matrix size to 8300 x 68, the function proces

[R] I need your help about graph. I would like to give color to text in the graph.

Hi. I want to show easily classified samples in the one biplot. My samples have 14 types and each type have 21 samples So, I cannot understand distribution of samples after plsda. I want to give a color each types text in the graph. How can I do that? [[alternative HTML version delete

Re: [R] How to properly shade the background panels of an xyplot?

On Thu, May 7, 2009 at 2:50 PM, Daniel Kornhauser wrote: > Hi: > > Here is simplified example of what I am having trouble with: > I want to set the gray shade of the background of each panel of a > xyplot with its mean. > My aim is to be able to compare at a glance which panel has the highest mean

Re: [R] by-group processing

how about: d=data[order(data$ID,-data$Type),] d[!duplicated(d$ID),] Max Webber wrote: > > Given a dataframe like > > > data > ID Type N > 1 45900A 1 > 2 45900B 2 > 3 45900C 3 > 4 45900D 4 > 5 45900E 5 > 6 45900F 6 > 7 45900I 7 > 8

Re: [R] RSPerl and Statistics::R

Being a Perl hacker for some time, and wanting to leverage what R provides, I've been trying to work with Statistics::R and RSPerl. The former has a race condition that breeds some unreliability and the latter seems to have issues all around, and neither has been updated in some time. I h

Re: [R] paste with apply, spaces and NA

> From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Sarah Goslee > Sent: Thursday, May 07, 2009 3:00 PM > To: r-help > Subject: [R] paste with apply, spaces and NA > > Hello everyone, > > I've come up with a problem with using paste() inside apply() that I >

[R] Hmisc::latex() help

Dear R-helpers, I have a data frame (info[[j]]) that looks like this: necklace originalPosition same run surprise1 pairsurprise2 triplesurprise3 1110 2 2.709511 11 4.81884167930642110 2.70951129135145 2120 2

Re: [R] Using lme() for split plot

Crossed Random effects are difficult using lme (wasn't designed for it). Try lmer in the lme4 package if you need this. Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Joshua Stults

Re: [R] Using lme() for split plot

I should say I'm using Google books to look at 'Mixed effects models...' so I can't see pp 49 - 50. On Thu, May 7, 2009 at 6:02 PM, Joshua Stults wrote: > That's a good example with a couple levels of nesting (similar to the > examples in the other book), but they still only have one factor, > 'V

Re: [R] PGF Device

baptiste auguie-2 wrote: > > I think the pgfSweave project on R-forge is working on this (as far as > i know it currently relies on eps2pgf) > Hi- I am one of the developers of pgfSweave. This is correct, pgfSweave currently relies on eps2pgf. Our goal for the summer is to develop a pgf d

[R] paste with apply, spaces and NA

Hello everyone, I've come up with a problem with using paste() inside apply() that I can't seem to solve. Briefly, if I'm using paste to collapse the rows of a data frame, AND the data frame contains strings with spaces, AND there are NA values in subsequent columns, then paste() introduces spaces

Re: [R] Using lme() for split plot

That's a good example with a couple levels of nesting (similar to the examples in the other book), but they still only have one factor, 'Variety', nested in each block. Am I missing something? Should I make up a psuedofactor with four levels to code my two two-level factors? On Thu, May 7, 2009

[R] How to properly shade the background panels of an xyplot?

Hi: Here is simplified example of what I am having trouble with: I want to set the gray shade of the background of each panel of a xyplot with its mean. My aim is to be able to compare at a glance which panel has the highest mean. But, in order to achieve this I have to normalize the means between

Re: [R] PGF Device

In fact, eps2pgf worked fine in my example after building it from svn. Am Donnerstag, den 07.05.2009, 13:28 -0600 schrieb Greg Snow: > Have you contacted the author of eps2pgf? > > A while back I had a plot generated by R that produced an error in eps2pgf, I > sent a copy of it to the author of

Re: [R] PGF Device

No I haven't done that, and of course, you're right. I'll do that. But I still wish there was a way to generate PGF code directly. Lasse Am Donnerstag, den 07.05.2009, 13:28 -0600 schrieb Greg Snow: > Have you contacted the author of eps2pgf? > > A while back I had a plot generated by R that prod

[R] making an Affybatch object for expresso

Hi All, I have microarray data that does not come in a CEL file. Currently it is in the form of columns = individual samples and rows = individual probes. There are about 79 columns and it is in a tab delimited text file. Is there a way to convert this file into an AffyBatch so that I can ru

[R] Using lme() for split plot

Hi, I'm trying to figure out how to use lme() for analyzing a split-plot experiment. I've been looking at the examples from the 'R Book', those are nested but with only one factor at the whole-plot level, my test is 2^2 at the whole-plot level, with a single many level factor at the sub-plot leve

Re: [R] Do you use R for data manipulation?

+1. I worked with Matthew for a while and saw in practice just how powerful that package is. I'm surprised it isn't more widely used. Martin Tom Short wrote: Another tool I find useful is Matthew Dowle's data.table package. It has very fast indexing, can have much lower memory requirements th

Re: [R] for loop vectorization

this 'ifelse' usage looks promising. thank you very much. On Thu, May 7, 2009 at 3:12 PM, Patrick Burns wrote: > If you haven't seen it yet, > 'The R Inferno' may be of use > to you. > > > Patrick Burns > patr...@burns-stat.com > +44 (0)20 8525 0696 > http://www.burns-stat.com > (home of "The R I

Re: [R] by-group processing

Max, Since the dataset is sorted by ID, with ties broken by N, the following should do it and do it quickly. It grabs the rows just before ID changes. > with(data, data[ c(ID[-1] != ID[-length(ID)], TRUE),, drop=FALSE]) ID Type N 7 45900I 7 10 49270E 3 24 46550I 7 Bill Dunlap

Re: [R] PGF Device

Have you contacted the author of eps2pgf? A while back I had a plot generated by R that produced an error in eps2pgf, I sent a copy of it to the author of eps2pgf and within a couple of days he sent me a new version that worked with my plot. But he cannot fix problems that he does not know abo

Re: [R] levelplot question

On Thu, May 7, 2009 at 12:02 PM, Greg Snow wrote: > From the help for panel.levelplot in the section on 'col.regions': > > " the exact number being one >          more than the length of 'at'" > > was what I based my assumption on (also the "not longer than at" phrase in > the last sentence on re

Re: [R] levelplot question

>From the help for panel.levelplot in the section on 'col.regions': " the exact number being one more than the length of 'at'" was what I based my assumption on (also the "not longer than at" phrase in the last sentence on recycling). The second part was somewhat speculation on my par

Re: [R] dependency ‘Matrix’ is not availabl e when installing Rcmdr

Uwe Ligges wrote: Argh, it disappeared. Probably a bug in my scripts that removed the old version of Matrix from the Windows binary repository for R-2.8.x (new versions of Matrix won't pass checks under the outdated version of R). I will compile an old version tomorrow and make it available a

Re: [R] for loop vectorization

uhoh, missed two lines on the top.Sorry about that. the whole code looks like this. nchars = 4 maxvals = c(2,2,2,2) tmp = 0:(maxvals[1]) for(i in 2:nchars) { tmp <- outer(tmp, 0:(maxvals[i]), FUN="paste", sep=".") } states = tmp stateidx = array(1:length(states), dim=dim(states)) transition <-

[R] GAM ordered probit

Dear All, Anyone know if there is a package that fits Generalized Linear Models(GAM) to data with ordered dependent variable(response) ? Simon Wood's mgcv has probit, logit,... other links, however, I could not find a way to do GAM *ordered *probit. Yee's VGAM claims to fit ordinal proportional o

[R] for loop vectorization

Hi, I am still not familiar with vectorization. Could you help with making this for loop more efficient? The code is trying to make a Q matrix for a multidimensional state space with specific conditions. thanks Mira tmp = 0:(maxvals[1]) for(i in 2:nchars) { tmp <- outer(tmp, 0:(maxvals[i]), FU

[R] CHM error

Hi, I would need some help in solving a problem with my R. Whenever I try to type " ?paste " or in general any command for the html help in R, I receive this error message. The strange thing is that I have never had this kind of problem in the past with R, and now I really don't understand why it

[R] increasing memory for R bg job

Hi, Is the following command used to increase the memory or any other command when a background R job is run? R --min-vsize=vl --max-vsize=vu --min-nsize=nl --max-nsize=nu --max-ppsize=N source: http://stat.ethz.ch/R-manual/R-patched/library/base/html/Memory.html Thx Carol

[R] Weighted multinomial logistic regression using the mlogit package

I have been trying to use the mlogit package to do a multinomial logistic regression, including both alternative-specific and individual-specific variables. I used the mlogit.data function to turn my dataframe into the correct format for the mlogit function, and have been able to run the re

Re: [R] Help with lme4 model specification

Thanks for your response. >This is suspicious : Spar_loc doesn't appear in your statement... Yes Spar_loc should look suspicious. I posted the real error messages but hid the variable names in the model. Spar_loc is the original form of S. >The precedence of modelling operators is not very c

[R] help~ how can I drow colorful text in biplot.

Hi. I want to show easily classified samples in the one biplot. My samples have 14 types and each type have 21 samples So, I cannot understand distribution of samples after plsda. I want to give a color each types text in graph. How can I do that? [[alternative HTML version deleted]]

Re: [R] How to do Naive Bayes in R?

On Wednesday 06 May 2009, spencerg wrote: > help.search('bayes') only searches installed packages. > > To go beyond that, you might try the following: Thanks for the clarification. Dylan -- Dylan Beaudette Soil Resource Laboratory http://casoilresource.lawr.ucdavis.edu/ University

Re: [R] 64 bit R double precision

Paul Emberson wrote: Hi, If R is compiled for a 64 bit platform, does the double data type make use of the additional precision or is it set at a fixed length? No. Double is double, 64 bits, although sometimes 80 bits when things are kept in registers. However, even though the FPU is identic

Re: [R] Step and stepAIC

Jun Shen gmail.com> writes: > >lme(biomarker~Cmax+AGE,data=cov,method='ML',random=~1+AGE|ID,keep.data=T) > Warning message: > In lme.formula(biomarker ~ Cmax + AGE, data = cov, method = "ML", : > Fewer observations than random effects in all level 1 groups That's why the posting guide says:

Re: [R] 64 bit R double precision

On 5/7/2009 12:56 PM, Paul Emberson wrote: Hi, If R is compiled for a 64 bit platform, does the double data type make use of the additional precision or is it set at a fixed length? It is still a 64 bit floating point value. R may use "long double" internally for some calculations, but the

[R] .First and loading a package in a workspace

Hi, I'm writing a package called rNMR and I'd like to be able to save a workspace with rNMR loaded, exit R, and then load the workspace, loading rNMR at the same time. This is no problem on Windows but on Linux and Macs I can't get it to work. Here's the deal. rNMR has includes a .onLoad functio

Re: [R] RSPerl and Statistics::R

RSPerl is not abandonded, but I haven't had much time to work on it. If you can send me an example to reproduce the race condition, we can probably think about a solution. D. Daniel Klein wrote: Greetings! Being a Perl hacker for some time, and wanting to leverage what R provides, I've be

Re: [R] Step and stepAIC

Hi, Dieter, > sessionInfo() R version 2.9.0 (2009-04-17) i386-pc-mingw32 > str(cov) 'data.frame': 40 obs. of 13 variables: $ ID : int 2 3 4 5 6 7 8 9 10 11 ... $ AMT : num 50 50 50 50 50 50 50 50 25 25 ... $ Cmax : num [1:40(1d)] 8901 5574 10496 10882 8130 ... ..- attr(*,

[R] 64 bit R double precision

Hi, If R is compiled for a 64 bit platform, does the double data type make use of the additional precision or is it set at a fixed length? Thanks, Paul __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Step and stepAIC

Jun Shen gmail.com> writes: > I changed the "method" to ML for my lme call. Now I get a different error > when running the stepAIC > > Error in terms.formula(formula, data = data) : > 'data' argument is of the wrong type > > What does that mean? > The same dataset was used for an lm call and

Re: [R] A potential bug for paste() ?

Thanks for everyone's attention on this issue! On Thu, May 7, 2009 at 11:20 AM, Phil Spector wrote: > Theoph is in the datasets package, so it should be available in all > versions of R: > > class(Theoph) >> > [1] "nfnGroupedData" "nfGroupedData" "groupedData""data.frame" > > The key to wha

Re: [R] A potential bug for paste() ?

Theoph is in the datasets package, so it should be available in all versions of R: class(Theoph) [1] "nfnGroupedData" "nfGroupedData" "groupedData""data.frame" The key to what's happening is that there are some methods for groupedData objects in the nlme library, which I suspect the orig

Re: [R] dotplot: labeling coordinates for each point

On Thu, May 7, 2009 at 7:42 AM, Qifei Zhu wrote: > Hi Deepayan, > > > > I guess I’m quite confusing when asking the question. I was actually looking > for labels beside the dot in the graph (shown below), rather than the labels > on the x-axis. Is that possible? Thanks a lot! In that case dotplo

Re: [R] A potential bug for paste() ?

On Thu, May 7, 2009 at 9:56 AM, Sarah Goslee wrote: > It probably has less to do with paste() than with Theoph, but since > we have no idea what that might be, it's hard to tell. > See the bit about "reproducible example", please. Well, actually Theoph is one of the datasets in the (required) da

Re: [R] How to call time series functions from C ?

On 5/7/2009 11:17 AM, goon83 wrote: Hi Duncan, Sorry, i miss the informatin in the previous message, you said that i can write a "front end" to call , does the front end can be implemented by c, and i can compile the c code to a standalone program. i will appreiciate for more detail informa

Re: [R] Linear least squares fit with errors in both x and y values.

James, look up "errors in variables" models or "instrumental variable" models in econometrics. The statistics alternative is a "random effects" or "mixed effects" model which plugs the variation in the x's into a randomly varying parameter - these are available in R (?lmer or glmm - I think). Som

Re: [R] Step and stepAIC

Thanks a lot, Dieter, I changed the "method" to ML for my lme call. Now I get a different error when running the stepAIC Error in terms.formula(formula, data = data) : 'data' argument is of the wrong type What does that mean? The same dataset was used for an lm call and everything is fine to r

Re: [R] I updated/reinstalled ggplot2 and the trouble started...

Hi, > What device are you using to plot the graph? I am not explicitly specifying any. I am on redhat linux and as I understand it R uses DISPLAY environment variable, which is set for an X11 server (X Darwin, v1.4.0, X11 6.8.2). Other qplot()'s I have used (straight out of the book) work fine.

Re: [R] A potential bug for paste() ?

On Thu, 2009-05-07 at 09:38 -0500, Jun Shen wrote: > Hi, everyone, > > Try the following command to see if you get TRUE or FALSE. I get FALSE on a > unix platform but TRUE on Windows. Any comment? > > all(paste(Theoph[1],Theoph[2])==paste(Theoph[[1]],Theoph[[2]])) And what is Theoph? Please do r

[R] lasso based selection for mixed model

Dear useRs (called Frank Harrell, most likely), after having preached for years to my medical colleagues to be cautious with stepwise selection procedures, they chanted back asking for an alternative when using mixed models. There is a half dozen laXXX packages around for all types of linear mod

Re: [R] Step and stepAIC

Jun Shen gmail.com> writes: > > I also wonder if stepAIC can work with an lme model. I tried. it doesn't for > now. Appreciate any comment. > It can do it, and it displays a nice error message, as I remember, why it won't do it in the default setting. Check parameter method of your lme call.

[R] visualization of sub-parts in 3D arrays

Hello all; Let's say we have a block of 3D data belonging to cubic cells, and increasing indexes are ordered x: from west to east y: from north to south z: from top to bottom --- nx <- 25 ny <- 12 nz <- 4 my.array <- array(NA,dim=c(nx,ny,nz)) for(z in 1:nz){ my.array[,,z] <- runif(nx*ny,z-1

Re: [R] Step and stepAIC

Jun Shen gmail.com> writes: > > I also wonder if stepAIC can work with an lme model. I tried. it doesn't for > now. Appreciate any comment. > It can, and it displays a nice error message telling you exactly why it did not work if you try out of the box. Check parameter "method" in your lme cal

Re: [R] How to call time series functions from C ?

Do you have to use the entire xxxGB at once? If not store it in a database and fetch the portions for the current analysis as you need them. On Thu, May 7, 2009 at 10:55 AM, goon83 wrote: > hi Duncan, > > tks for you response, it's depressed to know that i can not use time series > of R to fini

[R] compiling 64 bit JRI.dll with 64bit XP

Hi all, I was able to run my Java-R program using a 32-bit jdk on my 64-bit machine; however I wish to use the 64-bit jdk so my Java program can take advantage of the larger memory capabilities etc. I'm currently trying to compile JRI_0.4-1 on a 64-bit XP I got to a point that I was able to gener

Re: [R] Help with lme4 model specification

Thanks! I made this correction ONKELINX, Thierry wrote: > > > Please note that the case of characters is important in R. Hence "True" > differs from "TRUE". You need to use REML = TRUE instead of REML = True! > The only valid logical values are TRUE, FALSE, T and F. But I recommend > to use on

[R] How to save independent R sessions in BATCH mode???

Hi guys, I already have put this question some time ago however replies were not straight to the problem. Even so thanks a lot for the answers! I reformulate the problem. I have a script file with the following commands: for (i in 1:300) - open session --> load(“X_i.rdata”) - perform some c

[R] Linear least squares fit with errors in both x and y values.

HI, I'd like to perform a weighted linear least squares fit with R on data with varying errors on both vectors. I can do this with one axis using lm, but have no idea where to go from here. I've tried googling, but no idea. Any suggestions? Thanks, James _

[R] Draw a random number about the fitting of a glm model using family=quasi

Hi, I desire to draw a random number about the fitting of a glm model using a quasi linear error distribution. For example, if my family is poisson, then this is what I want: model1 = glm(amount~days,family=poisson(link="log"),data=dataset1) mean1 = predict(model1,newdata =data.frame(days=time),

Re: [R] How to call time series functions from C ?

you have a couple of options. If you require specific R functions to do what you want, then you will need to call R from C. I believe that Dirk has been working on an "RInside" package that does this. Alternatively, you can use my tslib package, which is a general time series library written in

Re: [R] Bumps chart in R

Thank you Thierry! That was very helpfull I have added "geom_point" which makes the lineending look more round. Also - since the chart becomes verye clustered I also added text to the right. like so: ggplot(df, aes(x = tal, y = value, label = text, group = text, size = size, col = Change))+ geom

Re: [R] How to call time series functions from C ?

Hi Duncan, Sorry, i miss the informatin in the previous message, you said that i can write a "front end" to call , does the front end can be implemented by c, and i can compile the c code to a standalone program. i will appreiciate for more detail information tks? ?

Re: [R] Step and stepAIC

I also wonder if stepAIC can work with an lme model. I tried. it doesn't for now. Appreciate any comment. Jun On Thu, May 7, 2009 at 9:15 AM, Dieter Menne wrote: > Peter-Heinz Fox yahoo.de> writes: > > > I’m using "step" and "stepAIC" for stepwise regression. After each > step, > > I would li

[R] Linear least squares fit with errors in both x and y values.

Hi, I'd like to perform a weighted linear least squares fit with R on data with varying errors on both vectors. I can do this with one axis using lm, but have no idea where to go from here. I've tried googling, but no idea. Any suggestions? I'm also trying to plot these errors on a graph. Again I

Re: [R] How to call time series functions from C ?

On 5/7/2009 10:55 AM, goon83 wrote: hi Duncan, tks for you response, it's depressed to know that i can not use time series of R to finish my study. I really think R is powerfull, i used it to finish all my statistic homework last year^_^. but i still wonder can i use the function of time se

Re: [R] PGF Device

I can only second Dieter and Baptiste. Am Donnerstag, den 07.05.2009, 14:57 +0200 schrieb baptiste auguie: > If one was to implement a Tikz output for R graphics, I'm not sure > what route to take: > > 1) focus exclusively on Grid graphics and implement a conversion at > the R level à la gri

Re: [R] A potential bug for paste() ?

It probably has less to do with paste() than with Theoph, but since we have no idea what that might be, it's hard to tell. See the bit about "reproducible example", please. Sarah On Thu, May 7, 2009 at 10:38 AM, Jun Shen wrote: > Hi, everyone, > > Try the following command to see if you get TRU

Re: [R] How to call time series functions from C ?

hi Duncan, tks for you response, it's depressed to know that i can not use time series of R to finish my study. I really think R is powerfull, i used it to finish all my statistic homework last year^_^. but i still wonder can i use the function of time series in the embedded way of C code.

[R] A potential bug for paste() ?

Hi, everyone, Try the following command to see if you get TRUE or FALSE. I get FALSE on a unix platform but TRUE on Windows. Any comment? all(paste(Theoph[1],Theoph[2])==paste(Theoph[[1]],Theoph[[2]])) -- Jun Shen PhD PK/PD Scientist BioPharma Services Millipore Corporation 15 Research Park Dr.

Re: [R] Bubble Plot Over a Map

Hi Rodrigo, did you consider having a look at package "ggplot2" (thanks to Hadley Wickham). I guess you'll find something very close to your need (functions ggplot, geom_point,...) Hope this helps. regards. Olivier Rodrigo Aluizio wrote: > > Hi everybody, > > Well here is another doubt! I’m el

Re: [R] error using lapply with oneway_test (coin package)

That's a scoping problem, I think a bug in oneway_test. Because the formula var ~ group is created with the anonymous function within lapply, its environment should be the evaluation frame of that function call and var should be visible. If I replace oneway_test() with lm() it works. I thi

Re: [R] Step and stepAIC

Peter-Heinz Fox yahoo.de> writes: > I’m using "step" and "stepAIC" for stepwise regression. After each step, > I would like to make an additional calculation based on the independent > variables that have been selected until this step and > their corresponding weights. Where do I have to add t

Re: [R] Generating plots with two lines, each line plotted with a different color

On Thu, 7 May 2009, John Sorkin wrote: R 2.9.0 windows xp I would like to have the two lines generated by my plot have different colors, in addition to the two different line types specified by the plot command below. I would appreciate advice on how I can specify line colors. plot(survfit(S

Re: [R] data transformation using gamma

Roslina, this code performs what you need: dt = matrix((1:(58*12))/58/12,58) # some numbers # if dt is a data.frame use dt = as.matrix(dt) a = (1:12)/12 # some a coef b = (12:1)/12 # some b coef dtgam = matrix(pgamma(dt,a,b),58) # dtgam is the transformation you're looking for

[R] Plotting counts vs. intervals for a shingle

Hello! Suppose I have a set of values: a <- c(1:10, 5:10) Suppose I also have a set of intervals: b <- cbind(c(0,2.5,4.5,6.5), c(5.5,7.5,9.5,11)) I can create a shingle that counts how many values are in each interval: c <- shingle(a, b) I can display the shingle to see the counts: summary(

Re: [R] How to call time series functions from C ?

On 5/7/2009 9:24 AM, goon83 wrote: Hi everyone, I meet the problem to call R from C. Can someone tell me is it possible to call time series functions(such as ar and arima)from C program on linux platform, and then compile the C program to a standalone executable file? I know we can ca

Re: [R] Print to File Formatting

Brigid Mooney wrote: Hello, I am writing out to a file and have two quick questions that I can't seem to track down the correct answers for. Luckily, I *think* they are both simple enough that someone might be able to point me in the right direction on them without too much trouble. Both ques

[R] Step and stepAIC

Hi all,   I’m using "step" and "stepAIC" for stepwise regression. After each step, I would like to make an additional calculation based on the independent variables that have been selected until this step and their corresponding weights. Where do I have to add this calculation?   And a secon

Re: [R] Show name of dataset in graph

The usual way of doing this in R functions is via deparse, as in the examples in ?deparse: fd<-function(x) plot(x, main=deparse(substitute(x))) q<-1:10 fd(q) as.character works for simple cases but is probably not the best option. Compare fd(log(q)) with fc<-function(x) plot(x, main=as.chara

[R] How to call time series functions from C ?

Hi everyone, I meet the problem to call R from C. Can someone tell me is it possible to call time series functions(such as ar and arima)from C program on linux platform, and then compile the C program to a standalone executable file? I know we can call some function defined in Rmath.h

[R] extending strsplit to handle missing text that doesn't have the target on which to split

I am sure there is an obvious answer to this that I'm missing but I can't find it. I'm parsing headers of Emails and most have a date like this: "Wed, 16 Nov 2005 05:28:00 -0800" and I can parse that using: tmp.dat.data <- matrix(unlist(strsplit(headers$Date.line,",")), ncol = 2, byrow = T

[R] RSPerl and Statistics::R

Greetings! Being a Perl hacker for some time, and wanting to leverage what R provides, I've been trying to work with Statistics::R and RSPerl. The former has a race condition that breeds some unreliability and the latter seems to have issues all around, and neither has been updated in some ti

Re: [R] When was the package relaimpo first introduced in R?

Vielen Dank, Uwe! Dimitri 2009/5/7 Uwe Ligges : > > > Dimitri Liakhovitski wrote: >> >> Thanks a lot, Uwe. >> Do you think you could expain how one looks this type of information up? >> What's obvious to some of us might not be obvious to others. >> Dimitri, der, obviously, deutsch spricht. > > >

[R] Bubble Plot Over a Map

Hi everybody, Well here is another doubt! I’m elaborating some maps, and I pretend to plot over it the abundance of species. I could do that using points(), but I’d like to use bubble() from sp package. When trying to use bubble() I face the fallowing issue: When I call the bubble function it ove

Re: [R] PGF Device

If one was to implement a Tikz output for R graphics, I'm not sure what route to take: 1) focus exclusively on Grid graphics and implement a conversion at the R level à la gridSVG (see Paul Murrell's webpage) 2) create a new device ( I have no idea where to find guidelines for this, and I

Re: [R] When was the package relaimpo first introduced in R?

Dimitri Liakhovitski wrote: Thanks a lot, Uwe. Do you think you could expain how one looks this type of information up? What's obvious to some of us might not be obvious to others. Dimitri, der, obviously, deutsch spricht. Dimitri, on CRAN there is a package page and the link to the archive

Re: [R] convert large integers to hex

Thanks for both answers. In the end I decided to use Gabor's bc package. Thanks, --sundar On Thu, May 7, 2009 at 5:10 AM, Gabor Grothendieck wrote: > There is an interface between R and bc -- not on CRAN but available > from its home page here: > http://r-bc.googlecode.com > >> source("http://r

Re: [R] When was the package relaimpo first introduced in R?

Thanks a lot, Uwe. Do you think you could expain how one looks this type of information up? What's obvious to some of us might not be obvious to others. Dimitri, der, obviously, deutsch spricht. 2009/5/7 Uwe Ligges : > > > Dimitri Liakhovitski wrote: >> >> Hello! >> I am wondering when the package

  1   2   >