Displaying 3 results from an estimated 3 matches for "logq".
Did you mean:
log
2013 Mar 21
0
"[[i]]$" <- "" indexing and lapply
...ian = c(55060, 55061, 55062, 55063, 55064,
55065), Month = c(10, 10, 10, 10, 10, 10), Day = c(275, 276,
277, 278, 279, 280), DecYear = c(2000.75, 2000.75273224044,
2000.75546448087,
2000.75819672131, 2000.76092896175, 2000.76366120219), MonthSeq =
c(1810,
1810, 1810, 1810, 1810, 1810), i = 1:6, LogQ = c(-0.120434728772853,
-0.153224551595844, -0.187126103271526, -0.222217423082796,
-0.258585067253671,
-0.296325395236517)), .Names = c("Date", "Q", "Julian", "Month",
"Day", "DecYear", "MonthSeq", "i", "LogQ&q...
2009 Sep 06
2
question about ... passed to two different functions
...### Markov chain we want to run (log likelihood + log prior)
ludfun <- function(beta) {
stopifnot(is.numeric(beta))
stopifnot(length(beta) == ncol(modmat))
eta <- as.numeric(modmat %*% beta)
logp <- ifelse(eta < 0, eta - log1p(exp(eta)), - log1p(exp(- eta)))
logq <- ifelse(eta < 0, - log1p(exp(eta)), - eta - log1p(exp(- eta)))
logl <- sum(logp[y == 1]) + sum(logq[y == 0])
val <- logl - sum(beta^2) / 2
return(val)
}
beta.initial <- as.vector(out$coefficients)
out <- metrop(ludfun, initial = beta.initial, nbatch = 20,...
2001 Jun 11
2
par
Dear all, I have a question about "par" options.
When I choose par(pty="s",mfrow=c(2,2)), how can I change the horizontal distance between
the plots????
I have tried to use some options, like "mai", "mar", "oma", "omd", so on but I didn''t get the expected result.
Some suggest will be appreciated.
JM.