similar to: Line of Organic Correlation (help with writing function)

mac osx 10.5.4 R 2.7.1 I have fit a model d<-lm(y~x) with an R^2 of 0.963 but when I issue the command abline(d) the line is below where it ought to be. Looks like the right slope, but not the right intercept. thanks -- Stephen Sefick Research Scientist Southeastern Natural Sciences Academy Let's not spend our time and resources thinking about things that are so little or so large

Dear friends i am trying to fit an Ornstein-Uhlenbeck process by MAXIMUM LIKELYHOOD method. i found these formulas on http://www.sitmo.com/article/calibrating-the-ornstein-uhlenbeck-model/ this is the mean-reverting process http://r.789695.n4.nabble.com/file/n4637271/process.txt process.txt and this is the script that i am using....... ouFit.ML=function(spread) { n=length(spread)

I don't even know if this is the right place to ask this question. I would like to download USGS stream gauging data for a couple of gauges on a daily basis- save the files to .csv files and append the nest days time series to this. Is there a way to do this automatically in R? thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or

Dear list, I just tried to do the same thing, and did not find anything on a weighted qqplot. My weights are actually counts (positive integers). Here is a modification of qqplot, following Duncan Murdoch's suggestion. Any feedback would be welcome! Thanks, Jean-Christophe weighted.qqplot <- function (x, y, plot.it = TRUE, xlab = deparse(substitute(x)), ylab = deparse(substitute(y)),

Hello, I've written a small function that's supposed to save me some time, and it's ending up killing it- the intention is to iteratively subset a dataset fram on framevec, fit a model (either lm or nls depending on type) and return the r2 or AIC from the model, respectively. Although as far as I can tell in my code the plots are dependent on the fit of the model to the data and the

I have a function in R^2, say f <- function(x,y) { ...skipped } I want to plot this function using contour, persp. wireframe, etc. I know that the function has a global minimum at (x0, y0) The naive approach is to evaluate the function on the outer product of two arrays, like this: sx <- c(seq(-3, x0, len = 100), seq(x0, 3, len = 100)[-1]) sy <- c(seq(-3, y0, len = 100), seq(y0, 3,

I don't know if there is a way to use the scale relation free argument in ggplot2 like in lattice. I have a feeling that there is not, but I would like to make a plea for this feature. It would be nice to be able to plot Total Inorganic Nitrogen Total Phosphorus and the ratio of the two- the numbers on the axis are not related, but the previous two are surely related to the last (this ratio

Dear list members, I am working with a multilevel growth curve, that in its simplest form goes like follows: Yit = Ai + Bi t + eit (the error term is assumed to follow an AR(1) autorregressive process) One major topic in my research is the convergence in the values of Y over time. Thus, I am interested in the relationship between the random effects for the intercept and the slope, and I

Hello ppl, we are using two servers, with dovecot 2.0.5, Sx and Sy, with: cat dovecot-ldap-ext: ... pass_attrs = uid=user,userPassword=password,=y=proxy,mailHost=host ... When I access server Sx via IMAP I see all my folders. When I access server Sy via IMAP it do not proxy to Sx. I want to see my IMAP folders independently of the server. any help? -- :) cumprimentos

Hello, I have a compiler in LLVM 2.9 for the KCPM3 processor. I'd like to create a subtarget for the new cpu version called KCPSM6. Besides a couple of new instructions which are not important at the moment, the KCPSM6 cpu has different instruction opcodes. Semantically the instructions are the same, hence I'd like to keep all the lowering and pattern matching stuff unmodified For

Hi, I am trying to do a backend to a very simple microcontroller. I have some questions. 1) I have instruction which do "r1 <- r1 op r2", from what I have ssen I must declare them like: let isTwoAddress = 1 in def ADD : FopRR< 0b01010, (outs CPURegs:$sX), (ins CPURegs:$isX, CPURegs:$sY), "ADD $sX, $sY"), [(set

On Apr 27, 2008, at 7:37 AM, Cédric Venet wrote: > Hi, > > I am trying to do a backend to a very simple microcontroller. I have > some > questions. Ok. > 1) I have instruction which do "r1 <- r1 op r2", from what I have > ssen I > must declare them like: > > let isTwoAddress = 1 in > def ADD : FopRR< 0b01010, > (outs

Hi, I am doing a project on the simulation of glucose metabolism based on a pharmacokinetic modeling in which we have 4 differential equations. I did this in R by using the odesolve package. It works very well, but I have two questions: Here is the odemodel function _________________________________________________ Ogtt.Odemodel <- function(t, y, p) { absx <- c(-60, -45, -30,

Hi R Users, How to use the function polygon() together with the package scatterplot3d? I am trying to color below of the curves defined for the function density(). I tried to use the site: R GRAPH GALLERY as tutorial. I tried to adapt the example of this page: [figure]: http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=30 [code]:

I have a time series from 1933-2005 of precipitation at Fayetteville NC. I get the following error messages when I plot the zoo series. Any help would be appreciated. If you need the data I can dput it or send the csv. I didn't include it here because I didn't want to clog up anybodies email account. I know that this is not reproducible, and I will send along the file if needed.

I have two sets of interval data.Below are my two dataset. In these dataset, there is a missing values in each of the data. I want to find the non-overlapping interval values. Here is my code: mysetdiff=function(x,y){ m=length(x) n=length(y) bx = logical(m) by = logical(n) for(i in 1:m){ for(j in 1:n){ if(x[i]<=y[j]){ bx[i] = T by[j] = T }

I'm trying to do a simple (?) analysis of a 1D spatial data set, allowing for spatial autocorrelation. (Actually, I'm comparing expected vs. observed for a spatial model of a 1D spatial data set.) I'm using models like gls(obs~exp,correlation=corExp(form=~pos),data=data) or gls(obs~exp,correlation=corLin(form=~pos),data=data) This form is supposed to fit a linear model of

rendercheck -f a8r8g8b8 passes as much as on a GK208, and xv appears to work. Very lightly tested. Instead of sticking coordinates into pushbufs, the vertex shader is modified to read them from a constbuf, indexed by vertex id. This approach could be used for all nvc0 generations, but I didn't want to rock the boat. Signed-off-by: Ilia Mirkin <imirkin at alum.mit.edu> ---

Dear lists, i want to find the non-overlapping interval values with this code: mysetdiff=function(x,y){ m=length(x) n=length(y) bx = logical(m) by = logical(n) for(i in 1:m){ for(j in 1:n){ if(x[i]<=y[j+1]){ bx[i] = T by[j] = T NA = NA } } } sx = x[!bx] sy = y[!by] s=c(sx,sy) return(s) } Below is my

Looks reasonable, some minor comments below. On 10/16/2016 02:06 AM, Ilia Mirkin wrote: > rendercheck -f a8r8g8b8 passes as much as on a GK208, and xv appears to > work. Very lightly tested. > > Instead of sticking coordinates into pushbufs, the vertex shader is > modified to read them from a constbuf, indexed by vertex id. This > approach could be used for all nvc0 generations,