search for: pi1

Hello, I am writing codes for time series computation but encountering some problems Given the quarterly data from 1983Q1 to 1984Q2 PI1<-ts(c(2.747365190,2.791594762, -0.009953715, -0.015059485, -1.190061246, -0.553031799, 0.686874720, 0.953911035), start=c(1983,1), frequency=4) > PI1 Qtr1 Qtr2 Qtr3 Qtr4 1983 2.747365190 2.791594762 -0.009953715 -0.015059485 1984 -1.190061246 -0.55303...

...k, because I get such odd values especially for sigma, that the > likelihood of some observations is zero and therefore the log is -Inf. > > Where is my mistake? > > My code: > > # EM algorithm manually > # dat is the data > > > # initial values > pi1<-0.5 > pi2<-0.5 > mu1<--0.01 > mu2<-0.01 > sigma1<-0.01 > sigma2<-0.02 > loglik[1]<-0 > loglik[2]<-sum(pi1*(log(pi1)+ > log(dnorm(dat,mu1,sigma1))))+sum(pi2*(log(pi2)+log(dnorm(dat,mu2,sigma2)))) > > > > tau...

...matically calculate 'JostD' for all combinations of columns, perhaps returning a matrix of distances. Thanks for any help or suggestions. Cheers Colin Function: JostD <- function(DF, X, Y) { Ni1 <- DF[,X] Ni2 <- DF[,Y] N1 <- sum(Ni1) N2 <- sum(Ni2) pi1 <-Ni1/N1 pi2 <-Ni2/N2 pisqr <- ((pi1+pi2)/2)^2 H1 <- 1 - sum(pi1^2) H2 <- 1 - sum(pi2^2) Ha <- 0.5*(H1+H2) Da <- 1/(1-Ha) Dg <- 1/sum(pisqr) Db <- Dg/Da D <- -2*((1/Db) - 1) D } Sample data: e<-c(0,0,0,4,27...

Dear Chales, Thanks for your quick reply. I have Raspberry with last Kernel [Linux pi1 3.18.6+ #753 PREEMPT Sun Feb 8 14:47:22 GMT 2015 armv6l GNU/Linux] If <upsdrvctl>, <upsd> and <upsmon> are stopped <nutdrv_atcl_usb> Works well without any error as "Device or resource busy". If I run <upsc> command never receive any value from my ups &lt...

On Feb 9, 2015, at 6:48 AM, PERE (yahoo.es) <pererivers at yahoo.es> wrote: > 0.028188 Can't set USB configuration: could not set config 1: Device or resource busy > 0.029246 detached kernel driver from USB device... > 0.030559 USB device [0001:0000] opened > 0.866216 read: (8 bytes) => 03 00 00 00 00 00 00 00 > 0.866413 reply[0] = 0x03

...nd there are 63 sites. If I were to analyze the data for multiple outcomes per site, the number of events per category given the category response rates can be assumed to be independently distributed. They can be modeled by a multinomial distribution with parameters n=60 and category response rates pi1, . . . , piC for site i. The individual variation in category response rates can be modeled by a Dirichlet distribution. Just wondering if I am thinking through this correctly. If so, could someone share some thoughts on how this could be done in R? Really Appreciate any help. Thanks Anamika [[...

...ould not install the nut in my raspberry. All installation good, but with my ups " Ovislink Chrome 1000VA" and "nutdrv_atcl_usb" driver, it inform me the following error details. I can startup nut but nutmon don't receive anything only (lost conection). My setup: root at pi1:/nut# lsusb Bus 001 Device 002: ID 0424:9512 Standard Microsystems Corp. Bus 001 Device 001: ID 1d6b:0002 Linux Foundation 2.0 root hub Bus 001 Device 003: ID 0424:ec00 Standard Microsystems Corp. Bus 001 Device 006: ID 0001:0000 Fry's Electronics <------------- my UPS ups.conf [ups]...

...eply, but I can't install nut_2.7.2 in my raspberry. > > 1) I create a [nut] directory and use wget to download three files: > [nut_2.7.2-1.dsc] [nut_2.7.2.orig.tar.gz] [nut_2.7.2-1.debian.tar.xz] 2)I > unpack: > dpkg-source -x nut_2.7.2-1.dsc > 3) I use command > root at pi1:/nut/nut-2.7.2# debuild -us -uc . > .. > And I have the following error code: > --------------------------------------------- > root at pi1:/nut/nut-2.7.2# debuild -us -uc -d dpkg-buildpackage -rfakeroot -d > -us -uc > dpkg-buildpackage: aviso: est? usando una orden para converti...

...e/cell/nbds/default/chunks/00000000000000030728' nbdkit: chunks[1]: error: Unable to fallocate '/home/cell/nbds/default/chunks/00000000000000031232' Indeed, there is definitely a problem with fallocate, as some of the chunks are the correct size (256k), and some are zero length: cell@pi1$ pwd /home/cell/nbds/default/chunks cell@pi1$ ls -l | tail -rw------- 1 cell cell 262144 Nov 7 06:01 00000000000000032256 -rw------- 1 cell cell 262144 Nov 7 06:01 00000000000000032257 -rw------- 1 cell cell 262144 Nov 7 06:01 00000000000000032258 -rw------- 1 cell cell 262144 Nov 7 06:01 00000...

...harles, Thanks for your quick reply, but I can't install nut_2.7.2 in my raspberry. 1) I create a [nut] directory and use wget to download three files: [nut_2.7.2-1.dsc] [nut_2.7.2.orig.tar.gz] [nut_2.7.2-1.debian.tar.xz] 2)I unpack: dpkg-source -x nut_2.7.2-1.dsc 3) I use command root at pi1:/nut/nut-2.7.2# debuild -us -uc . .. And I have the following error code: --------------------------------------------- root at pi1:/nut/nut-2.7.2# debuild -us -uc -d dpkg-buildpackage -rfakeroot -d -us -uc dpkg-buildpackage: aviso: est? usando una orden para convertirse en administrador (?gain-ro...

On Jul 13, 2014, at 9:46 AM, pere at riusnebot.com wrote: > Dear Charles Lepple, > Sorry for my email intrusion. for next time: http://lists.alioth.debian.org/cgi-bin/mailman/listinfo/nut-upsuser > I have an UPS called Ovislink 1000VA and it use a UPSMART2000R (with ?ATCL FOR UPS?). > I?m a beginner Linux user and I ?will install your new driver <nutdrv_atcl_usb > to work

...0030728' > nbdkit: chunks[1]: error: Unable to fallocate > '/home/cell/nbds/default/chunks/00000000000000031232' > > > Indeed, there is definitely a problem with fallocate, as some of the > chunks are the correct size (256k), and some are zero length: > > cell@pi1$ pwd > /home/cell/nbds/default/chunks > cell@pi1$ ls -l | tail > -rw------- 1 cell cell 262144 Nov 7 06:01 00000000000000032256 > -rw------- 1 cell cell 262144 Nov 7 06:01 00000000000000032257 > -rw------- 1 cell cell 262144 Nov 7 06:01 00000000000000032258 > -rw------- 1 cell...

...wer than I thought or find some missing tests I find interesting as I advance. All in all, I do learn fast and I'm really motivated as I was going to work with LLVM during summer anyway. Thank you for your attention if you read until here. [0] http://oredakore.es/index.php?p=torii [1] http://pi1.informatik.uni-mannheim.de/filepool/presentations/01-ulix-hardware.pdf [2] https://pi1.informatik.uni-mannheim.de/index.php?pagecontent=site%252FTeaching.menu%252FCourse%2520Archives.menu%252FFSS_09.menu%252FTeleseminar.page&show=true

...similar, except the typeface for \code was identical to that of the rest. Sample excerpt from the .Rd file: \value{returns a list: \item{singles}{data frame, one row per simulation, with the following columns: \tabular{ll}{ coefficients \tab one column per variable\cr \code{pi1} \tab conditional probability of sampling cases \cr Similar behavior in \arguments section. It is the pi1 above, for example, that is aligned poorly. Alignment of the outer level \item's is better (in fact, the opposite problem, their baseline is above the baseline of the first line of the s...

...n earlier; h is the forecasting horizon; since dynlm can't work for zoo object, the variables need to be transformed to ts objects. Xlast<-as.matrix(X[(t-h):(t-h-i+1)]) zlast<-t(as.matrix(z[(t-h):(t-h-j+1)])) # Make it a column matrix beta<-as.matrix(model$coefficient) pi1fore[t]<-pi1[t-h]+c(1, Xlast, zIlast)%*%beta # pi1fore is defined earlier . I just use the inner product as a vehicle for forecasting. I don't know if there's a way that the codes can be written neatly? return(pi1fore) Then I attempt to program similarly for the multi-factor cas...

...tical model, there are unmeasured scale variables y1* and y2* that are determined by y1* = b0 + b1 * y2 + b2 * X1 + b3 * X2 + e1 and y2* = c0 + c1 * y1 + c2 * X1 + c3 * X2 + e2 y1* and y2* are not observed, we see only the categorical outputs Y1 and Y2 that correspond to Y1 = 0 if y1* < pi1 Y1 = 1 if pi1 <= y1* < pi2 Y1 = 2 if pi2 <= y1* < pi3 Y1 = 3 if pi3 <= y1* and similarly for Y2. Since e1 is "going into" y1*, and y1* "goes into" y2*, then there is the good chance that the error term e1 is correlated with y2*. Running p...

...se starting value it eventually reports convergence, for some others it goes forever. ---Example code--- revise.mix.binbb_function(x,loss,no.inform){ if (any(is.nan(x)) || any(x<0) || any(x>c(1,1,.5,1))){ print(paste("x=",x)) return(Inf) } eta_x[1] pi1_x[2] omega1_x[3] p0_x[4] log.term_0 for (i in 1:length(no.inform)){ if ( (loss[i]-1) >= 0) { first.prod_1 for (r in 0:(loss[i]-1)){ first.prod_first.prod*(pi1+r*omega1) };} else {first.prod_1;} if ( (no....

...19 753 No 288 347 Desire is my explanatory and Drugs Using is my response... How can I specify the null model? Namely logit(pi) = n where mle should give n to be logit((219+288)/(753+347)) Also, if I want a one factor model logit(pi1) = n logit(pi2) = n + alpha where mle should give me n to be logit(219 / 972) and alpha be logit(288/635) - logit(219/972) Thanks! -- View this message in context: http://www.nabble.com/logit-newbie-question...-tp16850195p16850195.html Sent from the R help mailing list archive at Nab...

In vorbis/lib/mdct.h the following are defined: for integer: #define TRIGBITS 14 #define cPI3_8 6270 #define cPI2_8 11585 #define cPI1_8 15137 #define FLOAT_CONV(x) ((int)((x)*(1<<TRIGBITS)+.5)) for floats: #define cPI3_8 .38268343236508977175F #define cPI2_8 .70710678118654752441F #define cPI1_8 .92387953251128675613F #define FLOAT_CONV(x) = x Could someone explain where these values come from? What's the significanc...

...ot;CO1", "CO10", "CO11", "CO12", "CO13", "CO14", "CO15", "CO16", "CO17", "CO2", "CO3", "CO4", "CO5", "CO6", "CO7", "CO8", "CO9", "PI1", "PI2", "PI3", "PI4", "PI5", "PI6", "PI7", "PI8", "PO1", "PO10", "PO11", "PO12", "PO13", "PO2", "PO3", "PO4", "PO5", "PO6",...