R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(288.125,331.5,493.025,378.75,458.625,431.25,399,473.5,442,177.425,346,450.25,509.625,509.625,301.375,374,374,311.875,300.375,352.75,454.25,208.7,470.85,490.125,551.125,522.8,577.875,384.25,384.25,552,661.5,457.1,647.5,174,781,277.1,653,435.75,613.775,509.75,509.75,314.5,486,212,503.825,435,563,457.05,451.25,500.75,437.75,470.5,0,313.25,314,454,570.5,485,243,310,421.752,494.5,253.5,417.5,182.826,339.25,199,412.25,438.25,356,266.25,235.25,323.775,305.25,383.527,515.25,496.15,115.25,170.5,154.25,170,534.05,193.75,564.5,213.63,308.25,437.05,410.275,149.75,154.75,240.1,127.525,222.25,85.525,427.75,63.5,118.3,99.5,182.25,401,119.5,450.25,147.5,237,80.025,10.5,176.75,234,282.5,320,167.5,163.25,238.15,325.125,126.3,154.875,327.25,336.25,188,277.25) > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(MASS) > (f<-fitdistr(x, 'exponential')) rate 0.0028924821 (0.0002640463) > xlab <- paste('Exponential(rate=',round(f$estimate[[1]],2)) > xlab <- paste(xlab,')') > postscript(file="/var/www/html/rcomp/tmp/181t61270549130.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqplot(qlnorm(ppoints(x), meanlog=f$estimate[[1]]), x, main='QQ plot (Exponential)', xlab=xlab ) > grid() > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Parameter',1,TRUE) > a<-table.element(a,'Estimated Value',1,TRUE) > a<-table.element(a,'Standard Deviation',1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'rate',header=TRUE) > a<-table.element(a,f$estimate[1]) > a<-table.element(a,f$sd[1]) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/2b1sc1270549130.tab") > > try(system("convert tmp/181t61270549130.ps tmp/181t61270549130.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.412 0.171 1.711