Displaying 7 results from an estimated 7 matches for "bbar".
Did you mean:
bar
2011 Jul 19
2
loops and simulation
...t; ll <- function(y, k, c){
> n*log(c*k)-(k+1)*( sum(log(1+c*Y))+sum( R*log(1+c*Y)+1/k ) )
> }
>
> # Compute the log-likelihood for the initial values, and ignoring the
> missing data mechanism
> lltm1 <- ll(Yobs, theta.t, beta.t)
>
> repeat{
> # E-step
>
> Bbar<- function(theta.t,beta.t){
> sum(R*(1+beta.t*(theta.t+1)*Yobs)/(beta.t*(theta.t+1)*(1+beta.t*Yobs)))
> }
>
> Abar<- function(theta.t,beta.t){
> sum( R* ( log(1+beta.t*Yobs)+ 1/theta.t ))
> }
>
> theta.t1<- n/( sum(log(1+beta.t*Yobs)) + Abar(theta.t,beta.t))
>...
2005 May 19
1
R 2.1.0 RH Linux Built from Source Segmentation Fault
...hy)
> accepts=rep(0,R/keep)
> parameters=list(R=1,keep=1,s=scale)
> for (j in 1:R)
+ {
+ itime=proc.time()[3]
+ for (i in 1:nphy){
+ # for each physician i draw a beta[i] using MH algorithm
+ #
+ Data=list(m=nalt,X=X[rowx1[i]:rowx2[i],],y=Y[rowy1[i]:rowy2[i],])
+ bbar=t(theta)%*%z[i,]
+ prior=list(A=delta,betabar=bbar)
+ a=rmnlRwMetrop1(Data,prior,Mcmc=parameters,beta[i,])
+ beta[i,]=a$betadraw
+ accept[i]=a$acceptr
+ }
+ mregout=rmultireg(beta,z,thetabar,A1,nu1,V)
+ delta=mregout$Sigma
+ theta=mregout$B
+ if (j%%keep==0){
+ thetas[...
2011 Feb 22
5
"aggregate" in R
Hi, R users,
I'm wondering how I can aggregate data in R with different functions for
different columns. For example:
x<-rep(1:5,3)
y<-cbind(x,a=1:15,b=21:35)
y<-data.frame(y)
I want to aggregate "a" and "b" in y by "x". With "a", I want to use
function "mean"; with "b", I want to use function "sum". I tried:
2004 Mar 10
5
unexpected softupdate inconsistency
Hi there,
System is FreeBSD 4.9 Stable
Had a problem during a portupgrade that when portupgrade did some
cleaning a directory could not be deleted because it were not empty.
However, it was empty. rm -rf didn't help either.
Because I didn't know what to do I decided to boot.
At system startup fsck was started which produced the message: unexpected
softupdate inconsistency and also that
2011 Jun 02
0
[PATCH] pci: Use pr_<level> and pr_fmt
...t;%s: io_start: %x\n",
+ __func__, hpc_ptr->u.isa_ctlr.io_start);
+ debug("%s: io_end: %x\n",
+ __func__, hpc_ptr->u.isa_ctlr.io_end);
+ debug("%s: irq: %x\n", __func__, hpc_ptr->irq);
break;
case 2:
case 4:
- debug ("%s - wpegbbar: %lx\n", __func__, hpc_ptr->u.wpeg_ctlr.wpegbbar);
- debug ("%s - i2c_addr: %x\n", __func__, hpc_ptr->u.wpeg_ctlr.i2c_addr);
- debug ("%s - irq: %x\n", __func__, hpc_ptr->irq);
+ debug("%s: wpegbbar: %lx\n",
+ __func__, hpc_ptr->u.wpeg_ct...
2011 Jun 02
0
[PATCH] pci: Use pr_<level> and pr_fmt
...t;%s: io_start: %x\n",
+ __func__, hpc_ptr->u.isa_ctlr.io_start);
+ debug("%s: io_end: %x\n",
+ __func__, hpc_ptr->u.isa_ctlr.io_end);
+ debug("%s: irq: %x\n", __func__, hpc_ptr->irq);
break;
case 2:
case 4:
- debug ("%s - wpegbbar: %lx\n", __func__, hpc_ptr->u.wpeg_ctlr.wpegbbar);
- debug ("%s - i2c_addr: %x\n", __func__, hpc_ptr->u.wpeg_ctlr.i2c_addr);
- debug ("%s - irq: %x\n", __func__, hpc_ptr->irq);
+ debug("%s: wpegbbar: %lx\n",
+ __func__, hpc_ptr->u.wpeg_ct...
2011 Jun 02
0
[PATCH] pci: Use pr_<level> and pr_fmt
...t;%s: io_start: %x\n",
+ __func__, hpc_ptr->u.isa_ctlr.io_start);
+ debug("%s: io_end: %x\n",
+ __func__, hpc_ptr->u.isa_ctlr.io_end);
+ debug("%s: irq: %x\n", __func__, hpc_ptr->irq);
break;
case 2:
case 4:
- debug ("%s - wpegbbar: %lx\n", __func__, hpc_ptr->u.wpeg_ctlr.wpegbbar);
- debug ("%s - i2c_addr: %x\n", __func__, hpc_ptr->u.wpeg_ctlr.i2c_addr);
- debug ("%s - irq: %x\n", __func__, hpc_ptr->irq);
+ debug("%s: wpegbbar: %lx\n",
+ __func__, hpc_ptr->u.wpeg_ct...