Antje yahoo.de> writes:
> I've just discovered that the following code leads to boxplot
> (surprisingly to me).
> Can anybody explain to me why? Is this documented somewhere? I've never
> consider this option before.
> x <- rnorm(300)
> l <- c(rep("label1",100), rep("label2",50), rep("label3",15
Hi,
thank you both for your response.
I don't want to do anything like this - I just got some code like this from
someone else and was wondering about the result.
I would have used another approach to create a boxplot like this...
Ciao,
Antje
[EMAIL PROTECTED] schrieb:
hi: i'm not well verse
ArunPrasad wrote:
Hi,
Can anyone help me to understand why I am getting the error message like
below?
Tree.result <- PP.Tree("LDA",iris[train,5],iris[train,1:4])
These are *Warning* messages rather than error messages. These warnings
come from coding problems in the functions you are us
Building a package using R.2.8.0 on Windows XP gives problems when compiling C
code: Trying to run 'rcmd build' on R.2.8.0 on Windows XP gives
...
-- Making package gRbase
adding build stamp to DESCRIPTION
making DLL ...
gcc -std=gnu99 -Ic:/Programs/R/current/include
Hi,
how can I do the following in R?:
vec <- c("a", "b", "c")
magic.
print(str)
"abc"
or even better:
"(a,b,c)"?
Cheers
--
View this message in context:
http://www.nabble.com/merge-character-strings-tp20790116p20790116.html
Sent from the R help mailing list archive at Nabble.com.
__
check at paste(), e.g.,
vec <- c("a", "b", "c")
paste(vec, collapse = "")
I hope it helps.
Best,
Dimitris
mentor_ wrote:
Hi,
how can I do the following in R?:
vec <- c("a", "b", "c")
magic.
print(str)
"abc"
Cheers
--
Dimitris Rizopoulos
Assistant Professor
Department of Biostat
Since SVG has a lot of elements, Firefox is still under development to
support different elements in SVG; see:
http://www.mozilla.org/projects/svg/status.html
Currently animations won't work under Firefox 3.0.4, and all the rest
can be viewed using Firefox 3.0.4.
Regards,
Yihui
--
Yihui Xie <[EM
I would suggest sowas (Douglas Maraun) or wmtsa (Percival and Walden).
I don't know for sure but both of these probably rely on Rwave. I
like both of the packages mentioned it just depends on what you want
to do- wmtsa will do DWT and CWT + a lot more, and sowas is for
testing hypotheses and wav
Hello everyone,
i have been searching the whole day, trying ANY solution offered by the
internet and before headbanging against the wall would like to try asking
for your help.
As many others, i am haveing problems loading the Rcmdr package with Macbook
running Tiger 10.4.11. X11 version updated
Dear all,
I am facing the same problem. I followed the all instructions given in help
file and some steps of previous discussion (from google search) but I could not
install package Rgraphiviz in windows. I tried even the instructions given to
install Bioconductor package then only able to inst
Since I'm a SAS programmer, I'm used to creating command files in an editor for
submission later. Is there a way to do this in R? I'd need to retain an ouput
listing and a log to check for errors.
_
Send e-mail faster without impr
John Sorkin wrote:
R 2.7
Windows XP
I have two model that have been run using exactly the same data, both fit using glm(). One model is
a linear regression (gaussian(link = "identity")) the other a quasipoisson(link =
"log"). I have log likelihoods from each model. Is there any way I can de
Hello,
I am trying to read a dataset with 100,000 rows and around 365 columns
into R, using read.table/read.csv.
In Windows XP, with R 32 bit, I am able to read only 15266 rows and
not more than that.
I tried the same in R running in Ubuntu and it does the same and reads
only 15266 rows.
Using the
Hello,
after loading package "fGarch" the function is.na.data.frame() is behaving
different:
> sessionInfo()
R version 2.8.0 (2008-10-20)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_
Try
> source('myFirstScript.R')
Where myFirstScript.R as the following line
x <- rnorm(100)
y <- rnorm(100)
plot(x,y)
You could also use a editor like emacs with the ess-mode where one buffer can
be your script with a live R session in a second buffer.
Good luck
On 12/2/08 7:21 AM, "b g"
L.S.,
I am facing the same problem. I followed the all instructions given in help
file and some steps of previous discussion (from google search) but I could not
install package Rgraphiviz in windows. I tried even the instructions given to
install Bioconductor package then only able to insta
Hi Uwe,
Thanks for the reply.
I already tried that a couple of weeks ago, with the exactly the same
arguments you posted. So far no reply from the author.
I going to try again.
Regards,
Marcelo.
On Tue, Dec 2, 2008 at 9:42 AM, Uwe Ligges
<[EMAIL PROTECTED]>wrote:
>
>
> Marcelo Perlin wrote:
>
Thank you Uwe and Prof. Ripley.
The problem was solved. The row in question indeed have garbage data,
which possibly was truncating the number of lines read. I apologise
for the oversight.
Thank you once again.
Regards
Harsh Singhal
Bangalore, India
On Tue, Dec 2, 2008 at 2:50 PM, Prof Brian Ri
Hi,
Antje yahoo.de> writes:
>
> Hi folks,
>
> I've just discovered that the following code leads to boxplot
> (surprisingly to me).
> Can anybody explain to me why? Is this documented somewhere? I've never
> consider this option before.
>
> x <- rnorm(300)
> l <- c(rep("label1",100), rep("lab
Take a look at your dataset at around that row. Perhaps you have an
unmatched quote?
The limit on the number of rows of a data frame is far larger than 100,000
(2^31-1, but you will run out of address space on a 32-bit platform before
that - see ?"Memory-limits").
On Tue, 2 Dec 2008, Harsh
Rgraphviz is a Bioconductor package, so please ask on the apporpriate
mailing list (see the posting guide).
On Tue, 2 Dec 2008, Daren Tan wrote:
Hi,
I have problem loading Rgraphviz. Following the instructions specified by the
README in Rgraphviz_1.20.3.tar.gz didn't help either.
o. set th
Dear Ben,
It seems clear that the problem is in loading the tcltk package. To confirm
that, does library(tcltk) fail?
I'm not very knowledgeable about Mac OS X. I do have a MacBook, with OS X
10.5, however, and encountered no problems installing and running the Rcmdr
package. Have you seen the R
They are obviously growing desperate - I have now been asked to write
5 articles!
Hadley
--
http://had.co.nz/
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/po
ppaarrkk wrote:
I have two vectors, a and b. b is a text file. I want to find in b those
elements of a which occur at the beginning of the line in b. I have the
following code, but it only returns a value for the first value in a, but I
want both. Any ideas please.
a = c(2,3)
b = NULL
b[1] =
> I need to solve a equation like this :
>
> a = b/(1+x) + c/(1+x)^2 + d/(1+x)^3
>
> where a,b,c,d are known constant. Is there any R-way to do that?
Multiplying this expression with (1+x)^3 leads to a polynomial equation.
I would certainly recommend the 'PolynomF' package here:
# instal
Hi.
I am looking for a function for left-truncated data.
I have one data set with 2 variables (Hours~Yrs_Ed).
I already left-censored the data at 200 and left-truncated it at the same
spot, so that I am able to make 2 estimations (one for censoring and one for
truncation).
I know how to make the
Hi,
I have problem loading Rgraphviz. Following the instructions specified by the
README in Rgraphviz_1.20.3.tar.gz didn't help either.
o. set the following Windows environment variables accordingly
(control panel -> systems -> Advanced -> Environment Variables ):
(a) create new user
Hi folks,
I've just discovered that the following code leads to boxplot (surprisingly to
me).
Can anybody explain to me why? Is this documented somewhere? I've never
consider this option before.
x <- rnorm(300)
l <- c(rep("label1",100), rep("label2",50), rep("label3",150))
df <- data.frame(a
Marcelo Perlin wrote:
Hi Guys,
Recently I wrote a package for dealing with Markov Switching Regressions in
R and it is included in the Rmetrics project.
https://r-forge.r-project.org/projects/rmetrics/
Everything works fine when I use it in computer.
But, the package depends on the use of o
Hi there,
I know, this question is not directly an R-help question but probably someone
can give me a hint how to deal with the following problem.
I have a vector with file/folder names and want to filter for all entries which
have 6 numbers in a row and nothing else.
folders <- c("folder1
"Hans W. Borchers" <[EMAIL PROTECTED]> wrote:
> Why not use one of the global optimizers in R, for instance 'DEoptim', and
> then apply optim() to find the last six decimals? I am relatively sure that
> the Differential Evolution operator has a better chance to come near a
> global optimum than a
You could also look at littler
Then you could use, under Linux,
#!PTHTOR/r
x <- read.csv(...)
lm(...)
.
.
.
Rainer
On Tue, Dec 2, 2008 at 4:07 PM, Blanchette, Marco
<[EMAIL PROTECTED]> wrote:
> Try
>
>> source('myFirstScript.R')
>
> Where myFirstScript.R as the following line
>
> x <- rnorm(10
Harsh wrote:
Hello,
I am trying to read a dataset with 100,000 rows and around 365 columns
into R, using read.table/read.csv.
In Windows XP, with R 32 bit, I am able to read only 15266 rows and
not more than that.
I tried the same in R running in Ubuntu and it does the same and reads
only 15266
Dear all,
I just started to use the snow package to send multiple jobs on our cluster
using MPI and the Rpmi package as the communication method.
However, the Rmpi package have been behaving strangely. When I try to detach
the Rmpi package I get the following error message:
> library(Rmpi)
> d
Try this:
> folders <- c("folder1", "f2", "F234562", "12345678", "234567", "912345",
> "333")
> grep("^[0-9]{6}$", folders, value = TRUE)
[1] "234567" "912345"
On Tue, Dec 2, 2008 at 10:32 AM, Antje <[EMAIL PROTECTED]> wrote:
> Hi there,
>
> I know, this question is not directly an R-help quest
Dear list,
I've written a small utility function to add arbitrary legend(s) to a
lattice graph (or a combination of them), much like the legend
function of base graphics. I though perhaps it could be useful to
someone else, or improved by suggestions. I understand this goes
against the l
Hi Gabor,
it works! Thank you very much! But I still don't understand the difference
between [0-9] and [:digit:]...
Ciao,
Antje
Gabor Grothendieck schrieb:
Try this:
folders <- c("folder1", "f2", "F234562", "12345678", "234567", "912345", "333")
grep("^[0-9]{6}$", folders, value = TRUE)
Try this:
"^[[:digit:]]{6}$"
On Tue, Dec 2, 2008 at 10:48 AM, Antje <[EMAIL PROTECTED]> wrote:
> Hi Gabor,
>
> it works! Thank you very much! But I still don't understand the difference
> between [0-9] and [:digit:]...
>
> Ciao,
> Antje
>
>
> Gabor Grothendieck schrieb:
>>
>> Try this:
>>
>>> fol
Others have shown ways to remove your sample from the population, but this may
be doing things the long way. If you just want several samples from the same
population without overlap between any of the samples, just take one sample of
size equal to the sum of the individual sample sizes, random
Try this:
setdiff(pop, sample(pop, 2))
On Mon, Dec 1, 2008 at 5:16 PM, Hamid Hamid <[EMAIL PROTECTED]> wrote:
> Dear All,
> I am trying to build a program which will take repeated samples (w/o
> replacement) from a population of values. The interesting catch is that I
> would like the sample va
On Tue, Dec 2, 2008 at 9:48 AM, Antje <[EMAIL PROTECTED]> wrote:
> Hi Gabor,
>
> it works! Thank you very much! But I still don't understand the difference
> between [0-9] and [:digit:]...
You might find this site helpful:
http://regexp.resource.googlepages.com/analyzer.html
Copy in your attempt
Dear useRs,
Here is a weird behavior of transform function:
mtcars1<-matcars
transform(mtcars1,t1=3,t2=4)
Error in data.frame(`_data`, e[!matched]) :
arguments imply differing number of rows: 32, 1
instead, this works:
mtcars1$t1<-0
transform(mtcars1,t1=3,t2=4)
also works if applied
2008/12/2 b g <[EMAIL PROTECTED]>:
>
> Since I'm a SAS programmer, I'm used to creating command files in an editor
> for submission later. Is there a way to do this in R? I'd need to retain an
> ouput listing and a log to check for errors.
You probably want R CMD BATCH from a command-line. Fo
Antje wrote:
Hi Gabor,
it works! Thank you very much! But I still don't understand the
difference between [0-9] and [:digit:]...
If all else fails, read the help, here ?regex. Both [0-9] and
[[:digit:]] are character classes of digits, but the first contains only
arabic numerals. In some
Try:
cbind(mtcars, t1 = 3, t2 = 4)
On Tue, Dec 2, 2008 at 11:14 AM, Vitalie Spinu <[EMAIL PROTECTED]> wrote:
> Dear useRs,
>
> Here is a weird behavior of transform function:
>
> mtcars1<-matcars
> transform(mtcars1,t1=3,t2=4)
> Error in data.frame(`_data`, e[!matched]) :
> arguments imply dif
On Tue, Dec 2, 2008 at 11:18 AM, Barry Rowlingson
<[EMAIL PROTECTED]> wrote:
> 2008/12/2 b g <[EMAIL PROTECTED]>:
>>
>> Since I'm a SAS programmer, I'm used to creating command files in an editor
>> for submission later. Is there a way to do this in R? I'd need to retain
>> an ouput listing and
As the help page says
If some of the values are not vectors of the appropriate length,
you deserve whatever you get!
So you can use
mtcars1 <- mtcars
mtcars1[c("t1", "t2")] <- cbind(rep(3,32), rep(4, 32))
or even
mtcars1 <- transform(mtcars, t1=rep(3, 32), t2=rep(4, 32))
Vitalie Sp
[apologies if this appears twice]
Hi,
I have a situation where I have a set of pairs of X & Y variables for
each of which I have a (fairly) well-defined PDF. The PDF(x_i) 's and
PDF(y_i)'s are unfortunately often rather non-Gaussian although most
of the time not multi-modal.
For these data (est
You can get LoadLibrary failures if a dependency is not found (e.g. a required
dll or shlib used by the dll that you are explicitly loading). So make sure
that all necessary dependencies (e.g. graphviz libs) are available in the
system path.
Rory
Rory Winston
RBS Global Banking & Markets
Offi
On Tue, 02 Dec 2008 17:37:44 +0100, Prof Brian Ripley
<[EMAIL PROTECTED]> wrote:
As the help page says
If some of the values are not vectors of the appropriate length,
you deserve whatever you get!
So you can use
mtcars1 <- mtcars
mtcars1[c("t1", "t2")] <- cbind(rep(3,32), rep(4
On 12/2/08, baptiste auguie <[EMAIL PROTECTED]> wrote:
> Dear list,
>
>
> I've written a small utility function to add arbitrary legend(s) to a
> lattice graph (or a combination of them), much like the legend function of
> base graphics. I though perhaps it could be useful to someone else, or
> im
I wonder if you are using this term in its correct technical sense.
A linear functional relationship is
V = a + bU
X = U + e
Y = V + f
e and f are random errors (often but not necessarily independent) with
distributions possibly depending on U and V respectively.
and pairs from (X,Y) are obse
I wonder if you are using this term in its correct technical sense.
A linear functional relationship is
V = a + bU
X = U + e
Y = V + f
e and f are random errors (often but not necessarily independent) with
distributions possibly depending on U and V respectively.
and pairs from (X,Y) are obse
At 11:04 AM -0500 12/1/08, John Fox wrote:
Dear Bill,
Thanks for pointing out that this functionality is already in the psych
package. Shouldn't factor.residuals() avoid this computation for oblique
rotations?
John,
Good suggestion. I will add that in the next revision.
Bill
Regards,
Jo
Hi, I was wondering if there was some sort of package or function that
calculated asymmetric confidence intervals for small proportions. I thought of
both the epicalc and epitools package, but I am hoping to find something where
you can just plug in a standard error and point estimate and it wi
Roger,
I had success installing RCurl per Patrick's instructions.
I also got DAVIDQuery to install and run correctly but I had to run R CMD
INSTALL DAVIDQuery from a cygwin terminal after 1) I replaced cygwin's
make program v3.81 with the same version patched for windows from the
mingw project. M
Thanks a lot again to all of you!!!
Antje
Antje schrieb:
Hi Gabor,
it works! Thank you very much! But I still don't understand the
difference between [0-9] and [:digit:]...
Ciao,
Antje
Gabor Grothendieck schrieb:
Try this:
folders <- c("folder1", "f2", "F234562", "12345678", "234567",
Dear R-users,
I'm using lmer to fit two-level logistic models and I'm interested in
predicted probabilities that I get in this way (using "fitted"):
glm1 = lmer(XY$T1~X1 + X2 + X3 + (1|Cind), family=binomial) #estimation of a
two-level logit model
fit1=fitted(glm1) # I get the fitted line
suppose something like probability(passing test) is driven by
1. fixed effects -- sex
2. district effects - district funding
3. school effects - neighborhood income, racial composition, % two parent
families, ...
4. class effects - teacher quality measurement,
5. individual rand
Look at the sand package, which is available at
http://www.csm.ornl.gov/esh/statoed/
and the NADA package, which is available from CRAN. One or both may have
items of interest.
Tom
Zita wrote:
>
> Hi.
>
> I am looking for a function for left-truncated data.
> I have one data set with 2
Thanks for the reply!
On Tue, Dec 2, 2008 at 6:34 PM, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> I wonder if you are using this term in its correct technical sense.
> A linear functional relationship is
>
> V = a + bU
> X = U + e
> Y = V + f
>
> e and f are random errors (often but not necessa
Sarah E. McCormick wrote:
Hi, I was wondering if there was some sort of package or function that
calculated asymmetric confidence intervals for small proportions. I thought of
both the epicalc and epitools package, but I am hoping to find something where
you can just plug in a standard error
Richard
There is much more information that you need to provide before a
thoughtful answer can be provided. Maybe you can describe the structure
of your data, your outcome variable, etc. There is a vignette in the
lmer package called 'Implementation' that will show you some methods for
model fitti
List,
I am using the 'tisPlot' function in Jeff Hallman's excellent tis package
and was hoping that someone could spare me from having to dig into the
code of his 'tisPlot' function. So far as I can tell, the preferred
method of controlling the plotting of the x-axis is using the 'xTickFreq'
Hi, I would like to check which rows of 'types.prev' matrix pop up in
'types', following R in-built procedure. I tried 'match' function but it
works in case of the one dimensional vectors.
Will appreciate any suggestions.
best, robert
> types
edateK
[1,] 20060819 12.5
[2,] 200
The development version of zoo has a tis to zoo conversion function
that would allow you to plot your data using plot.zoo and xplot.zoo
in which case you can use classic graphics axis function (for plot.zoo)
or lattice facilities (with xyplot.zoo):
source("http://r-forge.r-project.org/plugins/scms
Sorry that should have been:
library(zoo)
z <- as.zoo(dat, class = "yearmon")
plot(z)
xyplot(z)
On Tue, Dec 2, 2008 at 2:35 PM, Gabor Grothendieck
<[EMAIL PROTECTED]> wrote:
> The development version of zoo has a tis to zoo conversion function
> that would allow you to plot your data using plot.z
Prof Brian Ripley wrote:
As the help page says
If some of the values are not vectors of the appropriate length,
you deserve whatever you get!
Yes (did I write that?). It is a bit annoying with things that almost
work, though.
[snip]
I often need to use this for creating new var
Thanks Harold, I will review the lmer vignette again.
On Tue, Dec 2, 2008 at 2:12 PM, Doran, Harold <[EMAIL PROTECTED]> wrote:
> Richard
>
> There is much more information that you need to provide before a
> thoughtful answer can be provided. Maybe you can describe the structure
> of your data,
Tena koe Robert
You could first apply paste
apply(types, 1, paste, collapse=':')
and then match.
HTH
Peter Alspach
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of threshold
> Sent: Wednesday, 3 December 2008 8:32 a.m.
> To: r-help@r-proje
Michael Styer wrote:
>> So my next question is, does anyone have any thoughts about how significant
>> a project it would be to compile R for 64-bit windows (using, e.g., the
>> Portland Group
>> compiler)?
Based on our experiences at REvolution Computing, it's quite a
significant project. We're
I am trying to suppress the tick labels on the x-axis of the following:
barchart(richness[Wood=="V"]~Sample[Wood=="V"])
I have tried col.axis="white"
I have tried removing the axis all together with axes=FALSE
I have tried xaxt="n"
I have also tried labels=c"(label1", "label2") to replace the
Isn't this a special case of structural equation modeling, handled
by the 'sem' package?
Spencer
Jarle Brinchmann wrote:
Thanks for the reply!
On Tue, Dec 2, 2008 at 6:34 PM, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
I wonder if you are using this term in its correct technic
Yes I think so if the errors were normally distributed. Unfortunately
I'm far from that but the combination of sem & its bootstrap is a good
way to deal with it in the normal case.
I must admit as a non-statistician I'm a not 100% sure what the
difference (if there is one) between a linear functio
I have been using rpart to get trees. May someone help me with how I can do
some logistic regression modelling to derive odds ratios on the branches.
Also help on how to derive estimates on prediction error.
--
Kind Regards,
Jestinah
Home Address:
Brusselsepoortstraat 85
9000
Gent
Belgium
Hi all,
I'm dealing with dates in R (2.7.2), but some basic operations raise a
warning.
Incompatible methods ("+.Date", "Ops.difftime") for "+"
I saw this topic in this mailing list, but I do not understand what to
do...
cf. https://stat.ethz.ch/pipermail/r-help/2008-June/165842.html
Do I
In case anyone is still reading this thread, I want to add this:
In a current problem (a data-shy five-parameter nonlinear
optimization), I found "nlminb" markedly more reliable than
"optim" with method "L-BFGS-B". In reviewing the fit I made, I
found that "optim" only came close to its own minimum
Dear R People:
In the DASL library, there is a story about hot dogs.
Here are the data:
Beef186 495
Beef181 477
Beef176 425
Beef149 322
Beef184 482
Beef190 587
Beef158 370
Beef139 322
Beef175 479
Beef148 375
Beef1
Hadley,
I don't know if I am doing something wrong or if it is ggplot please
see the two graphs at the bottom of the page (code).
melt.nut <- (structure(list(RiverMile = c(119L, 119L, 119L, 119L, 119L, 119L,
119L, 119L, 119L, 148L, 148L, 148L, 148L, 148L, 148L, 148L, 179L,
179L, 179L, 179L, 179L,
> [R] Suppressing tick labels?
> Graham Smith myotisone at gmail.com
> Tue Dec 2 21:26:03 CET 2008
>
> I am trying to suppress the tick labels on the x-axis of the
following:
>
> barchart(richness[Wood=="V"]~Sample[Wood=="V"])
Use the scales= argument. E.g.
d<-data.frame(richness=log(1:12),
Dear Erin,
What you characterize as "snaky" might also be called bimodal and
short-tailed. Try, e.g., plot(density(residuals(dog1.aov))) to see the
bimodality more clearly. A Box-Cox transformation can correct skewness, but
won't help here.
I hope this helps,
John
--
[EMAIL PROTECTED] writes:
> List,
>
> I am using the 'tisPlot' function in Jeff Hallman's excellent tis package
> and was hoping that someone could spare me from having to dig into the
> code of his 'tisPlot' function. So far as I can tell, the preferred
> method of controlling the plotting of
Is there a good and concise way of making simultaneous plots that are
identical, but directed to different devices?
I'm writing an R-script that produces a pdf file. I would really like to
check visually whether the pdf file shows what I expect. So I would like the
same commands to produce a plot
Thanks for the report,
the problem boils down to the call of "methods:::bind_activation(TRUE)"
in one of the depended package.
I can reproduce the problem with
> methods:::bind_activation(TRUE)
> dfr <- data.frame(matrix(0, nrow = 1 , ncol = 1000))
> dfr2 <- is.na(dfr)
I will forward you remar
Hi all,
I have a data frame with "clustered" rows as follows:
Cu1 x1 y1 z1 ...
Cu1 x2 y2 z2 ...
Cu1 x3 y3 z3 ... # end of first cluster Cu1
Cu2 x4 y4 z4 ...
Cu2 x5 y5 z5
Cu2 ... # end of second cluster Cu2
Cu3 ...
...
"cluster"-size is 3 in the example above (rows making up a cl
I would like to rotate the axis labels 45 deg.
--
Stephen Sefick
Let's not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods. We are mammals, and have not exhausted the
annoying little proble
'str' is your friend. Look at the results:
> alpha2=as.Date("2008-12-21")
> alpha1=as.Date("2008-12-21")-as.Date("2008-10-26")
> alpha1
Time difference of 56 days
> alpha2
[1] "2008-12-21"
> alpha1+alpha2
Time difference of 14290 days
Warning message:
Incompatible methods ("Ops.difftime", "+.Date
If you want to do the addition, 'unclass' the variable:
> alpha2+4
[1] "2008-12-25"
> alpha2 + unclass(alpha1)
[1] "2009-02-15"
>
On Tue, Dec 2, 2008 at 4:10 PM, Christophe Dutang <[EMAIL PROTECTED]> wrote:
> Hi all,
>
> I'm dealing with dates in R (2.7.2), but some basic operations raise a
> wa
Hi Stephen,
I think you will need to clarify what your problem is with the second plot.
HTH,
Thierry
-Oorspronkelijk bericht-
Van: [EMAIL PROTECTED] namens stephen sefick
Verzonden: di 2-12-2008 22:52
Aan: hadley wickham; R-help
Onderwerp: [R] ggplot2 facet_wrap problem
Hadley,
I don
Hi Stephen,
Have a look at the ggplot2 book. You will find an example on p. 6 of
http://had.co.nz/ggplot2/book/polishing.pdf
HTH,
Thierry
-Oorspronkelijk bericht-
Van: [EMAIL PROTECTED] namens stephen sefick
Verzonden: wo 3-12-2008 0:31
Aan: R-help
Onderwerp: [R] ggplot2 45deg axis la
Hi Stephen,
Have a look at the "polishing your plot for publication" chapter of
the ggplot2 book.
Regards,
Hadley
On Tue, Dec 2, 2008 at 5:31 PM, stephen sefick <[EMAIL PROTECTED]> wrote:
> I would like to rotate the axis labels 45 deg.
>
> --
> Stephen Sefick
>
> Let's not spend our time and r
Not sure exactly what you mean by 'sample' since you did not provide
an example of the expected output, or input data that could be used.
Here is an example of taking one sample from each cluster:
> df <- data.frame(id=paste("C", rep(1:5, each=3), sep=''), data=1:15)
> # sample 1 from each cluster
If you look at the TSS graph in the faceted example and then look at
the plot of just the GPP vs. TSS. They are different graphs all
together. The one that is not faceted is correct.
On Tue, Dec 2, 2008 at 6:36 PM, ONKELINX, Thierry
<[EMAIL PROTECTED]> wrote:
> Hi Stephen,
>
> I think you will n
Hello, My question likely got buried so I am reposting it in the hopes that
someone has an answer. I have thought more about the question and modified my
question. I hope tha
my specific question is:
I am attempting to create a bootstrap procedure for a finite sample using the
theory of Rao a
A few questions about maps...
(1) How can I find a listing of the internal data sets that map() from the maps
library contains?
For example, "usa", "county", "state", "nz" all work. Are there any others?
(2) Is there an easier, more generalized way to produce this
(http://www.ai.rug.nl/~hedder
On Tue, Dec 2, 2008 at 6:21 PM, Avram Aelony <[EMAIL PROTECTED]> wrote:
> A few questions about maps...
>
> (1) How can I find a listing of the internal data sets that map() from the
> maps library contains?
> For example, "usa", "county", "state", "nz" all work. Are there any others?
help(packa
I have a set of modeled climate data recorded at irregular intervals.
The format of the data is such that there are monthly measurements for
the years 2000, 2020, 2050, 2080, etc. Therefore I have 12 regular
records, a skip of some number of years, then 12 more monthly records,
another skip
Have a look at the zoo package. There are three vignettes (pdf documents)
included with the package that give many examples of its use. Also
see ?read.zoo, ?plot.zoo and ?xyplot.zoo
You will need something like:
library(zoo)
z <- read.zoo("myfile", ...whatever...)
plot(z)
On Tue, Dec 2, 2008
I have a series of csv files in several folders. All begin with a 7 digit
number and end with the letter "E" (eg. 0726016E.csv).
I want to be able to read a file in to R, take some of the data out of it
and store it in a matrix, then move on to the next file and do the same
thing.
I was planning
To get the file names in the current directory try:
list.files(pattern="[[:digit:]]{7}E")
On Tue, Dec 2, 2008 at 4:11 PM, Steven Kennedy <[EMAIL PROTECTED]> wrote:
>
> I have a series of csv files in several folders. All begin with a 7 digit
> number and end with the letter "E" (eg. 0726016E.csv)
1 - 100 of 119 matches
Mail list logo
