Displaying 5 results from an estimated 5 matches for "nzero".
Did you mean:
zero
2012 Apr 26
0
Problem with constrOptim when hitting boundary
...Muts) stop("x and m must be of same length")
if(coverageCutoff<1) stop("coverageCutoff must be at least 1")
if(any(m<coverageCutoff)){
if(verbose) warning("some mutations have no reads, will be ignored")
nBelowCutoff<-length(which(m<coverageCutoff))
nZero<-length(which(m==0))
x<-x[which(m>=coverageCutoff)]
m<-m[which(m>=coverageCutoff)]
}
else{
nBelowCutoff<-0
nZero<-0
}
if(length(m)<minMutations){
if(length(normCont)>0){
if(length(possAlleles)!=K+1) stop("invalid length for possible Allele Frequency...
2003 Aug 11
0
Designing and incorporating a digital filter
...lter using a web-
site. The site returns the filter in "ANSI C" source code. It is:-
/* Digital filter designed by mkfilter/mkshape/gencode A.J. Fisher
Command line: /www/usr/fisher/helpers/mkfilter -Bu -Bp -o 5 -a
1.0000000000e-02 1.0000000000e-01 -Z 2.0000000000e-01 -l */
#define NZEROS 12
#define NPOLES 12
#define GAIN 1.548068051e+03
static float xv[NZEROS+1], yv[NPOLES+1];
static void filterloop()
{ for (;;)
{ xv[0] = xv[1]; xv[1] = xv[2]; xv[2] = xv[3]; xv[3] = xv[4]; xv[4] =
xv[5]; xv[5] = xv[6]; xv[6] = xv[7]; xv[7] = xv[8]; xv[8] = xv[9]; xv[9] =
xv[10]; xv[10]...
2008 May 30
1
robust mlm in R?
I'm looking for something in R to fit a multivariate linear model
robustly, using
an M-estimator or any of the myriad of other robust methods for linear
models
implemented in robustbase or methods based on MCD or MVE covariance
estimation (package rrcov).
E.g., one can fit an mlm for the iris data as:
iris.mod <- lm(cbind(Sepal.Length, Sepal.Width, Petal.Length,
Petal.Width) ~ Species,
2006 Nov 24
2
low-variance warning in lmer
...c(-0.02,0.55),axes=FALSE)
axis(side=2)
box()
boxplot(rfrac~rvec,at=unique(rvec),add=TRUE,pars=list(boxwex=0.03),
col="gray")
points(jitter(rvec),rfrac,col=X[,"warn"]+1)
lines(unique(rvec),fracwarn,col="blue",type="b",lwd=2)
text(0.4,0.1,"fraction\nzero",col="blue")
abline(a=0,b=1,lwd=2)
points(unique(rvec),est.mean,cex=1.5,col=5)
##
plot(rvec,rfrac,type="n",xlim=c(-0.02,0.55),axes=FALSE,log="y")
axis(side=2)
axis(side=1,at=unique(rvec))
box()
points(jitter(rvec),rfrac,col=X[,"warn"]+1)
curve(1*x,add=TR...
2005 Jan 27
2
Request for help
My name is Michela Marignani and I'm an ecologist trying to solve a problem
linked to knight' s tour algorithm.
I need a program to create random matrices with presence/absence (i.e. 1,0
values), with defined colums and rows sums, to create null models for
statistical comparison of species distribution phenomena.
I've seen on the web many solutions of the problem, but none provides