Re: [R] How to write text in bar plot in R?

2012-06-13 Thread manish gupta
This does not meet my need as i need to write text as well. On Wed, Jun 13, 2012 at 9:41 AM, Manish Gupta mandecent.gu...@gmail.comwrote: Hi, I am working on R plot but i need to write some text for the bars where the value is = 0 counts 3 4 5 0 12 * 0* 4 1 3 *0* 1

Re: [R] replacing NA for zero

2012-06-13 Thread Rui Barradas
Hello, There's a shorter one, using merge(). DF2 - merge(data.frame(x=1:8), DF1, all.x = TRUE) DF2$y[is.na(DF2$y)] - 0 Rui Barradas Em 13-06-2012 02:02, arun escreveu: Hi, I got confused by your heading when I sent the previous reply. There is another way to reach the same conclusion.

Re: [R] How to build a large identity matrix faster?

2012-06-13 Thread Ceci Tam
that makes sense, thanks for the reply On 13 June 2012 05:53, Uwe Ligges lig...@statistik.tu-dortmund.de wrote: On 11.06.2012 11:45, Ceci Tam wrote: diag(n) is alright when n = 5e3, it took 0.7 sec on my machine for diag(5e3). However, it's slow when n = 23000, diag(23000) took 15 sec

Re: [R] How to write text in bar plot in R?

2012-06-13 Thread Manish Gupta
Here Text message is kind of error message for no data available. regards -- View this message in context: http://r.789695.n4.nabble.com/How-to-write-text-in-bar-plot-in-R-tp4633190p4633211.html Sent from the R help mailing list archive at Nabble.com.

[R] add horizontal reference lines in lattice bwplot

2012-06-13 Thread Gareth.Williams
I'm plotting a lattice bwplot from a data.frame and would like to add some horizontal lines indicating some reference values (preferably with text annotation) as well as a grid of regular log tics. The following illustrates the plot: data(ToothGrowth) require(lattice) p = bwplot(exp(len) ~

Re: [R] adjust space between horizontal legend text in a barplot

2012-06-13 Thread arun
Hi, Sorry about that.  My mistake. Try this:  barplot(VADeaths, names=rep(, 4), col=gray(1:3/4+0.01)) legend(bottom, legend=expression(scriptstyle(Strongly disagree/disagree),scriptstyle(Neutral),scriptstyle(Strongly agree/agree)), fill=gray(1:3/4+0.01), inset=-0.2, bty='n', horiz=T,

Re: [R] How to write text in bar plot in R?

2012-06-13 Thread arun
Hi, Try the format from this example: mp - barplot(VADeaths)  tot - colSums(VADeaths)     text(mp, tot+6, format(tot), xpd = TRUE, col = red)  text(mp, c(A,B,C,D), xpd = TRUE, col = blue) I hope it helps. A.K. - Original Message - From: Manish Gupta mandecent.gu...@gmail.com To:

Re: [R] add horizontal reference lines in lattice bwplot

2012-06-13 Thread PtitBleu
Hello Gareth, Concerning the lines, you can try : panel.abline(h=13) and panel.grid(h=2, v=0, y=c(11, 12), col.line=blue) I have no idea for the annotation. Sorry. Ptit Bleu. -- View this message in context:

Re: [R] Error in sort(abs(diff(genomdat)))[1:n.keep] : only 0's may be mixed with negative subscripts

2012-06-13 Thread sanshine
Ok, I`m going to check the n.keep values. Thanks! Sander From: Michael Weylandt [via R] [mailto:ml-node+s789695n463311...@n4.nabble.com] Sent: dinsdag 12 juni 2012 17:08 To: Zeeuw, S.A.J. van der (PATH) Subject: Re: Error in sort(abs(diff(genomdat)))[1:n.keep] :

[R] lme: extract result-function

2012-06-13 Thread Christof Kluß
Hi, mod - lme(A ~ -1 + B+C+D+E+F+G, random = ~1 | ...) results in summary(mod)$coeff B C D E F G (Intercept) b c d e f g i Now I'm interested in the function f - function(B,C,D,E,F,G) - { return(i + b*B + c*C + d*D + e*E + f*F + g*G) } Is there a easier way to create such function with

Re: [R] lme: extract result-function

2012-06-13 Thread Rui Barradas
Hello, Try f - function(response, regressors) as.formula(paste(response, paste(regressors, collapse= + ), sep= ~ )) (resp - A) (regr - c(-1, LETTERS[2:7])) fmla - f(resp, regr) Hope this helps, Rui Barradas Em 13-06-2012 09:21, Christof Kluß escreveu: Hi, mod - lme(A ~ -1 +

Re: [R] lme: extract result-function

2012-06-13 Thread ONKELINX, Thierry
Dear Christof, You want the predict() function. See ?predict.lme for the details. Best regards, Thierry PS Questions on lme() can be asked at r-sig-mixed models. ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie

Re: [R] lme: extract result-function

2012-06-13 Thread Dieter Menne
Christof Kluß wrote mod - lme(A ~ -1 + B+C+D+E+F+G, random = ~1 | ...) ... f - function(B,C,D,E,F,G) - { return(i + b*B + c*C + d*D + e*E + f*F + g*G) } It looks like you are trying to compute contrasts the ugly way. Check estimable in package gmodels or the vignette of the multcomp

Re: [R] lme: extract result-function

2012-06-13 Thread Pascal Oettli
Hello, Is it to re-create what predict is able to do? Best Regards, Pascal Le 12/06/13 17:21, Christof Kluß a écrit : Hi, mod- lme(A ~ -1 + B+C+D+E+F+G, random = ~1 | ...) results in summary(mod)$coeff B C D E F G (Intercept) b c d e f g i Now I'm interested in the function f-

Re: [R] GEE with Inverse Probability Weights

2012-06-13 Thread Thomas Lumley
On Wed, Jun 13, 2012 at 9:25 AM, RFrank spark...@gmail.com wrote: Greetings, I have a very, very, simple research question.  I want to predict one dichotomous variable using another dichotomous variable.  Straightforward, right?  The issue is that the dataset has two issues causing some

Re: [R] Rotating characters in text

2012-06-13 Thread Jim Lemon
On 06/12/2012 11:49 PM, Stuart Rosen wrote: For labelling a plot, I am trying to rotate a character string using text() so that characters are upright and reading down, for example, ... L i k e t h i s . It appears that par crt does not work with text. Does anyone have any other suggestions.

Re: [R] How to Group Categorical data in R?

2012-06-13 Thread Manish Gupta
If there are some unique value in each column. How to look up in R? -- View this message in context: http://r.789695.n4.nabble.com/How-to-Group-Categorical-data-in-R-tp4477622p4633235.html Sent from the R help mailing list archive at Nabble.com. __

[R] finding duplicates in a data frame

2012-06-13 Thread sathya7priya
I have two data frames which has 3 columns each.My first data frame is large like this below new.col ppm.p. freq.p. 1_3_diaminopropane 3.13859 5.67516 1_3_diaminopropane 3.137 6.65388 1_3_diaminopropane 3.13541 8.0142 1_3_diaminopropane 3.13383 9.64184 1_3_diaminopropane 3.12075 298.243

Re: [R] problem about set operation and computation after split

2012-06-13 Thread Rui Barradas
Hello, It seems you haven't discovered the *apply functions yet. In this case, since you have a list, the function would be lapply. 1. Reduce(intersect, lapply(mySplit, function(x) x$product)) 2. lapply(mySplit, function(x) within(x, relative - diff(c(value[1], value 3. lapply(mySplit,

[R] need help

2012-06-13 Thread shilpa rai
hello could you help in solving the following problem I want to replace same consecutive words by a single word in a sentence.. for example --- my name name name is micky so I want the output like this--my name is micky I want this solution for a text file can you tell me the code for it??

[R] Building normal qq plot

2012-06-13 Thread juliane0212
I have tried to build a normal qq plot on my own, but it does not match with the implemented qqnorm in R. What am I missing? As You can see, the points do not match. qqnorm(c(-2.8,qnorm(seq(0,1,by=(1/200))[2:199],0,1),2.8)) points(qnorm(seq(0,1,by=(1/200))[2:199],0,1),

[R] r-java for simulation

2012-06-13 Thread Mohan Radhakrishnan
Hi, I came across the integration package that lets a Java program interact with R. Does any one use this model for simulation ? I am plannning to get random variates generated by R which could be the number of bytes according to various distributions. My client load generator can

Re: [R] add horizontal reference lines in lattice bwplot

2012-06-13 Thread Gareth.Williams
-Original Message- From: PtitBleu [mailto:ptit_b...@yahoo.fr] Sent: Wednesday, 13 June 2012 5:46 PM To: r-help@r-project.org Subject: Re: [R] add horizontal reference lines in lattice bwplot Hello Gareth, Concerning the lines, you can try : panel.abline(h=13) and

Re: [R] finding duplicates in a data frame

2012-06-13 Thread Rui Barradas
Hello, Try merge(df1, df2, by.x=c(ppm.p., freq.p.), by.y=c(ppm.p., freq.p.)) ppm.p. freq.p. new.col.x new.col.y 1 4.81412 105.1100 1_amino_1_phenylmethyl_phosphonic_acid unknown 2 7.44687 7.1684 1_amino_1_phenylmethyl_phosphonic_acid unknown where 'df1'

Re: [R] finding duplicates in a data frame

2012-06-13 Thread Rui Barradas
What are the names of your data.frames columns? Paste the output of the following commands in a post: dput(head(df1, 30)) dput(head(df2, 30)) Tui Barradas Em 13-06-2012 12:55, Sathya Priya escreveu: I got this error while trying your code.. Can you please sort it out Error in fix.by(by.x,

Re: [R] finding duplicates in a data frame

2012-06-13 Thread Petr Savicky
On Wed, Jun 13, 2012 at 03:16:57AM -0700, sathya7priya wrote: I have two data frames which has 3 columns each.My first data frame is large like this below new.col ppm.p. freq.p. 1_3_diaminopropane 3.13859 5.67516 1_3_diaminopropane 3.137 6.65388 1_3_diaminopropane 3.13541 8.0142

Re: [R] How to Group Categorical data in R?

2012-06-13 Thread John Kane
This is not really enough information. PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Have you looked at ?unique John Kane Kingston ON Canada -Original Message- From:

[R] reading xls files using read.xls and xlsReadWrite

2012-06-13 Thread C W
Dear R list, I am trying to read in .xls files. I have tried using package 'gdata', I get the following error. dataset - read.xls(June.xls) Wide character in print at /Library/Frameworks/R.framework/Versions/2.14/Resources/library/gdata/perl/ xls2csv.pl line 262. Seems like this is a problem

Re: [R] Building normal qq plot

2012-06-13 Thread S Ellison
I have tried to build a normal qq plot on my own, but it does not match with the implemented qqnorm in R. The main issue is, that I would like to change the backgorund color auf my plot ... One thing you could try is setting the background colour using bar(bg=lightyellow) before calling

[R] regarding the function WebCorpus

2012-06-13 Thread raishilpa
hello, I have a question regarding the function WebCorpus..When I am using this function its showing error: couldnt find the function...I have downloaded the tm package...Can you help me in this regard? Also I want to know how we can data/information about a product or topic from google,facebook

[R] help

2012-06-13 Thread Natalia Martins
I'd like to check the adherence of some data will lognormal distribution, using ks.test (x, plognormal) However, not find the name of this lognormal function. Thank you -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e

Re: [R] templated use of aggregate

2012-06-13 Thread John Kane
I think the column is wrong in adf - aggregate(.~mm[,1], data=mm, sum) Try 2 rather than 1 but this is just a guess since the sample data is not easily useable.. Please supply your sample data using dput(() See ?dput for information. It makes llife much easier for the reader. John Kane

Re: [R] help

2012-06-13 Thread David L Carlson
?Lognormal -- David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77843-4352 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- project.org] On Behalf Of Natalia Martins

Re: [R] Building normal qq plot

2012-06-13 Thread juliane0212
Thank you very much. It was quite easy to adjust the code to get the background colour I have always used. I only tried to get the code using list(qqnorm). Thanks -- View this message in context: http://r.789695.n4.nabble.com/Building-normal-qq-plot-tp4633242p4633257.html Sent from the R

Re: [R] templated use of aggregate

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 9:38 AM, Matthew Johnson wrote: Sorry, i'll try and put more flesh on the bones. please note, i changed the data in the example, as fiddling has raised another question that's best illustrated with a slightly different data set. first of all, when i do as you suggest, i

[R] histogram fill lattice

2012-06-13 Thread Powell, Jeff
Dear all, I would like to change the fill pattern of a histogram using histogram() in the lattice package. I know how to do so using hist(), but would prefer to stay within lattice. dt1 - rnorm(100,0,1) hist(dt1, density=3, angle=45) library(lattice) histogram(dt1, xlab =

[R] separate the sentence after finding a particular word

2012-06-13 Thread raishilpa
hello, I want to know ..how we can separate the sentence after finding a particular word... for example I love to watch movies of Hollywood but should not be romantic...I want to join you school but due to bad financial condition I cant.. I want output in following format I love to watch movies

Re: [R] templated use of aggregate

2012-06-13 Thread David Winsemius
On Jun 12, 2012, at 11:32 PM, Matthew Johnson wrote: Dear R-help, I have an xts data set that i have subset by date. now it contains a date-time-stamp, and two columns (price and volume traded): my objective is to create tables of volume traded at a price - and i've been successfully

Re: [R] regarding the function WebCorpus

2012-06-13 Thread John Kane
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Kane Kingston ON Canada -Original Message- From: raishilpa@gmail.com Sent: Wed, 13 Jun 2012 05:11:29 -0700 (PDT) To:

Re: [R] Rotating characters in text

2012-06-13 Thread David L Carlson
Or for the y label and using gsub: plot(1:5, ylab=) text(.35,3, gsub((.), \\1\n, Like this), xpd=TRUE) Sticking it in ylab= does not work. -- David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77843-4352

Re: [R] templated use of aggregate

2012-06-13 Thread Matthew Johnson
Sorry, i'll try and put more flesh on the bones. please note, i changed the data in the example, as fiddling has raised another question that's best illustrated with a slightly different data set. first of all, when i do as you suggest, i obtain the following error: PxMat - aggregate(mm[,-1] ~

Re: [R] templated use of aggregate

2012-06-13 Thread Matthew Johnson
Sorry about the cross posting - i didn't realise it was bad etiquette. my sessioninfo was as follows: sessionInfo() R version 2.14.1 (2011-12-22) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8 attached base

Re: [R] templated use of aggregate

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 10:09 AM, Matthew Johnson wrote: my sessioninfo was as follows: sessionInfo() R version 2.14.1 (2011-12-22) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8 attached base packages: [1] stats

Re: [R] templated use of aggregate

2012-06-13 Thread Matthew Johnson
thank you for your patience. i assure you i will get better with the appropriate etiquette - and hopefully eventually contribute. On 13 June 2012 16:18, David Winsemius dwinsem...@comcast.net wrote: On Jun 13, 2012, at 10:09 AM, Matthew Johnson wrote: my sessioninfo was as follows:

Re: [R] reading xls files using read.xls and xlsReadWrite

2012-06-13 Thread Ivan Calandra
I don't know much about gdata, but I can tell you that there is still no version of xlsReadWrite on Mac (unfortunately). xlsx works well, though, on both Mac and Windows. HTH, Ivan -- Ivan CALANDRA Université de Bourgogne UMR CNRS/uB 6282 Biogéosciences 6 Boulevard Gabriel 21000 Dijon, FRANCE

[R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread Kristi Glover
Hi R experts, Could you please help me to fit a linear, cubic and quadratic curve in a figure? I was trying to show all these three fitting curves with different colour in one figure. I spent substantial time to figure it out, but I could not. I have given here a example and what I did for

Re: [R] histogram fill lattice

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 9:55 AM, Powell, Jeff wrote: Dear all, I would like to change the fill pattern of a histogram using histogram() in the lattice package. I know how to do so using hist(), but would prefer to stay within lattice. dt1 - rnorm(100,0,1) hist(dt1, density=3, angle=45)

Re: [R] Indexing Grouped Data

2012-06-13 Thread Peter Maclean
I need help in indexing grouped data. In this excample (df1 data), the first child had a first immunization at age 2. The second child had the first, second and third immunization at age 5,10, and 12, the third child had first and second immunization at age 4 and 6 and the fourth child had the

[R] Splitting Large Data Frame into Two

2012-06-13 Thread Joshua Budman
Hi I have a large data frame of the form: a 1 b 2 c 3 And I would like to split this data frame into two separate data frames based on the values in the first column, e.g. a 1 b 2 and c 3 Is there any way of doing this without having to write a different which statement for each

Re: [R] separate the sentence after finding a particular word

2012-06-13 Thread Rui Barradas
Hello, Sorry for my earlier reply, it was meant for another question you posted. In what follows, 'x' is your text. f - function(x, pattern) unlist(strsplit(x, pattern)) x - I love to watch movies of Hollywood but should not be romantic...I want to join you school but due to bad financial

Re: [R] finding duplicates in a data frame

2012-06-13 Thread Rui Barradas
Hello, Thanks! Now we can see that your df1 only has one column, your dataset was messed up in some previous step. So, it needs to be transformed into a 3 columns data.frame first. Then use merge(). x - levels(df1$new.col.ppm.p..freq.p.)[df1$new.col.ppm.p..freq.p.] x - t(sapply(x,

Re: [R] separate the sentence after finding a particular word

2012-06-13 Thread Rui Barradas
Hello, Try x - scan(what=character, text=my name name name is micky) rle(x) paste(rle(x)$values, collapse= ) Hope this helps, Rui Barradas Em 13-06-2012 13:06, raishilpa escreveu: hello, I want to know ..how we can separate the sentence after finding a particular word... for example I love

[R] asign variables in a for loop

2012-06-13 Thread Nympha Nymphaea
Dear R-helpers, I'm stuck with a little problem that surely has an easy solution but I can't think of a way to solve it. I'd really appreciate any help you can offer me! I'll provide a small example. Given a dataframe data.txt that looks like this: IDfreqVarVar_meanRatio_mean

Re: [R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread Kenneth Frost
Hi, Kristi- Here is an option. ?curve...to add the polynomials to your graph ?textfor adding the R2 to you plot Ken On 06/13/12, Kristi Glover wrote: Hi R experts, Could you please help me to fit a linear, cubic and quadratic curve in a figure? I was trying to show all these three

Re: [R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread Robert Baer
dput(test) structure(list(sp = c(4L, 5L, 9L, 12L, 14L), env = c(12L, 18L, 20L, 17L, 15L)), .Names = c(sp, env), class = data.frame, row.names = c(NA, -5L)) plot(test$sp~test$env, main = S vs. temp, xlim=c(0,20), ylim=c(0,14), ylab=S,xlab=env) linear-lm(test$sp~test$env)

[R] Plotted circle does not go through desired points - very long email with code

2012-06-13 Thread Monica Pisica
Hi,   I am trying to solve a problem related to Poincare circles ( for more info http://www.ms.uky.edu/~droyster/courses/spring08/math6118/Classnotes/Chapter09.pdf). In a nutshell i am trying to replicate the method in the above pdf section 9.2.1. that explains in broad terms how to draw the

Re: [R] Splitting Large Data Frame into Two

2012-06-13 Thread Sarah Goslee
How are you deciding which values in the first column go into which subset? If you have a vector containing those values, you could use %in% or if they're determined logically you could use that criterion. A reproducible example and a bit more information would get you more concrete answers.

Re: [R] Indexing Grouped Data

2012-06-13 Thread William Dunlap
df1$ind - ave(integer(nrow(df1)), df1$id, FUN=seq_along) There are faster ways to do this if you know that id is sorted. Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of

[R] logistic Regression with SSlogis but y == 0 ?!

2012-06-13 Thread pluto
Hello there! I got some data with x and y values. there are some y == 0. This is a problem for the selfstarting regression model SSlogis. The regression works if I use a non selfstarting model. The formula is the same. But this needs very detailed information of the starting list. I dont have

[R] phyloclim help

2012-06-13 Thread Brent Hendrixson
Hi all – I am a newbie to R and have had a terrible time trying to figure out how to set up an analysis.  I have installed and loaded the appropriate package (phyloclim), but I do not even know where to begin defining one of the arguments. The usage for the command I’m trying to perform is:

[R] Median line with stripchart

2012-06-13 Thread Luigi
Dear all, I would like to ask if it possible to draw a median line in a Stripchart, either using the built-in STRIPCHART function or with another package. I was expecting stripchart() to work more or less like boxplot(), but I couldn't find an easy way to draw a mean or median line, neither I

Re: [R] asign variables in a for loop

2012-06-13 Thread Sarah Goslee
Part of the answer is that you're overwriting your loop variable x and the variable Mean in each loop iteration. I have no idea what you think with(data, c()) is doing, either, and there's at least one comma missing. Also, data is the name of a function. That said, apply(mydata[ ,

Re: [R] Splitting Large Data Frame into Two

2012-06-13 Thread Joshua Budman
This is a sample of the data: ID chrom loc.start loc.end 1MA1B1Cy5 chr1 197008581 197026781 18 MA1B3Cy5 chr1 197079541 197080381 55 MA5B2Cy5 chr1 197088651 197118071 70 MA5B2Cy5 chr1 197172341 197189641 72 MA5B3Cy5 chr1 197008581 197010601 89 MA5B4Cy5 chr1

Re: [R] Splitting Large Data Frame into Two

2012-06-13 Thread Sarah Goslee
Then: Set1.names - c(MA1B1Cy5, MA5B4Cy5) Set1.data - mydata[mydata$ID %in% Set1.names, ] On Wed, Jun 13, 2012 at 12:15 PM, Joshua Budman josh.bud...@gmail.com wrote: This is a sample of the data:            ID chrom loc.start   loc.end 1    MA1B1Cy5  chr1 197008581 197026781 18   MA1B3Cy5 

Re: [R] Median line with stripchart

2012-06-13 Thread Uwe Ligges
On 13.06.2012 17:57, Luigi wrote: Dear all, I would like to ask if it possible to draw a median line in a Stripchart, either using the built-in STRIPCHART function or with another package. I was expecting stripchart() to work more or less like boxplot(), but I couldn't find an easy way to

Re: [R] Median line with stripchart

2012-06-13 Thread David L Carlson
It would help if you would give us a copy of at least part of the data you are using and show us what you have tried so far. This is a simple example. In this case the median is for the entire data set not for a single group. But you really didn't give us enough information. x - rnorm(75) g -

[R] Re: Not able to write to PostgreSQL database using dbWriteTable

2012-06-13 Thread Prakash Thomas
Dear R User's, Thank you, Mark. The following code suggested by you worked for me. dbGetQuery(connAE1, sprintf(INSERT INTO test1 (id) VALUES ( %d );, i)) But I have a issue in passing date and time data as variable.If I hard code the value like bellow it workings.

[R] Determining Legend for smoothScatter

2012-06-13 Thread Kurinji Pandiyan
Dear all, I am using the smoothScatter function in base R for a plot - Lab.palette.both - colorRampPalette(c(darkblue,lightblue,red,yellow), space = Lab) smoothScatter(X24fresh.sorted[,c(13,10)], colramp = Lab.palette.both) I understand that my yellow colored points represent highest density,

Re: [R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread William Dunlap
It is more direct to use predict() instead of reconstructing by hand the prediction expression from the formula given to lm(). E.g., x - seq(1,6,by=1/4) y - sin(x) + rnorm(length(x), 0, 1/4) plot(x, y) fits - lapply(1:3, function(degree)lm(y~poly(x, deg=degree))) xpred - pretty(x,

[R] How to calculate the statistcs for extracted region?

2012-06-13 Thread Jonsson
I have a binary file(a) with size of (360 720 )for the globe.I wrote the code given below to read and extract an area (south america)from that file. when I use summary for the whole file I got: summary(a, na.rm=FALSE) Min. 1st Qu. MedianMean 3rd Qu.Max.NA's 0.001.003.00

Re: [R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread Frank Harrell
See also the scripts in http://biostat.mc.vanderbilt.edu/BioMod Frank William Dunlap wrote It is more direct to use predict() instead of reconstructing by hand the prediction expression from the formula given to lm(). E.g., x - seq(1,6,by=1/4) y - sin(x) + rnorm(length(x), 0,

Re: [R] How to plot linear, cubic and quadratic fitting curve in a figure?

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 10:29 AM, Kristi Glover wrote: Hi R experts, Could you please help me to fit a linear, cubic and quadratic curve in a figure? I was trying to show all these three fitting curves with different colour in one figure. I spent substantial time to figure it out, but I

[R] gWidgets - gtable returning multiple values

2012-06-13 Thread michaelyb
Hello, I am having some trouble finding a solution for a probleIm I am facing with the package gWidgets. Here is the code: flavors-c(vanilla, chocolate, strawberry) f-function(h,...)print(svalue(h$obj)) w - gwindow(checkbox example) gp - ggroup(container=w) glabel(Favorite flavors:,cont=gp)

Re: [R] separate the sentence after finding a particular word

2012-06-13 Thread raishilpa
thank you so much for reply and this code is working... Can you tell me whether both this code will work for text file or not ..I mean if I am uploading a text file and then using the given code ,will this give output or there is need to do something else...Can you help me? thanking you in

[R] niche.width functions

2012-06-13 Thread Carmen Tattersfield
Hello, I want to calculate niche width of a species in a community using Smith's selectivity measure (Smith 1982), as described in Ecological Methodology (Krebs 1999). Is there a package/function available that can do this? I have searched and found the niche.width function in SPAA (below).

Re: [R] Splitting Large Data Frame into Two

2012-06-13 Thread arun
Hi, Try this:dat2-read.table(text= ID chrom loc.start  loc.end     MA1B1Cy5  chr1 197008581 197026781   MA1B3Cy5  chr1 197079541 197080381   MA5B2Cy5  chr1 197088651 197118071   MA5B2Cy5  chr1 197172341 197189641  MA5B3Cy5  chr1 197008581 197010601   MA5B4Cy5  chr1 197025421 197025701   MA5B4Cy5 

Re: [R] How to calculate the statistcs for extracted region?

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 12:55 PM, Jonsson wrote: I have a binary file(a) with size of (360 720 )for the globe.I wrote the code given below to read and extract an area (south america)from that file. when I use summary for the whole file I got: summary(a, na.rm=FALSE) Min. 1st Qu. Median

Re: [R] separate the sentence after finding a particular word

2012-06-13 Thread Rui Barradas
Hello, I believe it will work, but only you can say. Try it. And if it doesn't it's probably easy to adapt to a real dataset problem. Almost always the point in R-help is to understand what to do or how to do it. Rui Barradas Em 13-06-2012 17:13, raishilpa escreveu: thank you so much for

Re: [R] separate the sentence after finding a particular word

2012-06-13 Thread raishilpa
thanks a lot !!! -- View this message in context: http://r.789695.n4.nabble.com/separate-the-sentence-after-finding-a-particular-word-tp4633240p4633303.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]]

[R] side by side symmetric heatmap

2012-06-13 Thread Tanu Soni
I am trying to plot two heatmaps on same page . I am using heatmap() so that i can draw symmetrical heatmap but par() does not work with heatmap . Is there a way to plot two symmetrical heatmaps on same page. -- Tanu [[alternative HTML version deleted]]

[R] mask and reverse the selection

2012-06-13 Thread Celine
Hi I need some help with the function mask I have a polygon and I want to extract data for all the world except for the polygon I write this code but it doesn't work, do you have any idea of the problem ? how can I do that ? worldmask=mask(invHS,hotspot,inverse=TRUE)

[R] Reading several tables from stdin

2012-06-13 Thread Julio Sergio
I'm trying to write a Rscript program capable of reading several tables from the standard input. The problem is that the tables aren't in files because they are coming from another process that is generating them. From the R-console the following works pretty well: | f - stdin() | t -

Re: [R] How do I connect dots in the RGL package?

2012-06-13 Thread Maaike_020
David Winsemius wrote open3d() [1] 2 points3d(rnorm(10),rnorm(10),rnorm(10)) # First add standard axes axes3d() lines3d( x=c(-1,1), y=c(-1,1), z=c(-1,1) ,col=red) Thank you so much, it works! I have set a photo as a background by rgl.bg(sphere=TRUE,

Re: [R] side by side symmetric heatmap

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 4:37 PM, Tanu Soni wrote: I am trying to plot two heatmaps on same page . I am using heatmap() so that i can draw symmetrical heatmap but par() does not work with heatmap . Not exactly true or, at the very least, misleading. Is there a way to plot two symmetrical

Re: [R] How do I connect dots in the RGL package?

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 4:28 PM, Maaike_020 wrote: David Winsemius wrote open3d() [1] 2 points3d(rnorm(10),rnorm(10),rnorm(10)) # First add standard axes axes3d() lines3d( x=c(-1,1), y=c(-1,1), z=c(-1,1) ,col=red) Thank you so much, it works! I have set a photo as a background by

[R] phyloclim help

2012-06-13 Thread Brent Hendrixson
Hi all – My apologies if this is a redundant message.  I am a newbie to R and have had a terrible time trying to figure out how to set up an analysis. I have installed and loaded the appropriate package (phyloclim), but I'm not sure I have defined the arguments correctly. The usage for the

Re: [R] mask and reverse the selection

2012-06-13 Thread Rolf Turner
This question apparently relates to the mask() function from the raster package. This should have been made clear. You might find the function complement.owin() in the spatstat package helpful. cheers, Rolf Turner On 14/06/12 08:42, Celine wrote: Hi I need some help with the

Re: [R] Median line with stripchart

2012-06-13 Thread Luigi
Sorry for the inconvenience, I have added an EXAMPLE. Here I am plotting few data with the relative median values. The median are added ad hoc using the points() function, which is not really professional and is time consuming. With boxplot() the median values are automatically drawn, with

[R] Trouble with .bat files upon installation of 15.0

2012-06-13 Thread Sarah Henderson
Hello to all -- I'm hoping that someone more knowledgeable than me can shed some light on a problem I have been having. In point form: - I am running XP on a 64-bit processor - I run 5 automated R tasks every morning using .bat files and a little utility called System Scheduler by Splinterware.

Re: [R] Reading several tables from stdin

2012-06-13 Thread Julio Sergio
Julio Sergio juliosergio at gmail.com writes: I'm trying to write a Rscript program capable of reading several tables from the standard input. The problem is that the tables aren't in files because they ... Do you have any comments on this? Thanks, Well, maybe I'm not posting this

Re: [R] phyloclim help

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 6:29 PM, Brent Hendrixson wrote: Hi all – My apologies if this is a redundant message. I am a newbie to R and have had a terrible time trying to figure out how to set up an analysis. I have installed and loaded the appropriate package (phyloclim), but I'm not sure

Re: [R] Reading several tables from stdin

2012-06-13 Thread Jeff Newmiller
Your expectations of response time are out of line for a mailing list. The subject of your question might or might not be appropriate here depending on how portable you want the solution to be. Your dependence on operating-system-specific features (ctrl-d as eof) and expectation that such

Re: [R] phyloclim help

2012-06-13 Thread Brent Hendrixson
Hi David - Thanks for the response. On Jun 13, 2012, at 6:29 PM, Brent Hendrixson wrote: Hi all – My apologies if this is a redundant message. I am a newbie to R and have had a terrible time trying to figure out how to set up an analysis. I have installed and loaded the

Re: [R] phyloclim help

2012-06-13 Thread David Winsemius
On Jun 13, 2012, at 11:00 PM, Brent Hendrixson wrote: Hi David - Thanks for the response. On Jun 13, 2012, at 6:29 PM, Brent Hendrixson wrote: Hi all – My apologies if this is a redundant message. I am a newbie to R and have had a terrible time trying to figure out how to set up an

Re: [R] Not able to write to PostgreSQL database using dbWriteTable

2012-06-13 Thread Mark Dalphin
I still don't have any idea about your schema [eg CREATE TABLE (blah, blah, ...);], but I guess you don't have the right database type for id when you are storing a date. In PostgreSQL: CREATE TABLE myTable ( id INTEGER PRIMARY KEY, aDate TIMESTAMP ); In R: dbGetQuery(conn, paste(INSERT

[R] Help for boxplot

2012-06-13 Thread adamshun
Hi guys, I am doing a project to evaluate the 7 individual fund performance from a portfolio compared with all fund performance utilizing the same strategy. Lets say in total there are 10 strategies and in the portfolio there exist 5 strategies. First, i use the BOXPLOT() and SUBSET() to

[R] Tukey Kramer with ANOVA (glm)

2012-06-13 Thread Alaska_Man
Hello, I am performing a BACI analysis with ANOVA using the following glm: fit1-glm(log(Cucs_m+1)~(BA*Otter)+BA+Otter+ID+Primary, data=b1) The summary(aov(fit1)) shows significance in the interaction; however, now I would like to determine what combinations of BA and Otter are significantly

Re: [R] How do I connect dots in the RGL package?

2012-06-13 Thread Maaike_020
David Winsemius wrote And that is the only problem you see? The background is completely distracting the viewer from the data Yuck. This is a system for statistics. If you want to make movies go to Pixar. Wow... well ok.. -- View this message in context:

Re: [R] finding duplicates in a data frame

2012-06-13 Thread arun
Hi, Try this:  dat1-data.frame(x=letters[1:5],y1=4:8,y2=sample(1:4,5,replace=T))  dat2-data.frame(x=letters[6:10],y1=4:8,y2=sample(1:4,5,replace=T))  merge(dat1,dat2,by=c(y1,y2)) A.K. - Original Message - From: sathya7priya sathya7pr...@gmail.com To: r-help@r-project.org Cc: Sent:

Re: [R] Mapping one vector to another

2012-06-13 Thread Geoff Leyland
On 13/06/2012, at 12:35 AM, Rui Barradas wrote: Hello, ...merge... On 13/06/2012, at 12:36 AM, Jean V Adams wrote: Geoff, ...merge... Thanks Rui and Jean, merge appears to be what I was looking for, and much better than my workaround for loop.

  1   2   >