similar to: R crash on Mac

AA> In 2008, I have spotted some errors in a package, one which is AA> likely to have many users (I am not one myself). The more serious AA> errors are in the documentation, since they lead to a completely AA> distorted interpretation of the outcome; in addition, there is (at AA> least) one programming error which produces some wrong AA> computations. A few weeks later, the

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

Dear List-members, to download a file from the net, the function download.file(..) does the job. However, before embarking on the download, I would like to find out how large the file is. Is there a way to know it? Most easily, this question has been asked before, but I am new to the list. Regards, with thanks in advance, Adelchi Azzalini ---- Adelchi Azzalini <azzalini at

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

There seems to exist peculiar cases where optim does not take care of constraints on the parameters to be optimized over. The call to optim is of the form opt <- optim(cp, fn=sn.dev, gr=sn.dev.gh, method="L-BFGS-B", lower=c(-Inf, 1e-10, -0.99527), upper=c( Inf, Inf, 0.99527), control=control, X=X, y=y, hessian=FALSE) The code has worked fine

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@r-project.org # ###################################################### Today I have installed R 1.02.1 on my MSW-95 laptop; it essentially worked, but thre are two

Dear R helpers, I need to find out maximized log likelihoods, parameters estimates and standard errors (in parentheses) of r largest-order statistics model, with different values of r by using the function rlarg.fit. I want to specify required number of order statistics to the model. I attached my data file with this mail.please help me. Ruposh --- r-help-request at stat.math.ethz.ch wrote:

Dear Prof Azzalini, You have an interesting example of recursive use of uniroot(). [re-cited at the end] However, note that R currently has the same problem as S-plus: Uniroot() doesn't work reliably, recursively. When you found it to be better, you were just lucky. The relevant file in R's source, src/main/optimize.c says /* WARNING : As things stand, these routines should

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

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

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

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

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

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

I am using the MASS library function fitdistr(x, dpois, list(lambda=2)) but I get Error in optim(start, mylogfn, x = x, hessian = TRUE, ...) : Function cannot be evaluated at initial parameters In addition: There were 50 or more warnings (use warnings() to see the first 50) and all the first 50 warnings say 1: non-integer x = 1.452222 etc Can anyone tell me what I am doing

On 17/11/2017 10:47 AM, Adelchi Azzalini wrote: > Since I have updated my Linux system, and consequently re-installed R, I > am unable to see HTML documentation which I used to access via a web > browser. To be more specific, my R installation is declared to be here: > > > R RHOME > /usr/local/lib/R > > However, if I look at a location such as >

Since I have updated my Linux system, and consequently re-installed R, I am unable to see HTML documentation which I used to access via a web browser. To be more specific, my R installation is declared to be here: > R RHOME /usr/local/lib/R However, if I look at a location such as /usr/local/lib/R/library/base/html this contains only two entries File: 00Index.html File: R.css