search for: nstates

Hi, I have a question about running multiple in regressions in R and then storing the coefficients. I have a large dataset with several variables, one of which is a state variable, coded 1-50 for each state. I'd like to run a regression of 28 select variables on the remaining 27 variables of the dataset (there are 55 variables total), and specific for each state, ie run a regression of

...TRUE) emission_matrix=list(meanval_matrix,sigmaval_matrix,coeffval_matrix) powerdf is a column from my dataset which look like something = 19.0, 18.0, 24.0...............it has some 30 thousand rows I tried using the code below and got an error mod= depmix( response = power~1, data = powerdf, nstates=3, instart=intial_prob_matrix, trstart=transition_matrix, respstart=emission_matrix) Error in makeResponseModels(response = response, data = data, nstates = nstates, : 'respstart' has incorrect length, it should be 6 I cannot change my emission matrix, it wo...

Hi all. I am testing 0.25.4 to upgrade from 0.23.2. In testing, I found two strange messages on 0.25.4 which I have never found on 0.24.3. So could you tell me what these messages mean and how to fix them. [Environment] Puppetmasterd -version : 0.25.4 -OS : Solaris10 10/09 Puppetd -version : same as Puppetmasterd -OS : same as Puppetmasterd [Strange messages] 1.When Puppetd

Hi all, I am using ACS micro data (PUMS) with one of the columns as a factor for the place of birth (POBPF). I would like to create a column (POBR) containing a rank corresponding to the place of the observation in the POBPF rankings. For example, if a person is from Oregon, Oregon is the most popular Place of Birth, so each person who came from Oregon should have a 1 in their POBR.

Nice! This looks good to me but probably Chris or someone else should sign off on it. There were two minor warnings on Darwin: -- In file included from /Volumes/Data/Users/ddunbar/llvm/lib/Support/regexec.c:81: /Volumes/Data/Users/ddunbar/llvm/lib/Support/regengine.inc: In function 'sbackref': /Volumes/Data/Users/ddunbar/llvm/lib/Support/regengine.inc:665: warning: control reaches end of

...TTR) n <- 200 data <- getYahooData("GOOG", strftime(as.POSIXlt(Sys.time() - n*24*3600),format="%Y%m%d")) keep <- c("Open","High","Low","Close") data2 <- data[,keep] #Now fit the Hidden Markov Model library(RHmm); HMMFit(data2, nStates=4, nMixt=3,dis="MIXTURE") The HMMFit function returns usually an error (every once in a while it actually works, but if I run the command again it'll give the error) Rgui.exe-Application Error The instruction at "0x00530057" referenced memory at "0x00530057". The...

I am applying a hidden markov model on joint multivariate gaussian distribution for 2 vectors. I am using the depmixS4 package in R. Specifically, I am using the following code: mod<-depmix(list(response = mom ~ mkt + p0 + p1, mkt~1), data = regvar, nstates = 2, family = list(gaussian(), gaussian()),instart = delta, trstart=Pi) It seems that depmixS4 doesnt output the covariance estimates of the parameters. What is the best package to obtain them? -- View this message in context: http://r.789695.n4.nabble.com/covariance-matrix-of-model-parameters-tp...

Dear helpers, is there any possible that transition (in depmixS4) is in scale of two variable, e.g transition=~scale(x1,x2)? If it can be, how transition of two variable (covariate time) can be worked in depmixS4-hidden markov model for time series. Many thanks, nglthu -- View this message in context: http://r.789695.n4.nabble.com/depmixS4-transition-tp4572726p4572726.html Sent from the R help

Hi, I made a simple patch which makes switcher's window a splash window (_NET_WM_WINDOW_TYPE_SPLASH). This way I can easily add effects for it. The window is a splash window, so why not use this type instead of Unknown? Beryl's switcher window is a Splash window too. Is there a reason why Compiz uses Uknown for it? Patch: --- switcher.c 2006-12-28 10:14:52.000000000 +0100 +++

Hey all, I am trying to work under a State Space form, but I didn't get the help exactly. Have anyone eles used this functions? I was used to work with S-PLUS, but I have some codes I need to adpt. Thanks alot, Bernardo [[alternative HTML version deleted]]

Hello. I am trying to fit an asymptotic relationship (nonlinear) to some ecological data, and am having problems. I am interested in the upper bound on the data (i.e. if there is an upper limit to 'y' across a range of 'x'). As such, I am using the nonlinear quantile regression package (nlrq) to fit a michaelis mention type model. The errors I get (which are dependant on

Have an interesting one here, originally found on xen 4.1.0 but just upgraded to xen 4.1.1 and it''s still here. Creating a VM with a tap2 device pointed at /dev/mapper/something, when that device doesn''t exist, causes the tapdisk2 process to go into D mode and also manages to take out any process that queries it. For example, I have /dev/mapper/nathanxenuk1 as a valid disk

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >

I have a new problem: Register RBP is used in a function foo. (I am not allocating RBP to any virtual register, the instances of RBP in function foo are in the machine code when my register allocator starts.) Function foo calls function bar. Register RBP is not saved across the call, though it is live after the call. Function bar includes a virtual register. The code that I'm using to

On 11/30/2012 6:36 PM, Lang Hames wrote: > > > RBP is used as the frame pointer on x86 (hence its automatic > appearance in your code), and shouldn't be allocated to any vreg in > function bar. Loading/saving RBP should be managed by the stack frame > setup/teardown code. > If it doesn't already, your allocator should filter out reserved > registers (See

From: Tom Rix <trix at redhat.com> This is a upcoming change to clean up a new warning treewide. I am wondering if the change could be one mega patch (see below) or normal patch per file about 100 patches or somewhere half way by collecting early acks. clang has a number of useful, new warnings see

From: Tom Rix <trix at redhat.com> This is a upcoming change to clean up a new warning treewide. I am wondering if the change could be one mega patch (see below) or normal patch per file about 100 patches or somewhere half way by collecting early acks. clang has a number of useful, new warnings see