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 problem is in my definition of the prior. it wants a
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 factorized. I wanted to try the package polynom
command
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, original
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
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
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
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 the
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
allan clark allan at 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
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
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 list
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*pi)
x -
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,
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
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 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),
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
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:
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:
R
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
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
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
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
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
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 misfeatures
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
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 =
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, or
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
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 names
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 know
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.
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
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 definition
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 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 abbreviation of the above
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
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' returns a value with the same shape
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 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 coordinates
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
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
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
bnone3
d
Dieter Menne dieter.menne at 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
Yen Ngo yen_h_ngo at 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
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 to be the
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
On Mon, 2006-10-09 at 15:43 +, Gregor Gorjanc wrote:
Dieter Menne dieter.menne at 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
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 somewhere
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
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
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 result.
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
Subject:
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, Jonathan
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
Thomas Wutzler twutz at 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
Marc Schwartz (via MN mschwartz at 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
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
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
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
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 work'
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]
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
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)
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 to
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
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 use
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
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, 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
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.
+ seth
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
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 hexbin
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
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
png
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
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,
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 example happens
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
79 matches
Mail list logo