Re: [R] Hello

2022-03-08 Thread 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.


[R] Hello

2022-03-08 Thread Shrinivas Dharma

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

2020-12-14 Thread 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

2020-12-14 Thread 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

2020-12-14 Thread 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.


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

2020-12-14 Thread Nikolaos Anagnostakis
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

2016-10-29 Thread terry minton
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

2016-09-08 Thread 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

2016-09-08 Thread 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.


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

2016-09-07 Thread darth brando
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

2016-01-06 Thread 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.


[R] hello everyone and happy new year

2016-01-06 Thread Matthias Worni
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!

2013-12-17 Thread 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 R user!

2013-12-17 Thread 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  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!

2013-12-17 Thread Sarah Goslee
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!

2013-12-17 Thread 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  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!

2013-12-17 Thread 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  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!

2013-12-17 Thread 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.


[R] Hello R user!

2013-12-17 Thread bibek sharma
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

2013-05-01 Thread 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.


[R] Hello

2013-05-01 Thread Eric Mintah
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

2013-04-28 Thread 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,   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

2013-04-28 Thread ampemdo

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

2013-02-15 Thread Tanushree Tunstall
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

2012-12-14 Thread 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 
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

2012-12-14 Thread 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 
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

2012-12-14 Thread 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

2012-12-14 Thread 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

2012-12-14 Thread 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 
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

2012-12-14 Thread 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

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

2011-06-19 Thread Duncan Murdoch

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

2011-06-19 Thread p_connolly

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

2011-06-19 Thread 王海生
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!

2011-05-30 Thread 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
 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!

2011-05-29 Thread Amanda Zeqiri

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!

2011-02-28 Thread 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!

2011-02-27 Thread Joshua Wiley
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!

2011-02-27 Thread 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!

2011-02-27 Thread 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" 
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!

2011-02-27 Thread 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 
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!

2011-02-27 Thread Mark Knecht
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!

2011-02-27 Thread 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  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!

2011-02-27 Thread Jim Lemon

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!

2011-02-26 Thread 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  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!

2011-02-26 Thread tsvi sabo
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?

2009-12-21 Thread 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?

2009-12-20 Thread 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]"  
写道: -原始邮件-
发件人:"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?

2009-12-19 Thread 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  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?

2009-12-18 Thread 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  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?

2009-12-18 Thread bbslover

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

2009-08-20 Thread 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

2009-08-19 Thread 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

2009-08-19 Thread 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 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

2009-08-19 Thread dimple thyagarajan
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,

2008-08-12 Thread 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.


[R] Hello,

2008-07-30 Thread Kiran Annaiah
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

2008-06-11 Thread 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

2008-06-10 Thread 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" 
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

2008-06-10 Thread sumit gupta
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

2008-03-12 Thread bartjoosen

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

2008-03-12 Thread 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! i need help for a specific graphic issue

2008-03-12 Thread 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

2008-03-12 Thread 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.


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

2008-03-12 Thread Giacomo Prodi
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

2008-01-24 Thread 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.


Re: [R] hello

2008-01-24 Thread Jared O'Connell
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

2008-01-24 Thread elyakhlifi mustapha
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.