search for: yh

thanks to all for the 2d scatter plot. i have one more. how do i plot 'hist(y1, col="red") and hist(y2,col="blue") in the same window? thanks again.

...ight to compute the size of the box which is plotted with rect: z <- rnorm(10) # horizontal text works plot(rnorm(10)) x1 <- 5 y1 <- 0 label <- "Label" cha <- paste(" ", label, " ", sep = "") xh <- strwidth(cha, cex = par("cex")) yh <- strheight(cha, cex = par("cex")) * 5/3 rect(x1 - xh/2, y1 - yh/2, x1 + xh/2, y1 + yh/2, col = "white") text(x1, y1, cha, cex = par("cex")) Now I would like to do the same with rotated text. I found two solutions: - modify the par("usr") - compute rat...

Jean Marc & list, Following up on this- there seems to be an issue with this calculation in mdf.c (Smooth echo energy estimate over time) Pey += Eh*Yh; Pyy += Yh*Yh; it goes too large for floating point, which seems to be caused by extremely large values in st->Rf[j], st->Yf[j] and smaller values in st->Eh[j], st->Yh[j]- things were relatively quiet environmentally speaking. Anyhow, leak_estimate would then go to undefine...

...antities for one sample (y* is the response predicted from the model derived; i used one to test my R traduction compared to the R pls module ) Calculate the covariance matrix C=XXT and c*=Xx* for prediction y is centered and become y1 y*1=0 For h =1,2,3...hmax s=Cyh center s working scalar for prediction sample s*=c*Tyh orthogonalize s to previous t: for g=1,...(h-1), s=s-(tgTs/tgTtg)tg orthogonalize s* to previous t*: for g=1,...(h-1), s*=s*-(tgTs/tgTtg)t*g t*h=s* th=s th2=tTt betah=(tTyh)/t...

Did a bit of tuning in SVN. Could you try that? What frame size, filter length, sampling rate are you using? All the process you describe is actually the tip of the iceberg I think. There is something making the adaptation diverge, which causes Eh and Yh to be large and so on. The part I changed is probably the one causing divergence. If it doesn't work, reducing the coefficients even more should. Jean-Marc Le mercredi 16 novembre 2005 ? 10:25 -0800, Tom Harper a ?crit : > Jean Marc & list, > > Following up on this- there seems...

...frow=c(2,2)) > plot(fm2) Cheers, Berwin *** plot.lm.R.orig Fri Aug 15 22:13:08 2003 --- plot.lm.R Wed Mar 3 18:05:34 2004 *************** *** 16,32 **** show <- rep(FALSE, 4) show[which] <- TRUE r <- residuals(x) - n <- length(r) yh <- predict(x) # != fitted() for glm if (any(show[2:4])) { s <- if(inherits(x, "rlm")) x$s else sqrt(deviance(x)/df.residual(x)) hii <- lm.influence(x, do.coef=FALSE)$hat } if (any(show[2:3])) { ylab23 <- if(isGlm) "Std. devi...

...anual para incorporar funciones u propiedades adicionales a programas desarrollados?? Gracias Atte Diego Anexo la programación de interes ".fecov" <- function(x, n.ahead) { n.par<-sapply(x$varresult, function(x) summary(x)$df[2]) sigma.u <- crossprod(resid(x))/n.par Sigma.yh <- array(NA, dim = c(x$K, x$K, n.ahead)) Sigma.yh[, , 1] <- sigma.u Phi <- Phi(x, nstep = n.ahead) if (n.ahead > 1) { for (i in 2:n.ahead) { temp <- matrix(0, nrow = x$K, ncol = x$K) for (j in 2:i) { temp <- temp + Phi[, , j] %*% sigma.u %*% t(Phi[, ,...

I append below a suggested update for p.adjust(). 1. A new method "yh" for control of FDR is included which is valid for any dependency structure. Reference is Benjamini, Y., and Yekutieli, D. (2001). The control of the false discovery rate in multiple testing under dependency. Annals of Statistics 29, 1165-1188. 2. I've re-named the "fdr" meth...

Hi, Is it possible to program ssh tunneling in C or C++ rather than to use shell command? Are there any APIs for ssh tunneling available? I was told by friends that the openssh is working on it, please point me which version of openssh support API ssh tunneling. Thank you. Kind regards, J. yh

An embedded and charset-unspecified text was scrubbed... Name: cpumask:update-irq_desc-to-use-cpumask_var_t.patch Url: http://lists.linux-foundation.org/pipermail/virtualization/attachments/20090107/8ebe5843/attachment.txt

An embedded and charset-unspecified text was scrubbed... Name: cpumask:update-irq_desc-to-use-cpumask_var_t.patch Url: http://lists.linux-foundation.org/pipermail/virtualization/attachments/20090107/8ebe5843/attachment.txt

An embedded and charset-unspecified text was scrubbed... Name: cpumask:update-irq_desc-to-use-cpumask_var_t.patch Url: http://lists.linux-foundation.org/pipermail/virtualization/attachments/20090107/8ebe5843/attachment.txt

hi, i am trying to generate nice barplots with std-errors. do i really have to generate the std-errors myself by the segments() command ? thanks for help, jan -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

...stop("Use only with 'lm' objects") if (!is.numeric(which) || any(which < 1) || any(which > 5)) stop("`which' must be in 1:5") isGlm <- inherits(x, "glm") show <- rep(FALSE, 5) show[which] <- TRUE r <- residuals(x) yh <- predict(x) w <- weights(x) if (!is.null(w)) { wind <- w != 0 r <- r[wind] yh <- yh[wind] w <- w[wind] labels.id <- labels.id[wind] } n <- length(r) if (any(show[2:5])) { s <- if (inherits(x, "rlm")) x$s...

...--------------- # Plain Vanilla NADARAYA-WATSON estimator (or Local Constant regression, e.g. deg=0) # with gaussian kernel & fixed bandwidth mkern<-function(y,x,h){ Mx <- matrix(x,nrow=length(y),ncol=length(y),byrow=TRUE) Mxh <- (1/h)*dnorm((x-Mx)/h) Myxh<- (1/h)*y*dnorm((x-Mx)/h) yh <- rowMeans(Myxh)/rowMeans(Mxh) return(yh) } # Generating the design Y=m(x)+e n <- 10 h <- 0.5 x <- rnorm(n) y <- x + rnorm(n,mean=0,sd=0.5) # This is what I really want! mhat <- mkern(y,x,h) library(locfit) yhl.raw <- locfit(y~x,alpha=c(0,h),kern="gauss",ev="...

Dear list, I have a series of data points which I want to approximate with exactly two linear functions. I would like to choose the intervals so that the total deviation from my fitted lines is minimal. How do I best do this? Thanks! Kamila The information transmitted in this electronic communication...{{dropped}}

...me 2835ms rtt min/avg/max/mdev = 44.475/44.791/44.976/0.332 ms [igor at localhost ~]$ curl -I www.iquaid.org curl: (6) Couldn't resolve host 'www.iquaid.org' working URL [igor at localhost ~]$ ping www.google.com PING www.google.com (74.125.137.99) 56(84) bytes of data. 64 bytes from yh-in-f99.1e100.net (74.125.137.99): icmp_seq=1 ttl=40 time=24.0 ms 64 bytes from yh-in-f99.1e100.net (74.125.137.99): icmp_seq=2 ttl=40 time=24.4 ms 64 bytes from yh-in-f99.1e100.net (74.125.137.99): icmp_seq=3 ttl=40 time=24.1 ms C64 bytes from yh-in-f99.1e100.net (74.125.137.99): icmp_seq=4 ttl=40...

_mulx_u64 only exists when the target is x86_64. That's still not very portable. I'm not opposed to removing the bmi2 check, but gcc also has the same check so it doesn't improve portability much. ~Craig On Sat, Dec 29, 2018 at 4:44 PM Arthur O'Dwyer via llvm-dev < llvm-dev at lists.llvm.org> wrote: > Hi Pawel, > > There is the _mulx_u64 intrinsic, but it

> On Tue, Jul 15, 2008 at 10:33 AM, Suresh Siddha > <suresh.b.siddha at intel.com> wrote: > > On Sun, Jul 13, 2008 at 10:19:35PM -0700, Yinghai Lu wrote: > >> > >> fix for pv. > >> > >> Signed-off-by: Yinghai Lu <yhlu.kernel at gmail.com> > >> > >> --- > >> arch/x86/kernel/paravirt.c | 5 ---- > >> arch/x86/kernel/vmi_32.c | 51 ++++++++++++++++++++++++++++++++++++++++++--- > >> arch/x86/xen/enlighten.c | 19 +++++++--------- > >> include/asm...

HW: Lenovo T520 laptop with MS WIN 7 OR dedicated PC (DELL GX 280) SW: CentOS 6.5 (64-bit) installed using VMware Player on MS WIN 7 host CASE A) -- BRIDGED networking selected CASE B) -- NAT networking OR CentOS 6.5 i(32-bit) nstalled on real HW (DELL OptiPlex GX 280) ___issue____ at hand: The default Firefox ESR browser (17.0.10) does not find many web pages.