I assume this is (or was) a specification issue. I think write.dbf uses the
shapefile library (C not R library) so it applies to the use of shapefiles and
just happens to have been included in the foreign package because it has a
generic usefullness. (Is that a word?)
Since I very rarely care
Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Mulholland, Tom
Sent: Friday, 19 May 2006 11:48 AM
To: R-Help (E-mail)
Subject: [R] Converting character strings to numeric
I assume that I have missed something fundamental and that it
is there in front of me
.
Tom
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Friday, 19 May 2006 3:03 PM
To: Mulholland, Tom
Cc: R-Help (E-mail)
Subject: Re: [R] Converting character strings to numeric
On Fri, 19 May 2006, Mulholland, Tom wrote:
After replies off the list
I assume that I have missed something fundamental and that it is there in front
of me in An Introduction to R, but I need someone to point me in the right
direction.
x1 - 1159 1129 1124 -5 -0.44 -1.52
x2 - c(1159,1129,1124,-5,-0.44,-1.52)
x3 - unlist(strsplit(x1, ))
str(x2)
chr [1:6]
You might also consider looking at the R-sig-Geo list which has lots of
discussion about issues relating to file formats and the best ways to get data
in and out the various packages that are used.
Tom
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of
Well I downloaded the data using the link in your message which suggests that
the code is right. I don't have its loaded (I assume it's from the irregular
time series package) so I can't test the code as you have it.
Have you tried checking to see it it is an issue with the internet (your
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of
[EMAIL PROTECTED]
Sent: Wednesday, 14 September 2005 8:13 AM
To: vittorio
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Reading data from a serial port
snip
Now all this of course is written in
I searched the help for cosine distance and this was the first hit
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/3946.html
Tom
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Raymond K Pon
Sent: Tuesday, 13 September 2005 3:48 AM
To:
I cannot state this with the certainty that others might, but the Rd format is
a text format. If you want to produce something else then you need to choose an
alternative method. For instance, 1.4 of Writing R Extensions notes that
Documents in 'inst/doc' can be in arbitrary format, however we
For some reason (probably that our organisation has blocked the site) I could
not see the original articles that prompted the post. I however immediately
assumed that this was precipitated by Tufte and his comments about PowerPoint
(I recall seeing a good example of PowerPoint on his site)
I incorrectly relied upon my memory
...
and that
John Fox did something
http://ils.unc.edu/~jfox/powerpoint/introduction.html that I
enjoyed reading.
The work is that of Jackson Fox
Tom
__
R-help@stat.math.ethz.ch mailing list
I note that the axis help seems to refer to padj. After playing around it is
obvious that I don't know what is meant by this argument, so maybe I'm doing
something wrong. My practical soultion is
plot(1:50,axes = FALSE,ylab = )
axis(2,at = 1:50,labels = rep(,50),las = 2,padj = 0)
require( car )
set.seed(12345)
nn - sample( c( 2, 4 ), size=50, replace=TRUE )
rr - recode( nn, 2='TWO';4='FOUR' )
table( rr, exclude=NULL )
ss - recode( nn, 2='Num2';4='Num4' ) # Doesn't work as expected
table( ss, exclude=NULL )
ss - recode( nn, 2='Num2';4='Num4', TRUE ) #?
table( ss,
If this is a question about plotting 2 series on a plot then in general terms
you can plot 2 series as per this example in the list
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/5463.html
If it is about the use of rgev you should tell us the name of the package.
If you read the posting
Search the archives for zoom and you will find plenty of answers on this
question.
RSiteSearch(zoom)
Tom
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Henrik Andersson
Sent: Tuesday, 26 July 2005 4:16 PM
To: r-help@stat.math.ethz.ch
Subject: [R]
doubleEm - function(p1,p2,p3,p4){return(p1 *p1,p2 * p2, p3 * p3, p4 * p4)}
suppressWarnings(doubleEm(1,2,3,4))
[[1]]
[1] 1
[[2]]
[1] 4
[[3]]
[1] 9
[[4]]
[1] 16
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Paul Roebuck
Sent: Wednesday, 20 July 2005
You will find previous discussion about pairs in this list. There are limits on
what has been included within the functionality, but you can write your own
panel functions. Here's a starting point.
x - runif(100)
dim(x) - c(20,5)
panel.cor1 - function(x,y,...){
mycor- cor(x,y)
if (mycor
It's not so much a problem, as not working the way you expected.
cluster:::plot.partition is used to do the plotting. If you look at the code
for this you can see the difficulty in putting every possible permutation into
the code. If for example you want the silhouette plot to be red using col
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Henrik Bengtsson
Sent: Tuesday, 28 June 2005 2:54 AM
To: Spencer Graves
Cc: r-help@stat.math.ethz.ch; Dirk Eddelbuettel
Subject: Re: [R] How to convert c:\a\b to c:/a/b?
... snipped
Thus, you
Have you found the file Using Rmetasim? In windows you can access this file
by using the help and selecting browse directory. There appears to be a
reasonable amount of information here. It looks to me as if you need to work
your way through these files until you understand what is going on.
Well I skipped to the end so pardon me if I've missed something. My first
reaction was to go and look at the excellent article by Gabor in RNews 2004-1
on dates (p.32 in particular)
as.POSIXct(strptime(7/12/2001 10:32,format = %d/%m/%Y %H:%M))
[1] 2001-12-07 10:32:00 W. Australia Standard Time
It would have been helpful if you had written the code that actually showed the
issue you have rather than leaving it to us to try and reproduce.
That leaves generic options. If the screen is too small for you then try
plotting it to postscript or PDF device and setting the paper size to A3 or
What makes you think that there is a ylab parameter?
args(map)
function (database = world, regions = ., exact = FALSE, boundary = TRUE,
interior = TRUE, projection = , parameters = NULL, orientation = NULL,
fill = FALSE, col = 1, plot = TRUE, add = FALSE, namesonly = FALSE,
xlim =
I don't think Spencer replied to your message. It wasn't addressed to you. He
was replying to a specific post on nomograms. However you might try being less
specific as a search on spatial cluster gave back
[R] mclust - clustering by spatial patterns
Without a small example to see what you are doing it is hard to respond. There
are plenty of examples in the help for legend showing placement all over the
place. So I am guessing that this might help (The only thing added from the
help was the par(xpd = TRUE)
x - 0:64/64
y -
I think you may wish to look at the plotrix package, assumming that you have
taken care of the issues involved in breaking an axis and that your plots don't
result in misleading information.
I think to use the axis break you would have to calculate your own labels and
rescale the data, as it
Well since I know nothing about this topic I have lurked so far, but here's my
two bob's worth.
Firstly I tried to make sense of Brian's initial reply. I have got no idea who
Bellman is and you have not referenced (his/her) work in a way I can access the
issues you refer to. So I assumed
Just a little bit of trivia.
The 2001 Census in Australia had a significant group of people who responded to
the question of religion with the answer Jedi or Jedi Knight. Unfortunately the
Australian Bureau of Statistics is a bit fuddy duddy about the issue as they
see it as a trivialisation
I have snippets of code that I have either taken from examples or off of the
list. I apologise to those I have stolen it from but I didn't keep the proper
references.
n - 100
xx - c(0:n, n:0)
yy - c(c(0,cumsum(rnorm(n))), rev(c(0,cumsum(rnorm(n)
plot (xx, yy, type=n, xlab=Time,
Lattice is not my forte but here goes.
?panel.bwplot specifically notes that pch, col, cex: graphical parameters
controlling the dot. If you look at the code for panel.bwplot you will see
where the colours come from in which case you can probably set up your own
colour scheme using
I recall having this problem. I think I had a version that didn't work. Did you
download the htmlhelp.exe from http://www.murdoch-sutherland.com/Rtools/ . The
path does not matter as I have it in a folder on it's own. I think that I also
had a path to the version that didn't work and I had to
There's built in FTP in windowsXP, which I assume is what they are likely to
have. If they open file explorer there's an entry called My Network Places.
The problem I had seeing if it would work, was that
ftp://cran.r-project.org/pub/R/ doesn't want to know me. Not that I was at all
sure,
Fristly when you are using a package (in this case date) put it in your email.
Dates are stored as numbers and if I recall correctly as.date will be the
number of days since sometime in 1960. As with other objects there are
generally methods that will ensure that the correct printed format will
Of course Andy meant hist.factor(f)
In particular you should note that Andy uses the table function to transform
... the data from characters to numbers
Tom
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Liaw, Andy
Sent: Thursday, 21 April 2005 12:14
Well you have not given us anything to go on really. Are there more than 94
columns? Does each column have a valid fieldname? RODBC is not guaranteed to
work in every possible scenario. If you have a look through the list you will
find there are specific limitations which are not immediately
, 21 April 2005 10:06 AM
To: Mulholland, Tom
Subject: RE: [R] Histogram
Have you tried it? hist.factor() as defined would be the
hist method for
the factor class, so hist(f) would work if f is a factor.
Andy
From: Mulholland, Tom
Of course Andy meant hist.factor(f
I guess one of the reasons that you have not had a reply is that you have not
followed the posting guide. If you give the list something to work with (a
small reproducible example)
You use the word package which in R is very precise, but which I think you are
using to describe the file you are
from and what attempts you
have made to solve your problem.
Tom
* The posting guide explains all of this much better than I do.
-Original Message-
From: Owen Buchner [mailto:[EMAIL PROTECTED]
Sent: Tuesday, 19 April 2005 7:51 PM
To: Mulholland, Tom
Subject: RE: [R] using imported
Should you be using rsiTA(tsx[,2],14). If you look at the function you will see
it is expecting just the values you want in the calculation.
If you give it a matrix it treats the whole matrix as being price data.
Tom
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL
What does str(tsx) give?
So you are feeding in something the function has no idea about.
Try rsiTA([EMAIL PROTECTED],14)
Tom
-Original Message-
From: Neuro LeSuperHéros [mailto:[EMAIL PROTECTED]
Sent: Wednesday, 20 April 2005 9:42 AM
To: Mulholland, Tom; r-help@stat.math.ethz.ch
Well I don't know how I can live with myself. I guess I can't wait for the site
to mirror itself in case someone thinks I'm yesterday's man. ;-)
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Achim Zeileis
Sent: Monday, 18 April 2005 10:56 PM
...
Since it is not clear exactly what you require, I have assumed that what you
are looking for is the minimum value for each row or column depending upon
which is airports and which is the towers.
Is this what you are looking for
x - runif(100)
dim(x) - c(5,20)
apply(x,1,function(y) which(y ==
I think this might have been my code
mapply(paste,strwrap(levels(ncdata$Chapter),18,simplify = FALSE),collapse =
\n)
Tom
-Original Message-
From: Jan P. Smit [mailto:[EMAIL PROTECTED]
Sent: Thursday, 14 April 2005 5:15 PM
To: Mulholland, Tom
Cc: r-help@stat.math.ethz.ch
Subject
This may not be the best way but in the past I think I have done something like
levels(x) - paste(strwrap(levels(x),20,prefix = ),collapse = \n)
Tom
-Original Message-
From: Jan P. Smit [mailto:[EMAIL PROTECTED]
Sent: Thursday, 14 April 2005 11:48 AM
To: r-help@stat.math.ethz.ch
I'm not sure that this answers your questions but maybe they partly help.
p. 7 in An introduction to R notes
For most purposes the user will not be concerned if the numbers in a numeric
vector
are integers, reals or even complex. Internally calculations are done as double
precision real
, say col
= cl$cluster =
a set of group numbers say 1...10 should produce 10 distinct
colours points
grouped by colour.
how to do this when you have more than 8 group colours to plot
From: Mulholland, Tom [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch,Asha Jayanthi
[EMAIL PROTECTED
See upper.tri and lower.tri.
I think that you might also look for specific packages that function using
matrices, from what I have seen these often have the capacity to ignore the
diagonal or use just the upper or lower triangle. This is not an area that I
use very much, but I have seen
Does ?na.omit help
x - kmeans(na.omit(data),centres)
of course if you have too many NAs you need to be sure that their removal does
not unduly influence the results.
Although I am a bit confused as I thought that agnes did not allow NAs. I
assume that you are running an alternative clustering
PM
To: Mulholland, Tom; r-help@stat.math.ethz.ch
Subject: RE: [R] 2d plotting and colours
Thank you.
mycols - c(brown,orange,tomato)
plot(x,col = mycols[cl$cluster])
this works. I can define distinct colours and check the graph.
the rest of the examples does not give a wide palette
I think you need to read the posting guide (see the bottom of each post made)
and once you have done this take some time to compose your message.
The issue is that I have too little information about what you have done. It
looks to me as if you are using postscript, but I am not sure if you
The error message states that you are passing a parameter called cex which has
not been used. If you look at ?plclust more closely you will see it does not
have cex parameter. However the S3 method for class hclust, plot, does?
So does this help?
hc - hclust(dist(USArrests), ave)
plot(hc,cex =
There are two issues here identifying the outliers and highlighting them.
I have only a basic grasp of both of these concepts but will give what I have
in case it helps. There appears to have been a move in the last 2 decades to
improve the concepts of what actually constitutes an outlier,
Are you talking about something other than 'update.packages'?
as the help notes
Description:
These functions can be used to automatically compare the version
numbers of installed packages with the newest available version on
CRAN and update outdated packages on the fly.
I'm afraid you have lost me. What is it that you want that reordering the
formula does not achieve.
bwplot(yield ~ year | site, data = barley) has sites next to each other. If the
lattice structure is your issue (it appears you wish to remove the structure
and replace it with a wider space)
Do you have a really good reason to be using 1.9.1. If not then just keep using
1.9.0. Did you check what changes were made in the 1.9 upgrade. Often you will
find useful information about this type of issue in the change log.
You have not told us anything about the machine you are using and
I've used rimage to read in graphics files (jpeg.) If I recall correctly, I
think I had to install some libraries. What I can't recall is if it was any
faster than pixmap, as I was mainly concerned with the file format and the
forensic image processing possibilities.
Tom
-Original
It seems to me that you are trying to do too much at a time. Firstly I think it
would be a good idea to get you code working before you try and make a package.
Some possibilities are that you write somethin meaningful rather than the first
thing that pops into your head. What sort of output are
I'm not sure how to best explain what I am after but here goes. I have a data
frame with 2 geographical factors. One is the major region the other is the
component regions.
I am trying to process all the regions at the same time without using for. So
I need (think, I do) a list of matrices
reuseable. I guess I could pass a formula and work out a better
method of subsetting dimensions (where certain factor levels are not used. But
maybe someone has an elegant method they could share.
Tom
-Original Message-
From: Mulholland, Tom
Sent: Tuesday, 22 March 2005 1:35 PM
To: R
type ?NaN and the help will tell you.
-Original Message-
From: Brett Stansfield [mailto:[EMAIL PROTECTED]
Sent: Monday, 21 March 2005 12:10 PM
To: R help (E-mail)
Subject: [R] NaN
Dear R
What does NaN mean?
I recently did a correlation on a batch of data for some
reason it
I'm not sure I am answereing your question, but here goes
# Create a vector with 2400 items
x - runif(2400)
# Create a 600 by 4 matrix
y - matrix(x,ncol = 4)
#If you needed the matrix to be done row by row
#y - matrix(x,ncol = 4,byrow = T)
# Extract values from a particular part of the matrix
type ?cor and read the help file
In particular read it till you find the entry that tells you how to deal with
missing observations. You need to do this because as you get more competent you
will undoubdtedly come across other issues, so learning the format of the help
and realising why you
If it's specifically R that I am looking for I use the following link
http://finzi.psych.upenn.edu/search.html
If I were going to use google I would go to the advanced page and insert
finzi.psych.upenn.edu into the Domain field so as to restrict the search. As
you have found out the letter R
Sorry about the wrong link. The archives are at
http://sourceforge.net/mailarchive/forum.php?forum=rpy-list
-Original Message-
From: Mike R [mailto:[EMAIL PROTECTED]
Sent: Friday, 11 March 2005 2:44 PM
To: r-help@stat.math.ethz.ch
Subject: [R] howto: plot location, web search
I think I recall seeing a limited capability in the PBSmapping package.
Tom
-Original Message-
From: yyan liu [mailto:[EMAIL PROTECTED]
Sent: Wednesday, 9 March 2005 1:20 PM
To: r-help@stat.math.ethz.ch
Subject: [R] from long/lat to UTM
Hi:
Is there any function in R which
?axis
where you will find
See Also:
'axTicks' returns the axis tick locations corresponding to
'at=NULL'; 'pretty' is more flexible for computing pretty tick
coordinates and does _not_ depend on (nor adapt to) the coordinate
system in use.
Tom
-Original Message-
Does strsplit fit what you want?
Tom
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
Sent: Thursday, 24 February 2005 11:46 AM
To: r-help@stat.math.ethz.ch
Subject: [R] a simple question
Dear All,
I need to separate one column which is in a format like,
x - c(1,2,3,4,5,8,9,10,11,12,15,16,17,18,19,22,23,24,33,34,35)
require(cluster)
pam(x,5)
Medoids:
[,1]
[1,]3
[2,] 10
[3,] 17
[4,] 23
[5,] 34
Clustering vector:
[1] 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 5 5 5
Objective function:
build swap
1.285714 1.047619
Available
barplot(matrix(c(x,y),ncol = 2),beside=T)
Does this help
?barplot notes
height: either a vector or matrix of values describing the bars which
make up the plot. If 'height' is a vector, the plot consists
of a sequence of rectangular bars with heights given by the
I think you may need to check the list for posts relating to memory, large
datasets, etc
http://cran.r-project.org/search.html
It would help those on the list, if you gave more details about the machine you
are using. Given the vintage of R you areusing is considered quite old, there
is the
x - rep(seq(1,10),10)
y - rep(seq(1,10),each = 10)
z - matrix(100 + x*0.5 + y*0.5,ncol=10,byrow=T)
contour(z)
or back to where you were originally
x - 1:10
y - 1:10
# z - 100:110
z - 1:100
contour(matrix(z,ncol=10))
read the help on contour it states that
z: a matrix containing the values to
My first thought was that all it looked a bit complicated for something that
should be straightforward.
I created a file called t.txt. I worked out the way I would have done it and
then I tested to see which was fastest. One little hiccup is that the two
objects are not identical and I though
I'm sure others with more experience will answer this, but for what it is worth
my experience suggests that memory issues are more often with the user and not
the machine. I don't use Linux so I can't make specific comments about the
capacity of your machine. However it appears that there is
I am immediately reminded of something I read which goes
A sufficiently trained statistician can read the vagaries of a Q-Q plot like a
sharman can read a chicken's entrails, with a similar recourse to scientific
principles. Interpreting Q-Q plots is more a visceral than an intellectual
I am probably just displaying my ignorance, but I have obviously managed to
miss exactly what you are referring to.
Firstly I have to thank you for making me look closer at the article. I had
done so but I had obviously skipped over the Comparison table, which I would
have found useful in the
In general if you are using a package, it helps to identify the package. I
assume it is the tree package.
Is this what you are after (from the example on the ?tree page)
attr(ir.tr,ylevels)
[1] setosa versicolor virginica
Tom
-Original Message-
From: NICOLAS DEIG
I think Gabor gave you the answer, even if you didn't see it
x - rep(c(02/27/92, 02/27/92, 01/14/92, 02/28/92, 02/01/92),each = 5)
z - data.frame(Date = strptime(x, %m/%d/%y))
z$Value - trunc(runif(25) * 5) + 1
z$State - c(A,B,C,D,E)[z$Value]
z - z[order(z$Date),]
z
z[aggregate(1:nrow(z),
This seems to work
toPOSIX - function(x){
y - x - as.numeric(ISOdate(2005,1,1))
z - ISOdate(2005,1,1) + y
return(z)
}
test - as.numeric(ISOdate(2005,3,1) )
toPOSIX(test)
But whether one should be doing this I don't know. There are certainly
functions that play aorund with the POSIX
Well I am not sure that can call a single figure a cluster. Sure it's not near
the others but how can you conceptually measure it's cluster properties. It
seems reasonable that there has to be some form of doubt about it.
Back to that Google search hit number 3 www.stat.ncu.edu.tw/teacher/
What makes you trust any software?
There are some obvious points. First of all the code is open so if you know
enough you can actually read the code and make sure it does what you want.
Secondly you can replicate a process using two pieces of software and compare
the results. You can check
The answers to your questions are in the text below. However these questions
are generally answered in the base documentation. You might try going through
the Keywords by Topic which on the windows system at least) is accessed by
going through the html help entry on the help menu. Once your
?points has this example
plot(-4:4, -4:4, type = n)# setting up coord. system
points(rnorm(200), rnorm(200), col = red)
points(rnorm(100)/2, rnorm(100)/2, col = blue, cex = 1.5)
In general you might want to check out the keyword section of the help, in
particular the Graphics section which has
I hesitate to add this comment since it either completely confuses people or
they take to it very quickly.
The data that you are using is mostly categorical. I expect that tables will
have been used in the past and that to acertain extent the graphics are
suppossed to help with getting a quick
temp - table(rpois(100,5),rpois(100,1))
temp
temp[temp[,1] 1,]
temp[temp[,1] == 1,]
temp[temp[,1] == 2,]
temp[temp[,1] == 3,]
Of course you have not made it obvious if you are using the term table to mean
any output that looks like a table, but in general terms you will find
information in the
I use windows, so I can't state that the same is true of other operating
systems.
There is a file which is in the etc folder of the R directory called Rconsole.
part of mine looks like
## Font.
# Please use only fixed width font.
# If font=FixedFont the system fixed font is used; in this case
Try searching this list for delphi. This was the first post I found. You may
find it helpful.
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/19154.html
Tom
-Original Message-
From: YiYao_Jiang [mailto:[EMAIL PROTECTED]
Sent: Wednesday, 19 January 2005 10:12 AM
To:
(-0.25,1.25),ylim = c(-0.25,1.25))
points(y1,z1)
rect(0,0.1,0.5,0.9)
axis(2,at = seq(-0.25,1.25,length = 5),labels = seq(0,20, length = 5))
Tom
-Original Message-
From: Costas Vorlow [mailto:[EMAIL PROTECTED]
Sent: Wednesday, 12 January 2005 6:22 PM
To: Mulholland, Tom
Subject: Re: [R
I have often noted that statistics can't prove a damn thing, but they can be
really useful in disproving something. Having spent most of 80s and half of
the 90s with the Australian Bureau of Statistics to find out how you collect
these numbers, I am disconcerted at the apparent disregard for
You might find breakpoints in strucchange helpful
Tom
-Original Message-
From: Dr Carbon [mailto:[EMAIL PROTECTED]
Sent: Thursday, 13 January 2005 6:19 AM
To: r-help@stat.math.ethz.ch
Subject: [R] Finding seasonal peaks in a time series
I have a seasonal time series. I
You stand more chance if you do it yourself
https://stat.ethz.ch/mailman/listinfo/r-help
-Original Message-
From: Kevin Ita [mailto:[EMAIL PROTECTED]
Sent: Thursday, 13 January 2005 2:25 AM
To: R-help@stat.math.ethz.ch
Subject: [R] Please unsubscribe me from you list
Please
Sorry I didn't read the question properly. Please disregard, my mind was
elsewhere.
Tom
-Original Message-
From: Mulholland, Tom
Sent: Thursday, 13 January 2005 10:52 AM
To: Dr Carbon; r-help@stat.math.ethz.ch
Subject: RE: [R] Finding seasonal peaks in a time series
You
There are indeed speed advantages in using colSums etc. However the
disadvantage is that the newbie doesn't always find the power inherent in the
apply, sapply, tapply and mapply. For many things that I do, the speed is the
least of my worries; although I take the point that using apply for
apply(y,2,mean)
Is this what you are after. If it is I would suggest that you look at the
examples not just for this but for what I call the apply family sapply, tapply,
mapply. Once you get the hang of these they are really helpful.
Tom.
-Original Message-
From: Thomas Hopper
Dalgaard [mailto:[EMAIL PROTECTED]
Sent: Friday, 7 January 2005 4:46 PM
To: Mulholland, Tom
Cc: R-Help (E-mail)
Subject: Re: [R] Basic Linear Algebra
...
I can't. 4th element should be just 12.
I did indeed transcribe my solution incorrectly
...
--
O__ Peter Dalgaard
I don't normally have to go anywhere near this stuff , but it seems to me that
this should be a straight-forward process in R.
For the purposes of this enquiry I thought I would use something I can work out
on my own.
So I have my matrix and the right hand results from that matrix
tdata -
I thought this would be simple enough and in some ways it is, but I don't have
an answer I would call easy
The manual process would be to use barplot instead of screeplot, so that you
could return the x values. The other choice would be to hack the screeplot
function to return the values. (see
I thought it would be quite simple to generate a number of sequences, but my
only answer is to loop and this doesn't seem in keeping with R's capabilities.
In short I have the start and finish positions of several (well thousands
probably) sequences. So I am looking at how do I take
x -
Since I'm not sure what you are trying to achieve would [EMAIL PROTECTED] get
what you wanted without unclassing the object. I only ask this because my first
instinct when I started using the Rmetrics package was to try and use
ts1$positions which only gives you a null value.
Tom
Thanks for the hint on qcc. As for bollinger bands a generic porcess was put up
by Gabor Grothendieck.
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/22367.html
You might also like to have a look at the Rmetrics (fBasics, fSeries etc),
in particular RollingAnalysis {fSeries}
Tom
This has come up before you might try
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/36694.html
The bottom line is that you either have to make sure that the SPSS datasets are
in a format where no warnings come up or find a way to verify that no
corruption is happening. As one of the posts
1 - 100 of 132 matches
Mail list logo
