search for: azzalini

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

...volved in a minimization problem. For this reason, I could not make use of the CRAN package 'mvtnorm'. Exactly to avoid superposition with the CRAN package, 'mnormt' is made available somehere else, in case other people want to use it. The package 'mnormt' is at http://azzalini.stat.unipd.it/SN/Pkg-mnormt As explained, this is not uploaded to CRAN just to avoid clash with the existing package. However, if it is felt appropriate, I have no objection to upload it to CRAN. Best wishes, Adelchi Azzalini -- Adelchi Azzalini <azzalini at stat.unipd.it> Dipart.Scienze...

...en from outside R: sed s/\\./,/ outfile then I get it right (of course), something like -1,3612 -0,9772 0,1524 2,4046 0,4741 0,6659 -0,8277 0,5071 0,7190 0,4088 This is fine, but it would be good to have it working in first form. environment: R 1.7.0 on Debian linux. regards, Adelchi Azzalini -- Adelchi Azzalini <azzalini at stat.unipd.it> Dipart.Scienze Statistiche, Universit? di Padova, Italia http://azzalini.stat.unipd.it/

...ch the above is the core part, asking advice on how to handle this situation. First of all I would like to thank the list members that have expressed their views on this question. This has helped me to make up my mind: I reproduce below my 2008 message mentioned above. With best regards, Adelchi Azzalini Begin forwarded message: =========================================== Date: Tue, 26 Aug 2008 10:44:16 +0200 From: Adelchi Azzalini <azzalini at stat.unipd.it> To: Rmetrics-core at r-project.org Subject: dmvsnorm & mvst in fMultivar Dear Colleagues, I am writing about some functions in...

...irected to the R-help list, rather to the developer(s) of optim. Please instruct me on this point. Also, I appreciate that the above is not a reproducible example; this would be longish in text, and I though to ask first to whom it is appropriate that I direct my question. Best wishes, Adelchi Azzalini -- Adelchi Azzalini <azzalini at stat.unipd.it> Dipart.Scienze Statistiche, Universit? di Padova, Italia http://azzalini.stat.unipd.it/

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

...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 stat.unipd.it> Dipart.Scienze Statistiche, Universit? di Padova, Italia http://azzalini.stat.unipd.it/

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",

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.

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 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'?

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

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

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

...) Warning messages: 1: argument `lib` is missing: using C:/STAT/R/RW1021/library in: install.packages("sn") 2: No package "sn" on CRAN. in: download.packages(pkgs, destdir = tmpd, available = available, > Thanks for your attention and consideration. Best wishes Adelchi Azzalini --please do not edit the information below-- Version: platform = i386-pc-mingw32 arch = x86 os = Win32 system = x86, Win32 status = major = 1 minor = 2.1 year = 2001 month = 01 day = 15 language = R Windows 9x 4.0 (build 950) Search Path: .GlobalEnv, package:sm, Autoloads, package...

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

I'd welcome some comments or advice regarding the situation described below. The following illustrates what seems to me to be an inconsistency in the behaviour of matrix subsetting: > Z<-matrix(c(1.1,2.1,3.1,1.2,2.2,3.2,1.3,2.3,3.3),nrow=3) > Z [,1] [,2] [,3] [1,] 1.1 1.2 1.3 [2,] 2.1 2.2 2.3 [3,] 3.1 3.2 3.3 > dim(Z) [1] 3 3 >