I am trying to obtain multiple legends on a densityplot, using the
legend parameter. I am trying the following code, with just one
legend at the moment, which doesn't work. I get
Error in valid.data(rep(units, length.out = length(x)), data) :
No 'grob' supplied for 'grobwidth' unit
Please,
send your question to the ade4 list. Be more precise on your problem.
http://pbil.univ-lyon1.fr/ADE-4/adelist.html
Miguel Ribeiro wrote:
Hi, i'm having trouble using Statis in ADE4 package.
I want to study a matrix with 18x414...could any body help me?!?!'
Urgent
Thanks
Miguel
Please note the message is about the protect stack, and not the `stack'.
The message you quote did not come from a current version of R (and you
have not told us what version you are using, as asked in the posting
guide).
On Mon, 17 Oct 2005 [EMAIL PROTECTED] wrote:
Hi,
I'm using step() for
What version of R is this? That error message does not exist in current
versions of R as far as I can see (nor has it since 2.0.1). So I believe
the advice in the posting guide applies:
If you are using an old version of R and think it does not work
properly, upgrade.
In earlier
Hi
You are looking for functions
panel.*
especially panel.lmline, but I wondered if you can use linear for one
panel and quadratic for other panels. You could use a structure
provided in examples in xyplot help page to try to achieve what
you want. I do not have instant solution to your
Hi Everyone.
I am rather new to R and I've been trying to implement a function to
carry out the above test. For a couple of days now I've been stuck on
how to generate average rank differences.
Say I have a vector of average ranks:
averank- c(2,5,9,12)
I would like to subtract averank[1]
the data file has such structure:
1992 6245 49 . . 20 1
0 0 8.739536 0 . . .
. . . . .alabama
. 0 .
I don't know Dunn's rank test, but the following substracts each of
the sums of averanks from the next rank.
cumsum(averank)[-length(averank)] - averank[-1]
Hank
On Oct 17, 2005, at 4:30 AM, Iain Gallagher wrote:
Hi Everyone.
I am rather new to R and I've been trying to implement a
On Sun, 16 Oct 2005, Weijie Cai wrote:
Can anybody give me a quick answer that whether step() can be directly
applied to a glm model including both continuous variables and categorical
variables (by factor())? Maybe my data is not suitable, I applied step() on
my glm object and did not get
The help page means exactly what it says, but the English is too subtle
and I have reworded it.
I have no idea why you are interested in pairlists (they are hardly used
at user-visible level these days). The point is that pairlist() is NULL
and so strictly not a pairlist at all (try typeof()
Dear friends,
Do you guys have any idea about the release of JGR for R 2.2? Thanks!
--
***
Weimin Wang
Adrress:
Studio 10, 9 Chemin du Routoir,
ILLKIRCH, 67400, France
Tel:(+33)03-90-24-45-24(O),
(+33)03-88-65-42-45(D)
***
Dear R helpers,
I have a largish matrix (1300 x 1300) and I wish to find the row and column
numbers that identify particular elements whose values I know in advance
(for example, the row and column numbers for the maximum value of the
matrix).
I have looked in the help manual and found the
How about the function which?
?which
On 10/17/05, Jonathan Williams
[EMAIL PROTECTED] wrote:
Dear R helpers,
I have a largish matrix (1300 x 1300) and I wish to find the row and column
numbers that identify particular elements whose values I know in advance
(for example, the row and column
Hi Jose,
I am just beginning to plumb the depths of lattice, but perhaps my
recent experience can help.
I recently figured out (with encouragement from the list) how to plot
predicted values from a model into the appropriate panel. I am
certain that what I have done can be done better, but
Dear ronggui,
I didn't find any attachments, but using the data lines in your
message, and assuming that . represents missing data, the following
appears to do what you want:
as.data.frame(scan(c:/temp/ronggui.txt,
list(year=1, apps=1, top25=1, ver500=1,
mth500=1, stufac=1, bowl=1,
Hi Jose,
I need to make a small correction in my code -
mod$data$B works for glm objects, but not lm or aov objects. For
those use, mod$model$B.
On Oct 17, 2005, at 6:26 AM, Martin Henry H. Stevens wrote:
Hi Jose,
I am just beginning to plumb the depths of lattice, but perhaps my
recent
# create a 4x4 matrix with random values
smallmatrix = matrix(runif(16)*10, ncol = 4, nrow = 4)
# get the row and column number of the item in the
matrix with the highest value
which(smallmatrix == max(smallmatrix), arr.ind = 1)
HTH,
Martin
--- Jonathan Williams
[EMAIL PROTECTED] wrote:
Dear
In my job I write custom computer programs for data analysis, which
are used in our company's consulting business. Whenever I've needed
statistical analyses I've coded the algorithms myself, but my boss
wants me to start learning and using R, to speed up development.
I am very reluctuant to
R Language Definition
Writing R Extensions
1988 S book (the Blue Book)
1992 S book (the White Book)
S programing
Programing with data （about S4）
=== 2005-10-17 19:36:48 您在来信中写道：===
In my job I write custom computer programs for data analysis, which
are used in our
Dear R-list
I have a dataset like below (points), how can I produce a max value for
each column. I need a result like (I hope my eye correct):
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
[1,] 211 10 99 14 19695 8 5 7
This is a
apply(points, 2, max)
Should do the job.
For details look at ?apply
Best,
Matthias
Dear R-list
I have a dataset like below (points), how can I produce a max
value for
each column. I need a result like (I hope my eye correct):
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
KJ,
In my job I write custom computer programs for data analysis, which
are used in our company's consulting business. Whenever I've needed
statistical analyses I've coded the algorithms myself, but my boss
wants me to start learning and using R, to speed up development.
I'm curious as to what
Acchhh very easy, time to drink a cup of coffe, but
Thank you for your all.
apply(points, 2, max)
Best regards, Muhammad Subianto
On this day 17/10/2005 02:34 PM, Muhammad Subianto wrote:
Dear R-list
I have a dataset like below (points), how can I produce a max value for
each
Dear list member,
I am using the packages RDCOMClient and SWinTypeLibs and try to import a COM
object (created in Delphi) in R that is of type 'early binding' instead of
late 'late binding'. Is there a possibility to do this in R?
Currently, the following returns an error message:
l1 =
Hi
I do not make programs in R to be available for others to use. But I
**do** analyse data and make some conclusions from such analysis. And
if I have to analyse the same data every (day, week, month, year) I
prepare a program (function or several functions) to be able in
Upgrading worked. Thanks!
- Original Message -
From: Prof Brian Ripley [EMAIL PROTECTED]
To: taskletter [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Sent: Monday, October 17, 2005 12:41 AM
Subject: Re: [R] Error in lazyLoadDBfetch
What version of R is this? That error message does
[A lot of polite and constructive critique deleted]
Is my impression correct that R is simply not well-documented enough
for serious programming?
No.
Have I missed a key reference to programming R?
Yes.
How about reading the text that R displays when it starts (and follow
its suggestions)?
Schlatter Christian wrote:
Dear list members
I'm very new to R but a little informed about geostatistics.
As I was looking for possibilities of geostatistical analysis in R I
encountered at least two very interesting packages:
Rgeo and gstat
And of course I'm wondering now about
Hello,
I write to know how can I modify the x axis : when I plot a survival object, R
plots a graph with x values = 0, 10, 20, 30 while I want a graph with values 0,
6, 12, 18, 24 in the x axis. How can I do this? In R 2.1.1 version there was
time.inc in survplot, but in version R 2.2.0 there
All,
In the help for memory.size, there is:
Details:
Command-line flag '--max-mem-size' sets the maximum value of
obtainable memory (including a very small amount of housekeeping
overhead).
How does one implement a command-line flag in order to set the max for memory?
cheers,
hi,
first shut down automatic axis in plot with
plot( ...,axes=F)
then
customize your axis with
?axis
axis(side, at = NULL, labels = TRUE, tick = TRUE, line = NA,
pos = NA, outer = FALSE, font = NA, vfont = NULL,
lty = solid, lwd = 1, col = NULL, padj = NA, ...)
you can
Hi all,
This is my first post, I hope you will help me.
I've some data to present with histograms. I have few values with almost 99% of
the frequencies (thousands) and some other values with low frequencies (below
one hundred) that I want to emphasize. I think if I could present the
Afshartous, David wrote:
All,
In the help for memory.size, there is:
Details:
Command-line flag '--max-mem-size' sets the maximum value of
obtainable memory (including a very small amount of housekeeping
overhead).
How does one implement a command-line flag in order to
On Mon, 2005-10-17 at 16:51 +0200, Jorge Gaspar Sanz Salinas wrote:
Hi all,
This is my first post, I hope you will help me.
I've some data to present with histograms. I have few values with almost 99%
of
the frequencies (thousands) and some other values with low frequencies (below
one
Dear all,
Despite searching the fullrefman file I have been unable to locate a
function to calculate the relative risk and confidence intervals from a
simple 2 by 2 contingency table. Can anyone tell me how this is done.
Thanks in advance
Malcolm
--
MJ Price, Social
--- MJ Price, Social Medicine [EMAIL PROTECTED] wrote:
Dear all,
Despite searching the fullrefman file I have been unable to locate a
function to calculate the relative risk and confidence intervals from
a simple 2 by 2 contingency table. Can anyone tell me how this is
done.
Thanks in
Can someone, please, explain the difference in as.POSIXct results
before 1970 and on and after 1970 as illustrated below.
After 1970, the use of 'EST' or EST+5EDT' as the timezone does not
affect the result of asPOSIXct, but before 1970 on 10/28 the results
are different.
Dear All,
Does anybody know the meaning of this error from using ginv() function please?
Error in La.svd(x, nu, nv) : error code 1 from Lapack routine 'dgesdd'
Many thanks in advance.
Regards,
Werner
-
[[alternative HTML version
Event: 2005 Insightful User Conference
Dates: Oct 26-27, 2005
Location: Princeton, NJ
URL: http://www.insightful.com/news_events/2005uc/ for details on pricing,
hotel accommodations and to register for this event.
The Insightful 2005 User Conference is being held October 26th-27th in
But who is on the panel?
On 10/17/05, Michael O'Connell [EMAIL PROTECTED] wrote:
Event: 2005 Insightful User Conference
Dates: Oct 26-27, 2005
Location: Princeton, NJ
URL: http://www.insightful.com/news_events/2005uc/ for details on pricing,
hotel accommodations and to register for this
Your OS is? I cannot reproduce this on any of mine.
Sone OSes do not support dates before 1970, which would explain the
difference. R then has to guess when summer time ended in 1969, and I
suspect something is happening inconsistently.
Also, EST is not a real-world POSIX timezone and may
Dear R users,
I have been using the pdMat class pdIdnot (from the mgcv
package)instead of pdIdent to avoid overflow in GLMM fits with
the MASS package function glmmPQL, of the following form:
fit1 - glmmPQL(fixed=y0~-1+xx0, random=list(gp=pdIdent(~-1+zz0)),
family=binomial)
Pretty neat.
Jarek
\
Jarek Tuszynski, PhD. o / \
Science Applications International Corporation \__,|
(703) 676-4192 \
[EMAIL PROTECTED] ` \
Thanks for your reply Hank. It's not really what I'm
after (though it's good to know).
For the test ( as described in Statistics for the
Biosciences by W. Gardiner. Prentice Hall, 1997) I
have to rank my groups, calculate the average rank,
then subtratc each average rank from every other. Any
Hi,
I am trying to fit a ordinal GEE model using ordgee {geepack}. In order to
check the validity of the function, I specified the correlation structure as
independence (i.e. constr = independence) and compared the result with that
using polr {MASS}.
Because a GEE model with an independent
I think Martin told you the basic approach to the indexing:
averank-sort(sample(1:100,25,replace=TRUE))
averank[-1] - averank[-length(averank)]
[1] 1 1 6 3 4 14 1 1 8 1 2 6 5 4 10 0 3 2 1 11 1 1 2 0
averank
[1] 4 5 6 12 15 19 33 34 35 43 44 46 52 57 61 71 71 74 76 77 88
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
[It is best to ask these questions on the [EMAIL PROTECTED]
mailing list or directly to me as the maintainer of the code.]
There are new versions of the RDCOM packages. They are not
officially released because I haven't had time yet.
But they are
I think I misunderstood your follow-up question. Try this:
averank-sort(sample(1:100,10,replace=TRUE))
x=matrix(nrow=length(averank),ncol=length(averank))
for (i in 1:length(averank)){
+ for (j in 1:length(averank)){
+ x[i,j] - averank[i] - averank[j]
+ }}
x
[,1] [,2] [,3] [,4] [,5]
Hello!
I am a R user. recently,I am using Bayesian Network to inference
posterior.but,I don't know how to input prior probability and inference
posterior.I tried prob function,but fail.Can you help me? Thanks!
ChuanXiang Ma
[[alternative HTML version deleted]]
Dear List,
I have a two-part question related to bpplot(), a box-percentile plot
function in the Hmisc package.
Take the example given in the Help for bpplot(), for instance.
(1) How does one set but not draw the y-axis? What I did was,
bpplot(... , yaxt=n), but that apparently does not work
Dear R list,
I am reading in text file data prepared in Access database by someone. One
of the field contains timestamp data, how can I separate the timestamp data
into two varaibles: date and time. Can I specify the field is in timestamp
format when I first reading in ?
My reading in data are
So the answer seems to be Insightful folks and industry folks.
On 10/17/05, A.J. Rossini [EMAIL PROTECTED] wrote:
But who is on the panel?
On 10/17/05, Michael O'Connell [EMAIL PROTECTED] wrote:
Event: 2005 Insightful User Conference
Dates: Oct 26-27, 2005
Location: Princeton, NJ
On 18 October 2005 at 16:32, [EMAIL PROTECTED] wrote:
| I am reading in text file data prepared in Access database by someone. One
| of the field contains timestamp data, how can I separate the timestamp data
| into two varaibles: date and time. Can I specify the field is in timestamp
| format
Hi,
I am running a generalized additive semi-parametric regression with five
non-parametric variables, and I get the error vector memory exhausted
(limit reached?). Do you know how to increase vector memory in R?
Kent
__
R-help@stat.math.ethz.ch
xpRt.wannabe wrote:
Dear List,
I have a two-part question related to bpplot(), a box-percentile plot
function in the Hmisc package.
Take the example given in the Help for bpplot(), for instance.
(1) How does one set but not draw the y-axis? What I did was,
bpplot(... , yaxt=n), but
Dear helpers,
I am a Dutch student from the Erasmus University. For my Bachelor thesis I
have written a script in R using boosting by means of classification and
regression trees. This script uses the function the predefined function
rpart. My input file consists of about 4000 vectors each having
On Thu, 13-Oct-2005 at 05:33PM -0500, Marc Schwartz (via MN) wrote:
| On Thu, 2005-10-13 at 15:20 -0600, Chris Buddenhagen wrote:
| Dear all
|
| I am using R to produce ordinations library(vegan) and the plot function
| produced looks great on the screen but when I send it to jpg or pdf or
Hi All,
Currently I'm working in FIGARCH process [Fractionally Integrated
Generalized Autoregressive Conditional Heteroscedasticity]. I've already
got the codes to do the process in S-Plus. Can anyone help me to do it
in R?
Thanks,
SUMANTA BASAK.
58 matches
Mail list logo
