similar to: Probability(Markov chain transition matrix)

Hi, Does anybody know a function to simulate a Markov chain given a probability transition matrix and an initial state ? Thanks. Philippe -- -------------------------------------------------- Philippe Hup? Institut Curie - Equipe Bioinformatique 26, rue d'Ulm - 75005 PARIS France +33 (0)1 42 34 65 29 Philippe.Hupe at curie.fr <mailto:Philippe.Hupe at curie.fr>

Does anyone have any "R" code for computing the entropy of a simple first or second order Markov chain, given a transition matrix something like the following (or the symbol vector from which it is computed)? AGRe ARIe CSRe DIRe DSCe eos HRMe SPTe TOBe AGRe 0.0000000 0.0000000 0.0000000 0.0000000 1.0000000 0.0000000 0.0000000 0.0000000

Dear R users, I have been trying to obtain the MLE of the following model state 0: y_t = 2 + 0.5 * y_{t-1} + e_t state 1: y_t = 0.5 + 0.9 * y_{t-1} + e_t where e_t ~ iidN(0,1) transition probability between states is 0.2 I've generated some fake data and tried to estimate the parameters using the constrOptim() function but I can't get sensible answers using it. I've tried using

Dear kind R helpers, I have a vector of runway names in rwy ("31R", "31L",... the number is user selectable) arrgnd is a data frame with data for all flights and all runways, with a Runway column. I am trying to subset arrgnd into a dat frame for each selected runway, and then combine them back together using the following code: for (j in 1:nr) { # nr = number of

Dear R gurus, I have data for which I want to estimate the markov transition matrix that generated the sequence, and preferably obtain some measure of confidence for that estimation. e.g., for a series such as 1 3 4 1 2 3 1 2 1 3 4 3 2 4 2 1 4 1 2 4 1 2 4 1 2 1 2 1 3 1 I would want to get an estimate of the matrix that generated it [[originally: [,1] [,2] [,3] [,4] [1,] 0.00 0.33 0.33

Hello all, I am currently attempting to replicate data from a political science article that utilized a Markov chain transition model to predict voter turnout intention at time *t*; the data was separated into two different models based on whether prior intent was to vote or not to vote. The details don't really matter. Mostly I am curious how to run a Markov chain transition model in R,

I solved the problem in one more (and more elegant) way. So here's the program again. Where does R stand on the Anderson-Goodman test of 1957? I hunted around and nobody seems to be doing this in R. Is it that there has been much progress after 1957 and nobody uses it anymore? # Problem statement: # # You are holding a dataset where firms are observed for a fixed # (and small) set of years.

Hi all, I am working for nest box occupancy data for birds and would like to construct a Markov transition matrix, to derive transition probabilities for ALL years of the study (not separate sets of transition probabilities for each time step). The actual dataset I'm working with is 125 boxes over 14 years that can be occupied by 7 different species, though I have provided a slimmed down

Hola a todos: Estoy intentado realizar un modelo senoidal de unos datos de radiación solar con el fin de afrontar el relleno de la serie y aplicar la prueba de Thom para verificar su homogeneidad [0]. De momento me encuentro con los siguientes problemas: 1- ¿Existe la prueba de Thom en R? ¿O debo crearme mi propia función? 2- Para la realización del modelo senoidal estoy siguiendo los pasos

Folks, I am holding a dataset where firms are observed for a fixed (and small) set of years. The data is in "long" format - one record for one firm for one point in time. A state variable is observed (a factor). I wish to make a markov transition matrix about the time-series evolution of that state variable. The code below does this. But it's hardcoded to the specific years that I

Hi, Is there a package that will estimate simple Markov models and hidden Markov models for discrete time processes in R? Thanks in advance, David -- =============================================================== David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room, 1061 1025 W. Johnson Street Madison, WI 53706

Estimados, Muchas gracias por sua ayuda, ya pude instalar data.table y esta funcionando. Respecto al paquete biodiversityR, no se porque ahora no me corre. El script esta correcto, pero no se porque me tira error. > fase<-diversitycomp(biobsp, y=biopcat,factor1="Fase",index="richness",method=?s") Error in diversitycomp(biobsp, y = biopcat, factor1 =

Dear Users, i have this time series, the tree lines means different level, i would use a Markov switching model with two states to modelling this time series. i would obtain the relative transition matrix (2X2) the first state is above the value of 23.65 (the higher line) the second state is below the value of 23.65 You can ignore the other two lines

Hi, I am looking for a package to estimate regime switching models (states following a markov chain). I found packages for Hidden Markov Models but I am looking for something a little different: In the HMM the conditional distribution of the observations (give the state) is a known distribution (normal or others), while the package I need should allow to set a conditional distribution (given the

Hello, can any one point me to R-packages (if available) which include spatial Markov Chain methods? My second question is more general but hopefully not OT: Currently we are using the software TPROGS, which let people simulate property distributions in space by some Markov Chain approaches. We face some problems due to the lack of information between distances of samples along borehole path

Dear helpers, Does anyone know about a package or a function that allows to estimate Multivariate Markov-Switching Models, like MS-VAR as introduced by Krolzig(1997) with R ? Thanks a lot!! Carlo

Hi all, I am new to R. Can someone please give me some hints in how to do the following things: 1- Get ONE permutation of a set. I have looked at the gregmisc package's permutations() method, but I just want to get one permutation at a time. 2- Simulate a Markov chain in R. For instance, I want to simulate the simple random walk problem, in which a person can walk randomly around 4 places.

Greetings fellow R entusiasts! We have some problems converting a computer routine written initially for Gauss to estimate a Markov Regime Switching analysis with Time Varying Transition Probability. The source code in Gauss is here: http://www.econ.washington.edu/user/cnelson/markov/programs/hmt_tvp.opt We have converted the code to R, and it's running without errors, but we have some

Dear help group, Just fyi a markov chain is a sequence of transitions between states (say A,T,G,C - on a gene) with a given probability for each transition. In this case there'd be 16 different kinds, each with a different weight. Given a transition matrix (4x4) filled with all transition probabilities of course, how can I generate a random sequence of a given length, say 200?

Hi all, I'm looking for a package that will estimate Markov models and provide transition probabilities. I'm not speaking of MCMM estimation packages. Thanks in advance, David -- ======================================================================= David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room 1061