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,
On Mon, 21 Mar 2005, Thomas Lumley wrote:
On Mon, 21 Mar 2005, Tim Cutts wrote:
On 21 Mar 2005, at 4:42 pm, [EMAIL PROTECTED] wrote:
Hi,
I have a problem:I need to use the maximum amount of memory in order to
perform a very tough analysis. By purchasing the suitable computer, what's
the maximum
Hi !
I am interrested in counting process analysis (as described in the text book of
D.R. Cox the statistical analysis of series of events) and in particular
in calculating index of dispersion of counts IDC(t)=Var(Nt)/E(Nt). is there
some useful functions in R for such analysis ?
Many thanks
hi, all.
I am a beginner. Could you tell me how to use read.table() to read
following beams.txt file.
thanks
file content is here:
beams -
structure(list(strength = c(11.14, 12.74, 13.13, 11.51, 12.38,
12.6, 11.13, 11.7, 11.02, 11.41), SpecificGravity = c(0.499,
0.558, 0.604, 0.441, 0.55,
Thanks to Uwe Ligges! The problem seems to have been my fault. Using his
advice b) was the solution for me (under R-2.0.0). (My use of postscript()
must have been a relict from my migration from S-PLUS to R.)
a) You should try R-2.0.1 with a recent version of lattice (or help to
try out and
Try source(beams.txt) which reads and executes the commands in
beams.txt.
It might make more sense in future to store and work with dataframes,
especially with big datasets.
Regards, Adai
On Tue, 2005-03-22 at 17:00 +0800, chin wei wrote:
hi, all.
I am a beginner. Could you tell me how to
On Mon, 21 Mar 2005, Tae-Young Goo wrote:
Hello.
I've tried to install R to IBM AIX(v.5.1) machine.
I've used compile options indicated by R-admin.
There are several sets there, so which exactly?
Then, I met following error messages.
/home/local/R_2.0.1/lib/R/bin/exec/R is unchanged
Hi Pascal,
One thing you can do is to work on indexes, count unique occurences and
assign them, as following:
index - as.data.frame(table(index[,1],index[,2]))
index - do.call(cbind,lapply(index, as.numeric)) # ensures numeric
coding (as table turns into factors)
occurence[index[,1:2]] -
On Tue, Mar 22, 2005 at 11:59:25AM +0100, Christoph Lehmann wrote:
I produce graphics with R under linux, but my collaborators often use
windows and cannot import eps pics e.g. in msword
what is the standard way to get e.g. bmp's with the same quality as eps.
going the way: creating eps,
Dear useRs,
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Jan T. Kim
Sent: Tuesday, March 22, 2005 1:18 PM
To: r-help@stat.math.ethz.ch
Subject: Re: [R] r under linux: creating high quality bmp's
for win users
On Tue, Mar 22, 2005 at
On Tue, 22 Mar 2005, Christoph Lehmann wrote:
Hi
I produce graphics with R under linux, but my collaborators often use windows
and cannot import eps pics e.g. in msword
Are they sure? One can import eps images into Word, but you need to print
them on a postscript printer.
what is the standard
Dear R gurus
I've a question concerning the transformation of a new S4 class
(that I defined) and the existing class POSIXlt. I did the
following:
## Definition of the new class:
setClass(dtime,
representation(data = POSIXlt),
prototype(data = as.POSIXlt(2004/06/01)))
##
Sorry -- I meant to say dataframe instead of matrix.
Anyway I see that my troubles are gone when I scan the data in as a vector
then convert to matrix. (I had troubles doing such manipulations when I
read in the data using read.table)
x-scan(/home/wsimpson/papers/face/max.dat)
xx-matrix(x,
Dear Anthony
I don't know how SAS calculates the p-value, but in R the
p-value is calculated under the assumption that the parameters
of the distribution (you want to compare with your samples) are
known and not estimated from the data.
In your example you estimate them from the data (by mean(w)
Pascal BLEUYARD p.bleuyard at opgc.univ-bpclermont.fr writes:
:
: Hi all,
:
: I need to compute some occurence matrix: given a zero matrix and a set
: of paired indexes, I want to store the number of occurences of each paired
: index in a matrix. The paired indexes are stores as an index
Hello,
I have a problem with nls() when I want to fit a model with more than two
parameters to be fitted which are written as start=list(a1=,a2=,a3=,...).
Everytime, it displays:'error... singular gradient'
it's a real pain!!!
Hope sb knows something about this
Thanks
GS
On Tuesday 22 March 2005 03:06, Gerrit Eichner wrote:
Thanks to Uwe Ligges! The problem seems to have been my fault. Using
his advice b) was the solution for me (under R-2.0.0). (My use of
postscript() must have been a relict from my migration from S-PLUS to
R.)
Actually, that advice is more
a data file (close to 60M) with 90,000 rows and 173 columns. Is it too
big for R?
Thanks.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hello Christoph!
In the past, I used an utility called eps2wmf.
It only works under Windows though (maybe under Linux with wine?).
I believe it is available on the CTAN (Tex archives).
The nice thing is that wmf files are not bitmap and scale well.
Christophe Pallier
Christoph Lehmann wrote:
Hi
I
It might be too big to R on the machine your are using, and it's also highly
dependent on what you want to do with it. E.g., if you try to cluster the
rows using something that requires the full distance matrix, you're most
likely out of luck. I've dealt with data with 20,000 rows but more than
I have also recently come across a conversion program that works with
ghostscript to convert ps and pdf images into other vector formats,
including .wmf and .emf for M$ software: http://www.pstoedit.net/. This
might be helpful if your colleagues do not have, as Dr. Ripley
indicated, a postscript
Hi everybody,
I would like to know if it is possible to find the RPM
of R 2.0.1 for SuSe 9.2 AMD64 (Athlon 64 bit).
Alternatively, where to find clear instructions how to
compile it on the same OS (I'm not an expert).
Many thanks in advance.
Tino
On Mon, Mar 21, 2005 at 04:04:50PM -0800, Spencer Graves wrote:
1. Have you looked at the survival package? Venables and Ripley
(2002) Modern Applied Statistics with S (Springer) has a chapter on
survival analysis that I found quite helpful. The survival package
includes hazard plots,
Hello !
I'm struggling again against lattice graprhics. ;) I'm trying to produce
a conditionnal xyplot with two conditionning factors (let's say A and
B). I want the levels of those factors (A1, A2, etc) to show in the
margins of the lattice plot, not in the strips between the panels.
A1
Many people could help you, but the question is too general. In
brief, it means that the algorithm has found a place where the
(estimated?) matrix of first or second partial derivatives is of reduced
rank, and it refuses to do more. For such problems, I often use optim.
If you
Or at least read the posting guide!
http://www.R-project.org/posting-guide.html. Working through that can
help you get a more informative reply. spencer graves
Uwe Ligges wrote:
Noel Anel wrote:
Hello!
Which class(model) is most appropriate for numeric predition?
This is a joke, isn't
d'habitude je mets juste xyplot(y~x | factor(z)...) ou xyplot(y~factor(x)... )
il faut mettre la variable dont tu veux voir les niveaux comme factor()
GS
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
Can you provide a link. I did a google search and found something
on a Japanese site but it turned out that the writer had made a
mistake and it linked to wmf2eps, not eps2wmf.
Christophe Pallier pallier at lscp.ehess.fr writes:
:
: Hello Christoph!
:
: In the past, I used an utility called
I regularly work with data frames with around 144000 rows and 23 columns.
In fact, I work with two of them in the same session, one in
.GlobalEnv and the other in search()[2], attached using attach() on a
file previously saved with save().
-Don
At 9:11 AM -0500 3/22/05, Wensui Liu wrote:
a data
Hi All,
Does someone have an idea of how to cleverly convert a categorical
timeseries into a transition matrix?
Ie, I have something like:
x- c(1,1,2,1,1,2,2,2,1,2),
And I want a matrix with counts and/or probabilities:
tr - matrix(c(2,3,2,2),2,2)
tr
[,1] [,2]
[1,]22
[2,]32
Hi All,
I have a quick question and any help is greatly appreciated. For the following
data when I try to produce the image using xyplot function in lattice package,
the key has 4 rows instead of 2. Can anyone tell me what I'm doing wrong and
what is the way to fix the problem. Here the code
The following, using table, seems to work:
x - sample(letters[1:2], 10, replace=T)
x
[1] b a a b b a a b b b
table(x[-1],x[-length(x)])
a b
a 2 2
b 2 3
Hope this helps,
Giovanni Petris
Date: Tue, 22 Mar 2005 11:20:40 -0500
From: Ingmar Visser [EMAIL PROTECTED]
Sender:
I second the use of png() to create portable network graphics format
pictures, as MS Word (XP and 2003 at least) can handle these without any
problems. I use this format as a matter of course in preparing material for
embedding in web pages.
Regards,
Mike
-Original Message-
From:
All,
I'm currently using R 2.0.1 on a Powerbook G4 with OS X 10.3.8. So far
the only way I've found to set my proxy is by doing
Sys.putenv(http_proxy=insert proxy url:proxy port) everytime I
start up R.
This works fine, but I'd like to find a solution that doesnt require
manual input everytime
Thanks to everyone for your quick and helpful answers!
ingmar
On 3/22/05 11:40 AM, Christos Hatzis [EMAIL PROTECTED] wrote:
Here is a way to do this:
x - c(1,1,2,1,1,2,2,2,1,2)
y - cbind(x,c(x[-1],NA)) # time-shifted by one
aggregate(y, by=list(y[,1],y[,2]), length)
-Christos Hatzis
On Tuesday 22 March 2005 10:26, Ghosh, Sandeep wrote:
Hi All,
I have a quick question and any help is greatly appreciated. For the
following data when I try to produce the image using xyplot function
in lattice package, the key has 4 rows instead of 2. Can anyone tell
me what I'm doing wrong
Dear R list,
How can confidence interval be derived for e.g. the Tau-a coefficient or the
c index (area under ROC curve) such that I can compare the fitted lrm
(logistic) models with each other. Is this possible?
The aim is to conclude that one model is significantly better than other
model (a,
Let me modify my question about user-defined covariance structures for LME
models: Can somebody tell me how I can see the code for the definition of the
correlation structures that come with the NLME package. Specifically I like to
see the code for the functions coef, corMatrix, and intialize
Hello,
I am writing a package called 'DLM' containing a vignette.
The vignette contains a chunck with the function call 'library(DLM)'.
This worked fine with 'R CMD check DLM', but when it comes to building
the package with 'R CMD build DLM' I get the following error message:
* creating
Drew Balazs wrote:
All,
I'm currently using R 2.0.1 on a Powerbook G4 with OS X 10.3.8. So far
the only way I've found to set my proxy is by doing
Sys.putenv(http_proxy=insert proxy url:proxy port) everytime I
start up R.
This works fine, but I'd like to find a solution that doesnt require
On Tue, 22 Mar 2005, Shelby wrote:
I have also recently come across a conversion program that works with
ghostscript to convert ps and pdf images into other vector formats,
including .wmf and .emf for M$ software: http://www.pstoedit.net/. This
AFAIK, it converts to those formats only under
Jan Verbesselt wrote:
Dear R list,
How can confidence interval be derived for e.g. the Tau-a coefficient or the
c index (area under ROC curve) such that I can compare the fitted lrm
(logistic) models with each other. Is this possible?
The aim is to conclude that one model is significantly better
All,
I have a windows application in which data are written
to MySQL and R analysis initiated by a VB frontend.
After upgrading to DCOM 1.35, the analysis has become
excruciatingly slow when processing a simple RODBC
query. I am using DCOM 1.35 - I didn't have this
difficulty prior to upgrading.
Hi,
According to help(sub), the ^ should match the zero-length string at the
beginning of a string:
sub('^','var',1:3) # 1 2 3
sub('$','var',1:3) # 1var 2var 3var
# This generates what I expected from the first case:
sub('^.','var',11:13) # var1 var2 var3
Dave
--
Dr. David Forrest
[EMAIL
Hi,
I am using R 2.0.1 Mac OS X 10.3.5. I was able to authorize R to run system
command as root using login and password. But, I forgot to turn this function
on before installing Bioconductor, so I terminated the process of the
installation. Then, I tried to authorize R to run system command
Dear R-users,
I have an outcome variable and I'm unsure about how to treat it. Any
advice?
I have spending data for each county in the state of California (N=58).
Each county has been allocated money to spend on any one of the
following four categories: A, B, C, and D.
Each county may spend the
Hi All,
My question is in regards to an error generated when using randomForest
in R. Is there a special way to format the data in order to avoid this
error, or am I completely confused on what the error implies?
Error in randomForest.default(m, y, ...) :
Can not handle categorical
When trying to install an older version of R on a
windows XP system and recieved the following error.
R for windows GUI front end
App Rgui.exe
App version 1.6.1.21101.0
Mod MSUCET.dll
mod version 7.02600.1106
offset 0003213b
I do not know what this means. Please help.
Let me modify my question about user-defined covariance structures
for LME models: Can somebody tell me how I can see the code for the
definition of the correlation structures that come with the NLME
package. Specifically I like to see the code for the functions coef,
corMatrix, and intialize
Hi:
I am doing a confirmatory factor analysis now. In
the analysis, I have null hypothesis test which
specify some special structure for the loading matrix.
And the alternative is there is no such special
structure. Then the log likelihood ratio test can be
used. The problem I have is my data
It always helps to check whether you got the data into R correctly. Hint:
What does str(credit) tell you?
Andy
From: Melanie Vida
Hi All,
My question is in regards to an error generated when using
randomForest
in R. Is there a special way to format the data in order to
avoid this
Dear Liu,
I'm not aware of any way of doing what you want in R at present. You can
base the analysis on polychoric correlations, which can be computed using
the polycor package, and you can get consistent estimates of the factor
loadings and factor correlations from the sem packages, but the
Hi,
Is there any tool to check if there is update version of a package
available? I look for things alike YUM for linux?
YD
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi,
Is there some alternative to the 'groups' argument in lattice's bwplot
function for boxplots? Say in the example below:
bwplot(yield ~ site | year, data = barley)
you want to have two side by side boxplots per site, corresponding to each
year in the barley data frame. Ideally, the space
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.
Hello
I am facing the following problem using the R-version 1.9.1
The PDF or PS none of these device drivers are opening while I am using
R-1.9.1, the following error message is coming
Error in PS(file, old$paper, old$family, old$encoding, old$bg, old$fg,
:
unable to start device
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
Senhua Yu wrote:
Hi,
I am using R 2.0.1 Mac OS X 10.3.5. I was able to authorize R to run system command as root using login and password. But, I forgot to turn this function on before installing Bioconductor, so I terminated the process of the installation. Then, I tried to authorize R to run
Melanie Vida wrote:
Hi All,
My question is in regards to an error generated when using randomForest
in R. Is there a special way to format the data in order to avoid this
error, or am I completely confused on what the error implies?
Error in randomForest.default(m, y, ...) :
Can not
Derek Margetts wrote:
When trying to install an older version of R on a
windows XP system and recieved the following error.
R for windows GUI front end
App Rgui.exe
App version 1.6.1.21101.0
Mod MSUCET.dll
mod version 7.02600.1106
offset 0003213b
I do not know what this means. Please help.
R has
David == David Forrest [EMAIL PROTECTED]
on Tue, 22 Mar 2005 15:02:20 -0600 (CST) writes:
David According to help(sub), the ^ should match the
David zero-length string at the beginning of a string:
yes, indeed.
David sub('^','var',1:3) # 1 2 3
David sub('$','var',1:3) #
62 matches
Mail list logo