similar to: Help with loess

How to add trendline (i.e. straight line passing through maximum points) in graph. I have worked on the data given below. Please tell me how to add trendline in the graph. The script is as follows =================================== start ==================================================== # The data is as follows data <- c( 0.01, 0.02, 0.04, 0.13, 0.17 , 0.19 , 0.21 , 0.27 , 0.27 ,

Hi >i have this data > X [1] 5.79 1579.52 2323.70 68.85 426.07 110.29 108.29 1067.60 17.05 22.66 [11] 21.02 175.88 139.07 144.12 20.46 43.40 194.90 47.30 7.74 0.40 [21] 82.85 9.88 89.29 215.10 1.75 0.79 15.93 3.91 0.27 0.69 [31] 100.58 27.80 13.95 53.24 0.96 4.15 0.19 0.78 8.01 31.75 [41] 7.35 6.50

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I

Hello, I would be very grateful if someone could give me a hand with my split plot design problems. So here is my design : I am studying the crossed-effects of water (wet/dry) and mowing (mowed/not-mowed = nm) on plant height (PH) within 2 types of plant communities (Xerobromion and Mesobromion) : - Within each type of communities, I have localised 4 blocks - In each block, I have defined

On Mon, Dec 16, 2013 at 04:16:29PM -0800, Michael Dalton wrote: > Commit 2613af0ed18a ("virtio_net: migrate mergeable rx buffers to page frag > allocators") changed the mergeable receive buffer size from PAGE_SIZE to > MTU-size, introducing a single-stream regression for benchmarks with large > average packet size. There is no single optimal buffer size for all >

On Mon, Dec 16, 2013 at 04:16:29PM -0800, Michael Dalton wrote: > Commit 2613af0ed18a ("virtio_net: migrate mergeable rx buffers to page frag > allocators") changed the mergeable receive buffer size from PAGE_SIZE to > MTU-size, introducing a single-stream regression for benchmarks with large > average packet size. There is no single optimal buffer size for all >

I tried toy problems and there doesn't seem to be a basic problem between lm() and ts objects: X = data.frame(x=c(1,2,7,9), y=c(7,2,3,1)) lm(y ~ x, X) X <- lapply(X, function(x) ts(x, frequency=12, start=c(1994,7))) lm(y ~ x, X) and this works fine - whether you do an lm() before or after making ts objects, it's okay. But I have a situation where things aren't okay.

is there a way to use the actual index value for plotting zoo objects this is the way that the index is set up and a sample range of what I would like to plot 01/01/06 00:00:00 - 01/01/06 23:45:00 { library(zoo) # chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) }} x <- structure(c(15.57, 15.5,

This one really has me stumped. I'm running 866MHz, 256MB RAM, 80GB Disk (7200RPM) on Linux-mandrake 8.2 icecast 1.3.12 & ices 0.2.3 & lame 3.91 all mp3s encoded with lame 3.91 or higher and --r3mix (VBR) Sorry for the long email but I'm hoping to get someone out there that's seen this before and give me some advice... I'm continually getting errors like this when running

Hi, I was trying to read in the following data set by using: gas <- read.table("gas.tab", header = T) when I got the error message: Error in "row.names<-.data.frame"(*tmp*, value = row.names) : duplicate row.names are not allowed Any help are appreciated. Thanks, Ko-Kang Wang

Hi there are some good tips in "The R Inferno" http://www.burns-stat.com/documents/books/the-r-inferno/ or connect C++ to R with Rcpp http://dirk.eddelbuettel.com/code/rcpp.html or byte code compiler (library(compiler)) or library(data.table) but do you have an idea to fasten standard R source code, with the following Rprof output self.time self.pct total.time

I am hereby forwarding the data & method use to calculate the Kolmogorov-Smirnov goodness of fit test made manually by me in R launguage which deffers with the actual inbuilt formula as shown below. Further I have plot the graph in R. In that graph how to add trendline (i.e. straight line passing through maximum points in plot) to a Plot. R script is as follows please run this script to see

Dear all:Request your indulgence. The econophysics gurus do this stuff all the time: all their PDFs are smooth, with neat log x axis. 1. The kernel density estimate (KDE) function returns the empirical probability density at 2^n points (min: 512). The big question is how do I scale back the x-values (say, density$x) to x-values in terms of the original dataset? 2. To give you a concrete idea, i

# Just when I thought I had the basic stuff mastered.... # This has been quite perplexing, thanks for any help ## Here's the example: db1=data.frame( olditems=c('soup','','','','nuts'), prices=c(4.45, 3.25, 4.42, 2.25, 3.98)) db2=data.frame( newitems=c('stew','crackers','tofu','goatsmilk','peanuts'))

I am trying to interpolate missing values using spline and am running into some strange problems. first, this works just fine: x<- c(1:10, rep(NA, 3), 14:20) y<- c(rnorm(10), rep(NA,3), rnorm(7)) s<- spline(x,y, n=length(x)) plot(x,y) lines(s, col="blue") but look at what happens with my real data (sorry for the long vectors here): x<- c(-9.816, -9.801, -9.801, -9.795,

Hi All , Could you please tell using which library ,Dickey-Fuller Test can be run? Thanks a lot __________________________________________________ [[alternative HTML version deleted]]

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

Hello! Some time back, I added the Speex protocol to my version of VOCAL (www.vovida.org, VOIP tool). Recently, I also added PESQ (automated voice quality testing algorithm) to my tool and have been running some tests on a clean network. The source file is a woman reading some phrases meant to test various aspects of codecs... Speex has a respectable result of 3.67 Some other codecs I've

I just got Linux Mint 10, and for some reason, Toribash 3.91 just doesn't start at all. I tried it on Winebottler on OSX, and that version functioned fine. Install goes along normally. When I click on the exe that starts the game, though, it simply doesn't start. Doing it from terminal gives me this: err:menubuilder:WinMain unknown option -a err:menubuilder:WinMain unknown option

On 2016-05-25 19:13, Kelly Lesperance wrote: > Hdparm didn?t get far: > > [root at r1k1 ~] # hdparm -tT /dev/sda > > /dev/sda: > Timing cached reads: Alarm clock > [root at r1k1 ~] # Hi Kelly, Try running 'iostat -xdmc 1'. Look for a single drive that has substantially greater await than ~10msec. If all the drives except one are taking 6-8msec, but one is very