Displaying 9 results from an estimated 9 matches for "npop".
Did you mean:
nop
2010 May 20
1
Geneland error on unix: Error in MCMC(........ :, unused argument(s) (ploidy = 2, genotypes = geno)
...R","FH","SPP","STP")
+ for (sp in SPP.CODES){
+ path.sp<- paste(theWd,sp,"/",sep="")
+ dir.create(path.sp)
+ GENO.TABLE<-paste(theWd,sp,"geno",sep="")
+ XY.TABLE<-paste(theWd,sp,"xy",sep="")
+ MINPOP=1
+ MAXPOP=25
+ INITPOP=1
+ NITS=500000
+ THIN=NITS/1000
+ nrun <- 10
+ burnin <- 200
+ geno<-noquote(read.table(GENO.TABLE))
+ coord<-read.table(XY.TABLE)
+
+ ## Loop for multiple runs
+
+ for(irun in 1:nrun)
+ {
+ ## define path to MCMC directory
+
+ path.mcmc <- paste(path.sp,irun...
2010 Mar 11
2
Comparing matrices
...#39;d
rather be sure.
------- code ------------
# create a binary vector of size "len"
create_bin_Chromosome <- function(len)
{
sample(0:1, len, replace=T)
}
# create popsize members, each of length len
create_pop_2 <- function(popsize, len)
{
datasize=len*popsize
npop <- matrix(0, popsize, len, byrow=T)
for(i in 1:popsize)
npop[i,] = create_bin_Chromosome(len)
npop
}
POP_SIZE = 3
LEN = 8
pop = create_pop_2(POP_SIZE, LEN)
pop2 = pop
print(pop==pop2)
if (pop==pop2) { cat('equal\n') } else { cat('NOT equal\n') }
print(pop)
print...
2004 Jun 20
1
Sweave and echoing R comments
Is there any way to echo comments from an R source file into an
SWeave->LaTeX document?
Example:
# Npop is population total
# Npoph0..Npoph2 are stratum totals
# Npoph is vector of stratum totals
Npop<-sum(to2000)
Npoph0<-sum(to2000[bg==0])
Npoph1<-sum(to2000[bg==1])
Npoph2<-sum(to2000[bg==2])
Npoph<-c(Npoph0,Npoph1,Npoph2)
In the final LaTeX document, I'd like the comments t...
2009 Apr 26
3
3 questions regarding matrix copy/shuffle/compares
...function(len)
{
sample(0:1, len, replace=T)
}
############## create_population ###################
# create population of chromosomes of length len
# the matrix contains twice as much space as popsize
#
create_pop_2 <- function(popsize, len)
{
datasize=len*popsize
print(datasize)
npop <- matrix(0, popsize*2, len, byrow=T)
for(i in 1:popsize)
npop[i,] = create_bin_Chromosome(len)
npop
}
2010 Dec 15
0
Problems with the borders (High difficulty)
...der? The queen2 ring should be 16 cells, why in the 'SR2[]' expression divides by 8 and not by 16?
Do anyone know if the dispersed seed are properly included in the SB for the next iteration?
Here I attach you the full model to be this message understandable:
width = 10
height= 10
endt = 6
nPop = 1
SBLL = 1
SBUL = 10
LrD <- matrix(0, width, height)
x <- as.integer(runif(nPop,min=1,max=10)) #uniform random x ccordinates for occupied cells
y <- as.integer(runif(nPop,min=1,max=10)) #uniform random y ccordinates for occupied cells
z <- cbind(x,y)
LrD[z] = runif(nPop,min=SBLL,max...
2009 Jan 13
3
problem whit Geneland
I do the these passages:
library(Geneland)
set.seed(1)
data <- simdata(nindiv=200,
coord.lim=c(0,1,0,1) ,
number.nuclei=5 ,
allele.numbers=rep(10,20),
IBD=FALSE,
npop=2,
give.tess.grid=FALSE)
geno <- data$genotypes
coord <- t(data$coord.indiv)
path.mcmc <- paste(tempdir(),"/",sep="")
set.seed(1)
mcmcFmodel(coordinates=coord,
genotypes=geno,
path.mcmc=path.mcmc,
rate.max=10,...
2005 Feb 18
1
Is this a bug or by design? Workaround?
...------------------------------------------------------------------------------
N Y H E T E R!
- IP-telefoni, spara tusenlappar om ?ret!
- Rikst?ckande ADSL 0,25-24Mbit
- Internetaccess (Modem/ISDN64+128 via Ymex - utan abonnemangskostnad!
- Eposttj?nster, ?ven UUCP, Uppringd SMTP, MX fallback, Dom?nPOP
- Surf24 - en billig bredbandstj?nst fr?n Ymex f?r kunder i H?rn?sand/?landsbro.
-------------------------------------------------------------------------------------
Get your emailed Web-forms into a database of your choice!!!
Checkout DBFORM V1.0, see details at http://www.ymex.se
-------------...
2005 Feb 21
0
bug? Unterminated comment detected beginning on line 0
...------------------------------------------------------------------------------
N Y H E T E R!
- IP-telefoni, spara tusenlappar om ?ret!
- Rikst?ckande ADSL 0,25-24Mbit
- Internetaccess (Modem/ISDN64+128 via Ymex - utan abonnemangskostnad!
- Eposttj?nster, ?ven UUCP, Uppringd SMTP, MX fallback, Dom?nPOP
- Surf24 - en billig bredbandstj?nst fr?n Ymex f?r kunder i H?rn?sand/?landsbro.
-------------------------------------------------------------------------------------
Get your emailed Web-forms into a database of your choice!!!
Checkout DBFORM V1.0, see details at http://www.ymex.se
-------------...
2005 May 08
3
1.0-test69
http://dovecot.org/test/
- Several mbox fixes, upgrade recommended for test68 mbox users
- Possibly fixes some IMAP hangs where Dovecot just stopped replying
- Fixed delay-newmail workaround. It was badly broken before.
And somewhat off topic advertisement:
I got a bit distracted from Dovecot a week ago when a guy started
mailing me about wanting to write an irssi2 client as a project to