Displaying 20 results from an estimated 1100 matches similar to: "bootstrap: definition of original statistic"
2014 Feb 03
1
seq range argument
Hello dear developers,
I find myself often having the result of "range" oder "extendrange", which I want to create a sequence with.
But "seq" needs two seperate arguments "from" and "two".
Could an argument "range" be added?
Otherwise I will have to create an object with the range (may come from a longer calculation), index twice from
2011 May 26
4
predictive accuracy
I am trying to develop a prognostic model using logistic regression. I
built a full , approximate models with the use of penalization - design
package. Also, I tried Chi-square criteria, step-down techniques. Used
BS for model validation.
The main purpose is to develop a predictive model for future patient
population. One of the strong predictor pertains to the study design
and would not
2012 Mar 05
5
cron job not running
hello list,
I am attempting to backup a centos 5.4 (x86_64) server running mysql
with a cron job. Here's how the cron job looks:
[root at cloud:/home/bluethundr/backupdb] #crontab -l
* 3 * * * /usr/bin/mysqldump jfwiki >
/home/bluethundr/backupdb/wiki-$(date +%Y%m%d).sql
However if I run the command from the command line it seems to work fine.
If I grep syslog for cron this is what
2006 Nov 25
3
Impossible to change Primary group
I have Samba PDC and use tdbsam backend. I use the version 3.0.11, but
recently I have exchanged it for the version 3.0.23c. Now I have a
number of problem. Here is one of then.
I wish to move one user "gad" from Domain Users (RID 513) to Domain
Guests (RID 514).
Now Primary group of user "gad" is "Domain Users":
pdbedit -Lv gad -d0
Unix username: gad
NT
2011 Mar 01
1
Pairwise T-Tests and Dunnett's Test (possibly using multcomp)
Hello Everyone,
I've been learning to use R in my spare time over the past several months. I've read about 7-8 books on the subject. Lately I've been testing what I've learned by trying to replicate the analyses from some of my SAS books. This helps me make sure I know how to use R properly and also helps me to understand how the two programs are similar and different.
Below is
2007 Mar 16
3
ARIMA standard error
Hi,
Can anyone explain how the standard error in arima() is calculated?
Also, how can I extract it from the Arima object? I don't see it in there.
> x <- rnorm(1000)
> a <- arima(x, order = c(4, 0, 0))
> a
Call:
arima(x = x, order = c(4, 0, 0))
Coefficients:
ar1 ar2 ar3 ar4 intercept
-0.0451 0.0448 0.0139 -0.0688 0.0010
s.e.
2007 Mar 20
1
truehist bug?
Hi,
Is this a bug in truehist()?
> library(MASS)
> x <- rep(1, 10)
> truehist(x)
Error in pretty(data, nbins) : invalid 'n' value
Thanks,
Gad
> R.version
platform i486-pc-linux-gnu
arch i486
os linux-gnu
system i486, linux-gnu
status
major 2
minor 4.1
year 2006
month 12
day 18
svn
2007 Mar 20
1
truehist bug?
Hi,
Is this a bug in truehist()?
> library(MASS)
> x <- rep(1, 10)
> truehist(x)
Error in pretty(data, nbins) : invalid 'n' value
Thanks,
Gad
> R.version
platform i486-pc-linux-gnu
arch i486
os linux-gnu
system i486, linux-gnu
status
major 2
minor 4.1
year 2006
month 12
day 18
svn
2007 Oct 16
2
survreg's algorithm
Hi,
I'm using survreg() from the survival package for parametric survival
regression (modelling inter-arrival times of patients to a waiting list
as exponentially distributed, with various regressors such as queue size
and season).
Does anyone know which algorithm survreg() uses for this?
Thanks,
Gad
--
Gad Abraham
Department of Mathematics and Statistics
The University of Melbourne
2006 Jun 13
3
Multiple lag.plots per page
Hi,
I'm trying to plot several lag.plots on a page, however the second plot
replaces the first one (although it only takes up the upper half as it
should):
par(mfrow=c(2,1))
a<-sin(1:100)
b<-cos(1:100)
lag.plot(a)
lag.plot(b)
What's the trick to this?
I'm using R 2.2.1 (2005-12-20 r36812) on Ubuntu Linux.
Thanks,
Gad
--
Gad Abraham
Department of Mathematics and
2006 Jun 23
2
Time series labeling with Zoo
Hi,
I'm using zoo because it can automatically label the months of a time
series composed of daily observations.
This works well for certain time series lengths, but not for others, e.g.:
While:
> library(zoo)
> plot(zoo(runif(10), as.Date("2005-06-01") + 0:50))
Shows up the months and day of month,
> plot(zoo(runif(10), as.Date("2005-06-01") + 0:380))
2008 Aug 12
1
Conflict between octave3.0-headers and r-base-dev
Hi,
I'm using Ubuntu Hardy i386, and the R 2.7.1 Ubuntu packages from CRAN.
The CRAN r-base-dev package depends on refblas3-dev or atlas3-base-dev,
but octave3.0-headers depends on libblas-dev. The two blas packages seem
to conflict, so r-base-dev and octave3.0-headers can't be installed
together.
Is this a known issue?
Thanks,
Gad
$ sudo apt-get install octave3.0-headers
Reading
2009 Nov 27
2
Symmetric Matrix classes
Hi,
I'd like to store large covariance matrices using Matrix classes.
dsyMatrix seems like the right one, but I want to specify just the
upper/lower triangle and diagonal and not have to instantiate a huge
n^2 vector just for the sake of having half of it ignored:
Dumb example:
M <- new("dsyMatrix", uplo="U", x=rnorm(1e4), Dim=as.integer(c(100, 100)))
diag(M) <- 1
2008 May 18
1
predict.prcomp: 'newdata' does not have the correct number of columns
Hi,
I'm doing PCA on wide matrices and I don't understand why calling
predict.prcomp on it throws an error:
> x1 <- matrix(rnorm(100), 5, 20)
> x2 <- matrix(rnorm(100), 5, 20)
> p <- prcomp(x1)
> predict(p, x2)
Error in predict.prcomp(p, x2) :
'newdata' does not have the correct number of columns
> dim(x2)
[1] 5 20
> dim(p$rotation)
[1] 20 5
2008 May 16
1
Dimensions of svd V matrix
Hi,
I'm trying to do PCA on a n by p wide matrix (n < p), and I'd like to
get more principal components than there are rows. However, svd() only
returns a V matrix of with n columns (instead of p) unless the argument
nv=p is set (prcomp calls svd without setting it). Moreover, the
eigenvalues returned are always min(n, p) instead of p, even if nv is set:
> x <-
2010 May 26
2
cran2deb Packages.bz2 Hash Sum mismatch
Hi,
I'm using Ubuntu 9.10 amd64 with cran2deb packages (I know it's meant
for Debian, it's worked fine for me for many months).
Recently apt-get update has to started to complain:
...
Hit http://debian.cran.r-project.org testing/ Release
Get:1 http://debian.cran.r-project.org testing/ Packages [515kB]
Fetched 1B in 3s (0B/s)
W: Failed to fetch
2009 Mar 10
1
S4 generic masking S3 generic when using namespace
Hi,
I have two example packages, test1 and test2, where the only code in
them is:
setGeneric("predict", function(object, ...) standardGeneric("predict"))
(get them from http://www.cs.mu.oz.au/~gabraham/test1.tar and
http://www.cs.mu.oz.au/~gabraham/test2.tar)
The difference between them is that first does not have a namespace, and
loads fine. The second has a namespace
2011 Nov 23
1
Measure of separation for survival data
Dear all,
I am using R 2.9.2 on Windows XP.
I am undertaking a simulation study to consider methods of external validation in the context of missing covariates in the validation data set. I would like to use Royston's measure of prognostic separation as a method of external validation. Although there is no R code for this, the author informs me that it should be easy:
1. Calculate the
2008 Apr 03
1
Lapack error in Design:::ols
Hi,
I'm trying to use Frank Harrell's Design:::ols function to do regression
of y (numeric) on the interaction of two factors (x1 and x2), but Lapack
throws an error:
> library(Design)
...
> load(url("http://www.csse.unimelb.edu.au/~gabraham/x"))
> ols(y ~ x1 * x2, data=x)
Error in chol2inv(fit$qr$qr) : 'size' cannot exceed nrow(x) = 20
> traceback()
2007 Dec 07
1
Make natural splines constant outside boundary
Hi,
I'm using natural cubic splines from splines::ns() in survival
regression (regressing inter-arrival times of patients to a queue on
queue size). The queue size fluctuates between 3600 and 3900.
I would like to be able to run predict.survreg() for sizes <3600 and
>3900 by assuming that the rate for <3600 is the same as for 3600 and
that for >4000 it's the same as for