close Warning: Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.

Opened 15 years ago

Closed 15 years ago

#810 closed defect (fixed)

Duplicate parameter at some primary blocks

Reported by: pgea@… Owned by: Víctor de Buen Remiro
Priority: normal Milestone: BSR API
Component: Math Version:
Severity: major Keywords:
Cc:

Description

En MMS se consideró la posibilidad de que un mismo término
explicativo (un parámetro por un input) participase en más
de un nodo observacional con el mismo parámetro.

Este tipo de parámetros que podríamos denominar parámetro internodal
(en analogía con internacional) se intenta utilizar en BSR como
otro parámetro cualquiera simplemente especificando el mismo identificador
en todos los nodos de tipo "primary" que lo usan.

Imaginemos que hay dos ouptuts con el mismo parámetro:
Common__Input.Coef
así pues en los archivos primary.bsr de cada submodelo aparecerá:

...
/////////////////////////////////////////////////////////////////////////////
// Defining header of Regression Variables
/////////////////////////////////////////////////////////////////////////////
...
Common__Input.Coef <- +0.1;
...

Al utilizar esto en BSR encontramos errores del tipo:

ERROR: [1] Nombre de índice duplicado Common__Input.Coef
[Call stack]
 [7] Set StdLib::BysMcmc::Bsr::ModularSchema::Monophasic::MasterToJoint (Set master)
 [6] NameBlock StdLib::BysMcmc::Bsr::Gibbs::GetCyclerFromParsed (Set parsed_, NameBlock mcmcPartialHandler, NameBlock config, NameBlock notifier)
 [5] NameBlock StdLib::BysMcmc::Bsr::Gibbs::GetCyclerFromAscii (Text ascii_, Text moduleType, NameBlock mcmcPartialHandler, NameBlock config, NameBlock notifier)
 [4] NameBlock StdLib::BysMcmc::Bsr::Gibbs::GetCycler (Anything modelDef, NameBlock mcmcPartialHandler, NameBlock config, NameBlock notifier)
 ...

y termina con un:

ERROR: [8] No es posible aplicar CholeskiFactor a una matriz virtual no definida positiva  Cholmod.R.Sparse(839x58)

No sé si se trata de un error en BSR o de un problema de diseño o definición del modelo. Adjunto un ejemplo.

Attachments (1)

Ejemplo_ParametrosInternodales.zip (58.0 KB) - added by pgea@… 15 years ago.

Download all attachments as: .zip

Change History (3)

Changed 15 years ago by pgea@…

comment:1 Changed 15 years ago by Víctor de Buen Remiro

Status: newaccepted

comment:2 Changed 15 years ago by Víctor de Buen Remiro

Resolution: fixed
Status: acceptedclosed

(In [1821]) Fixes #810

Note: See TracTickets for help on using tickets.