Re: [R] Hello
On Tue, 8 Mar 2022 15:55:19 +0530 Shrinivas Dharma wrote: > I am not sure if this is the right place to ask the question of my > type. > I am wirking on a social network analysis project with R igraph > software. > > My graph data has multiple edges and multiple self loops Have you tried the igraph forum? It's probably a better fit because it's more likely to have people with igraph experience: https://igraph.discourse.group/c/usage/11 -- Best regards, Ivan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
[R] Hello
Hello I am not sure if this is the right place to ask the question of my type. Nonetheless, being desparate , let me ask I am wirking on a social network analysis project with R igraph software. My graph data has multiple edges and multiple self loops Going by the instructions I used the simplify function to make it simple graph code chunk as follows >simplify(r1,remove.multiple = T,remove.loops = T, edge.attr.comb = igraph_opt("edge.attr.comb") ) >plot.igraph(r1) still I am getting loops and multiple edges. I am using the latest igraph edition Please tell me how to go about it Regards SHRINIVAS __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] Hello! Problem - function ''matord'' in R
... and if so, then you are on the wrong mailing list... On December 14, 2020 1:43:49 PM PST, "T. A. Milne via R-help" wrote: >Might the desired function "matord" be part of a Bioconductor package? >This link >https://rdrr.io/bioc/clusterSeq/src/R/associatePosteriors.R > >suggests that could be the case. > >- T. Arthur Milne > >__ >R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >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, minimal, self-contained, reproducible code. -- Sent from my phone. Please excuse my brevity. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] Hello! Problem - function ''matord'' in R
Might the desired function "matord" be part of a Bioconductor package? This link https://rdrr.io/bioc/clusterSeq/src/R/associatePosteriors.R suggests that could be the case. - T. Arthur Milne __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] Hello! Problem - function ''matord'' in R
Are you sure you're not thinking of a matlab or some other software's function? (Don't answer... but R has fairly extensive linear algebra capabilities built into its standard distro, and this doesn't sound familiar. I could well be wrong about this, of course). Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Mon, Dec 14, 2020 at 11:40 AM Nikolaos Anagnostakis < nikos.anagnosta...@hotmail.com> wrote: > Hello! > > I am developing an algorithm via the R program to create phylogenetic > trees and calculate values that interest me. On 2019 I had used a function > called ''matord'' but I can't find it anymore. Specifically I needed it for > calculation of two matrices for CI and RI. > > Is there anyone who could help me? > Is there a similar function? > > Best Wishes > > Nikos > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
[R] Hello! Problem - function ''matord'' in R
Hello! I am developing an algorithm via the R program to create phylogenetic trees and calculate values that interest me. On 2019 I had used a function called ''matord'' but I can't find it anymore. Specifically I needed it for calculation of two matrices for CI and RI. Is there anyone who could help me? Is there a similar function? Best Wishes Nikos __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
[R] hello r
Good afternoon r http://wsitest01.com/minerals.php?france=y18m1xba0w8z My best to you terry minton __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] hello i have a question on music analysis and mathematical synthesis related to r code
Search on the Internet! "Analyze music in R" had hits for several R packages that seemed like they might be relevant. Apologies if you've already done this and found mothing to meet your needs. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Sep 7, 2016 at 10:57 PM, darth brando wrote: > Apologies for the long title but it is semi specific a topic and yes I am a > noobs user to the system. I have read the guide and will attempt to adhere to > the guide in this process and I do apologize in advance if I fail to do so, > this is my first time here. > > To the point; firstly version: > > I have windows 7 64 bit OS, I'm going to be working with most current up to > date version of R code for that OS with optional plug ins as needed, I am > going to be using R code in conjunction with Fruity Loop Studio, > dosbox/visual studio, and both synthesizer plug ins for Fruity Loop Studio > and a physical soundblaster sound card in my alienware 15 2015 version as > well as possibly a few emulator or OS instances for retro hardware to > experiment with. Its a lot of software to take in yea, some I know by heart > others like R Code I am a noob at. > > Now that the toolbox of software is out of the way; my Question: > > Music is largely mathematically based, R Code is perfect for analyzing very > large sets of data, naturally I saw some potential and wondered about the > specifics of the type of application which follows. > > I wish to use R Code to analyze portions of my personal library of music in > order to discover the main sets of underlying patterns within that portion of > my music library. I then wish to utilize patterns discovered to create a > program which using proper algorithms set to those patterns to digitally > synthesize music which conforms to those patterns but that would create said > music automatically and potentially non stop; as in it does not stop writing > the music and playing it until you turn it off. However; I do not wish this > program to create a few tracks and play it on a loop, yes due to the patterns > and the algorithms, certain bars and phrases will inevitably repeat but that > doesn't mean just loop the same X amount of minutes. > > Before I fall off tangent and go into semantics, I am asking help as; > > The Music Genome Project > and Sony's musical mood auto playlist generator > > are similar to this underlying theme but have gone down the road of separate > applications and to the most important part: > > I do not wish to infringe or plagiarize or violate copyright or IP on others > said similar themed projects/products. > > my idea; > > use themes and patterns present in a selection of music to then create a > potentially infinite and unique mathematical auto play of algorithmically and > digitally created music > > what I need help with; > > inputting a large data set of audio files into an R Code application for > pattern and algorithm analysis with out infringing on open, finished or > ongoing projects. > > I do not need help with porting the algorithm and pattern generator to an > audio synthesis program--> that I am familiar with how to do. > > I do not need help with making the end result potentially infinitely > continued computation of algorithms within the found and set pattern > parameters---> that I also know how to do. > > I apologize for the long message and it's redundancies, it is simply my first > time here and I wanted to be thorough. > > Thank You for Reading! > > Any advice on this will be supremely appreciated! > > ---Darth Brando > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] hello i have a question on music analysis and mathematical synthesis related to r code
Hi Darth, Have a look at the tuneR package. Jim On Thu, Sep 8, 2016 at 3:57 PM, darth brando wrote: > Apologies for the long title but it is semi specific a topic and yes I am a > noobs user to the system. I have read the guide and will attempt to adhere to > the guide in this process and I do apologize in advance if I fail to do so, > this is my first time here. > > To the point; firstly version: > > I have windows 7 64 bit OS, I'm going to be working with most current up to > date version of R code for that OS with optional plug ins as needed, I am > going to be using R code in conjunction with Fruity Loop Studio, > dosbox/visual studio, and both synthesizer plug ins for Fruity Loop Studio > and a physical soundblaster sound card in my alienware 15 2015 version as > well as possibly a few emulator or OS instances for retro hardware to > experiment with. Its a lot of software to take in yea, some I know by heart > others like R Code I am a noob at. > > Now that the toolbox of software is out of the way; my Question: > > Music is largely mathematically based, R Code is perfect for analyzing very > large sets of data, naturally I saw some potential and wondered about the > specifics of the type of application which follows. > > I wish to use R Code to analyze portions of my personal library of music in > order to discover the main sets of underlying patterns within that portion of > my music library. I then wish to utilize patterns discovered to create a > program which using proper algorithms set to those patterns to digitally > synthesize music which conforms to those patterns but that would create said > music automatically and potentially non stop; as in it does not stop writing > the music and playing it until you turn it off. However; I do not wish this > program to create a few tracks and play it on a loop, yes due to the patterns > and the algorithms, certain bars and phrases will inevitably repeat but that > doesn't mean just loop the same X amount of minutes. > > Before I fall off tangent and go into semantics, I am asking help as; > > The Music Genome Project > and Sony's musical mood auto playlist generator > > are similar to this underlying theme but have gone down the road of separate > applications and to the most important part: > > I do not wish to infringe or plagiarize or violate copyright or IP on others > said similar themed projects/products. > > my idea; > > use themes and patterns present in a selection of music to then create a > potentially infinite and unique mathematical auto play of algorithmically and > digitally created music > > what I need help with; > > inputting a large data set of audio files into an R Code application for > pattern and algorithm analysis with out infringing on open, finished or > ongoing projects. > > I do not need help with porting the algorithm and pattern generator to an > audio synthesis program--> that I am familiar with how to do. > > I do not need help with making the end result potentially infinitely > continued computation of algorithms within the found and set pattern > parameters---> that I also know how to do. > > I apologize for the long message and it's redundancies, it is simply my first > time here and I wanted to be thorough. > > Thank You for Reading! > > Any advice on this will be supremely appreciated! > > ---Darth Brando > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
[R] hello i have a question on music analysis and mathematical synthesis related to r code
Apologies for the long title but it is semi specific a topic and yes I am a noobs user to the system. I have read the guide and will attempt to adhere to the guide in this process and I do apologize in advance if I fail to do so, this is my first time here. To the point; firstly version: I have windows 7 64 bit OS, I'm going to be working with most current up to date version of R code for that OS with optional plug ins as needed, I am going to be using R code in conjunction with Fruity Loop Studio, dosbox/visual studio, and both synthesizer plug ins for Fruity Loop Studio and a physical soundblaster sound card in my alienware 15 2015 version as well as possibly a few emulator or OS instances for retro hardware to experiment with. Its a lot of software to take in yea, some I know by heart others like R Code I am a noob at. Now that the toolbox of software is out of the way; my Question: Music is largely mathematically based, R Code is perfect for analyzing very large sets of data, naturally I saw some potential and wondered about the specifics of the type of application which follows. I wish to use R Code to analyze portions of my personal library of music in order to discover the main sets of underlying patterns within that portion of my music library. I then wish to utilize patterns discovered to create a program which using proper algorithms set to those patterns to digitally synthesize music which conforms to those patterns but that would create said music automatically and potentially non stop; as in it does not stop writing the music and playing it until you turn it off. However; I do not wish this program to create a few tracks and play it on a loop, yes due to the patterns and the algorithms, certain bars and phrases will inevitably repeat but that doesn't mean just loop the same X amount of minutes. Before I fall off tangent and go into semantics, I am asking help as; The Music Genome Project and Sony's musical mood auto playlist generator are similar to this underlying theme but have gone down the road of separate applications and to the most important part: I do not wish to infringe or plagiarize or violate copyright or IP on others said similar themed projects/products. my idea; use themes and patterns present in a selection of music to then create a potentially infinite and unique mathematical auto play of algorithmically and digitally created music what I need help with; inputting a large data set of audio files into an R Code application for pattern and algorithm analysis with out infringing on open, finished or ongoing projects. I do not need help with porting the algorithm and pattern generator to an audio synthesis program--> that I am familiar with how to do. I do not need help with making the end result potentially infinitely continued computation of algorithms within the found and set pattern parameters---> that I also know how to do. I apologize for the long message and it's redundancies, it is simply my first time here and I wanted to be thorough. Thank You for Reading! Any advice on this will be supremely appreciated! ---Darth Brando __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] hello everyone and happy new year
Melt your data.frame library(reshape2) s3.melted <- melt(s3) and plot using something along the lines of ggplot(s3.melted,aes(x=time2,y=value, colour=variable) + geom_line() (not tested) Hope this helps Ulrik On Wed, 6 Jan 2016 at 17:01 Matthias Worni wrote: > I got the following problem in R studio. Im trying to make a plot of > different time series using different colours, which so far worked finde > with the ggplot tool. However I did not manage to accually write a name to > the correct layers. Below you can see the code that I wrote and also woks. > Remember, the code works, meaning I get a plot with all 9 timeseires from > aa until ii. I just could not figure out how to label the different layers. > What should I do when I would like to add a name to the red line, the green > line and so on beside the plot (in al list). Thank you for your time. > > Best Matthias > > aa <- OHC_d[40:70] > bb <- OHC_d[50:80] > cc <- OHC_d[173:203] > dd <- OHC_d[205:235] > ee <- OHC_d[273:303] > ff <- OHC_d[292:322] > gg <- OHC_d[302:332] > hh <- OHC_d[370:400] > ii <- OHC_d[381:411] > > s3 <- data.frame(aa,bb,cc,dd,ee,ff,gg,hh,ii) > y2 <- ggplot(s3,aes(x=time2,y=aa),colour="black") + > geom_line() + > geom_line(data=s3,aes(y=bb),colour="red") + > geom_line() + > geom_line(data=s3,aes(y=cc),colour="green") + > geom_line() + > geom_line(data=s3,aes(y=dd),colour="blue") + > geom_line() + > geom_line(data=s3,aes(y=ee),colour="yellow") + > geom_line() + > geom_line(data=s3,aes(y=ff),colour="pink") + > geom_line() + > geom_line(data=s3,aes(y=gg),colour="purple")+ > geom_line()+ > geom_line(data=s3,aes(y=hh),colour="violet") > geom_line()+ > geom_line(data=s3,aes(y=ii),colour="grey") > y2 > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
[R] hello everyone and happy new year
I got the following problem in R studio. Im trying to make a plot of different time series using different colours, which so far worked finde with the ggplot tool. However I did not manage to accually write a name to the correct layers. Below you can see the code that I wrote and also woks. Remember, the code works, meaning I get a plot with all 9 timeseires from aa until ii. I just could not figure out how to label the different layers. What should I do when I would like to add a name to the red line, the green line and so on beside the plot (in al list). Thank you for your time. Best Matthias aa <- OHC_d[40:70] bb <- OHC_d[50:80] cc <- OHC_d[173:203] dd <- OHC_d[205:235] ee <- OHC_d[273:303] ff <- OHC_d[292:322] gg <- OHC_d[302:332] hh <- OHC_d[370:400] ii <- OHC_d[381:411] s3 <- data.frame(aa,bb,cc,dd,ee,ff,gg,hh,ii) y2 <- ggplot(s3,aes(x=time2,y=aa),colour="black") + geom_line() + geom_line(data=s3,aes(y=bb),colour="red") + geom_line() + geom_line(data=s3,aes(y=cc),colour="green") + geom_line() + geom_line(data=s3,aes(y=dd),colour="blue") + geom_line() + geom_line(data=s3,aes(y=ee),colour="yellow") + geom_line() + geom_line(data=s3,aes(y=ff),colour="pink") + geom_line() + geom_line(data=s3,aes(y=gg),colour="purple")+ geom_line()+ geom_line(data=s3,aes(y=hh),colour="violet") geom_line()+ geom_line(data=s3,aes(y=ii),colour="grey") y2 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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, minimal, self-contained, reproducible code.
Re: [R] Hello R user!
On 12/18/2013 04:33 AM, bibek sharma wrote: Hello R user, I have created two plots (attached!) using the codes below and would like to merge these figures in one. any suggestions are highly appreciated! Thanks, plot(graph1$yod,graph1$xod,data=graph1) dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) abline(h=0,v=0) plot(graph2$yod,graph2$xod,data=graph2) dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) abline(h=0,v=0) Hi Bibek, From what I can guess of the above, you might be interested in the size_n_color plot (plotrix). Jim __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R user!
Here is a simple example (without the proportional size bubbles--you've been given some references on that) using the lattice package: # one dataframe holds the data from both "sources" I call them. # they would be data from your two separate dataframes, # that you call graph1 and graph2 dd <- data.frame(x=rnorm(10), y=rnorm(10), source=sample(c("A","B"), 10, replace=TRUE)) dd library(lattice) xyplot(y~x | source, data=dd) # or another way xyplot(y~x, groups=source, data=dd, auto.key=TRUE) You'll want to think about how you are storing your data. Certain ways of doing it lend themselves to certain ways of graphing. Some ways make things difficult . . . --Chris Ryan On Tue, Dec 17, 2013 at 1:24 PM, Sarah Goslee wrote: > On Tue, Dec 17, 2013 at 1:04 PM, bibek sharma wrote: >> Hi Sarah, >> It is not about mfrow or mfcol. I would like to see both sets of data in >> one figure. >> All I want was combining these two plots to one. >> Any suggestions? >> Bibek > > Suggestions? Yes. Read the link I and others provided about > reproducible questions. > > Then there's the suggestion I already provided, using points() or > lines() to add more data, possibly with xlim or ylim specified (see > ?par for details). Without a reproducible example, I can't give > specific details. > >> Also, size of the circle in the plots represents rates and so should be >> shown in different sizes. >> I tried using plots and points but this did not give me different sizes. > > If you want to use base graphics, then cex is what you need (see, you > guessed it, ?par). You can pass a vector of sizes for your plotting > character. > > For more sophisticated approaches, you might google "bubble plot R" for ideas. > > Sarah > >> >> >> On Tue, Dec 17, 2013 at 9:45 AM, Sarah Goslee >> wrote: >>> >>> What do you mean by "merge these figures in one"? If you want two >>> figures on one page, see ?par - specifically mfrow and mfcol. >>> >>> If you want both sets of data in one figure, maybe ?points or ?lines >>> though I see you're already familiar with at least ?lines. >>> >>> The list doesn't take most attachments, and you might also take a look at: >>> >>> http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example >>> >>> Asking intelligible questions is the best strategy for receiving >>> intelligible answers. >>> >>> Sarah >>> >>> On Tue, Dec 17, 2013 at 12:33 PM, bibek sharma >>> wrote: >>> > Hello R user, >>> > >>> > I have created two plots (attached!) using the codes below >>> > and would like to merge these figures in one. any suggestions are highly >>> > appreciated! >>> > Thanks, >>> > >>> > plot(graph1$yod,graph1$xod,data=graph1) >>> > dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) >>> > symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, >>> > bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) >>> > abline(h=0,v=0) >>> > >>> > plot(graph2$yod,graph2$xod,data=graph2) >>> > dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) >>> > lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, >>> > bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) >>> > abline(h=0,v=0) >>> > >>> > Best, >>> > Bibek >>> >>> > > > > -- > Sarah Goslee > http://www.functionaldiversity.org > > __ > 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R user!
On Tue, Dec 17, 2013 at 1:04 PM, bibek sharma wrote: > Hi Sarah, > It is not about mfrow or mfcol. I would like to see both sets of data in > one figure. > All I want was combining these two plots to one. > Any suggestions? > Bibek Suggestions? Yes. Read the link I and others provided about reproducible questions. Then there's the suggestion I already provided, using points() or lines() to add more data, possibly with xlim or ylim specified (see ?par for details). Without a reproducible example, I can't give specific details. > Also, size of the circle in the plots represents rates and so should be > shown in different sizes. > I tried using plots and points but this did not give me different sizes. If you want to use base graphics, then cex is what you need (see, you guessed it, ?par). You can pass a vector of sizes for your plotting character. For more sophisticated approaches, you might google "bubble plot R" for ideas. Sarah > > > On Tue, Dec 17, 2013 at 9:45 AM, Sarah Goslee > wrote: >> >> What do you mean by "merge these figures in one"? If you want two >> figures on one page, see ?par - specifically mfrow and mfcol. >> >> If you want both sets of data in one figure, maybe ?points or ?lines >> though I see you're already familiar with at least ?lines. >> >> The list doesn't take most attachments, and you might also take a look at: >> >> http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example >> >> Asking intelligible questions is the best strategy for receiving >> intelligible answers. >> >> Sarah >> >> On Tue, Dec 17, 2013 at 12:33 PM, bibek sharma >> wrote: >> > Hello R user, >> > >> > I have created two plots (attached!) using the codes below >> > and would like to merge these figures in one. any suggestions are highly >> > appreciated! >> > Thanks, >> > >> > plot(graph1$yod,graph1$xod,data=graph1) >> > dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) >> > symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, >> > bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) >> > abline(h=0,v=0) >> > >> > plot(graph2$yod,graph2$xod,data=graph2) >> > dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) >> > lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, >> > bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) >> > abline(h=0,v=0) >> > >> > Best, >> > Bibek >> >> -- Sarah Goslee http://www.functionaldiversity.org __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R user!
What do you mean by "merge these figures in one"? If you want two figures on one page, see ?par - specifically mfrow and mfcol. If you want both sets of data in one figure, maybe ?points or ?lines though I see you're already familiar with at least ?lines. The list doesn't take most attachments, and you might also take a look at: http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example Asking intelligible questions is the best strategy for receiving intelligible answers. Sarah On Tue, Dec 17, 2013 at 12:33 PM, bibek sharma wrote: > Hello R user, > > I have created two plots (attached!) using the codes below > and would like to merge these figures in one. any suggestions are highly > appreciated! > Thanks, > > plot(graph1$yod,graph1$xod,data=graph1) > dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) > symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, > bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) > abline(h=0,v=0) > > plot(graph2$yod,graph2$xod,data=graph2) > dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) > lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, > bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) > abline(h=0,v=0) > > Best, > Bibek -- Sarah Goslee http://www.functionaldiversity.org __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R user!
What do you mean by "merge" them into one? Make both graphs appear on the same page of a document? Make a single figure containing both graphs? Plot data from both dataframes on the same set of axes? --Chris Ryan On Tue, Dec 17, 2013 at 12:33 PM, bibek sharma wrote: > Hello R user, > > I have created two plots (attached!) using the codes below > and would like to merge these figures in one. any suggestions are highly > appreciated! > Thanks, > > plot(graph1$yod,graph1$xod,data=graph1) > dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) > symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, > bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) > abline(h=0,v=0) > > plot(graph2$yod,graph2$xod,data=graph2) > dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) > lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, > bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) > abline(h=0,v=0) > > Best, > Bibek > > __ > 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, minimal, self-contained, reproducible code. > __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R user!
The plots did not arrive. The R-help list is fussy about what it allows to go through. Actually the best way of doing things is to use dput() to provide sample data. See https://github.com/hadley/devtools/wiki/Reproducibility or http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example To answer your question probably you want to have a look a "mfcol" or "mfrow" under ?par. John Kane Kingston ON Canada > -Original Message- > From: mbhpat...@gmail.com > Sent: Tue, 17 Dec 2013 09:33:30 -0800 > To: r-help@r-project.org > Subject: [R] Hello R user! > > Hello R user, > > I have created two plots (attached!) using the codes below > and would like to merge these figures in one. any suggestions are highly > appreciated! > Thanks, > > plot(graph1$yod,graph1$xod,data=graph1) > dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) > symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, > bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) > abline(h=0,v=0) > > plot(graph2$yod,graph2$xod,data=graph2) > dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) > lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, > bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) > abline(h=0,v=0) > > Best, > Bibek > __ > 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, minimal, self-contained, reproducible code. FREE ONLINE PHOTOSHARING - Share your photos online with your friends and family! Visit http://www.inbox.com/photosharing to find out more! __ 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, minimal, self-contained, reproducible code.
[R] Hello R user!
Hello R user, I have created two plots (attached!) using the codes below and would like to merge these figures in one. any suggestions are highly appreciated! Thanks, plot(graph1$yod,graph1$xod,data=graph1) dfx = data.frame(ev1=graph1$xod, ev2=graph1$yod, ev3=abs(graph1$dif)) symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, bg="black", fg=NULL,xlim=c(-35,35),ylim=c(-35,35)) abline(h=0,v=0) plot(graph2$yod,graph2$xod,data=graph2) dfx = data.frame(ev1=graph2$xod, ev2=graph2$yod, ev3=abs(graph2$dif)) lines(symbols(x=dfx$ev1, y=dfx$ev2, circles=dfx$ev3,inches=1/8, ann=F, bg="blue", fg=NULL,xlim=c(-35,35),ylim=c(-35,35))) abline(h=0,v=0) Best, Bibek __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello
Hello, To compute the difference between two date/time objects use the minus operator. See ?difftime ?`-.POSIXt` Hope this helps, Rui Barradas Em 01-05-2013 16:58, Eric Mintah escreveu: Dear sir/madam, I used as.Date,strptime, POXIct, and POXIlt functions to work on my data in R. Please, I would like to know the function to find the difference between the dates and times in two different columns and make a new column for that. Thank you [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
[R] Hello
Dear sir/madam, I used as.Date,strptime, POXIct, and POXIlt functions to work on my data in R. Please, I would like to know the function to find the difference between the dates and times in two different columns and make a new column for that. Thank you [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello. PLEASE HELP
No. 1. If this is homework, this list has a no homework policy. 2. If it is not, please read the posting guide (link at bottom) and post a coherent question in a form that can be answered. In particular, you need to provide the code that failed, the error messages, the packages you used, the R version,... -- Bert On Sun, Apr 28, 2013 at 11:36 AM, wrote: > I am a student of NTNU-Norway. > > I have some issues in fitting an AR-GARCH with seasonal components. > > My AR model contains the terms AR(1), AR(5), AR(20) and AR(60) and a simple > GARCH(1,1) model. > > I received lot of errors when trying to do this. > > Can you please provide me the source code of fitting such model to my data? > > Best Regards > Ampem. > > __ > 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, minimal, self-contained, reproducible code. -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact Info: Phone: 467-7374 Website: http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm __ 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, minimal, self-contained, reproducible code.
[R] Hello. PLEASE HELP
I am a student of NTNU-Norway. I have some issues in fitting an AR-GARCH with seasonal components. My AR model contains the terms AR(1), AR(5), AR(20) and AR(60) and a simple GARCH(1,1) model. I received lot of errors when trying to do this. Can you please provide me the source code of fitting such model to my data? Best Regards Ampem. __ 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, minimal, self-contained, reproducible code.
[R] Hello
Hello All, I am new to the list. I have been learning to use R recently and its been great to see so much help available. However I must admit, I have stumbled upon a problem with correlation matrices and was hoping if someone could help. I have an excel workbook with a sheet per drug. Each sheet has 40 patients' drug response measured over different time points (several days). I have already been able to create a cor() matrix of all the drugs for all the patients on each day. However I have been asked to prepare a matrix of all days of all people for all drugs. So since there are 26 drugs, I am expected to form a 26x26 matrix for all this data. I am not even sure if this is possible. Can anyone point me in the right direction? PS: If I make individual matrix and combine them using cbind(), it is not a 26x26 sqaure matrix. Any help would be greatly appreciated. Thanks [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R User
HI, Try: dat2[,Time1:=c(0,diff(Time)),by=ID] dat2[,CumSTime1:=cumsum(Time1),by=ID] head(dat2,4) # ID Time Time1 CumSTime1 #1: 1 3 0 0 #2: 1 6 3 3 #3: 1 7 1 4 #4: 1 10 3 7 A.K. - Original Message - From: bibek sharma To: arun Cc: Sent: Friday, December 14, 2012 12:56 PM Subject: Re: [R] Hello R User Hi Arun, Great! Once we get, time1, I again wanna add time to its previous value for example, wanna get 0,3 4 etc... can I have suggestion? On Fri, Dec 14, 2012 at 9:42 AM, arun wrote: > Hi, > > You could also use library(data.table) to do this faster. > dat1<-read.table(text=" > ID Time > 1 3 > 1 6 > 1 7 > 1 10 > 1 16 > 2 12 > 2 18 > 2 19 > 2 25 > 2 28 > 2 30 > ",sep="",header=TRUE) > library(data.table) > dat2<-data.table(dat1) > res<-dat2[,Time1:=c(0,diff(Time)),by=ID] > head(res,3) > # ID Time Time1 > #1: 1 3 0 > #2: 1 6 3 > #3: 1 7 1 > > #Comparing different approaches: > set.seed(55) > dat3<- > data.frame(ID=rep(1:1000,each=500),Value=sample(1:800,5e5,replace=TRUE)) > dat4<-data.table(dat3) > system.time(dat3$Value1<-ave(dat3$Value,dat3$ID,FUN=function(x) c(0,diff(x > # user system elapsed > # 0.312 0.000 0.313 > > ids<-unique(dat3$ID) > system.time({ > for(id in ids){ > dat3$Value[dat3$ID==id]<-c(0,diff(dat3$Value[dat3$ID==id])) > } }) > # user system elapsed > # 36.938 0.868 37.873 > > system.time(dat5<-dat4[,Value1:=c(0,diff(Value)),by=ID]) > # user system elapsed > # 0.036 0.000 0.037 > head(dat5) > # ID Value Value1 > #1: 1 439 0 > #2: 1 175 -264 > #3: 1 28 -147 > #4: 1 634 606 > #5: 1 449 -185 > #6: 1 60 -389 > head(dat3) > # ID Value Value1 > #1 1 0 0 > #2 1 -264 -264 > #3 1 -147 -147 > #4 1 606 606 > #5 1 -185 -185 > #6 1 -389 -389 > > A.K. > > > > > > > > > > > > > - Original Message - > From: bibek sharma > To: R help > Cc: > Sent: Friday, December 14, 2012 10:51 AM > Subject: [R] Hello R User > > Hello R User, > In the sample data given below, time is recorded for each id > subsequently. For the analysis, for each id, I would like to set 1st > recorded time to zero and thereafter find the difference from previous > time. I.e. for ID==1, I would like to see Time=0,3,1,3,6. This needs > to be implemented to big data set. > Any suggestions are much appreciated! > Thanks, > Bibek > > ID Time > 1 3 > 1 6 > 1 7 > 1 10 > 1 16 > 2 12 > 2 18 > 2 19 > 2 25 > 2 28 > 2 30 > > __ > 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, minimal, self-contained, reproducible code. > __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R User
Hi, You could also use library(data.table) to do this faster. dat1<-read.table(text=" ID Time 1 3 1 6 1 7 1 10 1 16 2 12 2 18 2 19 2 25 2 28 2 30 ",sep="",header=TRUE) library(data.table) dat2<-data.table(dat1) res<-dat2[,Time1:=c(0,diff(Time)),by=ID] head(res,3) # ID Time Time1 #1: 1 3 0 #2: 1 6 3 #3: 1 7 1 #Comparing different approaches: set.seed(55) dat3<- data.frame(ID=rep(1:1000,each=500),Value=sample(1:800,5e5,replace=TRUE)) dat4<-data.table(dat3) system.time(dat3$Value1<-ave(dat3$Value,dat3$ID,FUN=function(x) c(0,diff(x # user system elapsed # 0.312 0.000 0.313 ids<-unique(dat3$ID) system.time({ for(id in ids){ dat3$Value[dat3$ID==id]<-c(0,diff(dat3$Value[dat3$ID==id])) } }) # user system elapsed # 36.938 0.868 37.873 system.time(dat5<-dat4[,Value1:=c(0,diff(Value)),by=ID]) # user system elapsed # 0.036 0.000 0.037 head(dat5) # ID Value Value1 #1: 1 439 0 #2: 1 175 -264 #3: 1 28 -147 #4: 1 634 606 #5: 1 449 -185 #6: 1 60 -389 head(dat3) # ID Value Value1 #1 1 0 0 #2 1 -264 -264 #3 1 -147 -147 #4 1 606 606 #5 1 -185 -185 #6 1 -389 -389 A.K. - Original Message - From: bibek sharma To: R help Cc: Sent: Friday, December 14, 2012 10:51 AM Subject: [R] Hello R User Hello R User, In the sample data given below, time is recorded for each id subsequently. For the analysis, for each id, I would like to set 1st recorded time to zero and thereafter find the difference from previous time. I.e. for ID==1, I would like to see Time=0,3,1,3,6. This needs to be implemented to big data set. Any suggestions are much appreciated! Thanks, Bibek ID Time 1 3 1 6 1 7 1 10 1 16 2 12 2 18 2 19 2 25 2 28 2 30 __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R User
Hi Bibek, how about this? dta<-read.table(textConnection("ID Time 1 3 1 6 1 7 1 10 1 16 2 12 2 18 2 19 2 25 2 28 2 30"),header=T) dta$delta<-with(dta,ave(Time,ID,FUN=function(x)c(0,diff(x dta hth. Am 14.12.2012 16:51, schrieb bibek sharma: > Hello R User, > In the sample data given below, time is recorded for each id > subsequently. For the analysis, for each id, I would like to set 1st > recorded time to zero and thereafter find the difference from previous > time. I.e. for ID==1, I would like to see Time=0,3,1,3,6. This needs > to be implemented to big data set. > Any suggestions are much appreciated! > Thanks, > Bibek > > IDTime > 1 3 > 1 6 > 1 7 > 1 10 > 1 16 > 2 12 > 2 18 > 2 19 > 2 25 > 2 28 > 2 30 > > __ > 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, minimal, self-contained, reproducible code. > __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R User
dataset<-data.frame(id=c(1,1,2,3,3,3),time=c(3,5,1,2,4,6)) dataset id time 1 13 2 15 3 21 4 32 5 34 6 36 ids<-unique(dataset$id) for(id in ids){ + dataset$time[dataset$id==id]<-c(0,diff(dataset$time[dataset$id==id])) + } dataset id time 1 10 2 12 3 20 4 30 5 32 6 32 might not be the fastest though. On 14.12.2012, at 16:51, bibek sharma wrote: > Hello R User, > In the sample data given below, time is recorded for each id > subsequently. For the analysis, for each id, I would like to set 1st > recorded time to zero and thereafter find the difference from previous > time. I.e. for ID==1, I would like to see Time=0,3,1,3,6. This needs > to be implemented to big data set. > Any suggestions are much appreciated! > Thanks, > Bibek > > IDTime > 1 3 > 1 6 > 1 7 > 1 10 > 1 16 > 2 12 > 2 18 > 2 19 > 2 25 > 2 28 > 2 30 > > __ > 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello R User
HI, Try this: dat1<-read.table(text=" ID Time 1 3 1 6 1 7 1 10 1 16 2 12 2 18 2 19 2 25 2 28 2 30 ",sep="",header=TRUE) dat1$Time1<-ave(dat1$Time,dat1$ID,FUN=function(x) c(0,diff(x))) head(dat1,3) # ID Time Time1 #1 1 3 0 #2 1 6 3 #3 1 7 1 #or dat2<-unsplit(lapply(split(dat1,dat1$ID),function(x) {x$Time<-c(0,diff(x[,2])); return(x)}),dat1$ID) head(dat2,3) # ID Time #1 1 0 #2 1 3 #3 1 1 A.K. - Original Message - From: bibek sharma To: R help Cc: Sent: Friday, December 14, 2012 10:51 AM Subject: [R] Hello R User Hello R User, In the sample data given below, time is recorded for each id subsequently. For the analysis, for each id, I would like to set 1st recorded time to zero and thereafter find the difference from previous time. I.e. for ID==1, I would like to see Time=0,3,1,3,6. This needs to be implemented to big data set. Any suggestions are much appreciated! Thanks, Bibek ID Time 1 3 1 6 1 7 1 10 1 16 2 12 2 18 2 19 2 25 2 28 2 30 __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
[R] Hello R User
Hello R User, In the sample data given below, time is recorded for each id subsequently. For the analysis, for each id, I would like to set 1st recorded time to zero and thereafter find the difference from previous time. I.e. for ID==1, I would like to see Time=0,3,1,3,6. This needs to be implemented to big data set. Any suggestions are much appreciated! Thanks, Bibek ID Time 1 3 1 6 1 7 1 10 1 16 2 12 2 18 2 19 2 25 2 28 2 30 __ 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, minimal, self-contained, reproducible code.
Re: [R] hello about proxy configuration
On 19/06/2011 7:29 PM, p_conno...@slingshot.co.nz wrote: Quoting 王海生: dear everyone system:windows XP R2.13.0 I download the windows binary from website and successfully install it ,because I use a proxy server ,when I follow the instrction as follows: I set a system property R_HOME=C:\Program Files\R\R-2.13.0 " R_HOME\bin\i386\Rgui.exe http_proxy=http://211.83.105.140:808/"; It's more likely to be 8080, not 808. If that doesn't work, there is an easier way. When you install the Windows binary, there is an option called 'Internet access'. Change from the default to internet2, and R will then use your browser settings. Assuming your browser works, R will have the correct proxy setting. You can make this change any time, not only at install time. Temporarily in a session using the setInternet2 function, permanently by editing the shortcut to include the --internet2 option on the command line. Duncan Murdoch HTH error:: '\i' is an unrecognized escape in character string starting "R_HOME\bin\i" what is wrong with my command? I just want to know is there any other way such as mannually edit a configuration file to set the proxy? kind regards edwin_uestc Center of E-Health technology and Engineering in UESTC [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
Re: [R] hello about proxy configuration
Quoting 王海生 : dear everyone system:windows XP R2.13.0 I download the windows binary from website and successfully install it ,because I use a proxy server ,when I follow the instrction as follows: I set a system property R_HOME=C:\Program Files\R\R-2.13.0 " R_HOME\bin\i386\Rgui.exe http_proxy=http://211.83.105.140:808/"; It's more likely to be 8080, not 808. If that doesn't work, there is an easier way. When you install the Windows binary, there is an option called 'Internet access'. Change from the default to internet2, and R will then use your browser settings. Assuming your browser works, R will have the correct proxy setting. HTH error:: '\i' is an unrecognized escape in character string starting "R_HOME\bin\i" what is wrong with my command? I just want to know is there any other way such as mannually edit a configuration file to set the proxy? kind regards edwin_uestc Center of E-Health technology and Engineering in UESTC [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
[R] hello about proxy configuration
dear everyone system:windows XP R2.13.0 I download the windows binary from website and successfully install it ,because I use a proxy server ,when I follow the instrction as follows: I set a system property R_HOME=C:\Program Files\R\R-2.13.0 " R_HOME\bin\i386\Rgui.exe http_proxy=http://211.83.105.140:808/"; error:: '\i' is an unrecognized escape in character string starting "R_HOME\bin\i" what is wrong with my command? I just want to know is there any other way such as mannually edit a configuration file to set the proxy? kind regards edwin_uestc Center of E-Health technology and Engineering in UESTC [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
Hi: Use the PolynomF package. Once you've installed it from CRAN, try the following code: library(PolynomF) example(PolynomF) The example contains a block of code to construct a list of Hermite polynomials which are stored in a list object named H. After you've run the example, try deriv(H) integral(H) A careful reading of the documentation should help you solve a number of useful problems. HTH, Dennis On Sun, May 29, 2011 at 12:27 PM, Amanda Zeqiri wrote: > > Hi, I'm a student studing for Math. & Infor. Ing. in Tirana,Albania and a > have a final project in R to finish in a week, so I badly need your help... > > The topic of the project is "Hermite Interpolation" and a I have allready > done a script that finds the approximation in a certan point x, H(x) ~f(x) > but i want to find the value of his derivative H'(x)~f'(x) as well and to do > this I need the polynomial form of H(x)=a1*x^(2n+1) + a2*x^(2n) + ... + an > but R requires x to be a value not a declared variable. Then I want to find > the derivative of this polynomial and evaluate it. > > So my point is how can I built a function f(x)=x^+x+1 and do its derivate > dipending on x. > > I also tryed: > >>library(Ryacas) >>x <- Sym("x") >>Simplify(deriv(sqrt(1 - x^2),x,2)) but It didn't work > > please help me or I'll have to repeat the year and it's my last year. > > Thank you > > [[alternative HTML version deleted]] > > __ > 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, minimal, self-contained, reproducible code. > __ 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, minimal, self-contained, reproducible code.
[R] Hello!
Hi, I'm a student studing for Math. & Infor. Ing. in Tirana,Albania and a have a final project in R to finish in a week, so I badly need your help... The topic of the project is "Hermite Interpolation" and a I have allready done a script that finds the approximation in a certan point x, H(x) ~f(x) but i want to find the value of his derivative H'(x)~f'(x) as well and to do this I need the polynomial form of H(x)=a1*x^(2n+1) + a2*x^(2n) + ... + an but R requires x to be a value not a declared variable. Then I want to find the derivative of this polynomial and evaluate it. So my point is how can I built a function f(x)=x^+x+1 and do its derivate dipending on x. I also tryed: >library(Ryacas) >x <- Sym("x") >Simplify(deriv(sqrt(1 - x^2),x,2)) but It didn't work please help me or I'll have to repeat the year and it's my last year. Thank you [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
> Date: Sun, 27 Feb 2011 14:31:57 -0800 [[elided Hotmail spam]] > From: jwiley.ps...@gmail.com > To: marchy...@hotmail.com > CC: r-help@r-project.org > > On Sun, Feb 27, 2011 at 9:16 AM, Mike Marchywka wrote: > > I can probably find more "negative" ones if you are interested LOL :) > > I was hoping to find another interesting technical conversation > > to which I could contribute a few thoughts but subjective matters > > do come up with statistical analysis and confirmation bias, rationalization, > > the "power of positive thinking", believing in your data when you > > should be testing it, is often the biggest problem to getting > > accurate perspective especially with post hoc analysis questions > > as often come up here. Look at some ongoing discussions about scientific > > literature and you find all kinds of problems with failure to publish > > negative > > results, including clinical drug trials. The power of "positive" thinking, > > seeing what you want to see, wastes a lot of other people's time. You only > > need to spend > > a few days trying to replicate "positive" results, or just integrate > > them into your understanding of a complex system like a living organism, > > or buy a few hyped securities to understand how bug this problem really is. > > I realize Ulysses was riddled with elaborate enigmas and subtle > references (though Bert and Ravi must have hated the ending), so I > tried to examine this paragraph from multiple angles, but I keep > coming back to, "Huh?" Well, the OP, besides criticizing one post in particular, criticized the list for being "too negative" or similar nebulous words to that effect. Not really sure what caused him to formulate that conclusion but as I mentioned in my private reply( you didn't yet clarify the "huh" question LOL) , I've been pursing a jihad against words that have no meaning, usually adjectives where numbers would make more sense. In this case, by "positive" I assumed he should have said "accurate" as making something "positive" by some person's perspective usually just means someone else has to waste a lot of time figuring out the reality of a given conclusion. In one common usage, a "positive" unreproducible result is often published without the context of all the negative ones. Optimistic "positive" adjectives are usually inserted to boost a given perspective with no justification. Again, the ways in which insidious "positive" adjectives waste the time of others is hard to state clearly in a terse post as there are many, many,many, like humnguous ways that strings of adjectives are used instead of numbers and it can be a total distraction for sure :) Did you buy a house in the US a few years ago based on more "positive" valuations than careful examination of fundamental issues may have justified? Often simple questions about data, to the point of being snide maybe, could offer more compelling explanations than " I took 3 data points, shopped for a post hoc analysis method that someone on the internet said was proper that gave me the most positive result and I published that one." Questions come up here quite a bit that may "best" be answered with statements the OP would label as "negative." > > > > > I personally > > question the utility of the post you cite but I often post things > > in a hurry and make stupid statements myself ( and I don't let money > > change that LOL). You can only type "see the posting guidelines" so > > many times. Punctuation problems can make things hard to read etc. > > I really wouldn't take it that "negatively" whatever that means. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
On Sun, Feb 27, 2011 at 9:16 AM, Mike Marchywka wrote: > I can probably find more "negative" ones if you are interested LOL :) > I was hoping to find another interesting technical conversation > to which I could contribute a few thoughts but subjective matters > do come up with statistical analysis and confirmation bias, rationalization, > the "power of positive thinking", believing in your data when you > should be testing it, is often the biggest problem to getting > accurate perspective especially with post hoc analysis questions > as often come up here. Look at some ongoing discussions about scientific > literature and you find all kinds of problems with failure to publish negative > results, including clinical drug trials. The power of "positive" thinking, > seeing what you want to see, wastes a lot of other people's time. You only > need to spend > a few days trying to replicate "positive" results, or just integrate > them into your understanding of a complex system like a living organism, > or buy a few hyped securities to understand how bug this problem really is. I realize Ulysses was riddled with elaborate enigmas and subtle references (though Bert and Ravi must have hated the ending), so I tried to examine this paragraph from multiple angles, but I keep coming back to, "Huh?" > > I personally > question the utility of the post you cite but I often post things > in a hurry and make stupid statements myself ( and I don't let money > change that LOL). You can only type "see the posting guidelines" so > many times. Punctuation problems can make things hard to read etc. > I really wouldn't take it that "negatively" whatever that means. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
> Date: Sun, 27 Feb 2011 08:27:00 -0800 > From: markkne...@gmail.com > To: gunter.ber...@gene.com > CC: r-help@r-project.org > Subject: Re: [R] Hello! > > On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter wrote: > > Are you a fan of James Joyce? Is the Caps key on your keyboard broken? > > > > -- Bert > > Are your snide comments adding anything to the conversation? [ ... ] > > As for me, this is a farewell to this list which is far and away the > most negative list I've ever subscribed to. I'll get answers to R > questions elsewhere. I don't need people like Mr. Berton Gunter in my > life. (He's now blacklisted in my email tools anyway as he seldom has > much to offer other than snide comments.) I can probably find more "negative" ones if you are interested LOL :) I was hoping to find another interesting technical conversation to which I could contribute a few thoughts but subjective matters do come up with statistical analysis and confirmation bias, rationalization, the "power of positive thinking", believing in your data when you should be testing it, is often the biggest problem to getting accurate perspective especially with post hoc analysis questions as often come up here. Look at some ongoing discussions about scientific literature and you find all kinds of problems with failure to publish negative results, including clinical drug trials. The power of "positive" thinking, seeing what you want to see, wastes a lot of other people's time. You only need to spend a few days trying to replicate "positive" results, or just integrate them into your understanding of a complex system like a living organism, or buy a few hyped securities to understand how bug this problem really is. I personally question the utility of the post you cite but I often post things in a hurry and make stupid statements myself ( and I don't let money change that LOL). You can only type "see the posting guidelines" so many times. Punctuation problems can make things hard to read etc. I really wouldn't take it that "negatively" whatever that means. > > Good bye, > Mark > __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
I actually thought that Bert's comment was tastefully funny. I, too, do not enjoy reading all lowercase emails or emails that use ridiculous abbreviations. In any case, I am not sure from where Mark gets the idea that the list is way too negative. In my earlier days of learning S-Plus and R, I used to get chided by some of the R gurus, but that only fired me up to learn more. Even now when I give presentations or send manuscripts to journals, I look for comments that force me to think outside of my comfort zone. I think this list has a good balance of helpfulness and challenging criticisms. Mark, by reacting to the comments in a personal manner, I am afraid that you are the loser. Best, Ravi. Ravi Varadhan, Ph.D. Assistant Professor, Division of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: rvarad...@jhmi.edu - Original Message - From: "Rios,Alfredo Arturo" Date: Sunday, February 27, 2011 12:40 pm Subject: Re: [R] Hello! To: Mark Knecht Cc: r-help@r-project.org, Bert Gunter > Hi Mark, > > I agree that comment was completely unnecessary, but Berton Gunter does > contribute (not defending him) especially, in the "caution make sure > you > know what you are doing and consult your local statistician > comments." I > would look for the nuggets of wisdom in his comments and ignore the > "extra material." By leaving the list you will be doing yourself a > disservice as in addition to the few people who give good comments + > something else, there is also people who care and give excellent help. > >From this last people I would not only learn stats+R but also about > their generosity and kindness to the other list members. > > Alfredo > > > > > On Sun, 27 Feb 2011 08:27:00 -0800, Mark Knecht > wrote: > > On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter > wrote: > >> Are you a fan of James Joyce? Is the Caps key on your keyboard broken? > >> > >> -- Bert > > > > Are your snide comments adding anything to the conversation? > > > > Do you allow for the possibility of an ESL speaker writing an email > to > > this list? > > > > Do you know that a lot of kids, right or wrong, don't use the Caps > key > > in email or test messages. > > > > Or, far more likely, do you just feel so far superior that you need > to > > get yourself off putting people down? > > > > > > As for me, this is a farewell to this list which is far and away the > > most negative list I've ever subscribed to. I'll get answers to R > > questions elsewhere. I don't need people like Mr. Berton Gunter in > my > > life. (He's now blacklisted in my email tools anyway as he seldom has > > much to offer other than snide comments.) > > > > Good bye, > > Mark > > > > __ > > R-help@r-project.org mailing list > > > > PLEASE do read the posting guide > > and provide commented, minimal, self-contained, reproducible code. > > __ > R-help@r-project.org mailing list > > PLEASE do read the posting guide > and provide commented, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
Hi Mark, I agree that comment was completely unnecessary, but Berton Gunter does contribute (not defending him) especially, in the "caution make sure you know what you are doing and consult your local statistician comments." I would look for the nuggets of wisdom in his comments and ignore the "extra material." By leaving the list you will be doing yourself a disservice as in addition to the few people who give good comments + something else, there is also people who care and give excellent help. >From this last people I would not only learn stats+R but also about their generosity and kindness to the other list members. Alfredo On Sun, 27 Feb 2011 08:27:00 -0800, Mark Knecht wrote: > On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter wrote: >> Are you a fan of James Joyce? Is the Caps key on your keyboard broken? >> >> -- Bert > > Are your snide comments adding anything to the conversation? > > Do you allow for the possibility of an ESL speaker writing an email to > this list? > > Do you know that a lot of kids, right or wrong, don't use the Caps key > in email or test messages. > > Or, far more likely, do you just feel so far superior that you need to > get yourself off putting people down? > > > As for me, this is a farewell to this list which is far and away the > most negative list I've ever subscribed to. I'll get answers to R > questions elsewhere. I don't need people like Mr. Berton Gunter in my > life. (He's now blacklisted in my email tools anyway as he seldom has > much to offer other than snide comments.) > > Good bye, > Mark > > __ > 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter wrote: > Are you a fan of James Joyce? Is the Caps key on your keyboard broken? > > -- Bert Are your snide comments adding anything to the conversation? Do you allow for the possibility of an ESL speaker writing an email to this list? Do you know that a lot of kids, right or wrong, don't use the Caps key in email or test messages. Or, far more likely, do you just feel so far superior that you need to get yourself off putting people down? As for me, this is a farewell to this list which is far and away the most negative list I've ever subscribed to. I'll get answers to R questions elsewhere. I don't need people like Mr. Berton Gunter in my life. (He's now blacklisted in my email tools anyway as he seldom has much to offer other than snide comments.) Good bye, Mark __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
Look at ?readline for input from the user. For reading Excel files, a good place to start would be page 26 of the manual: http://cran.r-project.org/doc/manuals/R-data.pdf -- Robert Tirrell | r...@stanford.edu | (607) 437-6532 Program in Biomedical Informatics | Butte Lab | Stanford University On Sat, Feb 26, 2011 at 22:11, Bert Gunter wrote: > Are you a fan of James Joyce? Is the Caps key on your keyboard broken? > > -- Bert > > On Sat, Feb 26, 2011 at 10:09 AM, tsvi sabo wrote: > > Hi,i would like to know what is the best way to write a procedure in R,it > seems that when i run a script it doesn't wait for previouse code lines to > be excuted,i would like to know how to force R to hold on a code line until > it is processed,for example:i want to choose from a list the name of the > sheet in which to do an analisys,then i want to run a goodness of fit test > for the data on that chosen sheet,how can i do that?and if it is not too > much i have couple more questions,1. how to retrieve the names of excel's > file sheets to an array ?2. how to retrieve the names of excel's Sheet > columns names to an array ? > > thanks a lot,Tsvi Sabo > > > > > > > >[[alternative HTML version deleted]] > > > > > > __ > > 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, minimal, self-contained, reproducible code. > > > > > > > > -- > Bert Gunter > Genentech Nonclinical Biostatistics > > __ > 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, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
tsvi sabo wrote: Hi,i would like to know what is the best way to write a procedure in R,it seems that when i run a script it doesn't wait for previouse code lines to be excuted,i would like to know how to force R to hold on a code line until it is processed,for example:i want to choose from a list the name of the sheet in which to do an analisys,then i want to run a goodness of fit test for the data on that chosen sheet,how can i do that?and if it is not too much i have couple more questions,1. how to retrieve the names of excel's file sheets to an array ?2. how to retrieve the names of excel's Sheet columns names to an array ? thanks a lot,Tsvi Sabo Hi Tsvi, I think you are asking about halting the execution of a script so that the output doesn't disappear up the screen like smoke. You can use: readline() which will halt the script execution until the user presses the key. As for the Excel questions, there are several packages like "xlsx" that might do what you want there. Jim __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello!
Are you a fan of James Joyce? Is the Caps key on your keyboard broken? -- Bert On Sat, Feb 26, 2011 at 10:09 AM, tsvi sabo wrote: > Hi,i would like to know what is the best way to write a procedure in R,it > seems that when i run a script it doesn't wait for previouse code lines to be > excuted,i would like to know how to force R to hold on a code line until it > is processed,for example:i want to choose from a list the name of the sheet > in which to do an analisys,then i want to run a goodness of fit test for the > data on that chosen sheet,how can i do that?and if it is not too much i have > couple more questions,1. how to retrieve the names of excel's file sheets to > an array ?2. how to retrieve the names of excel's Sheet columns names to an > array ? > thanks a lot,Tsvi Sabo > > > > [[alternative HTML version deleted]] > > > __ > 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, minimal, self-contained, reproducible code. > > -- Bert Gunter Genentech Nonclinical Biostatistics __ 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, minimal, self-contained, reproducible code.
[R] Hello!
Hi,i would like to know what is the best way to write a procedure in R,it seems that when i run a script it doesn't wait for previouse code lines to be excuted,i would like to know how to force R to hold on a code line until it is processed,for example:i want to choose from a list the name of the sheet in which to do an analisys,then i want to run a goodness of fit test for the data on that chosen sheet,how can i do that?and if it is not too much i have couple more questions,1. how to retrieve the names of excel's file sheets to an array ?2. how to retrieve the names of excel's Sheet columns names to an array ? thanks a lot,Tsvi Sabo [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071?
> The caret package is so good, I am learning it, but one problem is that > nearZeroVar function can be used to identify near zero–variance variables and > it only identify, how can I remove those variables that were identify, > because I have many zero- or near zero- ones, it is not realistic to removel > it by hand, can this function can identify and removel those ones > automatically? > newx <- x[, -nearZeroVar(x), drop = FALSE] should do it __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071?
Hello, Max The caret package is so good, I am learning it, but one problem is that nearZeroVar function can be used to identify near zeroâvariance variables and it only identify, how can I remove those variables that were identify, because I have many zero- or near zero- ones, it is not realistic to removel it by hand, can this function can identify and removel those ones automatically? looking for your reply. kevin å¨2009-12-19ï¼"Max Kuhn [via R]" åéï¼ -åå§é®ä»¶- å件人:"Max Kuhn [via R]" åéæ¶é´:2009å¹´12æ19æ¥ ææå æ¶ä»¶äºº:bbslover 主é¢:Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071? You can get this using the caret package. There are a few package vignettes that come with the package and a JSS article http://www.jstatsoft.org/v28/i05/paper about the package. Max On Fri, Dec 18, 2009 at 12:26 PM, bbslover <[hidden email]> wrote: > > as known, svm need tune some parameters like cost,gamma and epsilon to get > better performance,but one question appear, how can i monitor the > performance . generally speaking ,we chose the cross-validation MSE in the > training set, but It seems svm can not return the cross-validation MSE > value, we just get it from "summary model.svm", if I write a loop, have no > idear call the cros-validate MSE, and no way to monitor this performance > ,how can I do? > > > -- > View this message in > context:http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p974942.html > Sent from the R help mailing list archive at Nabble.com. > > __ >[hidden email]mailing list >https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guidehttp://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Max __ [hidden email]mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guidehttp://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. View message @http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p975038.html To unsubscribe from Hello all, How can I get corss-validation MSE of SVM in e1071?,click here. -- View this message in context: http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p975955.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071?
thank you for your help, caret package is so powerful , it can do many things. I now, need learn how to apply to my problems. Max Kuhn wrote: > > You can get this using the caret package. There are a few package > vignettes that come with the package and a JSS article > > http://www.jstatsoft.org/v28/i05/paper > > about the package. > > Max > > On Fri, Dec 18, 2009 at 12:26 PM, bbslover wrote: >> >> as known, svm need tune some parameters like cost,gamma and epsilon to >> get >> better performance,but one question appear, how can i monitor the >> performance . generally speaking ,we chose the cross-validation MSE in >> the >> training set, but It seems svm can not return the cross-validation MSE >> value, we just get it from "summary model.svm", if I write a loop, have >> no >> idear call the cros-validate MSE, and no way to monitor this performance >> ,how can I do? >> >> >> -- >> View this message in context: >> http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p974942.html >> Sent from the R help mailing list archive at Nabble.com. >> >> __ >> 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, minimal, self-contained, reproducible code. >> > > > > -- > > Max > > __ > 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, minimal, self-contained, reproducible code. > > -- View this message in context: http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p975231.html Sent from the R help mailing list archive at Nabble.com. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071?
You can get this using the caret package. There are a few package vignettes that come with the package and a JSS article http://www.jstatsoft.org/v28/i05/paper about the package. Max On Fri, Dec 18, 2009 at 12:26 PM, bbslover wrote: > > as known, svm need tune some parameters like cost,gamma and epsilon to get > better performance,but one question appear, how can i monitor the > performance . generally speaking ,we chose the cross-validation MSE in the > training set, but It seems svm can not return the cross-validation MSE > value, we just get it from "summary model.svm", if I write a loop, have no > idear call the cros-validate MSE, and no way to monitor this performance > ,how can I do? > > > -- > View this message in context: > http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p974942.html > Sent from the R help mailing list archive at Nabble.com. > > __ > 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, minimal, self-contained, reproducible code. > -- Max __ 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, minimal, self-contained, reproducible code.
[R] Hello all, How can I get corss-validation MSE of SVM in e1071?
as known, svm need tune some parameters like cost,gamma and epsilon to get better performance,but one question appear, how can i monitor the performance . generally speaking ,we chose the cross-validation MSE in the training set, but It seems svm can not return the cross-validation MSE value, we just get it from "summary model.svm", if I write a loop, have no idear call the cros-validate MSE, and no way to monitor this performance ,how can I do? -- View this message in context: http://n4.nabble.com/Hello-all-How-can-I-get-corss-validation-MSE-of-SVM-in-e1071-tp974942p974942.html Sent from the R help mailing list archive at Nabble.com. __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello
dimple thyagarajan wrote: Hello! I am trying to merge two xy-plot with different ylimits. It seems that you are using lattice graphics (which you do not mention), if so, look at the documentation of xyplot, specifically to the panel argument. cheers, Paul Can someone please give me possible way of achieving it.. Thank you Regards Dimple [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code. -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello
On Aug 19, 2009, at 6:24 PM, dimple thyagarajan wrote: Hello! I am trying to merge two xy-plot with different ylimits. Can someone please give me possible way of achieving it.. You cannot possibly be using a function named xy-plo. The particular function you have chosen will be critical to an appropriate answer. In all likelihood, the approach will be to plot with ylimits set by the maximum and minimum of the combined set of values. It's possible that the plotting function will allow group arguments. Beyond that it is hard to be of further help without an example and specifics. Note: An informative subject line is requested by the Posting Guide. Thank you Regards Dimple David Winsemius, MD Heritage Laboratories West Hartford, CT __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello
generate the plot with the limit of both ranges and then add the second with 'lines' plot(x1,y1,ylim=range(y1, y2)) lines(x2,y2) On Wed, Aug 19, 2009 at 6:24 PM, dimple thyagarajan wrote: > Hello! > > I am trying to merge two xy-plot with different ylimits. > > Can someone please give me possible way of achieving it.. > > Thank you > Regards > Dimple > > > > [[alternative HTML version deleted]] > > > __ > 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, minimal, self-contained, reproducible code. > > -- Jim Holtman Cincinnati, OH +1 513 646 9390 What is the problem that you are trying to solve? __ 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, minimal, self-contained, reproducible code.
[R] Hello
Hello! I am trying to merge two xy-plot with different ylimits. Can someone please give me possible way of achieving it.. Thank you Regards Dimple [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello,
Kiran exonmap is a bioconductor package - I'd suggest trying the bioconductor mailing list. Also, you're more likely to get a response if you include an informative subject line. Regards Richard. Kiran Annaiah wrote: Hello, A newbie to R. I am trying to use the exonmap package in R. According to the docs, the xmapDatabase() command should read the config file with all the connection parameters and connect to the DB. But i get the error shown below xmapDatabase("Human") Error in mysqlNewConnection(dbDriver(drv), ...) : RS-DBI driver: (could not connect [EMAIL PROTECTED] on dbname "NA" Error:Access denied for user 'cagadmin'@'cdev.cag.chop.edu' to database 'NA' ) when i connect to db using dbConnect() it works fine. conn<-dbConnect(drv="MySQL",username="cagadmin" ) Anybody who can shed some light on whats going on here? thank you all kiran Kiran Annaiah Center for Applied Genomics Informatics Children's Hospital of Philadelphia 267-426-9355 __ 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, minimal, self-contained, reproducible code. -- Richard D. Pearson [EMAIL PROTECTED] School of Computer Science,http://www.cs.man.ac.uk/~pearsonr University of Manchester, Tel: +44 161 275 6178 Oxford Road, Mob: +44 7971 221181 Manchester M13 9PL, UK.Fax: +44 161 275 6204 __ 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, minimal, self-contained, reproducible code.
[R] Hello,
Hello, A newbie to R. I am trying to use the exonmap package in R. According to the docs, the xmapDatabase() command should read the config file with all the connection parameters and connect to the DB. But i get the error shown below > xmapDatabase("Human") Error in mysqlNewConnection(dbDriver(drv), ...) : RS-DBI driver: (could not connect [EMAIL PROTECTED] on dbname "NA" Error:Access denied for user 'cagadmin'@'cdev.cag.chop.edu' to database 'NA' ) when i connect to db using dbConnect() it works fine. conn<-dbConnect(drv="MySQL",username="cagadmin" ) Anybody who can shed some light on whats going on here? thank you all kiran Kiran Annaiah Center for Applied Genomics Informatics Children's Hospital of Philadelphia 267-426-9355 __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello
"sumit gupta" <[EMAIL PROTECTED]> napsal dne 11.06.2008 07:53:54: > Hello > thanx for the info.. > > I still have tthe doubt regarding the heat map. > I am attching a plot of a 70X3 data matrix.In which 3rd variable has been > shown as colour of different boxes.Could you help me to draw this type of graph using R Well, as you did not provided some data mat<-cbind(sample(1:5, 70, replace=T), sample(1:5, 70, replace=T), rnorm(70)) dim(mat) [1] 70 3 library(akima) filled.contour(interp(mat[,1], mat[,2], mat[,3], duplicate="mean")) image(interp(mat[,1], mat[,2], mat[,3], duplicate="mean")) is this what you want? Regards Petr > > Regards, > Sumit > > > On 6/10/08, Petr PIKAL <[EMAIL PROTECTED]> wrote: > Hi > > [EMAIL PROTECTED] napsal dne 10.06.2008 14:49:55: > > > Hello, > > > > I am facing a problem in drawing heat map using R. > > > > I have a 70X3 matrix and I want to draw a heat map with 1 coloumn on X > axis > > another on Y axis and want to show the value of 3rd coloumn using heat > > colours. > > Could you please help me with this . > > You can use function interp from akima package > > Petr > > > Thank you . > > > > Regards, > > Sumit > > > >[[alternative HTML version deleted]] > > > > __ > > 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, minimal, self-contained, reproducible code. > > > > -- > Sumit Kumar Gupta > 4thYear ,Biotechnology > IIT Madras [příloha rfm.bmp odstraněna uživatelem Petr PIKAL/CTCAP] __ 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, minimal, self-contained, reproducible code.
Re: [R] Hello
Look at the "image" function and the "levelplot" function in the lattice package. -Original Message- From: "sumit gupta" <[EMAIL PROTECTED]> To: "r-help@r-project.org" Sent: 6/10/08 6:55 AM Subject: [R] Hello Hello, I am facing a problem in drawing heat map using R. I have a 70X3 matrix and I want to draw a heat map with 1 coloumn on X axis another on Y axis and want to show the value of 3rd coloumn using heat colours. Could you please help me with this . Thank you . Regards, Sumit [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code. __ 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, minimal, self-contained, reproducible code.
[R] Hello
Hello, I am facing a problem in drawing heat map using R. I have a 70X3 matrix and I want to draw a heat map with 1 coloumn on X axis another on Y axis and want to show the value of 3rd coloumn using heat colours. Could you please help me with this . Thank you . Regards, Sumit [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
Re: [R] hello! i need help for a specific graphic issue
quick (and dirty) solution: y.up<- means+stand.error y.dwn<- means-stand.error plot(means,ylim=c(3.5,10)) for (i in 1:length(means)) arrows(i,means[i],i,y.up[i],length=0.1) for (i in 1:length(means)) arrows(i,means[i],i,y.dwn[i],length=0.1) Giacomo Prodi wrote: > > hello, ladyes and gentlemans. > > check this: > > means<-c(4,6,8) > stand.error<-c(0.1,0.3,0.5) > > now i've strongly tryed to scatterplot the > means(y-axis),by showing their sd with the > arrow(..,code=3,angle=90) function. > The problem is that my x-axis has categorical values > (say, factor(x)), and the arrows() can't recognize > them as right coordinates. > ? > thank you all in advance > B.F. insubria university (varese) > > > > La web mail più usata al mondo. > http://it.docs.yahoo.com/mail/overview/index.html > > __ > 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, minimal, self-contained, reproducible code. > > -- View this message in context: http://www.nabble.com/hello%21-i-need-help-for-a-specific-graphic-issue-tp16002183p16002532.html Sent from the R help mailing list archive at Nabble.com. __ 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, minimal, self-contained, reproducible code.
Re: [R] hello! i need help for a specific graphic issue
I don't think that you're actually creating a scatterplot if the x-axis is a factor. You're getting a boxplot. If you want a scatterplot then something like this might work. == means<-c(4,6,8) stand.error<-c(0.1,0.3,0.5) aa <- factor(c("A","B","C","A", "C", "C")) dd <- 1:length(levels(aa)) plot(dd, means, xaxt="n") arrows(dd, means + stand.error, dd, means-stand.error, code=3, angle=90) axis(1, at=dd, labels=levels(aa)) === --- Giacomo Prodi <[EMAIL PROTECTED]> wrote: > hello, ladyes and gentlemans. > > check this: > > means<-c(4,6,8) > stand.error<-c(0.1,0.3,0.5) > > now i've strongly tryed to scatterplot the > means(y-axis),by showing their sd with the > arrow(..,code=3,angle=90) function. > The problem is that my x-axis has categorical values > (say, factor(x)), and the arrows() can't recognize > them as right coordinates. > ? > thank you all in advance > B.F. insubria university (varese) > > > > La web mail più usata al mondo. > http://it.docs.yahoo.com/mail/overview/index.html > > __ > 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, minimal, self-contained, > reproducible code. > __ 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, minimal, self-contained, reproducible code.
Re: [R] hello! i need help for a specific graphic issue
On Wed, Mar 12, 2008 at 5:37 AM, Giacomo Prodi <[EMAIL PROTECTED]> wrote: > hello, ladyes and gentlemans. > > check this: > > means<-c(4,6,8) > stand.error<-c(0.1,0.3,0.5) > > now i've strongly tryed to scatterplot the > means(y-axis),by showing their sd with the > arrow(..,code=3,angle=90) function. > The problem is that my x-axis has categorical values > (say, factor(x)), and the arrows() can't recognize > them as right coordinates. > ? > thank you all in advance > B.F. insubria university (varese) This is a bit easier to do with ggplot2: df <- data.frame( trt = factor(c("a", "b", "c")), mean = c(4,6,8), se = c(0.1,0.3,0.5) ) install.packages("ggplot2") library(ggplot2) qplot(trt, mean, data=df) qplot(trt, mean, data=df, min=mean - se, max = mean + se, geom="pointrange") Hadley -- http://had.co.nz/ __ 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, minimal, self-contained, reproducible code.
Re: [R] hello! i need help for a specific graphic issue
hi Giacomo, consider this: means<-c(4,6,8) stand.error<-c(0.1,0.3,0.5) lowlim<-means-stand.error uplim<-means+stand.error plot(as.factor(c("a","b","c")),means,ylim=c(min(lowlim),max(uplim))) arrows(c(1:3),lowlim,c(1:3),uplim,code=3,angle=90) Still, its just a workaround and I don't think you can get rid of the horizontal mean-lines because this is actually a box-plot. Otherwise you could try the 'errbar()' from the Hmisc package. I never used that, but maybe it works. Generally, if you have problems with plot-dimensions, you can get instant information with 'locator()', which returns coordinates of mouse clicked locations in plots. cheers, /rené > -Ursprüngliche Nachricht- > Von: "Giacomo Prodi" <[EMAIL PROTECTED]> > Gesendet: 12.03.08 11:56:45 > An: r-help@r-project.org > Betreff: [R] hello! i need help for a specific graphic issue > > hello, ladyes and gentlemans. > > check this: > > means<-c(4,6,8) > stand.error<-c(0.1,0.3,0.5) > > now i've strongly tryed to scatterplot the > means(y-axis),by showing their sd with the > arrow(..,code=3,angle=90) function. > The problem is that my x-axis has categorical values > (say, factor(x)), and the arrows() can't recognize > them as right coordinates. > ? > thank you all in advance > B.F. insubria university (varese) > > > > La web mail più usata al mondo. > http://it.docs.yahoo.com/mail/overview/index.html > > __ > 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, minimal, self-contained, reproducible code. > ___ Schon gehört? Der neue WEB.DE MultiMessenger kann`s mit allen: http://www.produkte.web.de/messenger/?did=3015 __ 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, minimal, self-contained, reproducible code.
[R] hello! i need help for a specific graphic issue
hello, ladyes and gentlemans. check this: means<-c(4,6,8) stand.error<-c(0.1,0.3,0.5) now i've strongly tryed to scatterplot the means(y-axis),by showing their sd with the arrow(..,code=3,angle=90) function. The problem is that my x-axis has categorical values (say, factor(x)), and the arrows() can't recognize them as right coordinates. ? thank you all in advance B.F. insubria university (varese) La web mail più usata al mondo. http://it.docs.yahoo.com/mail/overview/index.html __ 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, minimal, self-contained, reproducible code.
Re: [R] hello
Hi, > I want to know if it's possible to send a bar chart from the > software R to an directory. Do you mean saving your figure to a file ? If so, you can take a look at the device functions, which allow you to generate your figure to a file in different formats, especially pdf or postscript. For example, to generate your plot into a pdf file, you can use : pdf(file="filename.pdf") plot(...) dev.off() Hope that helps, Julien -- Julien Barnier Groupe de recherche sur la socialisation ENS-LSH - Lyon, France __ 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, minimal, self-contained, reproducible code.
Re: [R] hello
See the examples under ?png ?pdf On Jan 24, 2008 10:05 AM, elyakhlifi mustapha <[EMAIL PROTECTED]> wrote: > hello, > happy new year. > I want to know if it's possible to send a bar chart from the software R to > an directory. > Thanks. > > > > _ > > >[[alternative HTML version deleted]] > > __ > 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, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.
[R] hello
hello, happy new year. I want to know if it's possible to send a bar chart from the software R to an directory. Thanks. _ [[alternative HTML version deleted]] __ 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, minimal, self-contained, reproducible code.