search for: nstate

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 w...

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-t...

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

...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 +++ switcher_new.c 2006-12-28 10:16:40.000000000 +0100 @@ -799,6 +799,13 @@ XA_ATOM, 32, PropModeReplace, (unsigned char *) state, nState); } + + Atom type; + type = XInternAtom(dpy, "_NET_WM_WINDOW_TYPE_SPLASH", 0); + XChangeProperty(dpy, ss->popupWindow, + XInternAtom (dpy, "_NET_WM_WINDOW_TYPE", 0), + XA_ATOM, 32, PropModeReplace, + unsigned cha...

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

...ver doing anything that queries that tapdisk2 pid in /proc/ will fully hang. Doing an "strace -f ps auxf &" (Backgrounding so I can keep my console), it ends here and I can find the pid causing it: open("/proc/11327/status", O_RDONLY) = 6 read(6, "Name:\ttapdisk2\nState:\tD (disk sl"..., 1023) = 706 close(6) = 0 open("/proc/11327/cmdline", O_RDONLY) = 6 read(6, Trying to do almost anything against /proc/11327/ results in a hang, but I can see the FD''s ok: ukxen2 ~ # cd /proc/11327 ukxen2 11327 # ls -al &a...

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

...) ); list_character_set( file, out_char_set ); _IO_putc ('\n', file); } int *epsclosure( t, ns_addr, accset, nacc_addr, hv_addr ) int *t, *ns_addr, accset[], *nacc_addr, *hv_addr; { register int stkpos, ns, tsp; int numstates = *ns_addr, nacc, hashval, transsym, nfaccnum; int stkend, nstate; static int did_stk_init = 0, *stk; if ( ! did_stk_init ) { stk = (int *) allocate_array( current_max_dfa_size, sizeof( int ) ); did_stk_init = 1; } nacc = stkend = hashval = 0; for ( nstate = 1; nstate <= numstates; ++nstate ) { ns = t[nstate]; if ( ! (trans1[ns] < 0) )...

...for (i = 0; i < freq->nr; i++) t += scnprintf(str + t, PRIV_STR_SIZE - t, "mhz[%u]=%u\n", i, freq->mhz[i]); return t; } - break; case OID_TYPE_MLME:{ struct obj_mlme *mlme = r->ptr; return scnprintf(str, PRIV_STR_SIZE, "id=0x%X\nstate=0x%X\ncode=0x%X\n", mlme->id, mlme->state, mlme->code); } - break; case OID_TYPE_MLMEEX:{ struct obj_mlmeex *mlme = r->ptr; return scnprintf(str, PRIV_STR_SIZE, "id=0x%X\nstate=0x%X\n" "code=0x%X\nsize=0x%X\n", mlme->id,...

...for (i = 0; i < freq->nr; i++) t += scnprintf(str + t, PRIV_STR_SIZE - t, "mhz[%u]=%u\n", i, freq->mhz[i]); return t; } - break; case OID_TYPE_MLME:{ struct obj_mlme *mlme = r->ptr; return scnprintf(str, PRIV_STR_SIZE, "id=0x%X\nstate=0x%X\ncode=0x%X\n", mlme->id, mlme->state, mlme->code); } - break; case OID_TYPE_MLMEEX:{ struct obj_mlmeex *mlme = r->ptr; return scnprintf(str, PRIV_STR_SIZE, "id=0x%X\nstate=0x%X\n" "code=0x%X\nsize=0x%X\n", mlme->id,...