There is no such code in MASS, so please don't make random assignments of
blame. A competent R author would have used rowMeans(parmat).
R has debugging tools (see 'Writing R Extensions') and you should use them
to locate in which function the error is occurring. E.g. traceback()
would tell
Hi,
Given names - c(a, b, c, d). I would like to generate a vs. b, a
vs. c, a vs. d, b vs. c, b vs. d, c vs. d
Thanks
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
try this:
nams - c(a, b, c, d)
combn(nams, 2, paste, collapse = vs. )
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
On Tue, 11 Mar 2008, [EMAIL PROTECTED] wrote:
Hi guys,
I've got tthe following question:
Why do my fonts look a bit grubby when I use the jpeg() function? I see
the pixels of the font. The quality is set to 100%.
If this is normal with jpeg(), which function would you propose me to
use
Hi Brian,
Thank you for your fast answer. It looks a little bit better now. Nevertheless
I see pixels. May it depend on the font used?
In the praeambel I use
windowsFonts(Frutiger = windowsFont(Frutiger 45 Light))
par(family=Frutiger)
And later, to create a wide table on the
Hi,
could anyone explain to me what this warning message
exactly means and what the consequences are?
Is it due to the fact that there are very extreme
observations / outliers included or what is the reason
for it?
Thanks so much,
Werner
Machen Sie Yahoo! zu Ihrer Startseite. Los
On Tue, 11 Mar 2008, Werner Wernersen wrote:
Hi,
could anyone explain to me what this warning message
exactly means and what the consequences are?
Is it due to the fact that there are very extreme
observations / outliers included or what is the reason
for it?
See MASS4 pp.197-8.
Thanks again,
but I see the pixels on my print... The width is set to 1560, the height to 566.
Is the assignment to use Frutiger correct?
Regards
-Ursprüngliche Nachricht-
Von: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Gesendet: Dienstag, 11. März 2008 10:08
An: Schwander, Thomas,
Karen
levels returns the levels attribute of a variable, and a vector has no
such attribute. This is usually used with a factor, e.g.
temp - c(3, 5, 5, NA)
levels(factor(temp))
[1] 3 5
Best wishes
Richard
Chang Liu wrote:
Hello:
I'm trying to use levels function, but I don't know
On Tue, 11 Mar 2008, [EMAIL PROTECTED] wrote:
Hi Brian,
Thank you for your fast answer. It looks a little bit better now.
Nevertheless I see pixels. May it depend on the font used?
These _are_ bitmap formats: of course you will see pixels if you look
closely enough!
But beware that many
On 11-Mar-08 08:58:55, Werner Wernersen wrote:
Hi,
could anyone explain to me what this warning message
exactly means and what the consequences are?
Is it due to the fact that there are very extreme
observations / outliers included or what is the reason
for it?
Thanks so much,
Werner
On Mon, Mar 10, 2008 at 10:19:01AM -, michael watson (IAH-C) wrote:
I'd recommend outputting either as pdf or as a windows metafile
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Karin Lagesen
Sent: 10 March 2008 09:54
To:
On Tue, 11 Mar 2008, [EMAIL PROTECTED] wrote:
Thanks again,
but I see the pixels on my print... The width is set to 1560, the height to 566.
Is the assignment to use Frutiger correct?
How would I know? If you doubt this, try it with a standard font and
compare.
Note what I said about
Prof Brian Ripley wrote:
How would I know? If you doubt this, try it with a standard font
and compare.
Note what I said about artefacts: something other than R is involved
when you print. This is the sort of thing which happens if bitmap
images are not printed at their native
Hallo,
i get a warning message that NAs are introduced by coercion,
so my idea is to write a function to see which values are turned into
NA
For this i need to write a function to go through (loop) the original
data and the transformed (with the introduced na) to see which data were
transformed
Dear All,
Thank you for helping me doing PCA.
Please show me how to read Individuals factor map (PCA) and Variables factor
map (PCA).
http://www.nabble.com/file/p15975893/pca_var.jpg
http://www.nabble.com/file/p15975893/pca_ind_q1.jpg
Nguyen Tran Phuong Thao.
Hue University, Viet Nam.
--
I would like to reproduce a figure found in HBN Hynes The Ecology of
Running Waters on page 79 reprint first edition copyright 2001. This
is a graph with graphs of insect abundance through time lined up in 3d
as you proceed down river.
any help is appreciated
Stephen
--
Let's not spend our
On Tue, Mar 11, 2008 at 08:57:49AM +0100, Hans-Peter wrote:
I looked at the package (it's
http://cran.r-project.org/web/packages/xlsReadWrite/index.html) but
it's windows only. So no joy yet.
which platform would you need it?
Mac OS X and Linux.
--
Ajay Shah
thank you for your help.
Prof Brian Ripley wrote:
There is no such code in MASS, so please don't make random assignments of
blame. A competent R author would have used rowMeans(parmat).
that is why i asked for help. i've been using R for 2 months now. i
never claimed to be competent.
R has
You'll do much better with e.g. PDF, a format designed for printing and
that used by the printing industry.
On Tue, 11 Mar 2008, [EMAIL PROTECTED] wrote:
I'll try it. thx
-Ursprüngliche Nachricht-
Von: Alberto Monteiro [mailto:[EMAIL PROTECTED]
Gesendet: Dienstag, 11. März 2008 12:24
http://www.nabble.com/file/p15976467/table.pdf table.pdf
I want to create 8 trellis (|experiment) scatterplots with regression lines
for prey and pred in each plot with different colouring. So far i am
creating trellis plots- but the separation in these groups is never
possible.
Yvonne,
With ggplot2 you could tackle this problem with the code below.
library(ggplot2)
graphs- read.table(table.txt, header=T)
ggplot(data = graphs, aes(x = bif, y = abund)) + geom_point() +
geom_smooth(method = lm) + facet_grid(. ~ experiment)
Have a look at the ggplot2 website for more info
Wim Bertels [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Hallo,
i get a warning message that NAs are introduced by coercion,
so my idea is to write a function to see which values are turned
into NA
For this i need to write a function to go through (loop) the
original data and
I have a data frame of around 4500 rows which contain file name
prefixes and other interesting variables. each prefix corresponds to
multiple files and a I only have 4300 complete sets of files. If any
one of the files is missing for a given prefix, then the function I
use fails. The effect is
Dear list
Is there anywhere I could find further information on how to interpret
the output for a logistic regression for validate() from Design
package?. I tried ?validate and google but I cannot find information
on what the rows and the columns represent.
Thanks
David
station month bas
190 5 0.000
190 7 1.563
190 10 0.000
190 11 0.000
202 4 18.750
202 5 18.750
202 7 6.250
202 10 4.800
202 11 3.125
198 4 18.750
198 5 31.250
198 7 3.125
198 10
Simon Knapp wrote:
I have a data frame of around 4500 rows which contain file name
prefixes and other interesting variables. each prefix corresponds to
multiple files and a I only have 4300 complete sets of files. If any
one of the files is missing for a given prefix, then the function I
[EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Dear list
Is there anywhere I could find further information on how to
interpret the output for a logistic regression for validate() from
Design package?. I tried ?validate and google but I cannot find
information on what the rows and the
The standard errors are actually at the heart of the matter, so it's
good to know that they may be wrong. All of this is coming up as I'm
trying to correct the standard errors for design effects in the survey
data I'm using. Perhaps survey() is the way to go... Thanks,
Ben
On Tue, Mar 11, 2008
[EMAIL PROTECTED] wrote:
Dear list
Is there anywhere I could find further information on how to interpret
the output for a logistic regression for validate() from Design
package?. I tried ?validate and google but I cannot find information
on what the rows and the columns represent.
I have the following problem: how appropriate is my aov model under the
violation of anova assumptions?
Example:
a-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3)
b-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40,
50,50,60)
z-data.frame(a, b)
fligner.test(z$b,
I don't quite see how you can do this as a
3-dimenional graph but will something like this do
what you want?
Note I used xx rather than c for the data.frame name.
plot(x$bas~x$month, type=b, col=1,
ylim=c(min(xx$bas),max(xx$bas)),
xlim=c(max(xx$month),min(xx$month) ))
lines(y$bas~y$month,
stephen sefick [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
fish - station month bas
190 5 0.000
190 7 1.563
190 10 0.000
190 11 0.000
202 4 18.750
202 5 18.750
202 7 6.250
202 10 4.800
202 11 3.125
198
Have a look at the posting guide (link at bottom of
the post). :)
It is much easier to see what is happening if you
submit an working example. At the moment we'd just be
guessing at what you're doing.
If there is much data it would be much better to
supply a small subset (that includes all
Rerun your code and as soon as the error occurs type:
traceback()
This will show more information on where the error is actually occuring. If
this does not clear it up for you, then send more of an example of what you are
doing along with the output of traceback (see the link at the bottom
Subject: Re: Graphing question (multiple line graphs arranged spatially)
Newsgroups: R-help:gmane.comp.lang.r.general
To: stephen sefick [EMAIL PROTECTED]
stephen sefick [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
fish - station month bas
190 5 0.000
190 7 1.563
190
Hi,
I am running R on a cluster. Basic packages are installed at
~/MyHome/lib64/R/library/ and add-ons packages are installed at
~/R/x86_64-redhat-linux-gnu-library/2.5/. Both directories were created by R
by default.
Problems arises when I use add-on packages, R complains that they can't be
I am not sure what you mean by 3d in your graph (most attempts at adding 3d to
a 2d medium (paper, computer monitor) results in distortion rather than
illumination.
If what you really want is to add graphs to a diagram of the river at specific
locations then look at the subplot function in
John Fox wrote:
Dear JRG, Rolf, Ben, and Peter,
Frequency weights, possibly even non-integer weights, are useful for
surveys where observations are sampled with unequal probabilities of
selection. The approach in SPSS gives correct point estimates in this
situation but incorrect standard
I am having a little trouble getting R to do something without writing
a couple of very awkward loops.
I have a matrix of probabilities and I want to generate a new matrix
of ones and zeros where each element in the new matrix is the result
of a draw from a binomial distribution where the
Well, this can be realized with the OpenOffice API
and the rJava package i would guess:
http://rosuda.org/rJava/
I successfully use the OpenOffice APi with Rserve to get and send
data from or to the OpenOffice calc tool (Windows, Linux). The OpenOffice
API itself
and its use is good documented
# PROBLEM WITH ABLINE
# I have a question about a seemingly imprecise abline figure which I have
come accross now a few times.
# After reading variables a and b...
I think the following does what you want:
probMatrix - matrix(runif(5 * 5), 5, 5)
binomialMatrix - matrix(rbinom(5 * 5, 1, probMatrix), 5, 5)
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S User)
Economics Guy
imprecise are our eyes :)
check where 4336.543 actually is and note how the range of the data
confused you.
plot(a~b)
abline(model)
abline(h=4336.543, v=10)
b
On Mar 11, 2008, at 11:27 AM, Aad Termorshuizen wrote:
# PROBLEM WITH ABLINE
# I have a question about a seemingly imprecise
On Tuesday 11 March 2008 04:27:12 pm Aad Termorshuizen wrote:
AT # So my QUESTION is: is abline imprecise or am I making errors?
Either you are not precise or I did not got your problem.
plot(a~b)
abline(model)
points(10,4336,col=red)
as you see the predicted point (red) is exactly at the
I tried your code and could not get it to run on my installation of R,
so I may be missing something. But if you have a matrix of
probabilities (call it probs) and want to simulate random binomial
draws, can you not simply create a matrix of the same size of uniform
random numbers (runif()) (call
Dear R community,
by default, text in the strips of a trellis plot is centered in the
strip.
Is there a way to have the text left-aligned?
for example:
libary(lattice)
test - data.frame(x=rnorm(100), y=rnorm(100), a=rep(c(A: centered
text,B: centered text),50))
xyplot(y ~ x | a, data = test) #
I'm getting warnings when I plot a zoo object:
dts - chron(rep(2007-09-10, 5), paste(00:0, 0:4, :00, sep=),
c(y-m-d, h:m:s))
dat.zoo - zoo(1:5, dts)
plot(dat.zoo)
Warning messages:
1: In v[[perm[1]]] : partial match of 'm' to 'month'
2: In v[[perm[2]]] : partial match of 'd' to 'day'
The plot
I am unable to replicate that. I don't get any warnings. Perhaps
they are coming from chron and you have an old version of chron.
Try upgrading chron to the latest version (2.3-22) and try it again.
On Tue, Mar 11, 2008 at 12:25 PM, [EMAIL PROTECTED] wrote:
I'm getting warnings when I plot a
On 3/10/08, Thomas Steiner [EMAIL PROTECTED] wrote:
I cannot access internet through R.
My IT-guys told me that I should set the proxy and the port and then
everything will be fine. Where can I set them?
Recently, I also had problems installing packages when under a
transparent proxy. What
That did it, Gabor.
Thanks! I didn't think chron would be involved.
-- David
-Original Message-
From: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Sent: Tuesday, March 11, 2008 11:39 AM
To: David Reiner [EMAIL PROTECTED]
Cc: r-help@r-project.org
Subject: Re: [R] plot.zoo warnings -
All,
What purpose does the getTree function have in Random Forest? Can you graph it
as you can in rpart and can it be interpreted in the same way?
Helen Mills Poulos
Yale School of Forestry
__
R-help@r-project.org mailing list
Hi, maybe my question is not typical. I want to make R table-kind output to
be readable immediately in Latex but I am struggling with symbol '\' needed
for instance in '\hline'. For example when 'Table' is a matrix class object:
cat(,Table[1,1], ,Table[1,2], \n);flush.console()
cat(\hline
threshold wrote:
Hi, maybe my question is not typical. I want to make R table-kind output to
be readable immediately in Latex but I am struggling with symbol '\' needed
for instance in '\hline'. For example when 'Table' is a matrix class object:
cat(,Table[1,1], ,Table[1,2],
Thank you for that quick responses
2008/3/11, [EMAIL PROTECTED] [EMAIL PROTECTED]:
From: Gustave Lefou [EMAIL PROTECTED]
Date: 2008/03/11 Tue PM 12:32:25 CDT
To: [EMAIL PROTECTED]
Subject: [R] Multi arguments - return
if you're returning a list whose components
are vectors, then to return
I think you need to provide an example. We have no
way of knowing what you are putting in the list.
--- Gustave Lefou [EMAIL PROTECTED] wrote:
Hello,
I wanted a function to have two vectors x,y as its
output. I wrote
return(x,y) inside its definition. It worked but I
got a warning : it
Sorry, another newbie question :-(
I loaded a data set with 10 rows and 30 columns. The first column is
characters for names of car manufacturers:
Jeep
Nissan
Toyota1
Toyota2
Etc.
How can I replace Toyota2 with Scion?
Thanks again
[[alternative HTML
I am having difficulty loading a points shapefile (1600 records) I generated
from ArcGIS as a PointsSpatialDataFrame using the maptools readShapePoints
function.
library(maptools)
surveypts - readShapePoints(system.file(C:/temp/survey.shp,
package=maptools)[1]
Error in getinfo.shape(filen) :
Yep ! It works perfectly.
2008/3/11, [EMAIL PROTECTED] [EMAIL PROTECTED]:
From: Gustave Lefou [EMAIL PROTECTED]
Date: 2008/03/11 Tue PM 12:57:59 CDT
To: [EMAIL PROTECTED] [EMAIL PROTECTED],
Phil Spector [EMAIL PROTECTED]
Cc: [EMAIL PROTECTED]
Subject: Re: [R] Multi arguments -
On Tue, 2008-03-11 at 10:23 -0700, threshold wrote:
cat(\hline \n);flush.console()
you need to double the backslashes, as in
cat(\\hline \n);flush.console()
as in
cat(\\hline \n);flush.console()
\hline
The warning when using just one kind of hints at what is going on:
cat(\hline
Hi Luciana,
Try this:
yourData[,1] - sapply(yourData[,1], function(x){
x=as.character(x)
x[which(x==Toyota2)]-Scion
x
}
)
# Example
set.seed(123)
x=c(Jeep,Nissan, Toyota1, Toyota2)
y=rnorm(4)
DATA=data.frame(x,y)
DATA
x y
1Jeep -0.56047565
2 Nissan -0.23017749
3
Graphing: Not really. You might do somw work and get Rgraphviz to
produce the tree, but I'm not sure how useful that is.
Interpretation: You can read the tree just like any other decision
trees. However, remember that the splits are chosen with some
randomness, so I'm not sure how much trust
someone sent in a question earlier about doing
something in 3D so i took a stab at it purely
for educational purposes ( i'm not even sure that I understood the question
actually ).
Unfortunately, persp gives me an error that I don't understand because it says
object y not found. I'm sending y
If the column is factor:
levels(DATA$x)[levels(DATA$x) == Toyota2] - Scion
or character:
DATA$x[DATA$x == Toyota2] - Scion
On 11/03/2008, Suran, Luciana @ Torto Wheaton Research
[EMAIL PROTECTED] wrote:
Sorry, another newbie question :-(
I loaded a data set with 10 rows and 30
Dear useRs,
Is there a way one can add x-axis values to a histogram? Example:
z = rgamma(n=1000, shape = 8, scale = 1000)
hist(z, breaks = 30)
Currently I see 4 x-axis values: 5000, 1, 15000, 2. Can I add more
values or modify the ones currently showing? I tried ?hist, yet couldn't
Try this:
hv - hist(z, breaks = 30, xaxt='n')
axis(1, hv$breaks )
On 11/03/2008, diegol [EMAIL PROTECTED] wrote:
Dear useRs,
Is there a way one can add x-axis values to a histogram? Example:
z = rgamma(n=1000, shape = 8, scale = 1000)
hist(z, breaks = 30)
Currently I see 4 x-axis
On Tue, 11-Mar-2008 at 10:51PM +0800, Ng Stanley wrote:
| Hi,
|
| I am running R on a cluster. Basic packages are installed at
| ~/MyHome/lib64/R/library/ and add-ons packages are installed at
| ~/R/x86_64-redhat-linux-gnu-library/2.5/. Both directories were created by R
| by default.
|
|
The problem is that AIC is only defined for ML fitting, and gls defaults
to REML.
I have always maintained that it is a bug that nlme's logLik function
returns a log restricted likelihood for the default fits, and that this
is converted to a mis-named AIC. If you use
m1 - gls(height ~ age,
Dears useRs,
I am using the mle function but this gives me the follow erros that I
don't understand. Perhaps there is someone that can help me.
thank you for you atention.
Bernardo.
erizo - read.csv(Datos_Stokes_1.csv, header = TRUE)
head(erizo)
EDAD TALLA
10 7.7
21 14.5
31
probMatrix - matrix(runif(5*5), 5, 5) ## sets up the probs
binoMatrix - rbinom(length(probMatrix), 1, probMatrix) ## simulates
dim(binoMatrix) - dim(probMatrix) ## makes it a matrix also
What doesn't work?
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On
[EMAIL PROTECTED] schreef:
I am having difficulty loading a points shapefile (1600 records) I generated
from ArcGIS as a PointsSpatialDataFrame using the maptools readShapePoints
function.
library(maptools)
surveypts - readShapePoints(system.file(C:/temp/survey.shp,
package=maptools)[1]
This is not reproducible, so I am reduced to guessing.
Your loglike() is only defined for positive f, but you have not imposed
that constraint on your optimization. Further, if f=5.91e-05 is a typical
value, your problem is badly conditioned and the finite difference code in
optim() is
Dear list,
when I write .jinit(), the following error message coming out.
Error occurred during initialization of VM
Could not reserve enough space for object heap
Error in .jinit() : Cannot create Java virtual machine (-4)
I have C:\Program Files\Java\jre1.6.0_04 and C:\Program
I need to get more digits in predicting a test sample with a random
forests object. Format or options(digits=) do nothing. Any ideas?
Thank you,
Nagu
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
The help page for readShapePoints says that the file name that you
supply should not include the extension.
The system.file part of the example in the help page is only to use
the example file that came with the package. Since you want to open
your own file you omit the system.file() bit.
Do
[EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
someone sent in a question earlier about doing
something in 3D so i took a stab at it purely
for educational purposes ( i'm not even sure that I understood the
question actually ).
Unfortunately, persp gives me an error that I don't
I have an R program I am attempting to use to generate some SVGs. I've been
using the cairoDevice library. When running a session not connected to an X
session (like if I'm sshed in), attempting library(cairoDevice) causes:
*** caught segfault ***
address 0x8, cause 'memory not mapped'
Hi guys,
I remember that we had some ways to save a chart as ps file, which can be
used by latex directly. I can't find it now. Can you guys give me some
hints?
Thanks very much.
--
Samuel Wu
http://webclipping.com.cn
[[alternative HTML version deleted]]
Thanks Oldrich and Max,
I have some more queries.
If I need to train svm() with only one instance I get the following error:
Error in if (any(co)) { : missing value where TRUE/FALSE needed
Will it be wiser if I duplicate the instance with minute changes in the values,
or there is some other way
I got it solved via reading this material:
http://freshmeat.net/articles/view/2237/
Thanks all the same.
On Wed, Mar 12, 2008 at 12:04 PM, Samuel [EMAIL PROTECTED] wrote:
Hi guys,
I remember that we had some ways to save a chart as ps file, which can be
used by latex directly. I can't find
80 matches
Mail list logo