Re: [R] Hello

[Ivan Krylov]

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.

Re: [R] Hello! Problem - function ''matord'' in R

[Jeff Newmiller]

... 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

[T. A. Milne via R-help]

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

[Bert Gunter]

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.

Re: [R] hello i have a question on music analysis and mathematical synthesis related to r code

[Bert Gunter]

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

[Jim Lemon]

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.

Re: [R] hello everyone and happy new year

[Ulrik Stervbo]

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.

Re: [R] Hello R user!

[John Kane]

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.

Re: [R] Hello R user!

[Christopher W Ryan]

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 mbhpat...@gmail.com 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!

[Sarah Goslee]

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 mbhpat...@gmail.com 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!

[Sarah Goslee]

On Tue, Dec 17, 2013 at 1:04 PM, bibek sharma mbhpat...@gmail.com 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 sarah.gos...@gmail.com
 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 mbhpat...@gmail.com
 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!

[Christopher W Ryan]

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 sarah.gos...@gmail.com wrote:
 On Tue, Dec 17, 2013 at 1:04 PM, bibek sharma mbhpat...@gmail.com 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 sarah.gos...@gmail.com
 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 mbhpat...@gmail.com
 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!

[Jim Lemon]

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

[Rui Barradas]

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.

Re: [R] Hello. PLEASE HELP

[Bert Gunter]

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,  ampe...@stud.ntnu.no 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.

Re: [R] Hello R User

[arun]

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 mbhpat...@gmail.com
To: R help r-help@r-project.org
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

[Jessica Streicher]

 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

[Eik Vettorazzi]

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

[arun]

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 mbhpat...@gmail.com
To: R help r-help@r-project.org
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

[arun]

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 mbhpat...@gmail.com
To: arun smartpink...@yahoo.com
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 smartpink...@yahoo.com 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 mbhpat...@gmail.com
 To: R help r-help@r-project.org
 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 about proxy configuration

[p_connolly]

Quoting 王海生 edwin_ue...@163.com:


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.

Re: [R] hello about proxy configuration

[Duncan Murdoch]

On 19/06/2011 7:29 PM, p_conno...@slingshot.co.nz wrote:

Quoting 王海生edwin_ue...@163.com:


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!

[Dennis Murphy]

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
amandazeq...@hotmail.com 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.

Re: [R] Hello!

[Mike Marchywka]

 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!

[Jim Lemon]

tsvi sabotsvis...@yahoo.com  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 Return 
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!

[Rob Tirrell]

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 gunter.ber...@gene.com 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 tsvis...@yahoo.com 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!

[Mark Knecht]

On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter gunter.ber...@gene.com 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!

[Rios,Alfredo Arturo]

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 markkne...@gmail.com
wrote:
 On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter gunter.ber...@gene.com 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!

[Ravi Varadhan]

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 ar...@ufl.edu
Date: Sunday, February 27, 2011 12:40 pm
Subject: Re: [R] Hello!
To: Mark Knecht markkne...@gmail.com
Cc: r-help@r-project.org, Bert Gunter gunter.ber...@gene.com


 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 markkne...@gmail.com
  wrote:
   On Sat, Feb 26, 2011 at 10:11 PM, Bert Gunter 
 gunter.ber...@gene.com 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!

[Mike Marchywka]

 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!

[Joshua Wiley]

On Sun, Feb 27, 2011 at 9:16 AM, Mike Marchywka marchy...@hotmail.com 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!

[Bert Gunter]

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 tsvis...@yahoo.com 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.

Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071?

[Max Kuhn]

  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?

[bbslover]

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] ml-node+975038-1351220...@n4.nabble.com 
写道： -原始邮件-
发件人:Max Kuhn [via R] ml-node+975038-1351220...@n4.nabble.com
发送时间:2009年12月19日 星期六
收件人:bbslover dlu...@yeah.net
主题: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?

[bbslover]

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 dlu...@yeah.net 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?

[Max Kuhn]

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 dlu...@yeah.net 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.

Re: [R] Hello

[Paul Hiemstra]

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

[jim holtman]

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 thyagarajancalldim...@yahoo.com 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.

Re: [R] Hello

[David Winsemius]

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,

[Richard Pearson]

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.

Re: [R] Hello

[Petr PIKAL]

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

[Greg Snow]

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 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.

Re: [R] hello! i need help for a specific graphic issue

[René Capell]

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.

Re: [R] hello! i need help for a specific graphic issue

[hadley wickham]

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

[John Kane]

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

[Julien Barnier]

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.