[R] How can I simulate Pareto distribution in R?

2005-01-09 Thread Ni Xiao
Hi, guys, I need to simulate Pareto distribution. But I found 'rpareto' didn't exist in R. And it seems that Pareto distribution don't have mathematical relationships with other distributions. What can I do? Thanks a lot. Ni -

Re: [R] Does R accumulate memory

2005-01-09 Thread Prof Brian Ripley
On Sat, 8 Jan 2005, Duncan Murdoch wrote: On Sat, 8 Jan 2005 16:38:31 -0500, Doran, Harold [EMAIL PROTECTED] wrote: I am running into a memory issue that I haven't noticed before. I am running a simulation with all of the code used below. I have increased my memory to 712mb and have a total of 1

Re: [R] How can I simulate Pareto distribution in R?

2005-01-09 Thread Remigijus Lapinskas
If you define Pareto density as p(x)=c*x^(-(c+1)) for x1, then dpareto - function(x,c){ if(c=0)stop(c must be positive) # Diagnostic step ifelse(x1,0,c/x^(c+1))} ppareto - function(q,c){ if(c=0)stop(c must be positive 0) ifelse(q1,0,1-1/q^c)} qpareto - function(p,c){ if(c=0) stop(c must be

Re: [R] How can I simulate Pareto distribution in R?

2005-01-09 Thread Göran Broström
On Sun, Jan 09, 2005 at 12:21:47AM -0800, Ni Xiao wrote: Hi, guys, I need to simulate Pareto distribution. But I found 'rpareto' didn't exist in R. And it seems that Pareto distribution don't have mathematical relationships with other distributions. What can I do? Well, it has a relation

[R] Re: labels attached to variable names

2005-01-09 Thread Denis Chabot
Thank you to the many helpful list members who answered. I am grateful to have been told about attr and about Hmisc. Actually Hmisc also brings me solutions to many things on my wish list! Cheers, Denis Le Vendredi, 7 janv 2005, à 12:11 Europe/Paris, [EMAIL PROTECTED] a écrit : De: Denis

Re: [R] R-etiquette

2005-01-09 Thread Frank E Harrell Jr
Anne wrote: I'm about to present a report (for internal use of governmental agency). I used extensively R , contibuted packages, as well as communications on the R-list As well as citing R, I would like to know how to cite the contributed packages (it is not so easy, as some have been used

Re: [R] R-etiquette

2005-01-09 Thread Prof Brian Ripley
On Sun, 9 Jan 2005, Anne wrote: I'm about to present a report (for internal use of governmental agency). I used extensively R , contibuted packages, as well as communications on the R-list As well as citing R, I would like to know how to cite the contributed packages (it is not so easy, as

Re: [R] R-etiquette

2005-01-09 Thread Frank E Harrell Jr
Prof Brian Ripley wrote: On Sun, 9 Jan 2005, Anne wrote: I'm about to present a report (for internal use of governmental agency). I used extensively R , contibuted packages, as well as communications on the R-list As well as citing R, I would like to know how to cite the contributed packages

Re: [R] autoscaling plot font size in Sweave output possible?

2005-01-09 Thread Patrick Drechsler
Friedrich Leisch wrote on 08 Jan 2005 11:52:41 MET: [...] To me it sounds much more like you're one of those guys who like it when others solve their problems for free. I gave you an IMO perfectly working cookbook recipe and you didn't even bother to tell me what did not work out for you,

Re: [R] R-etiquette

2005-01-09 Thread Prof Brian Ripley
On Sun, 9 Jan 2005, Frank E Harrell Jr wrote: Prof Brian Ripley wrote: On Sun, 9 Jan 2005, Anne wrote: I'm about to present a report (for internal use of governmental agency). I used extensively R , contibuted packages, as well as communications on the R-list As well as citing R, I would like

Re: [R] R-etiquette

2005-01-09 Thread Gabor Grothendieck
Prof Brian Ripley ripley at stats.ox.ac.uk writes: : : On Sun, 9 Jan 2005, Anne wrote: : : I'm about to present a report (for internal use of governmental agency). : I used extensively R , contibuted packages, as well as communications on : the R-list : : As well as citing R, I would

Re: [R] R-etiquette

2005-01-09 Thread Romain François
hello, package MASS has one that looks just like that : citHeader(To cite the VR bundle (MASS, class, nnet, spatial) in publications use:) citEntry(entry=Book, title = Modern Applied Statistics with S, author = personList(as.person(W. N. Venables),

[R] plot.default and open ended limits

2005-01-09 Thread Dan Bolser
Hi, I would like to bound the lower limit of my y scale to zero, and let R chose an upper limit. Something like plot(x,ylim=c(0,)) or plot(x,ylim=c(0,na)) but nither of these do the job. I searched the docs, but I can't see a way to do this. Naturally its nothing I can't do 'by hand', I

Re: [R] plot.default and open ended limits

2005-01-09 Thread Uwe Ligges
Dan Bolser wrote: Hi, I would like to bound the lower limit of my y scale to zero, and let R chose an upper limit. Something like plot(x,ylim=c(0,)) or plot(x,ylim=c(0,na)) but nither of these do the job. I searched the docs, but I can't see a way to do this. Naturally its nothing I can't do 'by

Re: [R] plot.default and open ended limits

2005-01-09 Thread Romain François
Hello Dan, Look at the code of the plot.default function, you'll see that's not possible to specify one limit, nevertheless, you can do : plot(x,y,ylim=c(0,max(y))) Does it work for you ? Cordialement. Romain. Dan Bolser a écrit : Hi, I would like to bound the lower limit of my y scale to zero,

Re: [R] R-etiquette

2005-01-09 Thread Brian D Ripley
On Sun, 9 Jan 2005, Gabor Grothendieck wrote: Prof Brian Ripley ripley at stats.ox.ac.uk writes: : : On Sun, 9 Jan 2005, Anne wrote: : : I'm about to present a report (for internal use of governmental agency). : I used extensively R , contibuted packages, as well as communications on :

Re: [R] R-etiquette

2005-01-09 Thread Dirk Eddelbuettel
On Sun, Jan 09, 2005 at 04:26:05PM +, Gabor Grothendieck wrote: Prof Brian Ripley ripley at stats.ox.ac.uk writes: : See the citation() function, which works for packages too and a few : authors have taken advantage of the means of customizing it. [...] Could someone point out a few

[R] dist{amap} error??

2005-01-09 Thread Radu Dobrin
Dear all, I have come across a very confusing matter regarding dist() supplied by the amap package: --- m is just a test matrix library(amap) Loading required package: mva Warning message: package 'mva' has been merged into 'stats' m a b c aa 0.1 0.2 0.3 bb 2.0 3.0 4.0 cc 2.0 4.0

Re: [R] plot.default and open ended limits

2005-01-09 Thread Dan Bolser
On Sun, 9 Jan 2005, [ISO-8859-1] Romain François wrote: Hello Dan, Look at the code of the plot.default function, you'll see that's not possible to specify one limit, nevertheless, you can do : Suppose I wanted to contribute a 'fix' to the code to allow one of my suggested syntax solutions

Re: [R] plot.default and open ended limits

2005-01-09 Thread Romain François
Dan Bolser a écrit : Hello Dan, Look at the code of the plot.default function, you'll see that's not possible to specify one limit, nevertheless, you can do : Suppose I wanted to contribute a 'fix' to the code to allow one of my suggested syntax solutions below, how would I go about it?

Re: [R] plot.default and open ended limits

2005-01-09 Thread Thomas Lumley
On Sun, 9 Jan 2005, Dan Bolser wrote: On Sun, 9 Jan 2005, [ISO-8859-1] Romain François wrote: Hello Dan, Look at the code of the plot.default function, you'll see that's not possible to specify one limit, nevertheless, you can do : Suppose I wanted to contribute a 'fix' to the code to allow one of

[R] how to do by-processing using weighted.mean

2005-01-09 Thread Denis Chabot
Hi, I'd like to compute the weighted mean of some variables, all using the same weight variable, for each combination of 3 factor variables. I found how I could use summarize (from Hmisc) to do normal means for combinations of 3 factors, but I cannot find a way of doing weighted means. Is it

Re: [R] plot.default and open ended limits

2005-01-09 Thread Dan Bolser
On Sun, 9 Jan 2005, Thomas Lumley wrote: On Sun, 9 Jan 2005, Dan Bolser wrote: On Sun, 9 Jan 2005, [ISO-8859-1] Romain François wrote: Hello Dan, Look at the code of the plot.default function, you'll see that's not possible to specify one limit, nevertheless, you can do : Suppose I

Re: [R] how to do by-processing using weighted.mean

2005-01-09 Thread Frank E Harrell Jr
Denis Chabot wrote: Hi, I'd like to compute the weighted mean of some variables, all using the same weight variable, for each combination of 3 factor variables. I found how I could use summarize (from Hmisc) to do normal means for combinations of 3 factors, but I cannot find a way of doing

[R] (no subject)

2005-01-09 Thread fuerjing . com
ÏîÄ¿Ñ°×Ê ÏîÄ¿Ñ°×Ê ¶À¼Ò²Í×ÀÃÀ棬Óû¿ª×¨Âôµê200¼Ò£¬ÎåÄêÀûÈó2000ÍòÔªÒÔÉÏ£¬¼±ÐèÇ°ÆÚͶ×Ê20Íò¡£ºÏ×÷·½Ê½Ïê¼ûºÏͬ¼°ÒâÏò¡££¨²úÆ·Ôݲ»¹«¿ª£©¿ÉË÷×Ê ÁÏ £¬ÎÞ³ÏÒâÎðÈÅ£¡ µç»°£º69827192 83476289 __ R-help@stat.math.ethz.ch mailing list

RE: [R] Basic Linear Algebra

2005-01-09 Thread Mulholland, Tom
Thanks. You've confirmed that I really don't understand the language of mathematics, even if I understand how numbers work. I will keep plodding with my books until I can understand the implications for singularity in what I am trying to achieve. Tom -Original Message- From: Peter