Displaying 20 results from an estimated 10000 matches similar to: "Anyone use pls package?"
2002 Apr 24
2
Changing the colour in boxplots using bwplot()
I've managed to find how to change
trellis.settings$box.rectangle$col
and
trellis.settings$box.umbrella$col
but I can't figure out how to change the colour of the point used
to indicate outliers. There doesn't seem to be a general colour.
What am I overlooking?
best
--
*************************************************************
___ Patrick Connolly
{~._.~}
2002 Apr 28
2
dropterm() in MASS
To compare two different models, I've compared the result of using
dropterm() on both.
Single term deletions
Model:
growth ~ days + I(days^0.5)
Df Sum of Sq RSS AIC
<none> 2.8750 -0.2290
days 1 4.8594 7.7344 4.6984
I(days^0.5) 1 0.0234 2.8984 -2.1722
AND
Single term deletions
Model:
growth ~ days + I(days^2)
Df Sum
2002 Apr 18
2
Background in lattice plots using dotplot()
platform i586-pc-linux-gnu
arch i586
os linux-gnu
system i586, linux-gnu
status
major 1
minor 4.1
year 2002
month 01
day 30
language R
I can't seem to get rid of the slightly green background when I use
dotplot(). It shows on the screen in the
2002 Jan 25
1
Better idea than Poisson?
My data has three factors, a discrete response and an offset column.
>From the summary of a glm object derived using three way interactions,
I get this deviance information:
(Dispersion parameter for poisson family taken to be 1)
Null deviance: 39244 on 896 degrees of freedom
Residual deviance: 11913 on 795 degrees of freedom
AIC: 13905
Number of Fisher Scoring iterations: 6
The
2001 Nov 27
1
What name-value pairs can be used in gpar()?
To use grid graphics with lattice, the help for gpar indicates that
any number of named arguments can be specified but where do we find
what names make any sense?
I have ascertained that fontsize is given in points such as one
would use in a postscript device argument and uses that same name.
It's not as I'd have thought, something corresponding to cex in
regular par() settings.
2002 Sep 11
1
one question about title
Dear,
Thanks for your help. My question is the title can not completely show in the drawing area when the title is too long.
So, I need to change the value of cex.main each time. Could you help me to control the title display area?
Regrad,
ken
-----Original Message-----
From: Patrick Connolly [mailto:p.connolly at hortresearch.co.nz]
Sent: Wednesday, September 11, 2002 4:41 AM
To: Ken
2002 Sep 18
2
More on list to data frame (was: Re: List to Data Frame
Hi,
Now suppose I have just one list called FOO, which has 25 objects, e.g.:
[[1]]
1 2 3 4 5
[[2]]
6 7 8 9 10
.
.
.
And I want to do something like:
FRED <- data.frame(cbind(unlist(FOO[[1]]),
unlist(FOO[[2]]),
# ... for all 25 subsets
))
Is it possible to do this, without doing unlist(FOO[[i]]) 25
2002 Mar 21
1
Underdispersion with anova testing methods
Using anova of a glm with test = "Chisq", I get this:
Analysis of Deviance Table
Model: poisson, link: log
Response: Days
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev P(>|Chi|)
NULL 373 370.56
Block 3 71.05 370 299.51 2.543e-15
Variety 1 94.04 369
2002 Aug 12
1
Level sets of factors are different (panel.superpose)
I investigated why I was getting this error message...
Error in Ops.factor(groups[subscripts], vals[1]) :
Level sets of factors are different
which led me to putting a browser in panel.superpose:
Called from: panel.superpose(x, y, subscripts, ...)
Browse[1]> vals
vals
[1] 1 2 3
Levels: 1 2 3
Browse[1]> groups[subscripts]
groups[subscripts]
[1] 1 2 3 1 2 1 2 3 1 2 3 1 2 3
Levels: 1 2
2002 Oct 09
1
grid graphics and par settings
I see that adj is not a graphical parameter that can be used in a key
list to use with lattice functions.
Is that because it's not yet implemented in gpar? If so, it could
also explain why it's not available in grid.text.
Is there a work around? Centering looks silly on my plot.
best
--
Patrick Connolly
HortResearch
Mt Albert
Auckland
New Zealand
Ph: +64-9 815 4200 x 7188
2002 Sep 24
1
Vanishing tick marks using bwplot (lattice)
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major 1
minor 5.1
year 2002
month 06
day 17
language R
I've been using bwplot to plot a dataframe that has a structure like
this:
str(plot.df)
`data.frame': 5952
2002 Oct 04
1
Does the perl language have an equivalent to browser?
I'm asking on this list because the question will be more easily
understood than on a Perl list. Lots of talented people using R also
know how to use Perl, so it's easier to ask them.
Running a Perl script with the -d switch can do some pretty neat
things, but as far as I can tell, every line has to be done
individually. Is there a way to specify where to stop in the way
browser() does
2002 Jul 29
2
How do we omit the coastline using map?
I've been able to do most things I'd like to using the map package,
but I can't figure out how to fill a region without drawing a black
outline. I can draw the outline in a specified colour, but as soon as
I set fill to TRUE, I get a black outline in addition.
My map has a lot of activity on the coastline so I don't want a black
line obscuring what I do there. A pale colour to
2002 Nov 06
1
What else do tab-delimited and csv do differently?
Until recently, I naively believed that a csv was the same as a
tab-delimited file except commas replaced tab characters (provided
of course, the characters being separated weren't either).
Evidently, that is not the case. I tried reading a tab-delimited file
(created in Excel [2000, W2K] and ftp'd as ASCII to linux) into R
platform i686-pc-linux-gnu
arch i686
os
2000 Aug 14
0
Postscript device bug when more than one page is needed?
It would appear to me that there is a bug in the postscript device
when multiple pages are required. Line widths can be only 1 after the
first page. I have not seen reference to this as a problem before.
test.device <-
function()
{
# looks like the problem is in the postscript device
postscript(file = paste("Test.ps"), horizontal = F)
on.exit(dev.off())
2007 Dec 20
0
smoothScatter and geneplotter
On Tue, 18-Dec-2007 at 11:21AM -0500, James W. MacDonald wrote:
|> Duncan Murdoch wrote:
|> > Yes, I agree. (As an aside, there's actually a capital S in
|> > smoothScatter(), and it's a bit of a pain to install, because
|> > geneplotter depends on something that depends on DBI, which is not so
|> > easily available these days.)
|>
|> Somehow I always
2002 Jun 26
1
Output from make check
These are probably very basic questions, but I don't know how one
reads output such as this:
[hrapgc at biomat1 R-1.5.1]$ make check
make[1]: Entering directory `/home/hrapgc/local/R-1.5.1/tests'
make[2]: Entering directory `/home/hrapgc/local/R-1.5.1/tests'
make[3]: Entering directory `/home/hrapgc/local/R-1.5.1/tests/Examples'
make[4]: Entering directory
2002 Jan 22
1
glm.predict?
I've been attempting to calculate the predictions from a poisson glm
object, along these lines:
predict(foo.glm, type = "response")
and
predict(foo.glm, type = "response", se.fit = TRUE)
foo.glm is arrived at this way:
foo.glm <- glm(Insects ~ Dad * Mum + Location, offset = log(MM), family
= "poisson", data = model.df)
There are two
2002 Aug 14
1
Why can't Anove (car package) see the data?
At the end of the 'Details:' section of the help on Anova (car
package), it states:
The standard R `anova' function calculates sequential (type-I)
tests. These rarely test meaningful hypotheses.
So I thought I'd try it. However, I was perplexed to get this
message:
Browse[1]> Anova.glm(leaf.glm1, type ="II")
Error in terms.formula(formula, data = data) :
2002 Oct 28
2
subsetting character vector into groups of numerics
I'm sure there's a simple way to do this, but I can only think of
complicated ones.
I have a number of character vectors that look something like this:
"12 78 23 9 76 43 2 15 41 81 92 5(92 12) (81 78 5 76 9 41) (23 2 15 43)"
I wish to get it into a list of numerical vectors like this:
$Group
[1] 12 78 23 9 76 43 2 15 41 81 92 5
$Subgroup1
[1] 92 12
$Subgroup2
[1] 81 78 5