Re: [R] Loading .Rdata within an R function

Thank you for your consideration of this question. I have tried both your suggestions. However, the data is not loaded within the function. When I specify load(mydata.Rdata,.globalEnv), the data is loaded into the top level environment, and the function does access the data in the top level

[R] Appropriate tests for logistic regression with a continuous predictor variable and Bernoulli response variable

I have a data with binary response variable, repcnd (pregnant or not) and one predictor continuous variable, svl (body size) as shown below. I did Hosmer-Lemeshow test as a goodness of fit (as suggested by a kind “R-helper” previously). To test whether the predictor (svl, or body size) has

[R] Aide pour faire ACM

Bonjour, Je travaille sur R et je veux faire MCA, je veux savoir comment calcules les valeurs de Contribution, cosinus carré, coordonnée et distance au centre des modalités des différentes variables. Je vous remercie beaucoup d'avance. Housseima

[R] Non-parametric regression

I have two data sets, each a vector of 1000 numbers, each vector representing a distribution (i.e. 1000 numbers each of which representing a frequency at one point on a scale between 1 and 1000). For similfication, here an short version with only 5 points. a - c(8,10,8,12,4) b - c(7,11,8,10,5)

Re: [R] Non-parametric regression

From reviewing the first google page result for Non-parametric regression R, I hope this link will prove useful: http://socserv.mcmaster.ca/jfox/Courses/Oxford-2005/R-nonparametric-regression.html Contact Details:--- Contact

Re: [R] Kite diagrams

On 07/09/2010 07:23 AM, Graham Smith wrote: I asked the same question on R-sig-eco, and Ben Bolker provided this solution, which as I assume this should show up in a search I copy here. However, if someone can come up with a single function, that would be good. Hi Graham, library(plotrix)

[R] R^2 in loess and predict?

Parametric regression produces R^2 as a measure of how well the model predicts the sample and adjusted R^2 as a measure of how well it models the population. What is the equalvalent for non-parametric regression (e.g. loess function) ? Ralf __

Re: [R] R^2 in loess and predict?

I do not agree with your interpretation of the adjusted R^2. The R^2 is no more than the ratio of the explained variance by the total variance, expressed in sums of squares. The adjusted R^2 is adjusted for the degrees of freedom, and can only be used for selection purposes. The interpretation

Re: [R] Kite diagrams

Jim, This is very good news, not so much for me, as I don't use them, but I have colleagues who do, and its an expected graphic in student assignments. So I have passed on the information. So many thanks for adding this, I a sure many people will find it useful. Graham On 9 July 2010 09:41,

[R] accessing return variables from a function

Hi, I am trying to figure out a short way to access two values output from the sort function. x - c(3,4,3,6,78,3,1,2) sort(x, index.return=T) $x [1] 1 2 3 3 3 4 6 78 $ix [1] 7 8 1 3 6 2 4 5 It would be great to do something like this (doesn't work.): c(y, indexes) - sort(x,

[R] R crashes with large vectors

Good afternoon, I have been experiencing a lot of crashes working with large vectors in R. Specifically, I am using XTS of length of minimum 120k elements. My problem is that I cannot display the vector (otherwise R crashes), I cannot plot it either (otherwise R crashes). That could be solved

[R] Current script name from R

Is there a way for a script to find out about its own name ? Ralf __ 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,

[R] Plotting text in existing plot?

I would like to plot some text in a existing plot graph. Is there a very simple way to do that. It does not need to be pretty at all (just maybe a way to center it or define a position within the plot). ( ? ) Ralf __ R-help@r-project.org mailing list

Re: [R] Plotting text in existing plot?

see ?text Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English)

Re: [R] Loading .Rdata within an R function

Giles Crane wrote: Thank you for your consideration of this question. I have tried both your suggestions. However, the data is not loaded within the function. When I specify load(mydata.Rdata,.globalEnv), the data is loaded into the top level environment, and the function does access the data

Re: [R] Plotting text in existing plot?

On 07/09/2010 07:52 PM, Ralf B wrote: I would like to plot some text in a existing plot graph. Is there a very simple way to do that. It does not need to be pretty at all (just maybe a way to center it or define a position within the plot). ( ? ) Hi Ralf, The text function in the graphics

[R] split with list

Dear List I would like to ask you something concenting a better print of the R output: I have a bit data frame which has the following structure: CFISCALE RAGSOCBANNO VAR1VAR2. 9853312 astra 2005 6

Re: [R] C module causing rounding errors?

Joseph N. Paulson wrote: Hi all! I am currently writing a C-module for a for loop in which I permute columns in a matrix (bootstrapping) and I send a couple of variables into C initially. All of it is working, except the initial values I send to R are rounded/truncated (I believe rounded). I

Re: [R] accessing return variables from a function

On Jul 9, 2010, at 5:20 AM, Noah Silverman wrote: Hi, I am trying to figure out a short way to access two values output from the sort function. x - c(3,4,3,6,78,3,1,2) sort(x, index.return=T) $x [1] 1 2 3 3 3 4 6 78 $ix [1] 7 8 1 3 6 2 4 5 It would be great to do something

Re: [R] Non-parametric regression

On Jul 9, 2010, at 4:01 AM, Ralf B wrote: I have two data sets, each a vector of 1000 numbers, each vector representing a distribution (i.e. 1000 numbers each of which representing a frequency at one point on a scale between 1 and 1000). For similfication, here an short version with only 5

[R] KLdiv produces NA. Why?

I am trying to calculate a Kullback-Leibler divergence from two vectors with integers but get NA as a result when trying to calulate the measure. Why? x - cbind(stuff$X, morestuff$X) x[1:5,] [,1] [,2] [1,] 293 938 [2,] 293 942 [3,] 297 949 [4,] 290 956 [5,] 294 959 KLdiv(x)

Re: [R] accessing return variables from a function

On Fri, Jul 9, 2010 at 10:20 AM, Noah Silverman n...@smartmediacorp.com wrote: Hi, I am trying to figure out a short way to access two values output from the sort function. x - c(3,4,3,6,78,3,1,2) sort(x, index.return=T) $x [1]  1  2  3  3  3  4  6 78 $ix [1] 7 8 1 3 6 2 4 5 It would

Re: [R] Plotting text in existing plot?

On Fri, Jul 9, 2010 at 11:22 AM, Jim Lemon j...@bitwrit.com.au wrote: The text function in the graphics package will place text on your plot. It centers the text by default, so that: text(3,4,my neat\ntwo liner) will place the two lines of text centered at x=3 and y=4. BUT it only works

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

With appropriate design matrix, I mean the X matrix in the mixed-effects model y = Xb + u + e, where y is the vector of outcomes, u is a vector of (possibly correlated) random effects, and e is a vector of (possibly) random errors. The X matrix is specified via the 'mods' argument in the rma()

Re: [R] accessing return variables from a function

which is just unlist(sort(x, index.return=T)) but I think, Noah would like to have something like srt-sort(x, index.return=T) names(srt)-c(y,indexes) attach(srt) but that is wasteful either. Am 09.07.2010 12:28, schrieb David Winsemius: On Jul 9, 2010, at 5:20 AM, Noah Silverman wrote: Hi,

Re: [R] Current script name from R

I'm assuming you are using Rscript (please provide self-contained examples when posting) in which case you could look for the element in (base|R.utils)::commandArgs() that begin with the string --file= - the rest is the file name. See the asValues= parameter in help(commandArgs,

Re: [R] Current script name from R

I am using RGUI, the command line or the StatET Eclipse environment. Should this not all be the same? Ralf On Fri, Jul 9, 2010 at 7:11 AM, Allan Engelhardt all...@cybaea.com wrote: I'm assuming you are using Rscript (please provide self-contained examples when posting) in which case you could

Re: [R] transformation of data.frame

Hello Petr, sorry for the mixed up. your example works perfectly fine. The one from Søren has shown the mentioned error. But even after reading the columns as character go - read.table(go.txt, header= TRUE, colClasses = c(character, character)) or go - read.table(go.txt, header= TRUE, as.is

[R] sarima.Sim function

Does anyone have some a nice simple example of how this function is used? Thanks, Jason __ 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

[R] Mirror axis on hist2d plot - how?

The following code produces a heatmap based on normalized data. I would like to mirror x and y axis for this plot. Any idea how to do that? require(gplots) x - rnorm(500) y - rnorm(500) hist2d(x, y, freq=TRUE, nbins=50, col = c(white,heat.colors(256))) Best, Ralf

Re: [R] How not to print '\\' as '\\'

On Tue, 2010-07-06 at 19:08 -0700, Joshua Wiley wrote: Try cat(), for instance: cat(\\\n) \ The extra \n is needed for a proper line break. Or writeLines(\\) Which handles the newline for you. ?writeLines. HTH G HTH, Josh On Tue, Jul 6, 2010 at 7:01 PM,

Re: [R] Appropriate tests for logistic regression with a continuous predictor variable and Bernoulli response variable

Have a look at ?anova or rather ?anova.glm hth Am 09.07.2010 04:46, schrieb Kiyoshi Sasaki: I have a data with binary response variable, repcnd (pregnant or not) and one predictor continuous variable, svl (body size) as shown below. I did Hosmer-Lemeshow test as a goodness of fit (as

[R] Data format question for triangle.plot package ade4

hello, I am trying to develop a triangle plot but am having difficultly assigning the row.names to the 3 columns in the data.frame Here is what I've done, attach(SoilVegHydro) dim(SoilVegHydro) 129239 # now take 3 variables from main data.frame for plotting dat - cbind.data.frame(TP,

Re: [R] random sample from arrays

Hi Joris, I guess i did it wrong again. but your example didn't work either. I still get the error massage. but replicate function just fine. I can even replicate the whole array lines. THX Assa On Thu, Jul 8, 2010 at 15:20, Joris Meys jorism...@gmail.com wrote: Don't know what exactly

Re: [R] C module causing rounding errors?

Sorry for not including enough information everyone. I have quite a bit of code, so I will just enter relevant pieces... This is how I call C from R: The tstats are tstatistics (difference of mean, divided by sqrt of S1+S2) from an unpermuted matrix. The c code is below... dyn.load(testp.so)

Re: [R] Data format question for triangle.plot package ade4

On Jul 8, 2010, at 4:41 PM, steve_fried...@nps.gov wrote: hello, I am trying to develop a triangle plot but am having difficultly assigning the row.names to the 3 columns in the data.frame Here is what I've done, attach(SoilVegHydro) dim(SoilVegHydro) 129239 # now take 3

[R] how to plot two histograms overlapped in the same plane coordinate

Dear R-help listers, I am new. I just want to get helps on how to plot two histograms overlapped in the same plane coordinate. What I did is very ugly. Could you please help me to improve it? I want to got a plot with semi- transparent overlapping region. And, I want to know how to specify the

Re: [R] KLdiv produces NA. Why?

On 2010-07-09 4:31, Ralf B wrote: I am trying to calculate a Kullback-Leibler divergence from two vectors with integers but get NA as a result when trying to calulate the measure. Why? x- cbind(stuff$X, morestuff$X) x[1:5,] [,1] [,2] [1,] 293 938 [2,] 293 942 [3,] 297 949 [4,]

[R] distributing a value for a given month across the number of weeks in that month

Hello! Any hint would be greatly appreciated. I have a data frame that contains (a) monthly dates and (b) a value that corresponds to each month - see the data frame monthly below: monthly-data.frame(month=c(20100301,20100401,20100501),monthly.value=c(100,200,300))

[R] Function on columns of a dataframe

Hi, I would like to assign the largest value of a column to a specific category and repeat this for each column (v1 - v4). x=c(1:12) cat=c(cat1,cat5,cat2,cat2,cat1,cat5,cat3,cat4,cat5,cat2,cat3,cat6) v1=rnorm(12,0.5,0.1) v2=rnorm(12,0.3,0.2) v3=rnorm(12,0.4,0.1) v4=rnorm(12,0.6,0.3)

Re: [R] Function on columns of a dataframe

Hi Nils, have a look at ?tapply hth. Am 09.07.2010 15:37, schrieb LogLord: Hi, I would like to assign the largest value of a column to a specific category and repeat this for each column (v1 - v4). x=c(1:12) cat=c(cat1,cat5,cat2,cat2,cat1,cat5,cat3,cat4,cat5,cat2,cat3,cat6)

[R] strange floor rounding

Dear all, might seem and easy question but I cannot figure it out. floor(100*(.58)) [1] 57 where is the trick here? And how can I end up with the right answer? Thanks a lot everybody for your help. Trafim [[alternative HTML version deleted]]

Re: [R] strange floor rounding

You too have fallen victim to floating point error (see FAQ 7.31). (100*.58) == 58 [1] FALSE (100*.58) 58 [1] TRUE On Fri, Jul 9, 2010 at 9:46 AM, Trafim Vanishek rdapam...@gmail.com wrote: Dear all, might seem and easy question but I cannot figure it out. floor(100*(.58)) [1] 57

Re: [R] strange floor rounding

On Jul 9, 2010, at 8:46 AM, Trafim Vanishek wrote: Dear all, might seem and easy question but I cannot figure it out. floor(100*(.58)) [1] 57 where is the trick here? And how can I end up with the right answer? Thanks a lot everybody for your help. Trafim sprintf(%.20f, 100 *

Re: [R] Function on columns of a dataframe

On Jul 9, 2010, at 9:46 AM, Eik Vettorazzi wrote: Hi Nils, have a look at ?tapply hth. Perhaps this will be part way there (I couldn't really figure out the desired structure of the final object): lapply( bla[, -(1:2)], function(x) tapply(x, bla$cat, max) ) $v1 cat1 cat2

Re: [R] strange floor rounding

On Jul 9, 2010, at 9:46 AM, Trafim Vanishek wrote: Dear all, might seem and easy question but I cannot figure it out. floor(100*(.58)) [1] 57 where is the trick here? FAQ 7.31 And how can I end up with the right answer? Define right, please. (There have been several questions in the

[R] Help for mca

Good eveninig, I work on R and i want to do an MCA but, i find difficulty in calculating the contribution (CTR), cosine(CO2), coord(CORD) and distancefrom the center of the modalities of variables. Please help me how to calculate these values. In fact, when i enter the variables it writes

Re: [R] Non-parametric regression

Just to be correct : gam is mentioned on the page Tal linked to, but is a semi-parametric approach using maximum likelihood. It stays valid though. Another thing : you detect non-normality. But can you use a Poisson distribution for example? The framework of generalized linear models and

Re: [R] split with list

One solution is to put these unwanted entries to repor$9853312 [1:2,2:3] - Cheers Joris On Fri, Jul 9, 2010 at 12:18 PM, n.via...@libero.it n.via...@libero.it wrote: Dear List I would like to ask you something concenting a better print of the R output: I have a bit data frame which has

[R] how can achive step by step execution of the script

-- View this message in context: http://r.789695.n4.nabble.com/how-can-achive-step-by-step-execution-of-the-script-tp2283207p2283207.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] How can i draw a graph with high and low data points

Hi Tal, Thanks for your help. I've had a look at the site, and what i wanted to do was to plot X and Y where X is a characters and Y is numeric. The problem I'm having now is that the X axis isn't characters but just numbers from 1 onwards and when i plot it, the data i have is in descending

[R] installing packages over ssh without X forwarding

Hi, Is it possible to install packages without the testing if installed package can be loaded? I need to install bunch of packages on multiple computers over ssh. Some packages witch interact with X11 display cannot be installed in this way. for example after:

Re: [R] distributing a value for a given month across the number of weeks in that month

On Fri, Jul 9, 2010 at 9:35 AM, Dimitri Liakhovitski dimitri.liakhovit...@gmail.com wrote: Hello! Any hint would be greatly appreciated. I have a data frame that contains (a) monthly dates and (b) a value that corresponds to each month - see the data frame monthly below:

[R] ttrTests issue in cReturns

cr-cReturns(spData,ttr=MACD) Error in ind[t - k] - pos[t - k + 1] - pos[t - k] : replacement has length zero cr-cReturns(spData,ttr=macd4) Why is the above error coming? macd4 works alright ( which is a custom function built by the ttrTests author) while MACD doesnt work. Thanks -- 'Raghu'

[R] how can achive step by step execution of the script

Hi R Experts, I have certain code ,i want to achive interactive execution . For ex: 1. as part of input ,it should ask file name or table name as input. 2.in script so many graphs i need to draw,it should wait till certain key is pressed . 3:i am using windows R,rscript scriptname is not

Re: [R] Function on columns of a dataframe

you are right. But maybe aggregate is close to the desired result? aggregate(bla, list(bla$cat), max) Am 09.07.2010 16:01, schrieb David Winsemius: On Jul 9, 2010, at 9:46 AM, Eik Vettorazzi wrote: Hi Nils, have a look at ?tapply hth. Perhaps this will be part way there (I couldn't

Re: [R] strange floor rounding

Thanks everybody for referring me to FAQ 7.31 but I don't see how to solve it. I am giving a concrete number and I need to get 58 not 57. Seems, there is no way? On Fri, Jul 9, 2010 at 4:05 PM, David Winsemius dwinsem...@comcast.netwrote: On Jul 9, 2010, at 9:46 AM, Trafim Vanishek wrote:

[R] print.trellis draw.in

I am attempting to plot a trellis object on a grid. vplayout = viewport(layout.pos.row=x, layout.pos.col=y) grid.newpage() pushViewport(viewport(layout=grid.layout(2,2))) g1 = ggplot() ... g2 = ggplot() ... g3 = ggplot() ... p = xyplot() ... # works as expected print(g1, vp=vplayout(1,1))

Re: [R] installing packages over ssh without X forwarding

On 09/07/2010 7:37 AM, p...@orbit.umbr.cas.cz wrote: Hi, Is it possible to install packages without the testing if installed package can be loaded? I need to install bunch of packages on multiple computers over ssh. Some packages witch interact with X11 display cannot be installed in this way.

Re: [R] strange floor rounding

On Jul 9, 2010, at 10:27 AM, Trafim Vanishek wrote: Thanks everybody for referring me to FAQ 7.31 but I don't see how to solve it. I am giving a concrete number and I need to get 58 not 57. Seems, there is no way? Building on an example on the help page for as.integer, this seems to

Re: [R] Current script name from R

On 09/07/10 12:18, Ralf B wrote: I am using RGUI, the command line or the StatET Eclipse environment. Should this not all be the same? No, there is no particular reason why they should. Allan Ralf On Fri, Jul 9, 2010 at 7:11 AM, Allan Engelhardtall...@cybaea.com wrote: I'm

Re: [R] strange floor rounding

On 2010-07-09 8:27, Trafim Vanishek wrote: Thanks everybody for referring me to FAQ 7.31 but I don't see how to solve it. I am giving a concrete number and I need to get 58 not 57. Seems, there is no way? Sure there is: round(100*0.58). -Peter Ehlers On Fri, Jul 9, 2010 at 4:05 PM, David

Re: [R] strange floor rounding

On 09/07/2010 10:27 AM, Trafim Vanishek wrote: Thanks everybody for referring me to FAQ 7.31 but I don't see how to solve it. I am giving a concrete number and I need to get 58 not 57. Seems, there is no way? You aren't giving a concrete number. You're giving a string of three characters,

Re: [R] Function on columns of a dataframe

On Jul 9, 2010, at 10:26 AM, Eik Vettorazzi wrote: you are right. But maybe aggregate is close to the desired result? aggregate(bla, list(bla$cat), max) Right. I couldn't get it to work until I removed the first two columns: aggregate(bla[,-(1:2)], list(bla$cat), max) Then I got pretty

[R] print.trellis draw.in - plaintext (gmail mishap)

I am attempting to plot a trellis object on a grid. vplayout = viewport(layout.pos.row=x, layout.pos.col=y) grid.newpage() pushViewport(viewport(layout=grid.layout(2,2))) g1 = ggplot() ... g2 = ggplot() ... g3 = ggplot() ... p = xyplot() ... # works as expected print(g1, vp=vplayout(1,1))

Re: [R] distributing a value for a given month across the number of weeks in that month

Wow, Gabor - that's amazing - thank you so much! Dimitri On Fri, Jul 9, 2010 at 10:22 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Fri, Jul 9, 2010 at 9:35 AM, Dimitri Liakhovitski dimitri.liakhovit...@gmail.com wrote: Hello! Any hint would be greatly appreciated. I have a data

Re: [R] Data Frame Manipulation using function

Hi, Thanks a lot. The Vectorize method worked and its much faster than looping through the data frame. Regards, Harsh Yadav On Thu, Jul 8, 2010 at 11:06 PM, David Winsemius dwinsem...@comcast.netwrote: On Jul 8, 2010, at 10:33 PM, Erik Iverson wrote: I have a data frame: id url

Re: [R] random sample from arrays

Could you elaborate? Both x - 1:4 set - matrix(nrow = 50, ncol = 11) for(i in c(1:11)){ set[,i] -sample(x,50) print(c(i,-, set), quote = FALSE) } and x - 1:4 set - matrix(nrow = 50, ncol = 11) for(i in c(1:50)){ set[i,] -sample(x,11) print(c(i,-, set), quote =

Re: [R] Function on columns of a dataframe

just to satisfy my curiousity, aggregate(bla, list(bla$cat), max) works for me and resulted in Group.1 x catv1v2v3v4 1cat1 5 cat1 0.6337076 0.2887081 0.3629962 0.5328683 2cat2 10 cat2 0.5519426 0.6076447 0.4593770 0.9632341 3cat3 11 cat3 0.6094089

[R] print.xtable suppress my row.names

Dear list, someone knows why the print.xtable doesnt print row.names? I dident do anything with the options.may depends on the size of my table??? This is my code: \documentclass[a4paper]{article} \title{SCHEMA DI BILANCIO PER SINGOLE AZIENDE} \begin{document} \maketitle \hline echo=F=

Re: [R] Plotting text in existing plot?

Original poster wanted a simple way to do it, but when R has three graphics systems, four OO systems, and a zillion helpful people there's never a simple way :) -- Rather, I'd say it has a zillion simple ways. :) Bert Barry -- blog: http://geospaced.blogspot.com/ web:

Re: [R] Data Frame Manipulation using function

Really? I don't usually think of Vectorize as a performance enhancement, probably because my use of with a complex function then gets applied to 4.5 million records. I need to go out, get a cup of coffee, and leave it alone for about half an hour. I tried recently to figure out how I can

Re: [R] how can achive step by step execution of the script

Sounds like you want devAskNewPage(TRUE) or the related options(device.ask.default). See help(devAskNewPage, package=grDevices). Hope this helps. Allan On 09/07/10 13:54, vijaysheegi wrote: Hi R Experts, I have certain code ,i want to achive interactive execution . For ex: 1. as part of

Re: [R] Plotting text in existing plot?

Possible fortune. :) Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com (English)

Re: [R] How can i draw a graph with high and low data points

Hi Nathaniel , Could you give us a simple example of your data using the ?dput Function? Basically you might want to draw the axis yourself, and connect the lines is possible through using points(..., type = l) But I'd rather try and answer this with simple example data to be sure I understand

[R] interpretation of svm models with the e1071 package

Dear all, after having calibrated a svm model through the svm() command of the e1071 package, is there a way to i) represent the modeled relationships between the y and X variables (response variable vs. predictors)? ii) rank the influence of the predictors used in the model? Right now I am

Re: [R] how can achive step by step execution of the script

No - devAskNewPage is not what the OP asked for. 1. The following are, but probably only work when R is in interactive mode (?interactive), e.g. in the GUI: 2. ?winDialog ?file.choose ?choose.files ?select.list ?readline and friends and for keyboard: ?getGraphicsEvent Bert Gunter

Re: [R] how to plot two histograms overlapped in the same plane coordinate

I'm not sure what you are trying to do. Do you want one histogram for males and one for females on the same graph? If so, the simplest way to put two histograms together is to simply use the add parameter: age.males=age[which(sex==M)] age.females=age[which(sex==F)] hist(age.males,

Re: [R] how to plot two histograms overlapped in the same plane coordinate

Empirical CDFs are much better for this purpose, and allow superpositioning (see e.g. the Ecdf function in the Hmisc package). Otherwise look at histbackback in Hmisc. Frank On 07/09/2010 11:40 AM, Andrew Miles wrote: I'm not sure what you are trying to do. Do you want one histogram for

Re: [R] how to plot two histograms overlapped in the same plane coordinate

On Fri, Jul 9, 2010 at 2:29 PM, Mao Jianfeng jianfeng@gmail.com wrote: Dear R-help listers, I am new. I just want to get helps on how to plot two histograms overlapped in the same plane coordinate. What I did is very ugly. Could you please help me to improve it? I want to got a plot with

Re: [R] Query about using timestamps returned by SQL as 'factor' forsplit

Hi Ted, Well since you mentioned data.table (!) ... If risk_input is a data.table consisting of 3 columns (m_id, sale_date, return_date) where the dates are of class IDate (recently added to data.table by Tom) then try : risk_input[, fitdistr(return_date-sale_date,normal), by=list(m_id,

Re: [R] Calling Gnuplot from R

There is a basic interface between R and gnuplot in the TeachingDemos package, see ?gp.open Not much interest has been shown in this, so it is still pretty alpha level, but you can send your R data to gnuplot and have it create a basic plot. -- Gregory (Greg) L. Snow Ph.D. Statistical Data

[R] Ctree Question

Hello, I've been using ctree and have developed a 55 node - 28 terminal solution. As can be imagined, the plot is difficult to travel down each of the major branches. I've read the help files for ctree I saw where terminal nodes can be color coded. plot(airct, type = simple) plot(airct,

Re: [R] Plotting text in existing plot?

Others pointed you to the text function for the base graphics system. But if what you want to do is use text, but have a simple way of specifying the center of the plot without computing the user coordinates by hand of the center, then look at the grconvertX and grconvertY functions. --

Re: [R] how to plot two histograms overlapped in the same planecoordinate

Don't do this. The overlapping will confuse. Plot them in a lattice display with one group above the other on the same horizontal scale. See ?histogram. Bert Gunter Genentech Nonclinical Biostatistics -Original Message- From: r-help-boun...@r-project.org

[R] nls error regarding numerics vs logicals

I am trying to perform an nls for a valid negative exponential function: zz=nls(y~constant+a.est*2.7183^(b.est*x),start=list(constant=4.0,a.est=-4,b.est = -.005),trace=T) and am getting a number of different error messages, the most problematic of which is Error in nls(ring.area ~ constant +

Re: [R] R crashes with large vectors

Hi Jeremie, Maybe you can take a look at the bigmemory package. If you have multi core or have access to clusters, you may want to use any parallel computing strategy. For plotting of large data, if you are using basic R graphics, first try to use line instead of using 'point', if this still

Re: [R] nls error regarding numerics vs logicals

On 09/07/2010 1:51 PM, Jim Bouldin wrote: I am trying to perform an nls for a valid negative exponential function: zz=nls(y~constant+a.est*2.7183^(b.est*x),start=list(constant=4.0,a.est=-4,b.est = -.005),trace=T) and am getting a number of different error messages, the most problematic of

[R] select columns from vector of column names

Hi I want to extract columns from a data frame using a vector with the desired column names. This short example uses the select argument in the subset function to accomplish what I am trying to do. Is there a better solution? #names of desired columns colnames - c(col1,col3) #my data data -

Re: [R] Calling Gnuplot from R

Dear Greg, I keep on being amazed at the abundance of functions you have packed into the TeachingDemos package - thank you! Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com

[R] select columns from dataframe

Hi, I would like to extract columns from a dataframe using a vector of desired column names. The following working example uses the select argument in the subset function to accomplish what I am trying to do. Is there a better solution? Thanks. #my data data -

Re: [R] Ctree Question

Steve, I'm not sure if your task could be accomplished with a ready-made function in party. But, if you could manage to convert your tree structure to a dendrogram, then it's straightforward using dendrapply. In fact, there is an example in dendrapply help page showing how leaves are colored.

Re: [R] select columns from vector of column names

How about data[,colnames(data)%in%colnames] -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Jonathan Flowers Sent: Friday, July 09, 2010 11:27 AM To: r-help@r-project.org Subject: [R] select columns from vector of column names

[R] Not nice behaviour of nlminb (windows 32 bit, version 2.11.1)

nlminb( obj = function(x) x, start=1, lower=-Inf, upper=Inf ) $par [1] 0 $objective [1] 0 $convergence [1] 0 $message [1] absolute function convergence (6) $iterations [1] 1 $evaluations function gradient 22 [[alternative HTML version deleted]]

Re: [R] split with list

n.via...@libero.it n.via...@libero.it writes: Hi, Dear List I would like to ask you something concenting a better print of the R output: I have a bit data frame which has the following structure: CFISCALE RAGSOCBANNO VAR1VAR2. 9853312

[R] eval and assign in loop problem

deaR useRs, I am trying to assign different values to different objects in a for loop. The following is a toy example of the part that has been giving me a hard time. The first for loop generates four objects, b0, b1, b2, b3 with random numbers. And, the second for loop is equivalent to b1 = b0

[R] July BaselR Meeting

BaselR meeting July 2010 Thank you to everyone who supported our inaugural BaselR meeting on April 28th; we were fortunate to have three excellent presentations which prompted lively discussion - our thanks go to Andreas Krause, Yann Abraham and Charles Roosen for their presentations, details

[R] syntax for highlighting table rows and columns using R plugin for spss

I am using the R plug-in for spss 18, and would like to know if there is any R code that will highlight and bold rows (particularly the Totals row) within an spss table. The current option that spss has is to use a python plug-in which doesn't seem to work on my windows 7 machine. Chris

[R] Compress string memCompress/Decompress

Hello, I would like to compress a long string (character vector), store the compressed string in the text field of a SQLite database (using RSQLite), and then load the text back into memory and decompress it back into the the original string. My character vector can be compressed considerably

[R] Help finding valid permutations

So I have an array A of length n with multiple attributes and for a selected attribute y I need to list all valid permutations where a valid permutation is of the form: A[1,y] != A[n,y] A[i,y] != A[i+1,y] I've tried using the 'combinat' package, but with the vector lengths I'm using the permn

  1   2   >