I don't read R-help these days so have just seen this.
Both generalized inverse Gaussian and normal inverse Gaussian are in
GeneralizedHyperbolic.
HyperbolicDist is no longer being maintained.
David Scott
On 12/06/2012 5:41 a.m., David L Carlson wrote:
Should have been
For the normal
/htmlcodes.htm
David Scott
On 31/07/2012 9:21 p.m., ramonovelar wrote:
Thanks Arun,
Yes, I have Windows 7. I have tried 2 versions of R, 2.14.1 and 2.15.x, but
it did not change anything.
Right now I can't try a different version of win.
Ramón
On Tuesday, July 31, 2012, arun kirshna [via
to the package.
David Scott
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on behalf of
Joshua Wiley [jwiley.ps...@gmail.com]
Sent: Monday, January 02, 2012 9:31 AM
To: Michael
Cc: r-help
Subject: Re: [R] R report generator (for Word
there is this functionality in Rob Hyndman's forecast package.
See this post by him:
http://robjhyndman.com/researchtips/longseasonality/
which has some code which should also be of use.
David Scott
_
David Scott Department
There is a free eBook on tinn-R available from Rmetrics:
https://www.rmetrics.org/ebooks-tinnr
Written by the authors of tinn-R.
Please consider a donation to the Rmetrics Association.
--
_
David Scott Department
be able to access a css file but I don't see any
mechanism for doing that. Perhaps someone can enlighten me.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
On 20/07/11 18:56, Dieter Menne wrote:
David Scott-6 wrote:
I have been trying to read some data from an Excel workbook without
success.
...
faults- sqlFetch(channel, sqtable = 'Data',
+colnames = FALSE, as.is = TRUE)
faults
[1] HY001 -1040 [Microsoft][ODBC Excel
graphics grDevices utils datasets methods base
other attached packages:
[1] RODBC_1.3-2 djsmisc_1.0-1
loaded via a namespace (and not attached):
[1] tools_2.13.0
David Scott
--
_
David Scott Department of Statistics
of any named data
ranges.
Sample code from a student exercise:
require(RODBC)
channel - odbcConnectExcel(bikesWithDate.xls)
sqlTables(channel)
David
--
_
David Scott Department of Statistics
The University
pointers to
useful code, relevant ideas, or similar analyses.
David Scott
__
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
.
__
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.
Does
Sys.info()[nodename]
give you what you want?
David Scott
graphics you do. I don't know
where the qplot() function is from, but it looks as though it must be
from a package that uses grid.
Duncan Murdoch
_
qplot is from ggplot2, and is a simpler form of ggplot with limited
features.
David Scott
will go as
you expect. (But I haven't tested this.)
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 923 5055
to programming style.
David Scott
On 19/05/11 10:26, bill.venab...@csiro.au wrote:
Hi Bert,
I think people should know about the Google Sytle Guide for R because, as I said,
it represents a thoughtful contribution to the debate. Most of its advice is very
good (meaning I agree
to produce html
as well.
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
in
contributing to this thread was to have something explicit in the
archives for future reference, and your reply is excellent in that regard.
And many thanks for your work on lattice.
David Scott
--
_
David Scott Department
.
As to why this is the case .
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
On 05/04/11 13:14, David Winsemius wrote:
On Apr 4, 2011, at 9:03 PM, David Winsemius wrote:
On Apr 4, 2011, at 8:42 PM, David Scott wrote:
On 05/04/11 05:58, David Winsemius wrote:
On Apr 4, 2011, at 1:27 PM, Marius Hofert wrote:
Dear David,
do you know how to get plotmath-like symbols
.
Provisos: there may be problems with distributions with very heavy
tails, and generally the routines could be slow.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
take the name of the
file, prepend the path to the directory where your code is located and
then call source. Two or three lines would do it.
Also, who knows, maybe Gabor has something clever in his batch files, he
usually seems to.
David Scott
, 1900). In addition, Excel considers 1900-01-01 as day 1,
not day 0.
Hope this is helpful,
Dan
An explanation which seems reasonably authoritative is given here:
http://www.cpearson.com/excel/datetime.htm
David Scott
Daniel J. Nordlund Washington State Department of Social and Health
.
--
_
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
Thanks Brian, I stand corrected.
David Scott
On 27/02/2011 12:32 a.m., Prof Brian Ripley wrote:
It is less clear what you are after, but the canonical way to decide
if your R session is on Windows is
.Platform$OS.type == windows
Unlike {R.}version$os and Sys.info()[sysname], the set
{These are not to be called by the user.}
\keyword{internal}
I only mention this possibility for the sake of intellectual
completeness of the discussion of this topic.
It of course totally subverts the whole philosophy of package
construction in R, so I wouldn't dream of using it myself.
David Scott
can take time to get right. It does offer a bit more
control over line spacing which could be a concern with atop.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
]] - testList
}
head(megaList)
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
, please let me know thanks.
typhoong
Your test is wrong for NA. You must use is.na:
df
x y
1 1 345
2 6 NA
3 8 123
4 32 123
5 12 NA
6 6 124
7 7 NA
subset(df, is.na(y))
x y
2 6 NA
5 12 NA
7 7 NA
--
_
David
==.(param)))
then in the title I get Parameter = alpha,
whereas I want the Greek character alpha.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland
Thanks. Exactly what I wanted.
As usual, I played around with all sorts of things to try and get the
expression right, but never thought of as.name.
David Scott
On 25/01/2011 4:32 a.m., Charles C. Berry wrote:
On Tue, 25 Jan 2011, David Scott wrote:
I have a problem with expressions. I am
orthogonal polynomials of degree 1 to
'degree' over the specified set of points 'x'. These are all
orthogonal to the constant polynomial of degree 0. Alternatively,
evaluate raw polynomials.
David Scott
On 21/01/2011 3:50 a.m., Joshua Wiley wrote:
Hi Barth,
Here is an option
is that those who enter data into spreadsheets enter all
sorts of variations of what a human would wish to treat as the same
(Open, Open , open, etc.), even when told not to.
David Scott
On 14/01/2011 4:03 p.m., Jim Holtman wrote:
try strip.white=TRUE to strip out white space
Sent from my iPad
was the conversion of New
Zealand health ID numbers to dates. They are three letters then 4
digits, so AUG2699 became a date.
David Scott
On 14/01/2011 10:58 p.m., bgr...@dyson.brisnet.org.au wrote:
Hello David,
Thanks for your e-mail. The data was a report derived from a statewide
database, saved
has the inverse Gaussian as a special case.
David Scott
Secondly, is there a way to shift the pdf (code below) to the right (rather
than the data to the left)? I tried:
par(new=T)
shift-1
hist(iniSal_US_forHist,breaks=seq(1.1,21,by=0.625),col=grey,freq=F,xlim=c(0-shift,21-shift),plot=F)
curve
:
http://rwiki.sciviews.org/doku.php?id=tips:data-io:ms_windows
David Scott
There are several parts to my problem.
1) When I try the read.table command, even if I include the directory name
in the filename, R still cannot read the file, even if it is in .csv format
2) I was able to copy my
files, that I
just avoid them whenever I can.
Absolutely agree with this Duncan. Playing around with .csv files is
like playing with some sort of unstable explosive. I also avoid them as
much as possible.
David Scott
Duncan Murdoch
__
R-help@r
guy to add the MySQL/ODBC connector.
The instructions may not be quite up to date but hopefully you can find
your way through. Our students generally manage to set it all up on
their own computers.
David Scott
On Wed, Nov 24, 2010 at 1:48 PM, PtitBleuptit_b...@yahoo.fr wrote:
Hello,
I
periodogram (as the name suggests). What did you
test? Only a white noise process should stay inside the confidence
bands. There is some information about the use of the cumulative
periodogram in Venables and Ripley's book for which cpgram was written
(but admittedly not a lot).
David Scott
On 23/11/10 14:20, Stephen Liu wrote:
Hi folks,
Win7 64 bit
IE 64 bit
How to start IE on R? TIA
B.R.
Stephen L
?browseURL
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
can use it to open a local file on your machine as well. I use this
all the time with hwriter which writes html reports.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
*(0:(n-1)))
specx - spec.pgram(x, plot = FALSE)
spec.pgram(x)
spec.pgram(x, xaxt = n, xlab = frequency (Hz),
sub = paste(bandwidth = , round(1000*specx$bandwidth,2)))
axis(side = 1, at = (0:5)/10, labels = 1000*(0:5)/10)
David Scott
for free. timeDate is an Rmetrics package.
Go to https://www.rmetrics.org/ebooks to see what is available from
Rmetrics.
David Scott
-Original Message-
From: Li, Jing Yi [mailto:jingyi...@credit-suisse.com]
Sent: Friday, October 22, 2010 9:38 AM
To: David Reiner; David Winsemius
Cc: r
in .csv format. If the .csv files are similar there
shouldn't be too much overhead in defining table formats for the database.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland
like this in to me to grade.
David Scott
On 21/09/2010 5:32 p.m., snes1...@hotmail.com wrote:
I created a EM algorithm for Generalized hyperbolic distribution.
I want to estimate mutheldaplus, sigmatheldaplus, betasigmaplus in my code.
After getting use these value , then my iteration have
in particular, and possibly odfWeave.
Does anyone have sample documents using any of these approaches which
they could let me have?
David Scott
_
David Scott Department of Statistics
The University of Auckland, PB 92019
2008 2008 2008 2008 2008 2008 2008
2009 2009 2009 2009 2009 2009 2009
David Scott
On 1/09/2010 11:04 a.m., Alistair Gray wrote:
Dear List,
I'm getting weird and unexpected behaviour using time and %%, or %/%. It's
likely I'm not appreciating the nuances of floating point arithmetic. Or it
could
function.
There is one (logHist) in my package DistributionUtils on CRAN. You
don't need the rest of the package to use it. You could just extract
that particular function.
David Scott
--
_
David Scott Department of Statistics
in
the density correspond to straight lines in the logged y-axis plot.
Hope this helps.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142
.
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
code.
--
_
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
to the subroutine) to make this work?
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
Samuel Okoye wrote:
I can't find his email and I have asked the same question to
bioconduc...@stat.math.ethz.ch
Regards,
Samuel
So you need the handy function maintainer() added to 2.11.0:
require(GeneralizedHyperbolic)
maintainer(GeneralizedHyperbolic)
[1] David Scott d.sc
nuncio m wrote:
Hi all,
Is there any way in R to select the order of an ARIMA model
automatically
nuncio
Rob Hyndman's package forecast has a function auto.arima which produces
an automatic arima fit, including for seasonal models.
David Scott
Ron Burns wrote:
Hi all-
Is there an R function that returns the directory in which a package has
been installed?
Thanks
Ron
?system.file
For example:
system.file(package=DistributionUtils)
[1] C:/Users/dsco036/R/win-library/2.10/DistributionUtils
David Scott
In-memory Compression and Decompression
utils::unzipExtract or List Zip Archives
utils::zip.file.extract
Extract File from a Zip Archive
utils::win.version Auxiliary Functions for the Windows Port
which seems to provide some leads
David Scott
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
system.file available, which I presume works across platforms:
system.file(RunTests.R,package=DistributionUtils)
[1] /usr/local/david/R/DistributionUtils
system.file(unitTests/runTests.R,package=DistributionUtils)
[1] /usr/local/david/R/DistributionUtils/unitTests/runTests.R
David Scott
David Scott wrote:
Sharpie wrote:
Johannes W. Dietrich wrote:
There is certainly a trivial solution for my question, but I can't
find the answer in the documentation.
I need a platform independent method to obtain the file path of the
current R script.
My working group uses R on several
rule these days, when using Windows, I prefer to avoid the
.xls or .xlsx to .csv route to getting data into R from Excel and use
either RODBC or xlsReadWrite. Getting Excel to write a .csv file is just
adding some random noise to the process.
David Scott
.
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
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
guessing my MySQL is getting a bit decrepit and needs updating, but
welcome any suggestions about sorting this out.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Marc Schwartz wrote:
On Apr 19, 2010, at 7:02 PM, David Scott wrote:
I have a problem with some missing headers when trying to install RODBC on my
linux box.
install.packages(RODBC,
+ lib=/usr/local/david/R)
trying URL 'http://cran.stat.auckland.ac.nz/src/contrib
specialized and advanced
texts also. Off the top of my head, Pfaff, Hyndman, ...
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142,NEW
within Excel. Use
.csv files created by other software, yes, but not one saved from within
Excel.
David Scott
Nevertheless I got some valuable insight in the topic. What I'm trying to do
here is make it easier for people to upload files that might contain dates
to my webapp at yeroon.net/ggplot2
_
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
_
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
is taken from the keyboard
(or whatever ‘stdin()’ reads if input is redirected or R is
embedded).
I think the required function is readline which prompts for user input.
?readline
David Scott
_
David Scott
more work on your behalf.
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
to generalize to
this case.
Any pointers?
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 923 5055, or +64
before testing:
vtest - function(x, lookfor){
any(apply(x, 1, function(v)
{identical(sort(v), sort(lookfor))}))}
-Peter Ehlers
On 2010-03-27 2:46, Berend Hasselman wrote:
David Scott-6 wrote:
I am sure someone can come up with a clever way of doing what I
to make sure the paragraph you want to wrap is separated from other
paragraphs by blank lines, i.e. it must be clear it is a paragraph.
David Scott
--
_
David Scott Department of Statistics
The University
as an
appropriate method.
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
:/Salmon_data.accdb)
channel - odbcConnectAccess(testdb)
Data - sqlFetch(channel,Table1)
Try
odbcConnectAccess2007
instead.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
in the
department who can assist. Someone set up the linux machine where R is
running, they will know this stuff.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
questions he should be able to find his way
around the help system and other documentation by now.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland
is an environment which rewards a bit
of exploration and experimentation.
Every package has a maintainer whose email is given on CRAN. In the case
of the contrast package that is at
http://cran.r-project.org/web/packages/contrast/index.html
It ain't that hard.
David Scott
Microsoft SQL Server.
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
, scientific=F), sep=)
title(the_expression)
grid()
When confronted with this problem, after many failed attempts, I have
had to resort to mtext.
David Scott
--
_
David Scott Department of Statistics
xlsReadWritePro at
http://www.swissr.org/ which can write native xls files with multiple
sheets.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland
=2.25,cex=1.15)
Best wishes
David Scott
- Original Message
From: David Scott d.sc...@auckland.ac.nz
To: Jason Rupert jasonkrup...@yahoo.com
Cc: R-help@r-project.org R-help@r-project.org
Sent: Wed, January 20, 2010 3:49:40 PM
Subject: Re: [R] Greek letters on a multi-line plot title
guide to the build process, have a look at Rob
Hyndman's instructions:
http://robjhyndman.com/researchtips/building-r-packages-for-windows/
David Scott
_
David Scott Department of Statistics
The University
0.994 ...
..@ title : chr Gumbel Parameter Estimation
..@ description: chr Fri Jan 08 11:12:32 2010
p...@fit$par.ests
mu beta
0.03515609 0.99286204
David Scott
_
David Scott Department of Statistics
estimator...:-(
Any help would be great:-)
You can see an implementation of mps for the normal inverse Gaussian in
fBasics on R-Forge.
Have a look at the function .nigFit.mps in the file dist-nigFit.R
David Scott
--
_
David Scott
on their website led me to believe Stata 11
couldn't write in Stata 10 format. However Stata 11 can write Excel
files in XML, so presumably .xlsx files.
Any advice, correction of my misconceptions, etc. gladly received.
David Scott
in Version 11 of Stata will save a Version 10/11 file in the format that was
used by Versions 8 and 9 of Stata.
---JRG
On Wed, 16 Dec 2009, David Scott wrote:
I have a client who uses Stata 11.
Can anyone advise me on ways of transferring data from this version of Stata
to R?
Reading
Gaussian distribution. Its support is the real line and it has 4
parameters. The inverse Gaussian is a special case of the generalized
inverse Gaussian. It has support the positive half line and 2 parameters.
David Scott
On Sun, Dec 6, 2009 at 12:48 PM, noclue_ wrote:
I have a one-variable data
David Scott wrote:
Jorge Ivan Velez wrote:
Hi Tim,
Take a look at
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html
HTH,
Jorge'
i don't think this is what the OP wanted. The nig is the normal inverse
Gaussian distribution, which is a special case
the cospectrum
looks. That is really the only infallible way of determining what the
software is doing.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland
with a Gamma function
of high order in the denominator.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 923 5055
This is a reply to my own question. I thought I had found an answer but
it seems not so (some analysis follows below). Maybe Martin Maechler or
Robin Hankin or Duncan Murdoch may have some ideas---I know the question
is a bit specialized.
David Scott wrote:
I am looking for a method
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
to fit a (perhaps just a simple linear) model with t-distributed
residuals (or residuals of a different distribution)?
Package sn has this facility I believe.
David Scott
--
_
David Scott Department of Statistics
Kjetil Halvorsen wrote:
On Wed, Oct 28, 2009 at 7:25 AM, David Scott d.sc...@auckland.ac.nz wrote:
Karl Ove Hufthammer wrote:
On Tue, 27 Oct 2009 18:06:02 -0400 Ben Bolker bol...@ufl.edu wrote:
If transforming your data brings you closer to satisfying
the assumptions of your analytic methods
Tabulation with Tests for Factor Independence
Description
An implementation of a cross-tabulation function with output similar to
S-Plus crosstabs() and SAS Proc Freq (or SPSS format) with Chi-square,
Fisher and McNemar tests of the independence of all table factors.
David Scott
of the violation of the independence assumption.
If Antonio does not understand these comments he needs to consult a
local statistician.
David Scott
--
_
David Scott Department of Statistics
The University of Auckland
that it will give you the counts you want is just wrong. As
others have pointed out if you make things more explicit and don't take
the defaults it will do so.
Using hist to obtain counts like this is like using a hammer to drive in
a screw.
David Scott
are available for download.
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
nonparametric stats and this one has me scratching my head.
Many thanks!
Tom
The continuity correction? It is true by default for wilcox.test and is
not apparent in the help for kruskal.test.
David Scott
--
_
David Scott
does suffer from very
inconsistent naming and as I think Duncan said it makes it very
difficult sometimes to remember names when you have variations in case
and separators as with things related to system.
David
_
David Scott
to you.
cheers,
Rolf Turner
I have used the dip test for testing if a distribution is *unimodal*.
Possibly that is what John needs. See package diptest.
David Scott
--
_
David Scott Department of Statistics
1 - 100 of 165 matches
Mail list logo