hi Thomas/All,
I went through the thread(
https://stat.ethz.ch/pipermail/r-sig-finance/2006q1/000682.html which
concerns with swaps). Yeah it is correct that i would like to quote both
David and Krishna that the curve interpolation may vary considerably (for
e.g. any polynomial/parametric fit
hi Stefan
thanks for the reference but i do not have membership for the same.
Moreover i wanted to do it through R i.e. making of spot yield curve
through bootstrapping.
can you provide some links in R for the same.
-gaurav.
Stefan Grosse <[EMAIL PROTECTED]>
09-10-06 07:47 PM
To
[EMAIL P
On 10/9/06, Deepayan Sarkar <[EMAIL PROTECTED]> wrote:
> Just out of curiosity, what's wrong with
>
> xyplot(y1 + y2 ~ x | interaction(g, h, sep = "-"), data = DF, type =
> "l", auto.key = TRUE)
The reason would be that I want to add information to the plot that is
not of the form "i-j". My exampl
On 10/9/2006 7:58 PM, Richard A. O'Keefe wrote:
> I wrote:
> > R is, and should remain, editor-neutral.
> > Amongst other things, it should NOT acquire misfeatures
> > in order to support editors that happen not to support comment-region.
> > Block comments are indeed editor-neutral
Dear R-helpers,
I want to smooth a scatter plot (or 2d histogram) to generate the density
map using a squared Euclidean kernel function, with each data point
contributing a density of 1/(r*r + k) to each cell on the heatmap, where
r*r was the squared Euclidean distance between the coordinates of t
Hi,
I try to use persp to draw perspective plots of surfaces in a region.
the data matrix is 606x354
here is my simple code:
r<-read.table (file='prod.asc.txt', header=F)
r[r==-]<-NA # no values in surrounding region
r<-as.matrix(r)
y<-seq(1:ncol(r)) #354 cols
x<-seq(1:nrow(r)) #606 rows
pn
I wrote:
> R is, and should remain, editor-neutral.
> Amongst other things, it should NOT acquire misfeatures
> in order to support editors that happen not to support comment-region.
> Block comments are indeed editor-neutral,
> but they do not solve any problem that R currently
A little algebra up front can often help.
Note that
X <- shiftMatrixL(X, 1)*3 + shiftMatrixL(X,2)*5
(and similarly with more terms on the r.h.s)
is just
X <- X %*% mat
where mat is is a matrix of zeroes except that
mat[ i+1, i ] == 3
mat[ i+2, i ] == 5
and d
Hi
Li, Sue wrote:
> Dear R-users,
>
> I try to plot multiple plots of hexbin in the same page. However, par
> does not work when hexbin is used. Neither do xlim and ylim. Any ideas?
hexbin is now based on the grid package, so does not work with par().
You might want to take a look at the hexb
If you're able to work with the transpose of your matrix, you might
consider the function 'filter()', e.g.:
> filter(diag(1:5), c(2,3), sides=1)
Time Series:
Start = 1
End = 5
Frequency = 1
[,1] [,2] [,3] [,4] [,5]
1 NA NA NA NA NA
234000
30660
"Your EPEC ICT Team - Ricardo Rodríguez"
<[EMAIL PROTECTED]> writes:
> Please, does R support JDBC access to MySQL databases? Could you point
> me in the right direction to find some documents about this issue?
Have you found the RMySQL package? It allows you to talk to MySQL
databases from R.
+
Hi, is there any package/function that can tell if a
numeric vector (continuous data) has a bimodal or
trimodal distribution and caluclate the location of
the corresponding modes?
Thanks
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mail
I'm wondering what's the best way to shift a huge matrix left or right.
My current implementation is the following:
shiftMatrixL <- function(X, shift, padding=0) {
cbind(X[, -1:-shift], matrix(padding, dim(X)[1], shift))
}
X <- shiftMatrixL(X, 1)*3 + shiftMatrixL(X,2)*5...
However, it's still
Hi all,
Completely new to this list and completely new to R so, please, accept my
apologies in advance if any question is too basic or out of place. Thanks.
After reading some docs, googling around and searched the mail lists I've only
found spare references to accessing databases. Among thes
On 10/7/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> I have noticed that dispatch on functions seems not to work
> in another case too. We define + on functions (I have ignored
> the niceties of sorting out the environments as we don't really
> need it for this example) but when we try to u
Hi;
I have a matrix of 154 elements by 66241 sub-elements. The elements
are chains of characters, sub elements are simply sub-chains of a
certain length. For each element, I computed a count of the ocurrence
of sub-elements (scan of strings). I thus have a matrix of numerical
values (between 0 and
Kamila Naxerova wrote:
>
> Is there a way to plot elliptical shapes? symbols()
> only provides circles...
>
Elipse centered at (x0, y0) with axes _a_ and _b_:
theta <- seq(0, 2 * pi, length=(big enough))
x <- x0 + a * cos(theta)
y <- y0 + b * sin(theta)
plot(x, y, type = "l")
:-)
Or, if you want
Hi All,
I have measured set of induced expression values in
tumor samples with respect to normal samples. I
consider the expression value in normal samples as 1.
So, my final set of values (in total 10 values) are:
1:2.98 (normal:tumor), 1:1.5, etc. Is it possible to
find significant value(s) amo
Dinner time with a nice bit of red is a great occasion to ponder over
this sort of thing.
It turns out that I had added a path to my DYLD_LIBRARY_PATH that held
my own stats.so and constants.so . The former must explain why R's
stats.so didn't load, the latter why I had similar problems with a
Pyt
The cnvrt.coords fnuction in the TeachingDemos package may simplify the
procces of getting the coordinates.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[
On 10/9/06, Ritwik Sinha <[EMAIL PROTECTED]> wrote:
> Hi Gabor,
>
> Thanks for your reply. The solution you gave me and the one that I had
> attempted have the same problem. The plots don't differentiate between
> y1 and y2, in other words, they loose the groups attribute. This is
> what 'does not
On 10/9/2006 1:28 PM, Kamila Naxerova wrote:
> Is there a way to plot elliptical shapes? symbols() only provides circles...
There are several, including the ones in the ellipse package.
Duncan Murdoch
__
R-help@stat.math.ethz.ch mailing list
https://st
Is there a way to plot elliptical shapes? symbols() only provides circles...
Thanks!
Kamila
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.htm
On 10/9/2006 12:50 PM, Michael Jerosch-Herold wrote:
> Dear readers:
>
> I like very much how one can hightlight R-code sections in Tinn-R and send
> them to R for execution. Is there an editor for the Mac that similarly allows
> you to highlight code and send it to the R-console?
There's the
Marc Schwartz (via MN mn.rr.com> writes:
> If the content of this particular transcript is likely to be static,
> consider an alternative of making it available in a PDF document that is
> linked on that page.
>
> Then others can perhaps contribute by providing other relevant content
> in that s
Thomas Wutzler bgc-jena.mpg.de> writes:
> without grouped data (athough I recommend using grouped data after
> reading the book), the call can be formulated as:
Well, this is probably the only point where I disagree with Douglas Bates. I
found using grouped data introduces a level of indirection
Dear readers:
I like very much how one can hightlight R-code sections in Tinn-R and send them
to R for execution. Is there an editor for the Mac that similarly allows you to
highlight code and send it to the R-console?
Thank you in advance for recommendations!
Michael Jerosch-Herold
___
tryCatch() is doing exactly what you need. I consider tryCatch() an
updated and more flexible version of try() making the latter history.
/H
On 10/9/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> See:
>
> http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg09925.html
>
> On 10/9/06, Jonat
Thanks! I believe I have all I need to solve my issue.
Thank you very, very much,
mihai
- Original Message
From: Roger Bivand <[EMAIL PROTECTED]>
To: Mihai Nica <[EMAIL PROTECTED]>
Cc: Mihai Nica <[EMAIL PROTECTED]>; r-help@stat.math.ethz.ch
Sent: Monday, October 9, 2006 9:20:36 AM
Subj
See:
http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg09925.html
On 10/9/06, Jonathan Williams
<[EMAIL PROTECTED]> wrote:
> Dear R Helpers,
>
> I want to test if a procedure within a loop has produced an error or not.
> If the procedure has produced an error, then I want to ignore its resu
> I just installed R 2.4 on a windows machine and installed the lme4
> package
> version 0.9975-1. Unfortunately this version requires the 0.9975-2
> version of
> the Matrix package which I cannot seem to find. It seems the newest
> version I
> can find is 0.9975-1. Is there somewhere I could
Dear R Helpers,
I want to test if a procedure within a loop has produced an error or not.
If the procedure has produced an error, then I want to ignore its result.
If it has not produced an error, then I want to use the result. The problem
In order to run the loop without crashing if the procedure
Mark Lyman wrote:
> I just installed R 2.4 on a windows machine and installed the lme4 package
> version 0.9975-1. Unfortunately this version requires the 0.9975-2 version of
> the Matrix package which I cannot seem to find. It seems the newest version I
> can find is 0.9975-1. Is there somewh
On Mon, 2006-10-09 at 15:43 +, Gregor Gorjanc wrote:
> Dieter Menne menne-biomed.de> writes:
>
> >
> > Dear friends of lmer,
> >
> > http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
> >
> > I have put a transcript of the long thread on lmer/lme4 statistical test
> > into the W
Dear Douglas Bates and Spencer Graves,
thank you for your replies.
In the meantime I got the book and will recommend it to my colleages :)
I learned the answer to my problem and want to post it here:
If you use ungrouped data (a plain dataframe) you must specifiy the
grouping variable in the ra
On Wed, 4 Oct 2006, paul sorenson wrote:
> What is the recommended way to convert/coerce and integer to a POSIXct
> please?
>
> d <- as.POSIXct(Sys.Date())
> i <- as.integer(d)
>
> as.POSIXct(i)
> Error in as.POSIXct.default(i) : do not know how to convert 'i' to class
> "POSIXlt"
>
> This appears
Yen Ngo yahoo.se> writes:
>
> Hi all R-helpers,
>
> i am a new R-user and have problem with deleting some components in a
column. I have a dataset like
>
> Name Idx
>empty 2
>empty 3
> anone2
> bnone3
> d none 2
Dieter Menne menne-biomed.de> writes:
>
> Dear friends of lmer,
>
> http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
>
> I have put a transcript of the long thread on lmer/lme4 statistical test
> into the Wiki. For all those who missed it life, and for those like me, who
> don't l
On Mon, 2006-10-09 at 17:15 +0200, Yen Ngo wrote:
> Hi all R-helpers,
>
> i am a new R-user and have problem with deleting some components in a
> column. I have a dataset like
>
> Name Idx
>empty 2
>empty 3
> anone2
> bnone
Well... this may be a red herring after all, but it's an interesting one!
It turns out I have got DYLD_LIBRARY_PATH set. This gets exported to R
(or rather, affects the loader that loads the R process?), and somehow
mixes up the way symbols are exported. Something that applies to
launching R from
I just installed R 2.4 on a windows machine and installed the lme4 package
version 0.9975-1. Unfortunately this version requires the 0.9975-2 version of
the Matrix package which I cannot seem to find. It seems the newest version I
can find is 0.9975-1. Is there somewhere I could find this newer
On Mon, 2006-10-09 at 10:56 -0400, Mike Wolfgang wrote:
> Hi list,
>
> I want to add some lines at margin area of one figure. mtext could add text
> to these margins, can I add lines with different lty parameters? Thanks,
>
> mike
You can do it, but it will take some fiddling to get the coordina
Hi all R-helpers,
i am a new R-user and have problem with deleting some components in a column.
I have a dataset like
Name Idx
empty 2
empty 3
anone2
bnone3
d none 2
ad cfh 4
bf cdt 5
On Sat, 7 Oct 2006, Rolf Turner wrote:
> Peter Dalgaard writes:
>
>> Alberto Vieira Ferreira Monteiro <[EMAIL PROTECTED]> writes:
>>
>>> Why this kind of assignment does not work?
>>>
>>> n <- 1
>>> f <- ifelse(n == 1, sin, cos)
>>> f(pi)
>>
>> It's not supposed to.
>>
>> 'ifelse' retur
On Fri, 6 Oct 2006, Deepayan Sarkar wrote:
> On 10/6/06, Ted Harding <[EMAIL PROTECTED]> wrote:
>> Hi Folks,
>>
>> Given a series of n independent Bernoulli trials with
>> outcomes Yi (i=1...n) and Prob[Yi = 1] = Pi, I want
>>
>> P = Prob[sum(Yi) = r] (r = 0,1,...,n)
>>
>> OK, some abbreviati
Hi list,
I want to add some lines at margin area of one figure. mtext could add text
to these margins, can I add lines with different lty parameters? Thanks,
mike
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
ht
Apparently misunderstanding me,
On 10/8/06, Richard A. O'Keefe <[EMAIL PROTECTED]> wrote:
> Amongst other things, it should NOT acquire misfeatures
> in order to support editors that happen not to support comment-region.
...
One can only comment from one's experience. We have found the concept
of
On 10/9/06, hadley wickham <[EMAIL PROTECTED]> wrote:
> > Current .Rd documentation has some obvious problems:
> >
> > - the parser strips comments out of examples when it runs them
> > - there's no way to put images into the documentation
> > - the keywords aren't much use
> > - there's isn't a de
Dear R-users
I am trying to analyse the interaction between men and
women participation in the labour market. So I have
too equations of participation like this
h1=h1(w1,X1,h2) and h2=h2(w2,X2,h1) where 1=men,
2=women, h=hours worked, w=wages and X=taste
variables.
A basic approach would be to us
Apologies for cross-posting.
We would like to announce a 3-day R course in Lisbon, Portugal.
Full details can be found at:
http://www.brodgar.com/statscourse.htm
Kind regards,
Alain Zuur
Dr. Alain F. Zuur
First author of:
Analysing Ecological Data (2007).
Zuur, AF, Ieno, EN and Smith, GM.
Spr
> Current .Rd documentation has some obvious problems:
>
> - the parser strips comments out of examples when it runs them
> - there's no way to put images into the documentation
> - the keywords aren't much use
> - there's isn't a definition anywhere of what the format really is, so
> it's hard to
On Mon, 9 Oct 2006, Mihai Nica wrote:
> The answer to all your questions is simple. By the time I get a little
> grasp on things, they become "deprecated" :-). But my programming skills
> are so low, that I find this normal.
> My problem comes from the last two #comment lines :
> #(if data row n
Although it does not directly answer your question you might be
interested in:
Michalis Ioannides, A comparison of yield curve estimation techniques using UK
data, Journal of Banking & Finance, Volume 27, Issue 1, , January 2003, Pages
1-26.
Abstract: I compare different methods of estimating t
The answer to all your questions is simple. By the time I get a little grasp on
things, they become "deprecated" :-). But my programming skills are so low,
that I find this normal.
My problem comes from the last two #comment lines :
#(if data row names do not match polygon IDs, will reorder, o
Use panel.superose instead of panel.xyplot. Also you might
want to use grid.text instead of panel.text so that you can use
npc coordinates (i.e. 0 to 1) making it independent of the values
on the panel:
library(lattice)
library(grid)
xyplot(y1 + y2 ~ x | g * h, data = DF, type = "l", auto.key = T
Hi Gabor,
Thanks for your reply. The solution you gave me and the one that I had
attempted have the same problem. The plots don't differentiate between
y1 and y2, in other words, they loose the groups attribute. This is
what 'does not work' means. I am attaching a ps file (I hope that does
not ge
On 10/8/2006 8:03 PM, Richard A. O'Keefe wrote:
> Apparently misunderstanding the argument,
> BBands <[EMAIL PROTECTED]> wrote:
> Should R be editor specific? Or should it be editor neutral?
>
> R is, and should remain, editor-neutral.
> Amongst other things, it should NOT acquire misf
Hi All,
Thanks for all your discussions which I read offline and enhance my
knowledge.
I just want to know how to make yield curve by bootstrapping using R.
Please give me some subject-matter links and code links.
Especially the Spot Yield Curve and YTM Curve.
Sayonara With Smile & With Wa
To make things even weirder, everything seems to work just fine when I
launch R via the v1.17 3878 build of the R GUI for 2.4 ...
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http:
On 10/9/2006 6:48 AM, Tong Wang wrote:
> Hi,
> I am using the debug package, is there any way I can tell the debuger to
> stop when a warning message come out ?
There are some instructions here:
http://www.stats.uwo.ca/faculty/murdoch/software/debuggingR/#tryCatch
Duncan Murdoch
_
Thank you for the pointer to the FAQ. Thought I had searched the FAQ
thoroughly, obviously I didn't.
Unfortunately, my stats aren't up to fully understanding the explanation
and the proposed solution in the FAQ.
>For the time being, I would recommend using a Markov Chain Monte Carlo
>sample (fu
Unfortunately, no. No red herring, no success either. I removed the
.RData file that was indeed lurking, even moved aside .Rprofile and
the .R directory, installed the latest patch from r.research.att.com
... and still the claim about R_loess_raw
This symbol is called upon in simpleLoess, by the w
On Mon, 2006-10-09 at 13:25 +0200, René J.V. Bertin wrote:
> Hello,
>
> Since I (finally) upgraded to 10.4 "Tiger", I've been experiencing
> some strange behaviour with dynamic libraries.
>
> When I now launch R, and try to use the loess() function from the
> stats package, I get the following:
>
the following seems to work:
test.summary$"Error: Within"[[1]]$"Pr(>F)"
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32
Dear friends of lmer,
http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
I have put a transcript of the long thread on lmer/lme4 statistical test
into the Wiki. For all those who missed it life, and for those like me, who
don't like the special style of the R-list to keep full length q
Hi all,
In assessing the concordance of two sets of ANOVA results, I would like
to extract the p values from the summaries. I had a look at previous
answers to this problem, and the suggestions do not seem to work:
# toy example
test.df<-data.frame(subno=c(1:50,1:50),
ndrinks=c(rpois(50,4.5),rp
Hello,
Since I (finally) upgraded to 10.4 "Tiger", I've been experiencing
some strange behaviour with dynamic libraries.
When I now launch R, and try to use the loess() function from the
stats package, I get the following:
> R
R version 2.4.0 (2006-10-03)
Copyright (C) 2006 The R Foundation for
Hi,
I am using the debug package, is there any way I can tell the debuger to
stop when a warning message come out ?
Thanks a lot.
best
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posti
Dear R-help,
Why can't lme cope with an incomplete whole plot when analysing a split-plot
experiment? For example:
R : Copyright 2006, The R Foundation for Statistical Computing
Version 2.3.1 (2006-06-01)
> library(nlme)
> attach(Oats)
> nitro <- ordered(nitro)
> fit <- lme(yield ~ Variety*nitro
"S.Q. WEN" <[EMAIL PROTECTED]> writes:
> Hi,
> I want to get random number which is uniformly distributed on the unit
> disc.
> How can I do that with R?
Most easily, although perhaps not most efficiently, by using spherical
coordinates:
N <- 1
r <- sqrt(runif(N))
theta <- runif(N,0,2*p
S.Q. WEN wrote:
> Hi,
> I want to get random number which is uniformly distributed on the unit
> disc.
> How can I do that with R?
>
>
> Best wishes,
> WAN WAN
>
> [[alternative HTML version deleted]]
>
> __
> R-help@stat.math.ethz.ch mailing
Hi,
I want to get random number which is uniformly distributed on the unit
disc.
How can I do that with R?
Best wishes,
WAN WAN
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listin
allan clark stats.uct.ac.za> writes:
>
> hi all
>
> i would like to use r and winbugs in order to undertake seemingly unrelated
regression. i am using the
> R2WINBUGS library. i just simulated a simple example (sample size is 25) in
order to get the correct code.
>
> i suspect that the proble
Hi All,
Thanks for all your discussions which I read offline and enhance my
knowledge.
I just want to know how to make yield curve by bootstrapping using R.
Please give me some subject-matter links and code links.
Especially the Spot Yield Curve and YTM Curve.
Sayonara With Smile & With Wa
On Sun, 8 Oct 2006, Thomas P. Colson wrote:
> I am able to convert the flow accumulation grid into an area (for each
> pixel) grid, then import this into R as an ASCII file. The plot(ecdf)
> function in R seems to plot the opposite: curve starts at probability 0, for
> drainage area 0, should be t
HelponR wrote:
> I am trying to use lm to do a simple regression but on a batch of
> different files.
>
> Each file has different column names.
>
> I know the first column is the dependent variable and all the rest are
> explanatory variables.
I believe
> lm(data = thedataframe)
(i.e. with no f
Serguei Kaniovski <[EMAIL PROTECTED]> writes:
> Hi all,
>
> starting from a vector "v[1:n]" I would like to compute the coefficients
> of the polynomial (1+x^v[1])*(1+x^v[2])*...*(1+x^v[n]). The following
> code works but is extremely slow for a large "n" due to, I believe, the
> polynomial being
Hi, All,
I have come across a problem building a DLL from .f90 source (R 2.3.1,
Windows XP). When using the R CMD SHLIB procedure, the DLL itself was
being built, but its export table was empty.
Among the output from R CMD SHLIB the following message appeared:
c:\mingw\bin\nm.exe: 'a.out': No
"Richard A. O'Keefe" <[EMAIL PROTECTED]> writes:
> We find
> - languages which have copied PL/I (/* ... */)
> - languages which have copied BCPL (//)
> - languages which have copied Pascal (* ... *)
Not to stick my head further into the discussion at this point, but
actually, origina
79 matches
Mail list logo
