similar to: stop criteria when "L-BFGS-B needs finite values of 'fn' " in optim

Hi all ! I would like to simulate spatial lattice/areal data with a conditional autoregressive (CAR) structure, for a given neighbouring matrix and for a autocorrelation "rho". Is there any package or function in R to perform it ? I found the function "CARsimu" in the hdeco library, but this is not what I'm looking for Thanks in advance Dae-Jin --

Hi, When I run the following code, r <- c(3,4,4,3,5,4,5,9,8,11,12,13) n <- rep(15,12) x <- c(0, 1.1, 1.3, 2.0, 2.2, 2.8, 3.7, 3.9, 4.4, 4.8, 5.9, 6.8) x <- log10(x) fr <- function(c, alpha, beta) { P <- c + (1-c) * pnorm(alpha + beta * x) P <- pmax(pmin(P,1),0) -(sum(log(choose(n,r))) + sum(r * log(P)) + sum((n -r)* log(1-P))) } fit <- mle((fr), start = list(c

Dear All, I am trying to solve the optimization problem below, but I am always getting the following error: Error in optim(rep(20, nvar), f, gr, method = "L-BFGS-B", lower = rep(0, : L-BFGS-B needs finite values of 'fn' Any ideas? Thanks in advance, Paul ----------------------------------------------- k <- 10000 b <- 0.3 f <- function(x) { n <- length(x)

Hi all ! I've been trying to maximize a likelihood using optim( ) function, but it seems that the function has several local maxima. I've tried in my algorithm with different starting values and depending on them "optim" obtains different results... I use the "L-BFGS-B" method setting the lower values as 1e-06, because my parameters must be strictly positive. Also

Hello everyone! I'm trying to get tinc running, but couldn't make it until now... My current situation is: all tincs are connected, but there ist (almost) no traffic: Only two PCs are able to ping to another, but not to itself. Let me explain, what I'm planning to do: I've got four networks: PROM: 10.69.0.0/16 PAN: 192.168.30.0/24 (I'm afraid, I can't change PAN to

Vince, When that happens, one possible reason is that your DAE is of index > 1, which cannot be solved by daspk. The solver radau, also from deSolve can handle DAEs up to index 3, but you need to rewrite the problem in the form M*y' = f(x,y), where M is a mass matrix. If you do that for your problem, and solve it with radau, then radau complains that the "matrix is repeatedly

Hola, De esta manera sí funciona: #----------------------------------------------- #................... vec1a <- as.vector(1:3) n1a <- 3 n2a <- 3 #P1a <- matrix(rnorm(30),nrow=3) #P2a <- matrix(rnorm(30),nrow=3) Ka <- 10 n.val <- 30 # En forma de funcion producto <- function(n1a,n2a,Ka,vec1a,n.val){ # A eliminar una vez compruebes que la vectorización funciona

Dear Sirs, We, Anaya MUltimedia, are a Spanish publisher of computing and IT books.We have a catalogue of around 600 titles, growing at a rate of 120 new titles per year and all our titles are classified in different series, each one addressed to a different sector of the market. (IF you would like to take a look at our catalogue, you are very wellcome to visit our web site at the following

Hello. I'm from a Linux User Group in Madrid, Spain. We are planning to organize a series of speeches on November 14th. We would like someone from Freedektop to hold one, about whatever you feel feel like (what is Freedesktop about, xdg, X servers...). Are any of you interested? Maybe this is not the best place to ask, but I've not found any other contact address. Is there is a better

Hola a tod en s, me gustaría usar solo un(os) comando(s) específico(s) de una librería pero sin cargarla entera a la manera: library(laquesea) ¿se puede hacer de alguna manera? Y saludos y gracias jm~ _______________________________ J. Miguel Marin http://www.est.uc3m.es/jmmarin Dep. of Statistics University Carlos III of Madrid Spain (E.U.)

I'm getting this error on the attached code and breaking my head but can't figure it out. Any help is much appreciated. Thanks, Vince CODE: library(deSolve) Res_DAE=function(t, y, dy, pars) { with(as.list(c(y, dy, pars)), { res1 = -dS -dES-k2*ES res2 = -dP + k2*ES eq1 = Eo-E -ES eq2 = So-S -ES -P return(list(c(res1, res2, eq1, eq2))) }) } pars <- c(Eo=0.02,

Hola, tengo que leer más de 2000 ficheros en un programa largo de estimación sobre imágenes médicas, con una estructura como ésta: desde "IM-0005-0001.dcm" hasta "IM-0005-2021.dcm" Si hago esto: ui <- NULL for (i in 1:8){ ui[i] <- paste("IM-0005-", i,".dcm", sep="") } ui Obtendría: [1] "IM-0005-1.dcm"

Hola, con vistas a docencia, me gustaría saber si hay alguna librería que haga el test de la chi cuadrado para variables continuas, que no sean la normal. Lo encontré para variables discretas y para la normal. También, obviamente, se puede hacer un pequeño programa que lo calcule "como a mano", pero a los alumnos les obliga a entrar en detalles de programación... :-/ La idea sería

Hola a tod en s, tengo que hacer una operación con matrices que lleva un doble bucle. He intentado vectorizarlo pero sin mucho éxito con la función "Vectorize". ¿sabríais de alguna manera de evitar los bucles o de que funcionase Vectorize? Adjunto un ejemplo #................... vec1a <- as.vector(1:3) n1a <- 3 n2a <- 3 P1a <- matrix(rnorm(30),nrow=3) P2a <-

David, Peter, Let me try it one more time, now with an example… Short of upstreaming a fake function clone pass, let me only illustrate the issue in raw debug output. Hopefully it should be sufficient. Try this: clang -Os -g -fno-strict-aliasing test_arg_del.ll -mllvm -debug -mllvm -print-after-all After DAE - Removing argument 1 (reg) from foo DAE - Removing argument 2

Hi everyone, I am interested in estimating this type of random effects panel: y_it = x'_it * beta + u_it + e_it u_it = rho * u_it-1 + d_it rho belongs to (-1, 1) where: u and e are independent and normally zero-mean distributed. d is also independently normally zero-mean distributed. So, I want random effects for group i to be correlated in t, following an AR(1) process. I am

I am trying to follow this example for multinomial logistic regression http://www.ats.ucla.edu/stat/r/dae/mlogit.htm However, I cannot get it to work properly. This is the output I get, and I get an error when I try to use the mlogit function. Any ideas as to why this happens? > mydata <- read.csv(url("http://www.ats.ucla.edu/stat/r/dae/mlogit.csv")) > attach(mydata) >

Hi, I am starting tu use R on a linux environment (Red Hat 6.0), and when I tried to use library MASS, writting library(MASS) I got the following message: > library(MASS) Error: dynamic library `MASS' not found > Is there something I am doing wrong? Thanks in advance, Alberto --------------------------------------------------------------------- Alberto Munoz

Hola, ¿qué tal? attach y detach son funciones cuyas páginas de ayuda están llenas de avisos y caveats además de sugerencias acerca de alternativas a su uso. Yo siempre recomiendo no utilizarlas y yo en particular, las uso, si acaso, de manera sumamente excepcional. Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 21 de noviembre de 2013 13:18, Juan Bautista <jbautista

Hi, all Does anybody know if it is possible to include TeX-style formulae in R plots? The aim is to export such plots to LaTeX documents using the postscript() device. Many thanks in advance, Alberto Munoz --------------------------------------------------------------------- Alberto Munoz Phone: +34-91- 624 95 79 Dpto. de Estadistica y Econometria Fax: +34-91- 624