search for: deltab

Displaying 5 results from an estimated 5 matches for "deltab".

Did you mean: delta
2012 Oct 23
1
repeat takes along time
Hi I want help for "repeat" because it takes a long time repeat{ for (i in 1:n){ probb[i]=sum(Wc[z[,j]>=yb[i]]) } deltab[,b]=rbinom(n,1,probb) if(length(which(is.na(deltab[,b])==T))==0){break} } This code works but takes much computation time. Anyone could help me to find a solution to take a shor time. Thanks Jorge -- View this message in context: http://r.789695.n4.nabble.com/repeat-takes-along-time-tp464716...
2012 Apr 25
1
Using apply() with a function involving ode()
...k. Here is an example. Consider an epidemiology model of West Nile Virus: wnv.model <- function(Time, State, Pars){ with(as.list(c(State, Pars)), { dSB <- -(alphaB*betaB*SB*IM)/(SB + IB) dIB <- (alphaB*betaB*SB*IM)/(SB + IB) - deltaB*IB dSM <- bM*(SM + EM + IM) - (alphaM*betaB*IB*SM)/(SB + IB) - dM*SM dEM <- (alphaM*betaB*SM*IB)/(SB + IB) - kappaM*EM - dM*EM dIM <- kappaM*EM - dM*IM return(list(c(dSB, dIB, dSM, dEM, dIM)))...
2012 Sep 02
1
glmulti runs indefinitely when using genetic algorithm with lme4
...odel: y~1 Crit= 161.038899734164 Mean crit= 164.13629335762 Change in best IC: 0 / Change in mean IC: 0 After 19570 generations: Best model: y~1 Crit= 161.038899734164 Mean crit= 164.13629335762 etc. I have tried using calls that tell glmulti when to stop (deltaB = 0, deltaM = 0.01, conseq = 6) but nothing seems to work. I think the problem must lie with setting the function (?). It may be something really obvious however I'm new to R and I can't work it out. I am using R v.2.15.0, glmulti v.1.0.6 and lme4 v.0.999999-0 on Windows. Any help w...
2012 Aug 09
4
debug vs regular mode
Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the
2012 Aug 09
4
debug vs regular mode
Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the