Dear people
I've learnt that by using the "boxplot.stats" command in the "grDevices" library
I can get the 5-number summaries of a boxplot, plus other important information,
like the confidence interval around the median.
I'm interested in knowing the actual formula to used in that package to
cal
Dear R colleagues,
a friend and I are trying to develop a modest workflow for the problem
of decomposing tests of higher-order terms into interpretable sets of
tests of lower order terms with conditioning.
For example, if the interaction between A (3 levels) and C (2 levels)
is significant, it ma
Hi all,
I've been trying to run Sweave with R code embedded - using the ggplot2
package and in particular the qplot command. There appears to be a problem
in Sweave not picking up that qplot is a function. Has anybody else tried
to use qplot in Sweave and have you been successful? Any help woul
On Sunday 03 August 2008 01:54:43 pm Andrew Ramsey wrote:
> Hello--
>
> I am a relatively new user to R and I cannot find the information I
> need. Please help.
>
> I have a very large data set with values including letters, numbers,
> and symbols (sometimes within the same vector value [ie X9-].
>
~ The beta distribution only applies to data that are bounded
between 0 and 1 (in some cases strictly bounded, i.e. 0 wrote:
|> charter.net> writes:
|>
|>> If I have some data that based on the historgram and other plots it
"looks"
|> like a beta distribution. Is there
|>> a function or functi
Cody Hamilton Edwards.com> writes:
>
> Is there a way to create a 'bubble plot' in R?
>
> For example, if we define the following data frame containing the level of y
observed for 5 patients at three
> time points:
>
> time<-c(rep('time 1',5),rep('time 2',5),rep('time 3',5))
> y<-c('a','b','c
Exactly how is the translation of "X9-" -> 00911 done? Are there
unique mappings of character sequences to numbers? How many different
ones might there be? Why do you have leading zeros on the result?
If they are changed to numeric, then the default printing results in
'911'. Can you provide a
Hello--
I am a relatively new user to R and I cannot find the information I
need. Please help.
I have a very large data set with values including letters, numbers,
and symbols (sometimes within the same vector value [ie X9-].
I've imported the data using read.fwp and it arrives in list fo
charter.net> writes:
> If I have some data that based on the historgram and other plots it "looks"
like a beta distribution. Is there
> a function or functions within R to help me determine the model parameters for
such a distirbution?
library(MASS)
?fitdistr
_
This may be a begining question. If so, please bear with me.
If I have some data that based on the historgram and other plots it "looks"
like a beta distribution. Is there a function or functions within R to help me
determine the model parameters for such a distirbution? Similarily for other
"
You're right, the problem was with the print function not displaying enough
precision.
Thanks!
--
View this message in context:
http://www.nabble.com/Convert-date-to-decimal-days-tp18800462p18802117.html
Sent from the R help mailing list archive at Nabble.com.
_
I think you will find that the 'catastrophic loss of significance' is in
the printing at the default level. Try print(x, digits=15).
On Sun, 3 Aug 2008, cls59 wrote:
Hello all,
I have a quick question about formatting date strings..
Datetime strings, methinks. If you do just want days, u
If length( levels( a1$cat ) ) is not very large, and the
structure suggested in the toy example holds in the actual case (probes in
a2 do not 'overlap' and are in order according to a2$st and a2$en, then
this might do:
unsplit(lapply(levels(a1$cat), function(x) {
+ tmp1 <- subset( a1, ca
Hello all,
I have a quick question about formatting date strings..
I am currently debugging some Matlab code someone else wrote and since it is
so bad that I have to go through it line by line I figured that I would just
rewrite the thing in R.
The code produces plots of wave spectra with decim
On Sun, 3 Aug 2008, Dirk Eddelbuettel wrote:
On Sun, Aug 03, 2008 at 05:38:01PM +0100, Prof Brian Ripley wrote:
The Pemtium D 940 supports EMT64 according to
http://www.pcstats.com/articleview.cfm?articleID=1988 (I googled it).
That is the key: EMT64 is Intel's name for x86_64 and so your box
Hello,
I have encountered results which I am not sure how to interpret when using
anova.gam to compare 2 different models. For certain tests the results do not
include an F- or associated p-statistic. This happens when comparing certain
models and not others, and I do not discern a patten exp
On 03/08/2008 12:13 PM, Jonas wrote:
Duncan, Patrick,
Sorry for not being very clear (and having non-working code in my
example). What I'm actually trying to achieve is creating a wrapper
around the NelsonSiegel function in the fBonds package in order to
(among other things) create my own plot f
On Sun, Aug 03, 2008 at 05:38:01PM +0100, Prof Brian Ripley wrote:
> The Pemtium D 940 supports EMT64 according to
> http://www.pcstats.com/articleview.cfm?articleID=1988 (I googled it).
>
> That is the key: EMT64 is Intel's name for x86_64 and so your box should
> run x86_64 Linux which I belie
The Pemtium D 940 supports EMT64 according to
http://www.pcstats.com/articleview.cfm?articleID=1988 (I googled it).
That is the key: EMT64 is Intel's name for x86_64 and so your box should
run x86_64 Linux which I believe Debian/Ubuntu refer to as 'amd64' (which
understandably Intel does not f
I believe
?symbols
is what you are looking for. Of course, you need to convert categories to
numeric values for plotting:
For example, adding some actual data to your data frame:
> D <- cbind(D, P = runif(15))
> symbols(as.numeric(D$time), as.numeric(D$y), circles=D$P, inches=0.2,
ann=F, xaxt="n
On Sun, Aug 03, 2008 at 06:19:54PM +0200, pir2.jv wrote:
>
> - Peut-on poser des questions en français?
Pas vraiment. R-help est une liste organisee en anglais.
> - Très simple. Je fais un programme "truc.r" (sous emacs)
> Ce programme réalise des graphiques.
> Je voudrais les enregistrer au fur
On Sun, 3 Aug 2008, Thomas P C Chu wrote:
Thanks to Mr Dalgaard for his advice and everyone else who has contributed.
Inclusion of an error term at the end of sim.set$y = ... line did cure my
problems with drop1() and step().
I suppose it is my own inexperience in carrying out simulations cau
- Peut-on poser des questions en français?
- Très simple. Je fais un programme "truc.r" (sous emacs)
Ce programme réalise des graphiques.
Je voudrais les enregistrer au fur et à mesure -- sinon, quand la
session est finie, ces graphiques sont perdus--
J'ai essayé:
pdf("truc.pdf")
Mais le fichi
Duncan, Patrick,
Sorry for not being very clear (and having non-working code in my
example). What I'm actually trying to achieve is creating a wrapper
around the NelsonSiegel function in the fBonds package in order to
(among other things) create my own plot function. However, all of the
objects i
Duncan Murdoch wrote:
On 03/08/2008 11:29 AM, Jonas wrote:
hi,
my apologies if this has been covered numerous times before and it's
only my lack of search and/or R skills that is stopping me from
finding the solution.
i'm trying to access an object defined and created inside a function
(the ob
Running R version 2.6.1 under Gentoo Linux, I'm trying to produce a
thematic map of the USA using the gmaps package. The result thus far
has two problems from my point of view.
First, the projection (Miller cylindrical?) elongates southern states
and flattens northern ones unattractively. I'd
Thanks to all for your replies. It is appreciated.
Prof Brian Ripley wrote:
>
> On Sat, 2 Aug 2008, Burgs wrote:
>
>>
>> Yes, I'm using Windows. I'm getting an error when I use the source()
>> command, it seems not to like the path to the .R file I'm putting in,
>> although it is the correc
On 03/08/2008 11:29 AM, Jonas wrote:
hi,
my apologies if this has been covered numerous times before and it's
only my lack of search and/or R skills that is stopping me from
finding the solution.
i'm trying to access an object defined and created inside a function
(the object is not returned by
hi,
my apologies if this has been covered numerous times before and it's
only my lack of search and/or R skills that is stopping me from
finding the solution.
i'm trying to access an object defined and created inside a function
(the object is not returned by the function), but i can't seem to get
Hi List,
I have run R from Windows Terminal server and it works fine. I now need to
install it on apache Web server / Amazon EC 2 server to test R (with a GUI)
for computing on a cloud (using a framework I wrote about).Intent is to
benchmark R 's superior algorithms on very big files > 1-2 gb (aga
After doing some reading about 64-bit systems and software I am still
somewhat uncertain about some things. I have a Dell Dimension XPS 400 with a
dual core Intel Pentium D 940 (3.2 GHz) and 4 Gb of memory. I currently dual
boot the system with XP Professional and Ubuntu 8.04 (32 bit). If I simply
Hi R users,
I wonder if there is any function which would render the value of the
neighbors of the given element [i,j] of a matrix.
Thanks,
Rostam
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/ma
Constructing a zoo object is well covered in ?zoo
and in the three vignettes:
vignette("zoo")
vignette("zoo-quickref")
vignette("zoo-faq")
On Sat, Aug 2, 2008 at 10:27 AM, rcoder <[EMAIL PROTECTED]> wrote:
>
> Thanks for your reply Gabor.
>
> As I already have a dates column in my dataframe (in
On Sun, Aug 3, 2008 at 3:51 AM, Arthur Roberts <[EMAIL PROTECTED]> wrote:
> How do I get an inset graph (i.e. graph within a graph)? Your input is
> greatly appreciated.
See ?subplot in the TeachingDemos package.
__
R-help@r-project.org mailing list
ht
Hi Arthur,
This can be done quite easily using the appropriate arguments listed under
?par; and there are other approaches. Ready-made functions exist in several
packages. I tend to use ?add.scatter from package ade4. It's a short
function, so it's easy to customize it, but it works well straight
Thanks to Mr Dalgaard for his advice and everyone else who has
contributed. Inclusion of an error term at the end of sim.set$y = ...
line did cure my problems with drop1() and step().
I suppose it is my own inexperience in carrying out simulations caused
such gaffe.
Thomas
__
If you are using zoo for this you can do:
library(zoo)
set.seed(1)
x <- zoo(rnorm(7), as.numeric(1:7))
y <- zoo(rnorm(4), c(1, 2, 6, 7))
xy <- merge(x,y)
xy$y - xy$x
# or if you just want the difference at times existing in both
y - x
On Sat, Aug 2, 2008 at 5:16 PM, <[EMAIL PROTECTED]> wrote:
is this gam::gam or mgcv::gam? If the latter, and if I understand what
you want, if gam.out contains your GAM model fitted using gam() from
mgcv, then:
out.term <- predict(gam.out, type = "terms")
gives you the contributions of each covariate on the response for the
fitted model. read :?predict.g
Not quite true that you can't type anything. What happens (for me) is
that you are still in reverse-i-search, so you can get this effect
from "^R l ^C d".
>
(reverse-i-search)`l': ls()
(reverse-i-search)`l': ls()
>
(reverse-i-search)`l': ls()
(reverse-i-search)`ld': levels(ftpain3) <-
list
Peter Dalgaard wrote:
Gad Abraham wrote:
Prof Brian Ripley wrote:
A patch to do this was posted on 2007-09-29 by Glenn Davis. Some
people not addicted to Matlab find the behaviour very inconvenient
and prefer the getline/readline behaviour (triggered by ^R/^S) of
Rterm and R on Unixen.
Now
Em Sex, 2008-08-01 às 14:16 -0700, Rajasekaramya escreveu:
> hi
>
> I have list of matrix of lenggth 61 containg the mean values..I want
to make
> a boxplot for each of the matrix.
> I used a for loop but i cant figure out the way to save in the
boxplots
>
> > all.the.mean
> [[1]]
> mean
> 0.5
>
Also, your request can easily be formulated as an SQL statement,
for example utilizing the 'sqldf' package:
library(sqldf)
a1 <- data.frame(id = 1:6,
cat = paste('cat', rep(1:3, c(2,3,1))),
st = c(1, 7, 30, 40, 59, 91),
en = c(5, 25, 39, 55
First of all, thank you for your reply and for the links.
I see that my problem might be overly detailed in its description and not
very clear in its outlines.
In fact the question was: If there is any existing module for R to evaluate
data, with parameters being only boolean values. (Which u ge
I would like to request help with the following:
I am trying to use a Generalized Additive Model (gam) to examine the density
distribution of fish as a function of latitude and longitude as continuous
variables, and year as a categorical variable. The model is written as:
gam.out <- gam(Dens
Thanks for your reply Gabor.
As I already have a dates column in my dataframe (in column row.names), is
it possible to preserve this whilst still making a data set suitable for
rollapply()?
Thanks,
rcoder
Gabor Grothendieck wrote:
>
> If its regular you can convert it to ts or zoo.
> If its
Hi everyone,
I'm experiencing difficulty getting the results I want when I use a nested
for loop. I have a data set to which I perform some calculations, and then
try to apply a regression over a rolling window. The code runs, but the
regression results I am getting (intercept and slope) are simp
Hi All,
the new release v1.4-0 of the "sensitivity" package is now available on CRAN
The "sensitivity" package is devoted to factor screening and global sensitivity
analysis
of numerical model output.
Here are the new features list:
New functionalities:
* sequential bifurcation method
Hi the list,
A new version (0.92) of my package 'noia' will be available soon on CRAN
mirrors, and I think it might be a good opportunity to introduce it
shortly to the R community.
In summary: 'noia' will be of absolutely no interest for 99.99% of you.
The 0.01% remaining are quantitative g
Hi, all,
How do I get an inset graph (i.e. graph within a graph)? Your input
is greatly appreciated.
Best wishes,
Art
University of Washington
Department of Medicinal Chemistry
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/list
49 matches
Mail list logo