similar to: Random walk

Hi everybody, I am trying to fit a GAM model without intercept using library mgcv. However, the result has nothing to do with the observed data. In fact the predicted points are far from the predicted points obtained from the model with intercept. For example: #First I generate some simulated data: library(mgcv) x<-seq(0,10,length=100) y<-x^2+rnorm(100) #then I fit a gam model with

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

Dear all, I generated a covariance matrix and I would like to generate a 1D plot of the data that auto-correlate. any suggestions? Thanks, Marcelo -- Marcelo Andrade de Lima UNIFESP - Universidade Federal de São Paulo Departamento de Bioquímica Disciplina de Biologia Molecular Rua Três de Maio 100, 4 andar - Vila Clementino, 04044-020 Lab +55 11 55764438 R.1188 Cell +55 11 92725274

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

I want to use var.test function of ctest library. I was able to generate the output from R on the input given by rnorm. But when I store the same data set provided by rnorm into a file, and then read a file into a dataframe, then on using the dataframe as a parameter to var.test function - it gives an error "not enough x observations". I'll explain the above mentioned problem of

Hello everybody! With a HTB qdisc: # tc qdisc add dev eth0 root handle 99: htb How can I change the "final qdisc" of the default walk-throught class? # tc qdisc add dev eth0 parent 99:0 handle 88: sfq results in: "RTNETLINK answers: Invalid argument" CBQ automatically creates a xx:0 root class where I can attach final qdiscs... - If there is a root class in htb, how is it

Dear expeRts, recently I asked for a nice way to re-program a problem without using control constructs such as "for" or "sapply(1:length(x), ...". Is there a way to program a random walk with a reflecting boundary without resorting to such constructs? A working solution is ranwalk <- function(length, bound) { k <- cumsum(sample(c(-1, 1), length, replace=TRUE))

Following is my code, can some one help on the error at the bottom? > mh<-function(iterations,alpha,beta){ + data<-read.table("epidemic.txt",header = TRUE) + attach(data, warn.conflicts = F) + k<-97 + d <- (sqrt((x-x[k])^2 + (y-y[k])^2)) + p <- 1-exp(-alpha*d^(-beta)) + p.alpha<-1 - exp(-3*d^(-beta)) + p.beta <- 1 - exp(alpha*d^(-2)) +

This patch adds useful APIs for walking over the private data area. It is a prerequisite for the new event API stuff (specifically for the language bindings for that). Rich. -- Richard Jones, Virtualization Group, Red Hat http://people.redhat.com/~rjones virt-p2v converts physical machines to virtual machines. Boot with a live CD or over the network (PXE) and turn machines into Xen guests.

I would like to either talk to over the phone a * admin a few question on the deploy and a few question on the setup of * at the very least if some doesn't mind or a walk thru of there site if there in the Houston area. I am looking at cutting over to * for some things but I have some setup question that could easily be answer in about 30-60 minutes. Please reply to lists efastfunding com

From: Mircea C?rjaliu <mcirjaliu at bitdefender.com> This change is needed because the introspection tool can write-protect guest page tables, exec-protect heap/stack pages and let KVM emulate the instruction that caused these violations. Signed-off-by: Mircea C?rjaliu <mcirjaliu at bitdefender.com> Signed-off-by: Adalbert Laz?r <alazar at bitdefender.com> ---

There are some corrupted registry files that have invalid hbin cells but are still readable. This patch makes the following changes: * hivex_open - do not abort with complete failure if we run across a block with invalid size (unless it's the root block). Instead just log the event, and move on. This will allow open hives that have apparent invalid blocks but the ones of potential

Is anyone mirroring Johnny's write up on this? As fate would have it, I've finally got some time to sit down and configure a spare box to be my new mail server and Johnny's site has been defaced again. :-( Sigh....I even checked to make sure it was up this morning, drove to the office so I could do everything from the terminal and my handy dandy roadmap is gone.... Any pointers

Hi all, I'm having issues coding a random walk on a fully connected, undirected graph G with probability transition matrix P = (p_ij). Here is what I have so far ... for(i in 1:n){ for(m in 1:M){ x <- as.vector(matrix(rep(0,N+1),nc=N+1)) x[1] <- i for(k in 2:N+1){ y <- as.vector(matrix(rep(0,n),nc=n)) r <- runif(1) c <- c(0,cumsum(P[x[k-1],])) for(j in

Hi, I'm trying to use the kalman filter to estimate the variable drift of a random walk, given that I have a vector of time series data. Anyone have any thoughts on how to do this in R? Thanks, Alex [[alternative HTML version deleted]]