similar to: Calculate mean slope of a function

MailMarshal (an automated content monitoring gateway) has not delivered the following message: Message: B0000c89bd.00000001.mml From: samba@samba.org To: lianaj@balgrammar.vic.edu.au Subject: Question This is due to automatic rules that have determined that the intended recipient is not authorized to receive messages with certain potentially dangerous filetypes attached.

Hi there, I have two series of data. plotting the density function of both gives me an idea about the difference of the data. But I would like to quantify the difference I see. a <- rnorm(100) b <- rnorm(100) da <- density(a) db <- density(b) The problem is that da$x and db$x are different and so I have difficulties to compare them... Is there any way to force the density

I recently tried to install the latest version of spatstat, from CRAN, using the install.packages() function. It proceeded to install version 1.17-5 of spatstat, although the current version is 1.18-4. Checking the CRAN mirror that I used (New Zealand) via Firefox, I found that version 1.18-4 is indeed present on CRAN. I was able to download and install version 1.18-4 ``manually''

Hi all, I am writing a function (fun.R), but I dont know how to code the function so that the Help Text will be shown up when one types ?fun (of course, after he loads it up). Anyone has any advice for me how to do that? Thanks, D.

I am trying to construct a function to which I pass an expression as an argument. >From that expression I want to create a somewhat more complicated expression and then differentiate it using D() or deriv(). To give a simple example, I'd like to be able to do something like e1 <- expression(x^2 + y^2) e2 <- some.magic(e1,"n") or maybe e2 <- some.magic(e1,expression(n))

Hello everybody, I don't understand the warning I get when executing the following code: x <- rnorm(100) h <- hist(x, plot= FALSE, freq = FALSE) If I don't use the plot-argument, I don't get a warning... Any explanations for me? Antje

Dear all, I'm trying to interpolate a dataset to give it twice as many values (I'm giving the dataset a finer resolution by interpolating from 1 degree to 0.5 degrees) to match that of a corresponding dataset. I have the data in both a data frame format (longitude column header values along the top with latitude row header values down the side) or column format (in the format latitude,

Hi guys... I have problem with this excersise... Consider the pressure data frame again. (a) Plot pressure against temperature, and use the following command to pass a curve through these data: > curve((0.168 + 0.007*x)?(20/3), from=0, to=400, add=TRUE) (b) Now, apply the power transformation y3/20 to the pressure data values. Plot these transformed values against temperature. Is a linear

I am working on a project where we start with start with 2 long, equal-length vectors, massage them in various ways, and end up with a function mapping one interval to another. I'll call that function "f1." The last step in R is to generate f1 as the value of the approxfun function. I would like to put f1 into a package, but without having the package redo the creation of

Dear all, I have run into a problem when running some code implemented in the Bioconductor panp-package (applied to my own expression data), whereby gene expression values of known true negative probesets (x) are interpolated onto present/absent p-values (y) between 0 and 1 using the *approxfun - function*{stats}; when I have used R version 2.8, everything had worked fine, however, after updating

Hi r-devels, we are seeing a new problem with our packages RobAStRDA (just new on CRAN, thanks to Uwe and Kurt!) and RobExtremes (to be submitted). It must be something recent with the way you internally treat/store byte-code compiled functions, as we have no problems with R-3.1.3, but do see an "Error in fct(x) : byte code version mismatch" with R-devel SVNrev r70532. Background:

Readers, According to the help '?approxfun', the function can be used to obtain the interpolated values. The following test was tried: > testinterpolation<-read.csv('test.csv',header=FALSE) > testinterpolation V1 V2 1 10 2 2 20 NA 3 30 5 4 40 7 5 50 NA 6 60 NA 7 70 2 8 80 6 9 90 9 10 100 NA >

I am (desperately) trying to get used to using a Mac here at my new location. (Why *anyone* would ever use anything other than Linux, except under duress as I am, totally escapes me, but that's another story.) Fortunately much of the Mac OSX is actually Unix, so a civilized person can manage to carry on ... But there are some things. (Like this <expletive deleted> mailer ... But

Dear R Users, I have a couple of proto objects like: wedge <- proto(expr={ start.year <- 2008 end.year <- 2050 }) star.rating <- wedge$proto( star = c(4, 5, 8, 10), gain = c(0, .3, .5, .7), cost = c(0, 2100, 4000, 7500), star.rating <- function(., year) 6.0, setup = function(.){ .$cost.for.star <- approxfun(.$star,

Thanks, Luke, for having a look to it. Sure, I can give you some reproducible example -- even in two degrees of completeness ;-): see below. Thanks again, Peter %----------------------------------- (I) first example %----------------------------------- Just to reproduce the error, on r-devel, try: install.packages("RobAStRDA") require(RobAStRDA) getFromNamespace(".RMXE", ns

how can I import/export data from ms excel? also is there a way to import/export data of the hdf format? is there a way to import/export hdf format for ms excel as well? thank you very much in advance. /js ###################################################################### This e-mail message has been scanned and cleared by MailMarshal http://www.marshalsoftware.com

Hi there, I am wondering if it is possible to do an xyplot with two y-axes. I'd like to print two parameters in a time series but they both have different scales. Which parameter in xyplot can I add to achieve this result? Thanks a lot for any help. Antje [[alternative HTML version deleted]]

Hello, I'm not very experienced with lattice and I was wondering whether I get get some hints from you how to create a pure heatmap (using levelplot), without any axis, title, legend, margin at all... I just want to see the coloured squares, nothing else. Any suggestions? Antje

Hello, I tried to fit a poisson distribution but looking at the function fitdistr() it does not optimize lambda but simply estimates the mean of the data and returns it as lambda. I'm a bit confused because I was expecting an optimization of this parameter to gain a good fit... If I would use mle() of stats4 package or mle2() of bbmle package, I would have to write the function by myself

I have several hundred Excel 2007 data files in a folder. I would like to read every file in a single given folder using a loop. I have searched the FAQ, the forum archives here, other or older R boards and the R Import / Export documentation, and have asked some very knowledgeable R users without learning of a solution. I hope someone here can help. I understand that the most common