Displaying 2 results from an estimated 2 matches for "nresample".
Did you mean:
resample
2007 May 25
1
Speeding up resampling of rows from a large matrix
...over correlated SNPs of a simple statistic.
The code:
#-------------------------------------------------------------------------------
nSNPs <- 1000
H <- matrix(sample(0:1, 120*nSNPs , replace=T), nrow=120)
G <- matrix(0, nrow=3, ncol=nSNPs)
# Keep in mind that the real H is 120 x 65000
nResamples <- 3000
pair <- replicate(nResamples, sample(1:120, 2))
gen <- function(x){g <- sum(x); c(g==0, g==1, g==2)}
for (i in 1:nResamples){
G <- G + apply(H[pair[,i],], 2, gen)
}
#-------------------------------------------------------------------------------
The problem is that the...
2009 Mar 12
1
zooreg and lmrob problem (bug?)
...:
4 weights are ~= 1. The remaining 54 ones are summarized as
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.5522 0.8858 0.9537 0.9212 0.9831 0.9990
Algorithmic parameters:
tuning.chi bb tuning.psi refine.tol rel.tol
1.5476400 0.5000000 4.6850610 0.0000001 0.0000001
nResample max.it groups n.group best.r.s k.fast.s
k.max
500 50 5 400 2 1
200
trace.lev compute.rd
0 0
seed : int(0)
>
> mod1<-lmrob(as.zooreg(dad)~seq(58))
> summary(mod1)
Call:
lmrob(formula = as.zooreg(dad...