Displaying 20 results from an estimated 22 matches for "pi1".
Did you mean:
i1
2012 Apr 11
2
What is a better way to deal with lag/difference and loops in time series using R?
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...
2013 Apr 09
0
[R-SIG-Finance] EM algorithm with R manually implemented?
...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...
2009 Mar 27
1
General help for a function I'm attempting to write
...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...
2015 Feb 09
0
About your driver NUTDRV_ATCL_USB(8) for install
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 <...
2015 Feb 09
2
About your driver NUTDRV_ATCL_USB(8) for install
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
2013 May 12
1
Multinomial-Dirichlet using R
...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
[[...
2015 Feb 09
0
RV: About your driver NUTDRV_ATCL_USB(8) for install
...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]...
2014 Jul 24
3
RV: About your driver NUTDRV_ATCL_USB(8) for install
...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...
2015 Nov 07
2
mkfs.ext2 succeeds despite nbd write errors?
...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...
2014 Jul 24
0
RV: About your driver NUTDRV_ATCL_USB(8) for install
...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...
2014 Jul 17
3
About your driver NUTDRV_ATCL_USB(8) for install
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
2015 Nov 07
0
Re: mkfs.ext2 succeeds despite nbd write errors?
...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...
2010 Apr 05
0
[LLVMdev] [GSoC] llvm-test
...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
2003 Nov 01
0
html glitches with help?
...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...
2012 Aug 27
0
How can I find the principal components and run regression/forecasting using dynlm
...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...
2008 Apr 09
0
Endogenous variables in ordinal logistic (or probit) regression
...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...
2000 Feb 04
2
constrained optimiser doesn't obey constraints (PR#411)
...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....
2008 Apr 24
0
logit newbie question...
...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...
2002 Aug 01
2
mdct.h - PI1_8, PI2_8 etc.
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...
2013 Mar 29
1
Create values based on a table of conditions
...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",...