Timur == Timur Elzhov [EMAIL PROTECTED]
on Wed, 21 Jan 2004 14:54:58 +0300 writes:
Timur Dear R experts. 'D()' function recognizes some of
Timur the analitical functions, such as sin, cos, etc. But
Timur I'd like to take analytical derivatives from asin,
Timur atan etc.
--- In reply to Deepayan Sarkar: ---
Date:20.01.04 10:40 (-0600)
Thank you again for your very helpful and inspiring answer!
Some additional questions will follow below.
On Tuesday 20 January 2004 04:14, Wolfram Fischer wrote:
How one can add a text (e.g. the labels of an axis)
in a
PD == Peter Dalgaard [EMAIL PROTECTED]
on 21 Jan 2004 19:08:38 +0100 writes:
PD Martyn Plummer [EMAIL PROTECTED] writes:
Calculating the matrix exponential is harder than it
looks (I'm sure Peter knows this). In fact there is a
classic paper by Moler and Van Loan from the
Murad == Murad Nayal [EMAIL PROTECTED]
on Wed, 21 Jan 2004 15:19:28 -0500 writes:
Murad This might have been fixed in later versions (I am
Murad using R1.7.0),
yes, the bug has been fixed long ago,
from my ChangeLog (!), it was 2003-07-18.
Murad r-help archive contains
Hi,
The article Debugging Without (Too Many) Tears (Rnews 3/3) describes
the packages debug and mvbutils but it has no reference to where one can
find those packages and they are not in CRAN.
I've searched google but I can not find it :(
A reference will be welcome.
Thanks
EJ
Hi,
I read about packages debug and mvbutils in the R news but I cannot find
them on CRAN in the Package Sources. Does anyone know where to look?
Thanks,
Henning
--
Henning Rust
Potsdam Institute for Climate Impact Research
Dept. Integrated Systems Analysis
Tel.: #49/331/288-2596
Fax.:
Hi,
have got anybody experience with queing theory
and R or exist some material. Searching in archive i found nothing.
For first i search a function (..or have mysellf to write) which calculates
the frequency of parallel events. I have the starting time and duartion_time
of the events.
Many
Hello,
I'm fitting linear mixed models to gene-expression data from
microarrays, in a data set where 4608 genes are studied.
For a sample of 5 subjects and for each gene we observe the expression
level (Intensity) in four different tissues: N, Tp, Tx and M.
I want to test whether the expression
I read about packages debug and mvbutils in the R news but I cannot find
them on CRAN in the Package Sources. Does anyone know where to look?
Since it's the middle of the night in Tasmania I asked Mark yesterday
and he said they'd be on CRAN in the next day or two.
Hi,
I've got a quick question about file permissions and packages...
I'm creating my own package, and am having problems with its vignette
not being seen when I install it into R...
As I understand it, the permissions of the source tree should be as
follows:
o Directories - drwxrwxr--
o
Dear Sir/Madam,
I could succesfully normalise my microarray data using marrayNorm package.
However, i have not been able to get normalised red and green channel
intensities through R package. Is there a possibility to write a formula
to calculate back the red and green channel intensities
Hi, here's the calculations
M = log2(R/G) = [log rules] = log2(R)-log2(G) (1)
A = 1/2*log2(R*G) = [log rules] = (log2(R)+log2(G))/2 (2)
where log2(x) is the logarithm of x with base 2. If R,G 0, there is
a one-to-one relationship between (A,M) and (G,R) as follows
R =
H i, all!
First of all, I'd like to apologize for my poor English. It's for years
I don't use it.
This is a R-version of a function I wrote a long ago for my HP48 calculator.
It works with the binary expression of the power and just need to
duplicate the mem used by X.
Hope this helps.
hello
I am returning to some libraries that I had previously customised by amongst
other things adding additional functions. I had simply typed these new
functions into the file in library/R/thelibrary. These however do not seem
to be loaded now as they previously were under older versions of R.
Have you executed first
detach(package:library)
before loading the library again?
Eryk
*** REPLY SEPARATOR ***
On 1/22/2004 at 1:47 PM Stephen Henderson wrote:
hello
I am returning to some libraries that I had previously customised by
amongst
other things adding additional
If you want to access these utilities before M.Bravington post them on CRAN
they are available on his site:
(ftp://ftp.marine.csiro.au/software/bravington/)
Gérald Jean
Analyste-conseil (statistiques), Actuariat
télephone: (418) 835-4900 poste (7639)
télecopieur : (418)
I had quit q() R, reloaded and loaded ipred anew. Trying
detach(package:ipred) makes no difference. I added a few pointless lines to
one existing function and can see these when I call the function, passing no
parameters. thats why i was assuming it was an R thing.
is this just some mysterious
There are the concept of private and public functions in the new version.
There are a file in the root directory of the lib where public functions can be
entered.
I by myself do not know much more. There are more info in the R-ext.pdf (R extensions).
Probably you have to add your new function
ha ha ha ha..hmmm.
Well although I don't follow the developments and deprecations of R
religiously that doesn't mean I'm innumerate. Thanks for the advice
nonetheless-- its working now.
-Original Message-
From: Prof Brian Ripley
To: Stephen Henderson
Cc: ''[EMAIL PROTECTED]' '
Sent:
The ipred package has a namespace. Do follow the instruction in the R-exts
manual if you want functions available via a package.
Andy
From: Stephen Henderson
I had quit q() R, reloaded and loaded ipred anew. Trying
detach(package:ipred) makes no difference. I added a few
pointless lines
ipred has a NAMESPACE: you cannot do this to packages that have
namespaces without changing the NAMESPACE file (and from the next release,
you will need to reinstall as well).
You can read about this in `Writing R Extensions', but is it a good idea
to be altering a package that uses concepts
I have added the differentiation rules for asin, acos and atan to the file
deriv.c.
Please let me know where to send it.
Date: Thu, 22 Jan 2004 09:05:39 +0100
From: Martin Maechler [EMAIL PROTECTED]
Subject: Re: [R] derivative of atan(x) and similar functions
To: Timur Elzhov [EMAIL
Dear all,
I try to submit a library to CRAN but can't overcome the last R CMD CHECK.
Can someone enlighten me how to put the \usage{} section for an S3-Method
extractor defined as
args(get([.refdata))
function (x, i = NULL, j = NULL, drop = FALSE, ref = FALSE)
NULL
I read the Writing R
Vicente == Vicente Canto Casasola [EMAIL PROTECTED]
on Thu, 22 Jan 2004 14:34:01 +0100 writes:
Vicente H i, all! First of all, I'd like to apologize for
Vicente my poor English. It's for years I don't use it.
no problem at all.
Vicente This is a R-version of a function I
I think you may be trying too hard. Look at Extract.data.frame.Rd and
Extract.factor.Rd. I would just have an Extract.refdata.Rd.
You only need the \methods{generic}{class} notation when you want to
document a generic and a method, or more than one method, in the same help
file. And at
Crispin == Crispin Miller [EMAIL PROTECTED]
on Thu, 22 Jan 2004 12:08:01 - writes:
Crispin Hi,
Crispin I've got a quick question about file permissions and packages...
Crispin I'm creating my own package, and am having problems
Crispin with its vignette not being seen
Dear R users
I have two questions about estimating the spectral power of a
time series:
1) I came across a funny thing with the following code:
data(co2)
par(mfrow=c(2,1))
co2.sp1-spectrum(co2,detrend=T,demean=T,span=3)
co2.sp2-spectrum(co2[1:468],detrend=T,demean=T,span=3)
The first plot
In R, storage.mode- is an interpreted function that calls as.double (in
your case) and copies attributes across so it is definitely worse than the
second.
Is x$one usually double? If so then
.Fortran(foo, if(is.double(x$one)) x$one else as.double(x$one), ...)
would save a copy.
Does your
Martin Maechler [EMAIL PROTECTED] writes:
The update is actually available online
from http://epubs.siam.org/sam-bin/dbq/article/41801
with the extended title , 25 Years Later .
To some, that is. I can download it to the work machine, but if I do
it from home, all I get is a
...
I am a new and unexperienced user of R and got so far as to know how to produce
boxplots. I have no experience of messing with function code, so presently I do
not know how to create a boxplot with group means instead of group medians. If
somebody could help me either replace the median with
Martin Maechler wrote:
Murad == Murad Nayal [EMAIL PROTECTED]
on Wed, 21 Jan 2004 15:19:28 -0500 writes:
Murad This might have been fixed in later versions (I am
Murad using R1.7.0),
yes, the bug has been fixed long ago,
from my ChangeLog (!), it was 2003-07-18.
sorry
This will add triangles at the mean value. To change the behavior of
boxplot() to draw means instead of medians would involve rewriting the
bxp() code I believe. You could change the points in the code below to
segments. See ?segments. The archives have quite a few examples of
people modifying the
Apologies, basic question on plot.
y - c(-4,3,-2,1);
x - c(time 1, time 2, time 3, time 4);
plot(x,y, type=b);
of course fails.
x - 1:4
makes it succeed, but then I have too many ticks on my X axis. I want
exactly 4 tickmarks. It would also be nicer if I could name the
On Thu, 2004-01-22 at 12:11, Johannes Ullrich wrote:
I am a new and unexperienced user of R and got so far as to know how to produce
boxplots. I have no experience of messing with function code, so presently I do
not know how to create a boxplot with group means instead of group medians. If
On Thu, 22 Jan 2004 19:11:08 +0100
Johannes Ullrich [EMAIL PROTECTED] wrote:
I am a new and unexperienced user of R and got so far as to know how to
produce boxplots. I have no experience of messing with function code, so
presently I do not know how to create a boxplot with group means instead
On Thu, 2004-01-22 at 12:32, ivo welch wrote:
Apologies, basic question on plot.
y - c(-4,3,-2,1);
x - c(time 1, time 2, time 3, time 4);
plot(x,y, type=b);
of course fails.
x - 1:4
makes it succeed, but then I have too many ticks on my X axis. I want
exactly
On Thu, Jan 22, 2004 at 01:32:29PM -0500, ivo welch wrote:
Apologies, basic question on plot.
y - c(-4,3,-2,1);
x - c(time 1, time 2, time 3, time 4);
plot(x,y, type=b);
of course fails.
x - 1:4
makes it succeed, but then I have too many ticks on my X axis. I
You need to look at ?axis. Try this:
y - c(-4,3,-2,1)
x - c(time 1, time 2, time 3, time 4)
plot(seq(1,length(x), by = 1), y, type = b, axes = FALSE, xlab =
Times, ylab = Stuff)
axis(1, at = seq(1,length(x), by = 1), labels = x)
axis(2)
HTH, Andy
__
Is there a wavelet toolbox in R available that can be compared to the
Matlab toolbox for wavelets. Does this toolbox offer the possibility to
calculated approximations with different wavelets such as Daubechies 2,
Daubechies 3, etc, Symmlet 3, Symmlets 5, etc, Coiflets and so on. I
have tried
Hi,
Is this a bug in termplot, or (once again) do I just not understand what R
is really doing?
I am using termplot to contruct partial residual plots,
1. For all terms at once
2. One term at a time
but I get different results from these two methods. To give a concrete
example, I
I'm reading a file into a list by:
PDF = scan(file,what=character,sep=\10)
\10 is the newline character in this file, also tried \n originally
On lines that are ended by \13\10, both are dropped from the list entry
I want scan to keep the \13 in the list entry.
Is this a bug or just a strange
Hi,
How can I extract the standard deviation of the coefficients when using the lm
function. I know $coefficient gives me the individual betas. thank you
-
[[alternative HTML version deleted]]
__
[EMAIL
A while back I started looking into using R to fit kinetic models
and I'm finally getting back to it. It was suggested that I use
lsoda (thanks Peter Dalgaard). So I've tried that, even though
Peter warned me that it'd be tricky. I've put the following
example together from the lsoda help
On Thu, 22 Jan 2004, jenny smith wrote:
Hi, How can I extract the standard deviation of the coefficients when
using the lm function. I know $coefficient gives me the individual
betas. thank you
$coefficient is not the recommended way to get the coefficients.
use coef(your.model) for
Hello,
Here is my problem. I am trying to run and learn R trough the
terminal application of my Mac (it runs on Panther).
So here it goes I just recently installed R, in unix mode (terminal)
by following the attached web page. The trouble is when I try to run
either demo(graphics),
If you are interested in study each response isolated and because you have a
temporal time-series it could be a possibility for your study to modelize
that
time-serie with ARMA or GARCH process. For that you have ts and tseries
packages.
Palhoto.
Dear All,
A few weeks ago I posted a
If the model is called mymodel, you could do
summary(mymodel)$coefficients[,2]
or
sqrt(diag(vcov(mymodel)))
This gives the estimated standard errors of the coefficients, which I
think is what you wanted.
David
On Thursday, Jan 22, 2004, at 20:35 Europe/London, jenny smith wrote:
Hi,
On Thu, 22 Jan 2004, Greg Riddick wrote:
I'm reading a file into a list by:
PDF = scan(file,what=character,sep=\10)
\10 is the newline character in this file, also tried \n originally
On lines that are ended by \13\10, both are dropped from the list entry
I want scan to keep the \13 in
[EMAIL PROTECTED] (Michael A. Miller) writes:
A while back I started looking into using R to fit kinetic models
and I'm finally getting back to it. It was suggested that I use
lsoda (thanks Peter Dalgaard). So I've tried that, even though
Peter warned me that it'd be tricky. I've put the
I started to use the rgdal package in order to work with DEM files. Since
these files are based on non-rectangular quads, there are lots of
NA-values in the corresponding datasets. Does anybody know how to remove
the NA-values from an object of class GDALDataset in order to get a
rectangular
Thanks,
Right, I can see why sep=\n might grab the entire \13\10 but it seems
like sep=\10 should not strip the \13 also.
I need to read in this file (PDF file) and create a list of lines defined by
the \10 delimiter.
Any suggestions how I could use ReadBin to do that?
Not a strange
On Thu, 22 Jan 2004, Greg Riddick wrote:
Thanks,
Right, I can see why sep=\n might grab the entire \13\10 but it seems
like sep=\10 should not strip the \13 also.
Why do you claim so? Text mode is *documented* to use any of CR, CRLF or
LF as the line delimiter. Please don't keep telling
Please use the CRAN versions, not the ftp versions which are now out-of-date. Both
packages are on base CRAN now, but may not have propagated to all mirrors yet. (But
thanks to Gerald for responding-- I had de-subscribed from R-help.)
For anyone using the HANDY package (a Windows-specific
Dear R People:
When I use data.entry for a new variable, the first entry is a NA
The function works all right, but then I get a Warning message:
data.entry(xb,Modes=numeric)
Warning message:
NAs introduced by coercion
Since the function works acceptably, should I just ignore the
warning
On Thu, 22 Jan 2004, Sebastien Durand wrote:
Hello,
Here is my problem. I am trying to run and learn R trough the
terminal application of my Mac (it runs on Panther).
So here it goes I just recently installed R, in unix mode (terminal)
by following the attached web page. The trouble is
Dear Colleagues,
Our group is also working on implementing the use of R for pharmacokinetic
compartmental analysis. Perhaps I have missed something, but
fit - nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=0.5, k2=0.5)),
+data=C1.lsoda,
+start=list(K1=0.3,
Hi,
After fitting with lm, I want to get the 95 % confidence intervals of the
expected predicted value.
I know in S-plus I can use pointwise() to get the confidence intervals .
But in R, I couldn't find such a function.
Is there any R package available for such functionalities.
Please advise.
Hi everyone,
I've been looking for a function that calculates the Turnbull estimate
for left, right and interval censored data. None of the data that I am
using has exact failure times. The only function I seem to find can handle
data that has no left censoring. It seems that I will have to
Dear Michael,
One key is adjustment of nls optimizer tolerance. I notice it has to be
higher than usual, but, I recovered your noisy known parameter values
with an error of K1 (-7%) and k1 (-6%):
Miller problem with Dalgaard modifications
## Linares 1/22/2004
## Solution 1
nls(noisy ~
Hi Yun Fan!
Have you looked at predict.lm?
It can give you confidence and prediction intervals.
Hope this helps!
Sincerely,
Erin Hodgess
mailto: [EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
On Thu, 22 Jan 2004, Simon Wotherspoon wrote:
Hi,
Is this a bug in termplot, or (once again) do I just not
understand what R is really doing?
Yes, it's a bug. The lower right plot has the partial residuals for x1
not for x2. The fitted line is correct.
-thomas
I am using
Hi,
I have estimated a step function and need to take the derivatives of
this function at all points in the range. Does anyone know of any clever
ways to do this?
(I have already tried to fit a polynomial through the points in order to
obtain a smooth representation and then take derivatives
It might be helpful to know just why you want to do this. Just because one
function is a smooth approximation to another doesn't imply anything about
the derivatives approximating each other, (well, not much).
The GENERALIZED derivative of a step function can be written down explicitly
as the
63 matches
Mail list logo