Subject: Re: [R] how to avoid NaN in optim()
Here is how you do it:
library(BB)
Amat - matrix(c(1,0,0,1,-1,-1), 3, 2, byrow=TRUE)
b - c(0, 0, -1)
p0 - c(0.5, 0.4)
spg(p0, lik ( 176,182 , 60 ,17) , project=projectLinear,
projectArgs=list(A=Amat, b=b, meq=0))
Hope this helps,
Ravi
hi ,
lik - function(nO, nA, nB, nAB){
loglik - function(par)
{
p=par[1]
q=par[2]
r - 1 - p - q
if (c(p,q,r) rep(0,3) c(p,q,r) rep(1,3) )
{
-(2 * nO * log (r) + nA * log (p^2 + 2 * p * r)
+ nB * log (q^2 + 2 * q * r)
+ nAB *
: [R] how to avoid NaN in optim()
To: r-help@r-project.org
hi ,
lik - function(nO, nA, nB, nAB){
loglik - function(par)
{
p=par[1]
q=par[2]
r - 1 - p - q
if (c(p,q,r) rep(0,3) c(p,q,r) rep(1,3) )
{
-(2 * nO * log (r) + nA
Subject: [R] how to avoid NaN in optim()
To: r-help@r-project.org
hi ,
lik - function(nO, nA, nB, nAB){
loglik - function(par)
{
p=par[1]
q=par[2]
r - 1 - p - q
if (c(p,q,r) rep(0,3) c(p,q,r) rep(1,3) )
{
-(2 * nO * log (r
of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins University
Ph. (410) 502-2619
email: rvarad...@jhmi.edu
- Original Message -
From: Ravi Varadhan rvarad...@jhmi.edu
Date: Thursday, September 30, 2010 2:54 pm
Subject: Re: [R] how to avoid NaN in optim()
To: arindam fadikar
- Original Message -
From: arindam fadikar arindam.fadi...@gmail.com
Date: Thursday, September 30, 2010 2:17 pm
Subject: [R] how to avoid NaN in optim()
To: r-help@r-project.org
hi ,
lik - function(nO, nA, nB, nAB){
loglik - function(par)
{
p=par[1]
q=par[2
3:04 pm
Subject: Re: [R] how to avoid NaN in optim()
To: Ravi Varadhan rvarad...@jhmi.edu
Cc: arindam fadikar arindam.fadi...@gmail.com, r-help@r-project.org
You also need the constrain that par[1] + par[2] 1 in order to avoid
NaNs.
You can do this using the `projectLinear' argument
) 502-2619
email: rvarad...@jhmi.edu
- Original Message -
From: arindam fadikar arindam.fadi...@gmail.com
Date: Thursday, September 30, 2010 2:17 pm
Subject: [R] how to avoid NaN in optim()
To: r-help@r-project.org
hi ,
lik - function(nO, nA, nB, nAB){
loglik - function
arindam fadikar wrote:
loglik - function(par)
{
p=par[1]
q=par[2]
r - 1 - p - q
if (c(p,q,r) rep(0,3) c(p,q,r) rep(1,3) )
{
-(2 * nO * log (r) + nA * log (p^2 + 2 * p * r)
+ nB * log (q^2 + 2 * q * r)
+ nAB *
9 matches
Mail list logo