similar to: plot empirical pdf

Dear R-people, are there any way to change the size of legend in barplot? I have tried various versions of cex, both as par(cex.*= ) and barplot(..., cex.*= ). So long without success. Sincerely, Ott > version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1

Hi, I have an RH9 box (gcc (GCC) 3.2.2, perl 5.8.0, make 3.79.1). My locale is LANG="et_EE.UTF-8" I do: $ ./configure --prefix=/usr/local ... $ make ... make[4]: Entering directory `/home/otoomet/a/R-1.9.0/src/library/stats4' dumping R code in package 'stats4' Error in structure(c(unlist(lapply(list(...), unclass))), class = c("POSIXt", : couldn't

Hi, I have a following vector: > smallch [1] 0.0652840 0.1181300 0.0319370 0.0155700 0.0464110 0.0107850 [7] 0.0158970 0.0375900 0.0603090 0.0310300 0.0105920 0.0540580 [13] -0.0177740 0.0039393 Pretty (R 1.5.1) has problems with zero: > pretty(smallch) [1] -2.000000e-02 -3.469447e-18 2.000000e-02 4.000000e-02 6.000000e-02 [6] 8.000000e-02 1.000000e-01 1.200000e-01

Hi, I have a following vector: > smallch [1] 0.0652840 0.1181300 0.0319370 0.0155700 0.0464110 0.0107850 [7] 0.0158970 0.0375900 0.0603090 0.0310300 0.0105920 0.0540580 [13] -0.0177740 0.0039393 Pretty (R 1.5.1) has problems with zero: > pretty(smallch) [1] -2.000000e-02 -3.469447e-18 2.000000e-02 4.000000e-02 6.000000e-02 [6] 8.000000e-02 1.000000e-01 1.200000e-01

Hi everybody, I am using gzipped files quite extensively and cannot figure out how I can close (or reuse) connections. Example: > for(i in 1:300) {cat(i, "\n"); a <- readLines(gzfile("file.gz"), n=1) } 1 2 ... 46 47 Error in gzfile("~/tyyq/andmebaasiq/ETU/ETU1997.tsv.gz") : all connections are in use To use just gzfile without any open/close should

Does R do cumulative frequency distribution plots? -- Tommy E. Cathey, Senior Scientific Application Consultant High Performance Computing & Scientific Visualization SAIC, Supporting the EPA Research Triangle Park, NC 919-541-1500 EMail: cathey.tommy at epa.gov My e-mail does not reflect the opinion of SAIC or the EPA. Federal Contact - John B. Smith 919-541-1087 - smith.johnb at epa.gov

Hi, I've estimated a simple kernel density of a univariate variable with density(), but after I would like to find out the CDF at specific values. How can I do it? thanks for your help, with it I am very close to finish my first little bit more serious work in R, Viktor

I have written the following function to calculate the weighted mean difference for univariate data (see http://www.xycoon.com/gini_mean_difference.htm for a related formula). Unsurprisingly, the function is slow (compared to sd or mad) for long vectors. I wonder if there's a way to make the function faster, short of creating an external C function. Thanks very much for your advice. gmd

In the document "R Data Import/Export", section "Output to connections", there is the following portion of code: ## convert decimal point to comma in output, using a pipe (Unix) zz <- pipe(paste("sed s/\\./,/ >", "outfile"), "w") cat(format(round(rnorm(100), 4)), sep = "\n", file = zz) close(zz) ## now look at the output

I have justed upgraded R v1.7.0 on Windows NT 4 and have installed the latest RMySQL (version 0.5-1)and DBI (version 0.1-5) packages. When I issue the following commands (tactfully adjusted) R just crashes and disappears, any ideas? require(RMySQL) m <- dbDriver("MySQL") con <- dbConnect(m, dbname="xxx", user="xxx", password="xxx",

Hi all, I'd like to quantize a variable to map it into a limited set of integers for use with a colormap. "image" and filled.contour" do this mapping inside somewhere, but I'd like to choose the colors for plotting a set of polygons. Is there a pre-existing function that does something like this well? i.e., is capable of using 'breaks'?

Hello, How can I obtain the likelihood ratio of a Poisson regression model? Regards. _____________________________________________ dr. Marziliano Ciro Facolta' di Economia Universita' degli Studi di L'Aquila p.zza del Santuario, 19 67040 Roio Poggio, L'Aquila tel.: 0862 434836 fax: 0862 434803

I am learning glm function, but how do you fit data using exponential distribution with glm? In the help file, under "Family Objects for Models", no ready made option seems available for the distribution as well as for other distributions satisfying GLM requirements not listed there.

Dear list, Does anyone know how to fit the power law distribution? I have the empirical distribution and would like to check whether it fits the power law (with the power estimated from the data). Any hints are appreciated. Tanks a lot! Galina [[alternative HTML version deleted]]

Hello... Is it possible to use "density" or another kernel density estimator to identify the mode of a distribution? When I use 'density', the resulting density plot of my data is much cleaner than the original noisy histogram, and I can clearly see the signal that I am interested in. E.g., suppose my data is actually drawn from two or more normal (or other)

Hi R-users,  I have a data from multi skew t and would like to transform each of the data to uniform data.  I tried using 'pmst' but only got one output:   > rr1 <- as.vector(r1);rr1  [1]  0.7207582  5.2250906  1.7422237  0.5677233  0.7473555 -0.6020626 -2.1947872 -1.1128313 -0.6587316 -1.1409261     > pmst(rr1, xi=rep(0,10), Omega=diag(10), alpha=rep(1,10), df=5) [1] 3.676525e-09

Dear R-help list members, I have created a package 'mnormt' with facilities for the multivariate normal and t distributions. The core part is simply an interface to Fortran routines by Alan Genz for computing the integral of two densities over rectangular regions, using an adaptive integration method. Other R functions compute densities and generate random numbers. The starting

> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Caio Lucidius > Naberezny Azevedo > Sent: 01 September 2005 12:09 > To: Help mailing list - R > Subject: [R] Multivariate Skew Normal distribution > > > Hi all, > > Could anyone tell me if there is any package (or function)

Is there code for bivariate kernel density estimation? For bivariate kernels there is kde2d in MASS kde2d.g in GRASS KernSur in GenKern (list probably incomplete) but none of them seems to accept a weight parameter (like density does since R 2.2.0) -- Erich Neuwirth, University of Vienna Faculty of Computer Science Computer Supported Didactics Working Group Visit our SunSITE at

Hello! I am writing some numbers and character vectors to an ascii file and would like to get rid of [1] ... as shown bellow (a dummy example) R> runif(20) [1] 0.653574 0.164053 0.036031 0.127208 0.134274 0.103252 0.506480 0.547759 [9] 0.912421 0.584382 0.987208 0.996846 0.666760 0.053637 0.327590 0.370737 [17] 0.505706 0.412316 0.887421 0.812151 I have managed to work up to remove quotes