search for: dy2

...esults. Can anyone tell my what I am doing wrongly? Many thanks, Werner set.seed(123456) e1 <- rnorm(100) e2 <- rnorm(100) y1 <- ts(cumsum(e1)) y2 <- ts(0.6*y1 + e2) lr.reg <- lm(y2 ~ y1) error <- ts(residuals(lr.reg)) error.lagged <- error[-c(99, 100)] dy1 <- diff(y1) dy2 <- diff(y2) diff.dat <- data.frame(embed(cbind(dy1, dy2), 2)) colnames(diff.dat) <- c('dy1', 'dy2', 'dy1.1', 'dy2.1') ecm.reg <- lm(dy2 ~ error.lagged + dy1.1 + dy2.1, data=diff.dat) ecm.dynreg <- dynlm(d(y2) ~ L(error) + L(d(y1),1) + L(d...

Dear List, Wonder if you have some thoughts on the following question using lsoda in desolve: I have the following data and function: require(deSolve) times <- c(0:24) tin  <- 0.5 D <- 400 V    <- 26.3 k <-0.056 k12  <- 0.197118 k21  <- 0.022665 yini <- c(dy1 = 0,dy2 = 0)  events <- data.frame(var = "dy1",time = c(10,15),value = c(200,100),method = "add") pkmod <- function(t, y, p) {   if (t < tin) R <- (D/tin) else R <- 0   dy1 <- R - (p["k12"]+p["k"])* y[1] + p["k21"]*y[2]   dy2 <- p[&quot...

...um proved to be very helpful, I got another question... I'd like to fit data points on which I have an error, dx and dy, on each x and y. What would be the common procedure to fit this data by a linear model taking into account uncertainty on each point? Would weighting each point by 1/sqrt(dx2+dy2) (and taking dx and dy as relative errors) in a lm() fit do the job? I would like to propagate uncertainty of the points into the uncertainty of the fit, would that be the case? Thanks for all the help -- View this message in context: http://r.789695.n4.nabble.com/Fitting-with-error-on-data-tp271...

...0.00000 22.33107 334.44880 4402.99293 562.91532 I wondered if you could give me any idea where I am gong wrong. Many thanks Angela Dr Angela Moss Public Health Information Analyst Dudley PCT St. John's House Union Street Dudley DY2 8PP Tel: 01384 366091 Fax: 01384 366485 [[alternative HTML version deleted]]

I need to solve a system of non-linear fourth-order differential equations. Is there a command which solves this system? Thanks in advance.

...55.47692 89.46138 20 246.18142 62.57087 73.32186 21 186.47726 277.25464 386.75818 Many thanks for your help Dr Angela Moss Senior Public Health Intelligence Specialist Dudley PCT St. John's House Union Street Dudley DY2 8PP Tel: 01384 366493 Fax: 01384 366485 [[alternative HTML version deleted]]

...hould I look at hit rates and if so how? Can you present the outcome of the model graphically? I would appreciate any help you could provide. Many thanks Yours faithfully Angela Moss Dr Angela Moss Public Health Information Analyst Dudley PCT St. John's House Union Street Dudley DY2 8PP Tel: 01384 366091 Fax: 01384 366485 [[alternative HTML version deleted]]

...if (method == 0) { GPolyline(ns, xxx, yyy, USER, dd); @@ -871,7 +871,7 @@ if (gotLabel) { /* find which plot edge we are closest to */ - int closest; // 0 = index, 1 = index+range + int closest; /* 0 = index, 1 = index+range */ double dx1, dx2, dy1, dy2, dmin; dx1 = fmin2((xxx[index] - dd->gp.usr[0]), (dd->gp.usr[1] - xxx[index])); @@ -971,8 +971,8 @@ */ GConvert(&ux, &uy, USER, INCHES, dd); GConvert(&vx, &vy, USER, INCHES, dd); - // 1, 1 => sans serif, basic font - // 0, .5 => left, c...

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking