On Tue, 25 Aug 2009 12:13:21 +0200, Uwe Ligges
<lig...@statistik.tu-dortmund.de> wrote:
Vitalie S. wrote:
Hello Uwe,
Just a related question, OpenBugs is using C and T for Censoring and
Truncation. But this does not seem to work with BRugs.
I am using this document
http://mathstat.helsinki.fi/openbugs/Manuals/ModelSpecification.html#TheBUGSLanguageStochasticNodes
What's the problem? Error message? Reproducible code?
Oh, sorry,I justs assumed it would be an easy answer, like "not
implemented yet".
Here is the code:
fModel<- function()
{
beta ~ dnorm(0,1) %_% T(0, 1)
}
writeModel(fModel,"Model.txt");
modelCheck("Model.txt");
#ERROR: this density cannot be truncated error pos 33 (error on line 3)
Note that C and I work fine in the code above.
My info:
sessionInfo("BRugs")
R version 2.9.0 (2009-04-17)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
character(0)
other attached packages:
[1] BRugs_0.5-1
Vitalie.
Specifying, e.g.
model <- function()
x ~ dnorm(mu, tau)%_%T(lower, upper)
writeModel(model)
works for me.
Best,
Uwe
Thanks,
Vitalie.
On Tue, 25 Aug 2009 10:00:44 +0200, Uwe Ligges
<lig...@statistik.tu-dortmund.de> wrote:
R Heberto Ghezzo, Dr wrote:
Hello, I am sorry, I have this problem before and Uwe send me the
answer but I misplaced it
Oh dear! But it is not lost, since the answer to the second part of
your problem was: Please read the documentation!
and can not find it.
> writing a model for BRugs
library(BRugs)
Loading required package: coda
Loading required package: lattice
Welcome to BRugs running on OpenBUGS version 3.0.3
setwd("c:/tmp")
Error in setwd("c:/tmp") : cannot change working directory
So c:/tmp does not exist or you don't have permissions there?
mo <- function(){
+ for (k in 1:p){
+ delta[1,k] ~ dnorm(0,0.1)I(,delta[2,k])
Error: unexpected symbol in:
The answer was: Please read the documentation!
?writeModel (which you want to use on this function) tells you:
"As a difference, BUGS syntax allows truncation specification like
this: dnorm(...) I(...) but this is illegal in R. To overcome this
incompatibility, use %_% before I(...): dnorm(...) %_% I(...). The
dummy operator %_% will be removed before the BUGS code is saved. "
" for (k in 1:p){
delta[1,k] ~ dnorm(0,0.1)I"
delta[2,k] ~ dnorm(0,0.1)I(delta[1,k],delta[3,k])
Error: unexpected symbol in " delta[2,k] ~ dnorm(0,0.1)I"
delta[3,k] ~ dnorm(0,0.1)I(delta[2,k],)}
Error: unexpected symbol in " delta[3,k] ~ dnorm(0,0.1)I"
}
Error: unexpected '}' in "}"
so R parser does not like the I(,) construct
which is *not* the problem, see above.
Best wishes,
Uwe
, What is the alternative way of propgramming the
constrain I(lower,upper)
Thanks
Heberto Ghezzo
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
--
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
