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
*(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
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
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
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
.
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
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
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
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
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 .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
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
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
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
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
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
_
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
_
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
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
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
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
.
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
, 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
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
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
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
, but there is an ODBC
driver for linux and I have used it to connect it quite successfully
(thank you Brian) to access MySQL on a linux box from both Windows and
Linux boxes.
The installation of the ODBC driver on linux is described in the README
file which is part of the RODBC package.
David Scott
Not sure about answering your original question, but why not write the
graphics output to a file, then fire up a viewer?
David Scott
On Wed, 24 Oct 2007, Luke Spadavecchia wrote:
Hi there,
I am trying to plot some output from a FORTRAN (ifort) program using
R (2.5.1) under batch mode
In regard to validation of R, there are test directories with test
code. See for example:
https://svn.r-project.org/R/trunk/tests/
and
https://svn.r-project.org/R/trunk/src/library/stats/tests/
David Scott
_
David Scott
There is a whole package to do this: nortest
David Scott
On Thu, 8 Nov 2007, [EMAIL PROTECTED] wrote:
Hi,
I would like to know if there is an algorithm in R for testing if a data set
as a normal destribution.
Thank you for your time,
Pedro Marques
If you are after the lognormal distribution, I would recommend you use the
lognormal functions from R itself. Check them out with
?dlnorm
David Scott
On Wed, 14 Nov 2007, Scionforbai wrote:
If what you want is a lognormal distribution of n values you can use
the following transformations
yvals - histRes$counts
length(xvals)
length(yvals)
xvals - c(xvals,xvals[length(xvals)])
yvals - c(0,yvals,0)
plot(xvals,yvals,type=S)
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University
use the uniroot function which is far more
efficient than bisection. See
?uniroot
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland
I am aware of one (unofficial) guide to style for R programming:
http://www1.maths.lth.se/help/R/RCC/
from Henrik Bengtsson.
Can anyone provide further pointers to good style?
Views on Bengtsson's ideas would interest me as well.
David Scott
the ' characters.
Overall, I was reasonably happy with the approach I took. If I have a
future need I will have some ready-made functions to work from which will
make life easier. I will be interested to try some of the other
suggestions in this thread.
David Scott
would welcome any advice. I had a look at the rtools section in
the Administration Manual and the rtools site without gaining sufficient
enlightenment to proceed.
I am running Windows Vista and R 2.6.2
David Scott
_
David Scott
. Pick up the latest version
from Hans-Peter Suter's website
http://treetron.googlepages.com/
David Scott
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142
including underlying C and Fortran are
explained in the article by Uwe Ligges:
Uwe Ligges. R Help Desk: Accessing the sources. R News, 6(4):43-45,
October 2006.
This article is essential reading for anyone contemplating much use of R.
David Scott
.
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 373 7599 ext 86830 Fax: +64 9 373 7000
Email
My guess is that you maybe have some leading or trailing spaces in the
direction_ variable.
Try a tabulation of the direction_ variable to see exactly what values are
in your data set.
David Scott
On Thu, 14 Aug 2008, Farley, Robert wrote:
I can't figure out the syntax I need to get
)
SurveyData$direction_ - temp
should do it. You could just replace the direction_ variable without using
temp but it is safer this way and you should check the values of temp
before replacing your original data.
David Scott
_
David
.
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 373 7599 ext 86830 Fax: +64 9 373 7000
probably be improved
by a regular expression whiz, but works:
strngs - c(06/05/08,23/11/2008)
sub(([0-9][0-9]/[0-9][0-9]/)([0-9][0-9]$),\\120\\2,strngs)
[1] 06/05/2008 23/11/2008
David Scott
_
David Scott Department of Statistics
really simple.
Thanks,
Edna Bell
getAnywhere(t.test.default)
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW
On Thu, 18 Sep 2008, Deepayan Sarkar wrote:
On 9/18/08, David Scott [EMAIL PROTECTED] wrote:
I have a data set concerning ferritin levels in blood. There are three
relevant columns for this question, ferritin (continuous), score (ordered,
from 0 to 8) and gender. There is a good linear
of a lot else being available.
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 373 7599 ext 86830
[,sort(names(df))]
sortdf
a y z
1 3 5 1
2 4 6 2
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9
k1 k2 z a y
1 1 3 5 7 9
2 2 4 6 8 10
revdf - df[,rev(names(df))]
revdf
y a z k2 k1
1 9 7 5 3 1
2 10 8 6 4 2
David Scott
milicic.marko wrote:
Hi,
I have the data.frame with 4 columns. I simply want to invert dataset
so that last row becomes first...
I tried with rev(my_data
On Fri, 26 Sep 2008, David Scott wrote:
On Thu, 25 Sep 2008, Julian Burgos wrote:
How about something like
my.data=my.data[,4:1]
Julian
This can be tackled in a similar way to the question by Mark Na which I just
answered. Using the same data frame construction we have
df
Googling MDI gives a wikipedia entry in the first place which should
explain things.
David Scott
On Thu, 24 Jul 2008, Duncan Murdoch wrote:
On 7/24/2008 8:24 AM, Megh Dal wrote:
Thanks for this. Can you please explain me what are those MDI and SDI and
the diff in laymans language
the single quote which starts at 'select. (That is you need a ' just
before the closing bracket.
David scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
think what you are missing is as.name.
I think you need (untested) something like:
names(as.name(paste(data,i,sep=.))) - y1
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University
://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
QRMlib has routines for fitting t distributions. Have a look at that
package. Also sn has routines for skew-t distributions
David Scott
On Thu, 8 May 2008, kate wrote:
I have a data with 236 observations. After plotting the histogram, I found that
it looks like non-central t
the argument plot=FALSE and hist will produce the
counts etc for you
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW
who uses TeX and R it is an excellent
option.
I haven't checked John Fox's advice concerning XEmacs lately but in the
past it has been the best source of information about XEmacs and R if you
do still want to stick with XEmacs.
David Scott
On Mon, 12 May 2008, Vincent Goulet wrote:
Le lun
gratitude for a nice piece of work.
David Scott
On Mon, 2 Jun 2008, Dan Montgomery wrote:
I would like to set up a Windows XP environment with R and MySQL.
This set-up was recommended by a friend who uses this set up, but I
have been unable to get RJDBC to work. I have installed all of the
software
:
http://www.stat.auckland.ac.nz/~dscott/779/UsingMySQLAndODBC.pdf
David Scott
Many thanks for your help!
Me.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
and in some cases getting an
empty vector.
x - numeric(0)
max(x)
[1] -Inf
Warning message:
In max(x) : no non-missing arguments to max; returning -Inf
David Scott
Best as I can tell the code is working properly and the output is as
expected. I would like some help in understanding why he
.
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
.
--
_
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
Graduate Officer, Department
The preface to Heyde and Seneta can be seen on Google books and possibly
on Amazon and gives a number of recent references.
David Scott
_
David Scott Department of Statistics
The University of Auckland, PB 92019
.
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
Auckland 1142,NEW ZEALAND
Phone: +64 9 373 7599 ext 86830 Fax: +64 9 373 7000
Email: [EMAIL PROTECTED
truly works in
misterious ways.
read.csv has an argument fill which should fix this problem. In my case I
was actually reading the .csv file into mySQL and the solution was to
select the whole of the .xls file and format it as text before writing the
.csv file.
David SCott
with the character encoding. On the
linux box I have
LANG=en_US.UTF-8
Any advice would be welcome
David Scott
_
David Scott Department of Statistics, Tamaki Campus
The University of Auckland, PB 92019
On Tue, 22 Jan 2008, Prof Brian Ripley wrote:
On Wed, 23 Jan 2008, David Scott wrote:
I have encountered a problem with reading a .csv file on a linux box. I
can read the file on my windows machine (under XP) but on the linux box it
gives :
patients - read.csv(../Patients.csv, header
am on. Painful but at
least it is system independent.
Thanks again
David
On Tue, 22 Jan 2008, Prof Brian Ripley wrote:
On Wed, 23 Jan 2008, David Scott wrote:
On Tue, 22 Jan 2008, Prof Brian Ripley wrote:
On Wed, 23 Jan 2008, David Scott wrote:
I have encountered a problem
I am looking for a recommendation for an open source competitor to Winzip.
I seem to recall Brian Ripley mentioning one in the last year or so, but
couldn't find it in the mail archives. (Searching on Ripley there is
somehow not terribly useful.)
Suggestions?
David Scott
and Visualizing Data.
David Scott
On Tue, 29 Jan 2008, Greg Snow wrote:
I would question if a stacked area chart is really the best way to display
information on 12 groups. You can put a lot of information into the plot,
but the viewer will probably experience information overload and not be able
on CRAN xlsReadWrite is not currently available. There is an
archived version available however.
David Scott
_
David Scott Department of Statistics
The University of Auckland, PB 92019
Auckland 1142
1 - 100 of 165 matches
Mail list logo