Re: [R] why package ZIGP is not there anymore?
Funnily, the bivpois package is no longer available too... -- View this message in context: http://r.789695.n4.nabble.com/why-package-ZIGP-is-not-there-anymore-tp4660663p4660726.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] horizontal yaxis label
Yes, it works, but I wonder if I can make the y title horizontal just setting any "hidden" parameter... -- View this message in context: http://r.789695.n4.nabble.com/horizontal-yaxis-label-tp802761p4660725.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] wireframe and margins
Thanks, I'm setting now the mai parameter, but it is not so intuitive how it works with persp... -- View this message in context: http://r.789695.n4.nabble.com/wireframe-and-margins-tp4653838p4653865.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] wireframe and margins
Hi everybody, I'm working with wireframe (and persp) functions, and I'm not able to set the margins between the 3d plot and the R window... It should be easy, innit?, but I can't find out which parameters I have to change... Any help? Thanks, ktk -- View this message in context: http://r.789695.n4.nabble.com/wireframe-and-margins-tp4653838.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Compiling Rgraphiz on Windows 7 64bit with R-2.13.0
Thanks, I'm going to update it! -- View this message in context: http://r.789695.n4.nabble.com/Compiling-Rgraphiz-on-Windows-7-64bit-with-R-2-13-0-tp3493750p3557431.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Compiling Rgraphiz on Windows 7 64bit with R-2.13.0
lol, I'm having problems with texi2dvi as well... The Rdlatex file says: Hmm ... looks like a package Error in texi2dvi("Rd2.tex", pdf = (out_ext == "pdf"), quiet = FALSE, : pdflatex is not available Error in running tools::texi2dvi You may want to clean up by 'rm -rf C:/DOCUME~1/AB/IMPOST~1/Temp/RtmpOJ8zBP/Rd2pdf2ff2d7a' With the last version of Rtools, R 2.13 and Miktex 2.6 running on Windows XP. Thanks in advance for your help KtK -- View this message in context: http://r.789695.n4.nabble.com/Compiling-Rgraphiz-on-Windows-7-64bit-with-R-2-13-0-tp3493750p3556483.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] dbetagen function in mc2d package
bbolker wrote: > > krusty the klown libero.it> writes: > >> I found this useful package for generalized beta, yet the function that >> calculates its density leaves me puzzled, especially when I plotted it: >> >> plot(function(y) dbetagen(y,4,1.2,min=0,max=40),xlim=c(0,40)) >> >> The area between the density curve and the x axis does not seem to >> measure >> 1!!! > > It would help to tell us that this is in the "mc2d" package ... > > * I wrote it on the title! :-) > >> integrate(dbetagen,shape1=4,shape2=1.2,min=0,max=40,0,40) > 40 with absolute error < 2e-06 > > I would suggest that you report this as a bug to the > package maintainers [maintainer("mc2d")]. Seems fairly easy to > fix for yourself, though -- just scale by the size of the domain > (max-min). > > * That's ok, thanks! > > Ben Bolker > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- View this message in context: http://r.789695.n4.nabble.com/dbetagen-function-in-mc2d-package-tp3535411p3538128.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] dbetagen function in mc2d package
Hi all, I found this useful package for generalized beta, yet the function that calculates its density leaves me puzzled, especially when I plotted it: plot(function(y) dbetagen(y,4,1.2,min=0,max=40),xlim=c(0,40)) The area between the density curve and the x axis does not seem to measure 1!!! -- View this message in context: http://r.789695.n4.nabble.com/dbetagen-function-in-mc2d-package-tp3535411p3535411.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Status of the "bs" Package
Perhaps there's a problem with the function rbs, that generates a random sample from the Birnbaum-Saunders distribution. The rgbs function (package gbs) does it better. Look at these code and results: > library(bs) > set.seed(1) > alpha<-0.17 > beta<-130 > x<-rbs(n=1000,alpha,beta) > # sample mean > mean(x) [1] 150.8288 > # expected value > beta*(1+alpha^2/2) [1] 131.8785 > est1bs(x) $beta.start [1] 149.9556 $alpha [1] 0.1079173 $beta [1] 149.9558 $converge [1] "TRUE" $iteration [1] 2 > #not fine > library(gbs) > # > set.seed(1) > x<-rgbs(n=1000,alpha,beta) > # sample mean > mean(x) [1] 131.7481 > # expected value > beta*(1+alpha^2/2) [1] 131.8785 > est1bs(x) $beta.start [1] 129.7422 $alpha [1] 0.1758475 $beta [1] 129.7421 $converge [1] "TRUE" $iteration [1] 2 > #fine -- View this message in context: http://r.789695.n4.nabble.com/Status-of-the-bs-Package-tp827806p3033734.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Status of the "bs" Package
Perhaps the problem stands on the rbs function which generates random samples from the Birnbaum-Saunders distribution: > library(bs) > set.seed(1) > x<-rbs(n=1000,alpha=0.5,beta=1.0) > # sample mean > mean(x) [1] 1.117749 > # expected value > beta*(1+alpha^2/2) [1] 2.125 > # so different! -- View this message in context: http://r.789695.n4.nabble.com/Status-of-the-bs-Package-tp827806p3033546.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Status of the "bs" Package
Hi, the package bs, version 1.0, is available. Yet, I've gote some problems with it, specifically with the maximum likelihood estimation of the parameters \alpha and \beta of Birnbaum-Saunders distribution. Let the following code run: library(bs) alpha<-1.5 beta<-0.5 n<-100 set.seed(1) x<-rbs(n,alpha,beta) est1bs(x) The MLE algorithm does not converge. Not only for this pair of values and for this sample... Please have a look... -- View this message in context: http://r.789695.n4.nabble.com/Status-of-the-bs-Package-tp827806p3033442.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] homals (level numerical, scoremat and catscores...)
Just run this code: library(homals) data(sleeping) sleeping res <- homals(sleeping[,2:4], level = c(rep("numerical",2),"nominal"), ndim=1) plot(sleeping[,2],res$scoremat[,,1][,1],axes=F) axis(1,sleeping[,2]) abline(h=res$scoremat[,,1][,1]) box() It looks strange, innit? There is something which does not come out right among catscores and scoremat... which are the right transformed values? :working: -- View this message in context: http://old.nabble.com/homals-%28level-numerical%2C-scoremat-and-catscores...%29-tp26526489p26526489.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] (exact) confidence bounds for lognormal parameters \mu and \sigma
Sorry, probbaly it is...:blush: Krusty the Klown wrote: > > Anyway, the results do not provide confidence intervals, but mle estimates > and the hessian matrix, whose extradiagonal elements are null... is it > correct? > http://www.weibull.com/LifeDataWeb/confidence_bounds_log.htm > http://www.weibull.com/LifeDataWeb/confidence_bounds_log.htm > > -- View this message in context: http://old.nabble.com/%28exact%29-confidence-bounds-for-lognormal-parameters-%5Cmu-and-%5Csigma-tp26394177p26404958.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] (exact) confidence bounds for lognormal parameters \mu and \sigma
Thank you, it is what I was looking for: library(fitdistrplus) x<-rlnorm(100) res<-mledist(x, distr="lnorm") res Anyway, the results do not provide confidence intervals, but mle estimates and the hessian matrix, whose extradiagonal elements are null... is it correct? http://www.weibull.com/LifeDataWeb/confidence_bounds_log.htm http://www.weibull.com/LifeDataWeb/confidence_bounds_log.htm I will investigate... David Scott-6 wrote: > > Krusty the Klown wrote: >> Dear all, >> a statistical question: how can I compute exact confidence intervals for >> the >> lognormal distribution parameters? I found something only on >> www.weibull.com www.weibull.com . Does exist a package in R which can >> compute them? >> Thanks in advance, >> KTK > > I think fitdistrplus can do this > > David Scott > -- > _ > David Scott Department of Statistics > The University of Auckland, PB 92019 > Auckland 1142,NEW ZEALAND > Phone: +64 9 923 5055, or +64 9 373 7599 ext 85055 > Email:d.sc...@auckland.ac.nz, Fax: +64 9 373 7018 > > Director of Consulting, Department of Statistics > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > -- View this message in context: http://old.nabble.com/%28exact%29-confidence-bounds-for-lognormal-parameters-%5Cmu-and-%5Csigma-tp26394177p26404934.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] (exact) confidence bounds for lognormal parameters \mu and \sigma
Dear all, a statistical question: how can I compute exact confidence intervals for the lognormal distribution parameters? I found something only on www.weibull.com www.weibull.com . Does exist a package in R which can compute them? Thanks in advance, KTK -- View this message in context: http://old.nabble.com/%28exact%29-confidence-bounds-for-lognormal-parameters-%5Cmu-and-%5Csigma-tp26394177p26394177.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] how to modify a function in a R package calling other (invisible?) functions
Hey hey kids! I'm facing this rough problem: I have to modify a function in a R package (namely, homals) but I'm not able to do it from the R interface, since this function recalls other functions which looks like invisible... I downloaded the source package from CRAN, http://cran.r-project.org/src/contrib/homals_0.9-10.tar.gz http://cran.r-project.org/src/contrib/homals_0.9-10.tar.gz how can I work it out? Yours, KtK:-D -- View this message in context: http://www.nabble.com/how-to-modify-a-function-in-a-R-package-calling-other-%28invisible-%29-functions-tp23390639p23390639.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.