search for: adelchi

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

...not be directed 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/

...) and then 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/

...e of which 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 fun...

...se 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 Statistiche, Universit?? di Padova, Italia tel. +39 049 8274147, http://azzalini.stat.unipd.it/

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

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

...R version 3.4.2 (2017-09-28) -- "Short Summer" but the same problem existed with 3.4.1 after the system update. When I (re-)install R, everything seems fine; no error messages are generated. Can anyone tell me how to put the HTML documentation files in the right place? Best regards, Adelchi Azzalini

...ks like: Examples x <- toeplitz(rev(1:4)) x.inv <- pd.solve(x) print(x.inv logDet <- attr(x.inv, "log.det") print(abs(logDet - determinant(x, logarithm=TRUE)$modulus)) It seems that the examples are parsed as ordinary text, not R commands. I am really puzzled. Any explanation? Adelchi Azzalini ---- [mac:aa:~/SW/R/Pkg-mnormt] R CMD Check mnormt running .Rprofile (Sun 16-Jan-2011, 11:17:59) * using log directory ?/Users/aa/SW/R/Pkg-mnormt/mnormt.Rcheck? * using R version 2.12.1 (2010-12-16) * using platform: i386-apple-darwin9.8.0 (32-bit) * using session charset: UTF-8 * check...

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

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

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.

Hi, are there any function to plot the empirical probability distribution function? I just don't want to reinvent the wheel... Best wishes, Ott -- Ott Toomet PhD Student Dept. of Economics ?rhus University Building 322 Universitetsparken 8000 ?rhus C Denmark otoomet (a) econ au dk ph: (+45) 89 42 20 27 ------------------------------------------- (o_ (*_ (O_

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

I am keen to look at Rcgi as I want to put up some simple bits of R to do prescribed tasks on HTML form input. Rweb is overkill and worryingly flexible for what I want and it sounds as if Rcgi is more what I need. However, I can't get any of the URLs I've found for it to work over the last few days. Does anyone have a recent copy they could Email me or a working URL for it? TIA,