1"lambda" <- c(535,NA)
2"theta" <- 5
3"tau" <-  0.1
4"T" <- c(1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,
5         1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2,1,2)
6