similar to: classification using zero-inflated negative binomial mixture model

>From copyVector in duplicate.c : void copyVector(SEXP s, SEXP t) { int i, ns, nt; nt = LENGTH(t); ns = LENGTH(s); switch (TYPEOF(s)) { ... case INTSXP: for (i = 0; i < ns; i++) INTEGER(s)[i] = INTEGER(t)[i % nt]; break; ... could that be replaced with : case INTSXP: for (i=0; i<ns/nt; i++) memcpy((char *)DATAPTR(s)+i*nt*sizeof(int),

Hello when I set qualify = yes on trunk I can't do outgoing call. Incoming is always working. [Feb 15 23:01:41] WARNING[12909][C-00000012]: app_dial.c:2525 dial_exec_full: Unable to create channel of type 'SIP' (cause 20 - Subscriber absent) but my linphone is registered all the time. when set qualify = no outgoing call is working (but i have problems when WAN IP is changed after

Hello, I am having a problem with the zero-inflated negative binomial (package pscl). I have 6 sites with plant populations, and I am trying to model the number of seeds produced as a function of their size and their site. There are a lot of zero's because many of my plants get eaten before flowering, thereby producing 0 seeds, and that varies by site. Because of that and because the

Hi all, I have to install 25 IP Phone in some building. I want just basic IP Phones like: Cisco-Linksys SPA922 u$s 146 Grandstream GXP-2000 u$s 105 Snom 300 u$s 119 The most valuables parameters for me are (in importance order from high to low): - Stability (device don't hang in any way) - Voice quality using G729 - Provisioning So what device do

I found a form of the if statement that works but it is very long. I'm attempting to check the value of of two different variables to see if they evaluate to either of two different values, which will result in a division by 0 in the final equation. This first if statement works for me _______ if ((x0.trial01 == 0.0)||(x0.trial01 == npt01)||(x1.trial01 == 0.0)||(x1.trial01 == npt01))

>> Amsterdam!! > I had my laptop and suitcase stolen in Amsterdam the one time I went > there, after hearing someone talk about how safe a city it was over > dinner. Most importantly, also stolen was my (apparently irreplacable) > copyleft shirt (yellow/gold with large blue backwards (C) symbol on front > and GPL preamble on back) which no amount of effort has managed to

Last thursday ( Feb 2nd ) Dag's place was broken into a significant amount of hardware was taken, along with other valuables. You can see the details posted by Dag at http://lists.freshrpms.net/pipermail/freshrpms-list/2006-February/013868.html We all know the excellent work he has done for the entire linux community over the years, and specifically his continued support for the CentOS

Dear All, Attached files is related to the changes made to add the Schedmodel for a AMD bulldozer target, Please note that , the model is incomplete but has some of the valuables features implemented. Request to the group or someone from AMD for the comments on the implementation. Thanks ~umesh On Wed, Nov 13, 2013 at 8:14 PM, Umesh Kalappa <umesh.kalappa0 at gmail.com>wrote: >

Hi, I am facing the problem that I want to plot a histogram chart set freq to true and overlay with normal or weibull or exponential distribution curve. The sample code is shown as below: >samp<-c(-8.2262,-8.2262,-8.2262,-8.20209,-8.09294,-8.07321,-8.07321, -8.07321,-8.07175,-8.04948,-8.04948,-8.04948,-8.03848,-8.03848,

I would like to use samp as a part of a filename that I can change. My source files are .csv files with date as the file name, and I would like to be able to type in the date (later perhaps automate this using list.files) and then read the csv and write the pdf automatically. I have tried different combinations with "" and () around samp, but I keep getting the error "object

"lpt1" through "lpt9" and "com1" through "com9" are in fact "special" file names under various versions of Microsoft Windows. However, "R Extensions" says only that "lpt1" through "lpt4" and "com1" through "com3" cannot be used as filenames, and "R CMD check" checks for just these

my reproducible example test<-structure(list(site = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L), .Label = c("A", "B", "C", "D", "E"), class = "factor"),

Hello guys, I would like to ask for help to understand what is going on in "func2". My plan is to generalize "func1", so that are expected same results in "func2" as in "func1". Executing "func1" returns... 0.25 with absolute error < 8.4e-05 But for "func2" I get... Error in dpois(1, 0.1, 23.3065168689948, 0.000429064542600244,

Hi! On my notebook I got a strage ext3-related problem: root at jura98 /usr/portage/app-misc/obexftp # ls -l total 4 dr-Sr-sr-t 2 8242 15720 4096 Dec 28 1993 metadata.xml [how this got created is beyond my knowledge, possibly because of some random memory corruption a bad kernel patch caused a few weeks ago] root at jura98 /usr/portage/app-misc/obexftp # rm metadata.xml/ rm: cannot remove

Dear List, We are using the mix() function in the mixdist package to fit mixture models to some of our data. The package provides a function to compare the fits of nested models using an ANOVA function. However, we were wondering whether there are methods that can be used to compare models that differ in the number of distributions fitted and/or the shape of the distributions. Is there a way,

I'm pretty new in R... try to learn ASAP for a project, I wonder if anyone can provide some advice on which package to use when fit Mixture regression model with both categorical(more than two levels) and continuous variables. I tried to use Flexmix... and transposed the multi-level variables as several 0-1 dummy variable and tried to put in the package and find the fit. But it always tells

Hi, I'm trying to fit a hidden Markov model (mixture of Gaussians) to stock prices (Open, High, Low, Close) data. Here's my code: #First download the data with package TTR library(TTR) n <- 200 data <- getYahooData("GOOG", strftime(as.POSIXlt(Sys.time() - n*24*3600),format="%Y%m%d")) keep <-

Hi, Does anyone know how to compute the quantile of a mixture of four bivariate normal distriutions? Many thanks! Hannah [[alternative HTML version deleted]]

I couldn't find a direct answer in CRAN to this question, so I'm asking with some trepidation. I have a multivariate dataset (data.frame) with columns that can be expressed as a set of mixed normals (at least I think) and need to impute values that have constraints (truncated mixture of normals where the values cannot be below zero). If there isn't a package that can do this, is there

hi i'm looking for a gibbs sampling algorithm for R for the case of mixture of K normals, and in particular for the case of bivariate normals. i'd be grateful if anyone could send its own R-routine, at least for the univariate case. thank you in advance matteo