search for: tbar

...roblem and solution below posted on r-help could have been a bit slicker if %*% worked with multidimensional arrays multiplying them so that if the first arg is a multidimensional array it is mulitplied along the last dimension (and first dimension for the second arg). Then one could have written: Tbar <- tarray %*% t(wt) / rep(wti, each = 9) which is a bit nicer than what had to be done, see below, given that %*% only works with matrices. I suggest that %*% be so extended to multidimensional arrays. Note that this is upwardly compatible and all existing cases would continue to work unchang...

...a faster construct. I tried things like "kronecker" and "outer" combined with apply, but couldn't get it to work. Here is a sample code: ########################## n <- 120 sigerr <- 5 covmat <- diag(c(8,6,3.5)) mu <- c(105,12,10) mcsamp <- 10000 Tbar <- array(0, dim=c(3,3,n)) # theta is a mcsamp x 3 matrix theta <- mvrnorm(mcsamp, mu = mu, Sigma = covmat) wt <- matrix(runif(n*mcsamp),n,mcsamp) wti <- apply(wt,1,sum) tarray <- array(apply(theta,1,function(x)outer(x,x)),dim=c(3,3,mcsamp)) for (i in 1:n) { for (k in...

...,m_j]>med[m_j]) med_ge[m_i,m_j]=0 else med_ge[m_i,m_j]=1 } y=c(1,1,1,1,1,1,0,0,0,0) n<-c(sum(y == 1),sum(y==0)) touse3 <- y==1 T1<- apply(med_ge[touse3,], 2, mean) T0<- apply(med_ge[!touse3,], 2, mean) T=rbind(T1,T0) Tbar=colMeans(T) Tdiff=T-Tbar cov_rank=cov(med_ge) inv_cov_rank=ginv(cov_rank) LN=0 for(m_i in 1:length(n)) { LN <- LN+((Tdiff[m_i,]%*%inv_cov_rank)%*%t(Tdiff)[,m_i])*n[m_i] } return(LN) }} func_LN(data) Now, I want to try this function on subgroups of data. S...

Hi, I have the following numbers as a result of bootstrap: Bootstrap Statistics : original bias std. error t1* 0.2700797 0.02168322 0.05843803 Intervals : Level Percentile 95% ( 0.2048, 0.4256 ) Can somebodly please show me how to plot a bar with confidence interval? I've searched the archive, but i'm a novice, so it didn't help

...+ break; + case 3: + tname = "ISR"; + break; + case 4: + tname = "DeviceCfg"; + break; + default: + tname = "<unknown>"; + break; + } + + printf("VirtIO: %s\n", tname); + + if (verbose < 2) + return; + + printf("\t\tBAR=%d offset=%08x size=%08x\n", + get_conf_byte(d, where + 4), + get_conf_long(d, where + 8), + get_conf_long(d, where + 12)); + + if (type != 2 || length < 20) + return; + + printf("\t\tmultiplier=%08x\n", + get_conf_long(d, where+16)); +} diff --git a/ls-caps.c b/ls-caps...

...+ break; + case 3: + tname = "ISR"; + break; + case 4: + tname = "DeviceCfg"; + break; + default: + tname = "<unknown>"; + break; + } + + printf("VirtIO: %s\n", tname); + + if (verbose < 2) + return; + + printf("\t\tBAR=%d offset=%08x size=%08x\n", + get_conf_byte(d, where + 4), + get_conf_long(d, where + 8), + get_conf_long(d, where + 12)); + + if (type != 2 || length < 20) + return; + + printf("\t\tmultiplier=%08x\n", + get_conf_long(d, where+16)); +} diff --git a/ls-caps.c b/ls-caps...

...e 3: + tname = "ISR"; + break; + case 4: + tname = "DeviceCfg"; + break; + default: + tname = "<unknown>"; + break; + } + + printf("VirtIO: %s\n", tname); + + if (verbose < 2) + return; + + printf("\t\tBAR=%d offset=%08x size=%08x\n", + get_conf_byte(d, where + 4), + get_conf_long(d, where + 8), + get_conf_long(d, where + 12)); + + if (type != 2 || length < 20) + return; + + printf("\t\tmultiplier=%08x\n", + get_conf_long(d, where+16)); +} + +void +show_vendor_caps(struct...

...e 3: + tname = "ISR"; + break; + case 4: + tname = "DeviceCfg"; + break; + default: + tname = "<unknown>"; + break; + } + + printf("VirtIO: %s\n", tname); + + if (verbose < 2) + return; + + printf("\t\tBAR=%d offset=%08x size=%08x\n", + get_conf_byte(d, where + 4), + get_conf_long(d, where + 8), + get_conf_long(d, where + 12)); + + if (type != 2 || length < 20) + return; + + printf("\t\tmultiplier=%08x\n", + get_conf_long(d, where+16)); +} + +void +show_vendor_caps(struct...

Hi, > BTW: is there a tool (or pciutils patch) which can decode the virtio > capabilities? Searched for a patch today, and all google found me was this mail asking for one :-o So I went ahead and coded one up. Attached. While hacking it up I've noticed spec doesn't match reality. The "Virtio Structure PCI Capabilities" section here ...

Hi, > BTW: is there a tool (or pciutils patch) which can decode the virtio > capabilities? Searched for a patch today, and all google found me was this mail asking for one :-o So I went ahead and coded one up. Attached. While hacking it up I've noticed spec doesn't match reality. The "Virtio Structure PCI Capabilities" section here ...