similar to: Training parameters for a HMM

Hi r-help, I have been using your RHmm package for some time and have recently had to try using the package for a new dataset. Basically I have a dataset with a number of discrete observation variables that change over time, and I would love to try modeling them using a HMM. Basically I was wondering if RHmm can be used to model a multivariate discrete HMM, i.e., the observations are a vector

Dear List, We are using RHmm to cluster data through HMM. We would like to restrict the transition matrix of HMM, to get hierarchical connections between clusters. But, RHmm doesn't seem to support these restrictions. Can any one suggest a library to do that. Thanks, Kishor

We are trying to build a human respiration model. Preliminary analysis of some breathing signals has shown that humans breathe through switching among a finite number of patterns. Hidden Markov seems to be the right approach. Since most of our code is written in R scripting language, finding an R package implementing an HMM that we can use for our prototype would be very helpful. I have been

Thank you for your prompt answer. The breathing signal observations are the amplitude values as a function of time and phase. According to our model the hidden states are the different breathing types. Subjects, whose respiratiion process is regular, are likely to breathe, keeping the same cycle pattern/type, for many consecutive cycles. therefore dwelling in the same hidden state. The more

I admire package-builders being proactive and having their facilities ready for upcoming R releases. However, if the publicly released version of R is 2.10.1 and a package is built for R 2.11.0, users get the embarrassing notice about a disconnect, and cannot do much about it. If it is tedious to roll out a new package on time, perhaps there's a possibility for automation here. Thanks. I

Hi all, I am using a discrete Hidden Markov Model with discrete observations in order to detect a sequence of integers. I am using the "hmm.discnp" package. I am using the following code: signature <- c(-89, -98, -90, -84, -77, -75, -64, -60, -58, -55, -56, -57, -57, -63, -77, -81, -82, -91, -85, -89, -93) quant <- length(-110:-6) # Initialize and train the hmm with the

Hi, I need to draw hmm's in R.I also need to include a small plot near each state of the hmm.How can I do this? -- Rajesh.J [[alternative HTML version deleted]]

Hi all, I'm extracting the name of the term in a regression model that dropterm specifies as the least significant one, and I'm assigning this name to an object. However, when I use update(), it ignores this object. Is there a way I can make it not ignore it? A reproducible example is below: > lm(x1~1+y1*y2+y3+y4,data=anscombe)->my.lm >

Dear R-users; Previously I posted a question about the problem of rank deficiency in summary.manova. As somebody suggested, I'm attaching a small part of the data set. #*************************************************** "test" <- structure(.Data = list(structure(.Data = c(rep(1,3),rep(2,18),rep(3,10)), levels = c("1", "2", "3"), class =

I wish to rearrange the matrix, df, such that all there are not repeated x values. Particularly, for each value of x that is reated, the corresponded y value should fall under the appropriate column. For example, the x value 3 appears 4 times under the different columns of y, i.e. y1,y2,y3,y4. The output should be such that for the lone value of 3 selected for x, the corresponding row entries

Hi R-experts, suppose I have a list with containing data frame elements: [[1]] (Intercept) y1 y2 y3 y4 -6.64 0.761 0.383 0.775 0.163 [[2]] (Intercept) y2 y3 -3.858 0.854 0.834 Now I want to put them into ONE dataframe like this: (Intercept) y1

Hi, I would like to create a plot of y1,y2,y3,y4 against x for several subjects such that y1 and y2 are plotted against x in one panel and y3 and y4 against x in another panel. Thus if there are 3 subjects I should end up with 6 panels. Is there a simple way of doing so (i.e. without calling xyplot() several times, and then padding the results together)?? Regards S?ren

RNG in R and Splus 3.4 Prof. Ripley asked the details of the example. We were doing parametric bootstrap, so it is similar to simulation. Anyway here is the details. We start with a sample of 19 positive numbers. We know the sample is from truncated exp(0.3)...only the truncation point, theta, is unknown. In other words, the sample can be generated from something like x1 <- rexp(100,

RNG in R and Splus 3.4 Prof. Ripley asked the details of the example. We were doing parametric bootstrap, so it is similar to simulation. Anyway here is the details. We start with a sample of 19 positive numbers. We know the sample is from truncated exp(0.3)...only the truncation point, theta, is unknown. In other words, the sample can be generated from something like x1 <- rexp(100,

Hi Yves lavaan looks like a very nice package. From the tutorial introduction I see you create path diagrams for some of the models you describe. How did you do this? I don't see a function for this in the package. I know there is a path.diagram function in the sem package that uses dot to draw the diagram, but I've always found the layouts from dot somewhat strange for path diagrams

Dear members, I want to create 8 graphs and write it into one page using mfrow=c(4,2). How to make all graphs (including the titles, legends, line types) to be scale down (resized proportionally). As an illustration, below is the code: pdf("testmfrow.pdf") par(mfrow=c(4,2)) x<-seq(1:10) y1<-rnorm(10) y2<-rnorm(10,mean=2,sd=1) y3<-rnorm(10,mean=3,sd=1)

Dear R- and Omega-list-members, I am trying to make statistical inference about the following design: A dependent variable y has been measured multiple times, i.e. 4 times (y1,y2, y3, y4), unfortunately suffering from some successive dropouts (i.e. the sample sizes varies for y1, y2, y3, and y4). For every y, two other variables (covariates) were also measured: x & z, and both do presumably

Hello all! I want to estimate parameters in a MIMIC model. I have one latent variable (ksi), four reflexive indicators (y1, y2, y3 and y4) and four formative indicators (x1, x2, x3, x4). Is there a way to do it in R? I know there is the SEM library, but it seems not to be possible to specify formative indicators, that is, observed exogenous variables which causes the latent variable. Thanks,

Hi   I have four groups   y1=c(1.214,1.180,1.199) y2=c(1.614,1.710,1.867,1.479) y3=c(1.361,1.270,1.375,1.299) y4=c(1.459,1.335) Is there a function that can give me the length for each, like the made up example below?   >function(length(y1:y2) [1] 3 4 4 2 [[alternative HTML version deleted]]

Full_Name: Benjamin Tyner Version: 2.6.0 (43063) OS: WinXP Submission from: (NULL) (171.161.224.10) At the end of the examples for isoreg, there is cat("R^2 =", formatC(sum(residuals(ir4)^2) / (9*var(y4)), digits=2),"\n") I think this should be cat("R^2 =", formatC(1 - sum(residuals(ir4)^2) / (9*var(y4)), digits=2),"\n" Thanks Ben platform