, and # brief margin notes
I find it useful, but all things change over time.
Steven McKinney
> -Original Message-
> From: ESS-help [mailto:ess-help-boun...@r-project.org] On Behalf Of Paul
> Johnson
> Sent: May-04-17 7:50 PM
> To: ess-help
> Subject: Re: [
/datasets/data/morley.tab
mm - read.table(filepath)
head(mm)
Expt Run Speed
0011 1 850
0021 2 740
0031 3 900
0041 4 1070
0051 5 930
0061 6 850
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Victor Gabillon
Sent: April-28-11 8:22 PM
To: r-help@r-project.org
Subject: [R] Change the text size of the title in a legend of a R plot.
Hello,
Is it possible to
of, though the multcomp tables and plots yield logically
equivalent results and plots.
Writing a few lines of code to accomplish your graph is fairly straightforward.
HTH
Steven McKinney
THANKS !
version
platform x86_64-apple-darwin9.8.0
arch x86_64
os
operating system you are using, how much RAM you have
or what sessionInfo() reports on your machine. That information will help to
figure this out.
Steven McKinney
From: Jim Silverton [jim.silver...@gmail.com]
Sent: April 9, 2011 9:21 AM
To: Steven McKinney
?
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Jim Silverton [jim.silver...@gmail.com]
Sent: April 8, 2011 9:43 AM
To: r-help@r-project.org
Subject: Re: [R] Fast version of Fisher's Exact Test
Is anyone aware
to the registry
If you install tcltk2, is tcltk then available?
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Rolf Turner [r.tur...@auckland.ac.nz]
Sent: April 7, 2011 5:15 PM
To: r-help@r-project.org
?
Steven McKinney
From: Rolf Turner [r.tur...@auckland.ac.nz]
Sent: April 7, 2011 5:38 PM
To: Steven McKinney
Cc: r-help@r-project.org
Subject: Re: [R] Where is the tcltk package?
On 08/04/11 12:31, Steven McKinney wrote:
Searching for tcltk on the BioConductor
-Original Message-
From: stephen sefick [mailto:ssef...@gmail.com]
Sent: April-03-11 5:35 PM
To: Steven McKinney
Cc: R help
Subject: Re: [R] Linear Model with curve fitting parameter?
Steven:
You are exactly right sorry I was confused
-Original Message-
From: stephen sefick [mailto:ssef...@gmail.com]
Sent: April-04-11 2:49 PM
To: Steven McKinney
Subject: Re: [R] Linear Model with curve fitting parameter?
Steven:
I am really sorry for my confusion. I hope this now makes sense.
b0 == y intercept == y
similarly. Don't
forget to adjust for multiple comparisons...
HTH
Steve
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org
-Original Message-
From: stephen sefick [mailto:ssef...@gmail.com]
Sent: April-01-11 5:44 AM
To: Steven McKinney
Cc: R help
Subject: Re: [R] Linear Model with curve fitting parameter?
Setting Z=Q-A would be the incorrect dimensions. I could Z=Q/A.
I suspect this is confusion
data.
Fitting lm(Z ~ R + S, data = x) should yield
intercept parameter estimate = estimate for log(K)
R coefficient parameter estimate = estimate for r
S coefficient parameter estimate = estimate for s
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia
Hi Cristoph,
Glad to hear it worked - replying to r-help in case others face this issue.
Cheers
Steven McKinney
From: Martin Knapp [mkna...@aucklanduni.ac.nz]
Sent: March 28, 2011 7:28 PM
To: Steven McKinney
Subject: Re: [R] xlsx problem
Thanks
anything, just about anywhere.
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Martin Knapp [mkna...@aucklanduni.ac.nz]
Sent: March 28, 2011 1:43 PM
To: r-help
Subject: [R] xlsx problem
Dear list,
I'm running
c4Z
56d4Z
57a5Z
58b5Z
59c5Z
60d5Z
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
BORGMANN,Kathi [kborgm...@audubon.org]
Sent: March 28
Here's a start...
require(maps)
foo - map(state, new york)
lines(x = range(foo$x, na.rm = TRUE), y = range(foo$y, na.rm = TRUE))
You can figure out from this how to specify the coordinates that you
want for dividing up the map, put them in a file etc.
Steven McKinney
How about this?
x-c(1:5,NA,NA,8:10)
y-1:10
plot(0,0,xlim=c(0,10), ylim=c(0,10),type=n,main=Dont show the bloody 0
values!!)
lines(x~y, col=blue, lwd=2, subset = !is.na(x))
NAs let you do lots of useful manipulations in R.
Steven McKinney
From: r
(Control Panel - Administrative
Tools
on Windows 7). Set up the R script you want to run daily, and invoke it
with the R CMD BATCH command.
In Unix, you can use the cron scheduler to do the same.
Steven McKinney
Thanks for any input.
--
View this message in context:
http://r
, linfct = mcp(Provenancef = Tukey)))
Error in `[.data.frame`(mf, nhypo[checknm]) : undefined columns selected)
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Steven McKinney [smckin...@bccrc.ca]
Sent
Comments in-line below
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Lilith
Sent: December-02-10 9:39 AM
To: r-help@r-project.org
Subject: [R] Tukey Test, lme, error: less than two groups
Dear R-group,
I am trying
data, use beside=FALSE to stack, etc.
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Par Leijonhufvud [...@hunter-gatherer.org]
Sent: October 28, 2010 8:37 PM
To: rhelp
Subject: [R] wilcox.test; data type
to the database, and running data queries.
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
lord12 [trexi...@yahoo.com]
Sent: October 19, 2010 11:42 AM
To: r-help@r-project.org
Subject: [R] How to write to sqlite
3.8 3.9 4
[4,]4 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9 5
[5,]5 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 6
[6,]6 6.1 6.2 6.3 6.4 6.5 6.6 6.7 6.8 6.9 7
HTH
Steven McKinney
From: r-help-boun...@r
This is better suited for R-help than R-devel, so I'm copying to the R-help
list:
-Original Message-
From: r-devel-boun...@r-project.org [mailto:r-devel-boun...@r-project.org] On
Behalf Of Martin Kerr
Sent: October-08-10 3:09 AM
To: r-de...@r-project.org
Subject: [Rd] Selecting
For those that want it all...
{cat(?); a-readLines(n=1)
+ print(hey)
+ print(b-paste(t,a,sep=))}
?ada
[1] hey
[1] tada
b
[1] tada
Steven McKinney
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Peter Dalgaard
Sent
()) ...
So see if you can get package lmtest.
HTH
Cheers
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
robm [r.malp...@ntlworld.com]
Sent: September 24, 2010 6:53 AM
To: r-help@r-project.org
Subject: [R
} for a-aigu).
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Jennifer Young [jennifer.yo
out
strategies to share bits of data and information across multiple R sessions.
HTH
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Olga Lyashevska [o...@herenstraat.nl]
Sent: August 3, 2010 8:04 AM
To: r
(edu)
[,1] [,2]
[1,] 0.26550871
[2,] 0.37212391
[3,] 0.57285340
[4,] 0.90820780
[5,] 0.20168191
[6,] 0.89838970
Steven McKinney
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Hey Sky
Does function
showNonASCII(x)
in package tools do what you need?
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org
Quebec
Quebec0 42
Mississippi 42 0
quartz() ## Your plot window type here (quartz() works on a Mac)
bwplot( conc ~ Type : Treatment, data = CO2 )
quartz()
bwplot( conc ~ Type2 : Treatment, data = CO2 )
Steven McKinney
Statistician
Molecular Oncology
ODBC Data Source Administrator widget in
Windows administrative tools.
Install RMySQL, configure and enjoy.
Good luck!
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
neatgadgets [neatgadg...@gmail.com
.
What OS are you running?
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
neatgadgets [neatgadg
I think this is a question for R-devel so I'm
cross-posting there with apologies.
I've just acquired a Windows 7 64-bit box and
also will need RMySQL eventually.
Is there any information about issues involved
with compiling RMySQL for Windows 64-bit?
Steven McKinney
Statistician
Molecular
?cast
A reproducible example would get you more feedback.
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
n.via...@libero.it [n.via...@libero.it]
Sent: June 18, 2010 6:39 AM
To: r-help@r-project.org
the matrix to a dataframe after the loop completes if needed.
HTH
Steven McKinney
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf Of
Dimitri Liakhovitski [ld7...@gmail.com]
Sent: March 26, 2010 6:40 PM
To: Bert Gunter
Cc: r-help
From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] On Behalf
Of duncandonutz [dwads...@unm.edu]
Sent: March 19, 2010 1:11 PM
To: r-help@r-project.org
Subject: [R] Merging Matrices
I have two symmetric matrices, but of different
) is a smooth function of X.
Hence you might want to rethink why you'd want a
factor variable as a predictor variable in a GAM.
This is why the gam machinery doesn't just do the
factor conversion to indicator variables as is done in
lm.
HTH
Steven McKinney
From: r
email: dmhul...@metstat.com
web: http://www.metstat.com
-
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney -at- bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular
row of coefficients in the parameter estimates matrix mA
You can spot check with e.g.
mA[1,] %*% mB[15,]
[,1]
[1,] 214.7864
mA[1,] %*% mB[5,]
[,1]
[1,] 27.31023
mA[1,]
[1] 0.2875775 0.9568333 0.8895393
mB[5,]
x
1 5 25
HTH
Steven McKinney, Ph.D.
Statistician
?NumericConstants
will bring up a help page that mentions
All other numeric constants start with a digit or period and are either a
decimal or hexadecimal constant optionally followed by L.
and
An numeric constant immediately followed by L is regarded as an integer number
when possible
-Original Message-
From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Sent: Monday, November 16, 2009 4:52 PM
To: Duncan Murdoch
Cc: Steven McKinney; R Help
Subject: Re: [R] Where are usages like == 2L documented?
On Mon, Nov 16, 2009 at 7:25 PM, Duncan Murdoch murd
, IV1, IV2, error.bars=se, ylim = c(27, 99))
or whatever your desired ylim values are.
HTH
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
___
×ñçóéìïðïéåßôå Yahoo!;
ÂáñåèÞêáôå ôá åíï÷ëçôéêÜ ìçíýìáôá (spam); Ôï Yahoo! Mail
äéáèÝôåé ôçí êáëýôåñç
# Same as lapply loop result
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
I researched on this, but I can't find a solution.
I would really appreciate any help.
P.S. Excuse my lack of terminology :).
Iurie
some key point? Any reasons why this
would be a Bad Idea?
Best
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney -at- bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor
)
(but this does not work)?
Any info appreciated
Best
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hi Mary,
Something such as
plot(1,1)
title(expression(paste(LBAuo = -, infinity)))
Is this what you're after?
Best
Steve McKinney
From: Mary A. Marion [mms...@comcast.net]
Sent: July 23, 2009 4:23 PM
To: Steven McKinney
Subject: Re: [R] mathematical
-Original Message-
From: Marc Schwartz [mailto:marc_schwa...@me.com]
Sent: Thursday, July 23, 2009 2:04 PM
To: Steven McKinney
Cc: R-help@r-project.org
Subject: Re: [R] Navigate to Index page of a package from R command
prompt
On Jul 23, 2009, at 3:30 PM, Steven McKinney wrote
Does this approach what you're looking for?
See
?is.finite
for more info
LBAuo- -
LBAuo
[1] -
if (LBAuo = -) LBAuo - -Inf
LBAuo
[1] -Inf
HTH
Steven McKinney
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf
foo.bar.baz myrbind
foo.bar.baz
X1 X2 X3
1 1 5 9
2 2 6 10
3 3 7 11
4 4 8 12
5 101 105 109
6 102 106 110
7 103 107 111
8 104 108 112
HTH
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Warning message:
In anova.lm(lm(vtot ~ fac)) :
ANOVA F-tests on an essentially perfect fit are unreliable
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research
Hi John,
I can't reproduce your case with built-in data sets.
Your data set is small. Can you show the data for the
variables involved in your example?
(Time30, Died, Rx, Age)
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer
graphics grDevices utils datasets methods
[8] base
other attached packages:
[1] survival_2.35-4
Steven McKinney
-Original Message-
From: John Sorkin [mailto:jsor...@grecc.umaryland.edu]
Sent: Monday, May 11, 2009 8:53 PM
To: Steven McKinney; r-h
graphics grDevices utils datasets methods base
other attached packages:
[1] nlme_3.1-90
loaded via a namespace (and not attached):
[1] grid_2.9.0 lattice_0.17-22 tools_2.9.0
Also occurs on Windows box with R 2.8.1
Steven McKinney
Statistician
Molecular Oncology
Thanks Duncan,
Comments and a proposed bug fix in-line below:
-Original Message-
From: Duncan Murdoch [mailto:murd...@stats.uwo.ca]
Sent: Wednesday, April 29, 2009 5:10 PM
To: Steven McKinney
Cc: R-help@r-project.org
Subject: Re: [R] Corrupt data frame construction - bug
-
guide.html
and provide commented, minimal, self-contained, reproducible code.
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckin...@bccrc.ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th
#Here my first try [error message on the line within the loop,
saying
something like:
# 'recursive indexing on level 2 failed']
sites_object_list - vector(list,99)
If there's a site number larger than 99 this will be problematic.
Sorry, this isn't correct - R will just increase
.
The functions as.single and single are identical to as.double and double
except they set
the attribute Csingle that is used in the .C and .Fortran interface, and they
are intended
only to be used in that context.
HTH
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
by the
following rule. The environment variables
TMPDIR, TMP and TEMP are checked in turn
and the first found which points to a
writable directory is used: if none
succeeds /tmp is used.
HTH
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer
$a - as.Date(c$a)
c
a
1 2008-07-27
2 2008-10-01
3 2008-08-15
4 2008-08-14
5 2008-08-14
6 2008-09-20
7 2008-07-27
8 2008-10-01
lapply(c, class)
$a
[1] Date
Steven McKinney, Ph.D.
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email
Hi Jason
Have you checked out Stat/Transfer?
www.stattransfer.com
They state they handle SYSTAT files, and
SigmaPlot is a SYSTAT product.
Stat/Transfer has a demo download you
could use to test.
I've used it to good effect (though for
other than SYSTAT datasets) in the past.
HTH
Steven
lme4_0.99875-9
[5] nlme_3.1-89
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C.
V5Z 1L3
$wt = 150.0, age]
Apologies and thanks for setting me straight.
Best
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave
Use the is.na() function to assign
NA values:
is.na(A) - !B
A
[,1] [,2] [,3]
[1,]3 NA NA
[2,]333
[3,]33 NA
C - matrix(c(3,3,3,NA,3,3,NA,3,NA),3,3)
all.equal(A, C)
[1] TRUE
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
A - list(c(3,2,3), c(1,2,3,4), c(5,6))
matrixFromList(A)
[,1] [,2] [,3] [,4]
[1,]323 NA
[2,]1234
[3,]56 NA NA
class(matrixFromList(A))
[1] matrix
Best
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia
See also the dist() function documentation.
If you use indexing as described in ?dist
it is straightforward to maintain and
use a vector of the distances.
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney
This may just be version incompatibilities.
A similar discussion recently transpired on
r-sig-mixed-models, see e.g.
https://stat.ethz.ch/pipermail/r-sig-mixed-models/2008q4/001526.html
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer
a chromosome graphic to the
bottom of data plots of genomic data.
You can see examples in the vignette
PDFs.
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: [EMAIL PROTECTED]
tel: 604-675-8000 x7561
BCCRC
Molecular
.
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: [EMAIL PROTECTED]
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C.
V5Z 1L3
Canada
__
R-help
- if(scale 0) RSS/scale - n else n * log(RSS/n)
c(edf, dev + k * edf)
}
extractAIC.aov - extractAIC.lm
extractAIC.negbin - function(fit, scale, k = 2, ...)
{
n - length(fit$residuals)
edf - n - fit$df.residual
c(edf, -fit$twologlik + k * edf)
}
HTH
Steven McKinney
Statistician
AIC(ll(object), k = k)
}
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C.
V5Z 1L3
Canada
(Apologies if this repost is a duplicate,
my first submission did not appear to
make it through.)
-Original Message-
From: Steven McKinney
Sent: Mon 11/24/2008 4:28 PM
To: 'Blanchette, Marco'; R-help
Subject: RE: [R] 64bit R for Mac
Dear Marco
Check out
http://r.research.att.com
- letters[4:6]
x - intersect(y,z)
class(x)
[1] character
length(x)
[1] 0
y - 1:3
z - 4:6
x - intersect(y,z)
class(x)
[1] integer
length(x)
[1] 0
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney
is working well.
You might also want to follow the
r-sig-mac list for mac-specific
information.
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: [EMAIL PROTECTED]
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West
as to the best way to
handle model assessment, not all classical summary
statistics and associated p-values map to lmer models
in a straightforward fashion.
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc
should be
done and the error thrown. Otherwise, the data objects
y and x that you set up should have been passed downwards
in some fashion for evaluation. R-core members who know
the rules better than I will have to determine how best to
handle this one.
HTH
Steven McKinney
Statistician
in the lower left, not top-down as in
a table. )
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C
greater than zero,
so that's not an issue in this calculation (this
'B' not being the same as the 'B' in Amy's original
description).
Thanks again to Duncan Murdoch, and my apologies
for my misinterpretation.
Steven McKinney
-Original Message-
From: Robin Hankin [mailto:[EMAIL PROTECTED
enough to give your labels enough room.
HTH
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C.
V5Z
then they need
to be set aside and investigated
separately, etc.
Steven McKinney
Statistician
Molecular Oncology and Breast Cancer Program
British Columbia Cancer Research Centre
email: smckinney +at+ bccrc +dot+ ca
tel: 604-675-8000 x7561
BCCRC
Molecular Oncology
675 West 10th Ave, Floor 4
Vancouver B.C
Your column of data has some
character data in it, perhaps
an Excel #VALUE! or a blank or
some such entry not strictly numeric.
When R reads in such a column, it
makes that column variable into
a 'factor' variable instead of a
numeric variable, because the values
are not all numeric.
You can
Here's the current method for printing a data frame:
print.data.frame
function (x, ..., digits = NULL, quote = FALSE, right = TRUE,
row.names = TRUE)
{
n - length(row.names(x))
if (length(x) == 0L) {
cat(NULL data frame with, n, rows\n)
}
else if (n == 0L) {
Is something missing in the melt()?
x-data.frame(x1=c(1,7),x2=c(4,6),x3=c(8,2))
require(reshape)
Loading required package: reshape
dfm - melt(x, id = c())
Error in if (!missing(id.var) !(id.var %in% varnames)) { :
missing value where TRUE/FALSE needed
dfm[order(dfm$value), ]
Error: object
One thing that R does very well that SAS does not is
graphics - graphical portrayal of data is important,
and you can keep up with R by supplementing your
SAS analyses with R graphics.
Steve McKinney
-Original Message-
From: [EMAIL PROTECTED] on behalf of Barry Rowlingson
Sent: Fri
Not likely that anyone can explain, as
there is not enough information in your
email.
Including the contents of the freqtest.txt file
was a good idea, as the posting guide suggests
(the posting guide is that clearly labeled bit
at the bottom that looks like this:
PLEASE do read the posting guide
Hello,
Your theta() function is returning different
sets of coefficients depending on the results of
step().
You'll need to add code to theta() to figure
out which variables were selected, and store
them into the right positions of a vector
of length 20 (the apparent number of covariates
you
This is close, but maybe not optimal lattice coding.
I haven't yet figured out how to suppress the x axis
labeling.
bwplot(yield ~ 1|year, panel = function(x, y, ...){panel.bwplot(x, y, ..., pch
= |); panel.points(x, mean(y), ..., pch=17)}, data = barley, horizontal =
FALSE, xlab = )
Steve
Hi Ed,
Here's a simple example showing your needs:
myfun - function(n1, n2, n3) {
mat1 - matrix(rep(1), nrow = n1, ncol = 3)
mat2 - matrix(rep(2), nrow = n2, ncol = 4)
mat3 - matrix(rep(3), nrow = n3, ncol = 5)
require(survival) ## make sure the package you need is loaded
mypkgfun
Did you try it with the vector '' and operator?
d-sapply(res,function(.df){(.df$TimesVisited[.df$Tick912 .df$Id0])})
(The '' operator is designed for use in
e.g. if() clauses where you want a scalar logical
answer)
HTH
Steve McKinney
-Original Message-
From: [EMAIL PROTECTED] on
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
On Behalf Of [EMAIL PROTECTED]
Sent: Thursday, August 14, 2008 2:00 PM
To: r-help@r-project.org
Subject: [R] Detecting duplicate values
This is another how do I do it type of question. It seems that for a
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
On Behalf Of Altaweel, Mark R.
Sent: Wednesday, August 13, 2008 3:03 PM
To: r-help@r-project.org
Subject: [R] Conditional statement used in sapply()
Hi,
I have data stored in a list that I would like to
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
On Behalf Of rcoder
Sent: Wednesday, August 13, 2008 3:17 PM
To: r-help@r-project.org
Subject: [R] missing TRUE/FALSE error in conditional construct
Hi everyone,
I posted something similar to this in reply
The help page for is.na()
is worth reading repeatedly.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
On Behalf Of stephen sefick
Sent: Tuesday, August 12, 2008 6:53 PM
To: Charles C. Berry
Cc: Mike Prager; [EMAIL PROTECTED]
Subject: Re: [R] ignoring
Hi useRs
I'm trying to figure out how to source
an R script file straight from a subversion
repository, without having to put a copy
of the script into the local working directory.
Has anyone done this?
Something such as
source(file = paste(svn://myrepo.xxx.org/opt/svn/repos/,
Steve McKinney
-Original Message-
From: Marc Schwartz [mailto:[EMAIL PROTECTED]
Sent: Fri 8/1/2008 5:42 PM
To: Steven McKinney
Cc: r-help@r-project.org
Subject: Re: [R] source a script file straight from a subversion repository
on 08/01/2008 06:49 PM Steven McKinney wrote:
Hi
I came up with the same solution as Mark Leeds.
a1roworders - t(apply(a1, 1, order))
a2ord - t(sapply(seq(nrow(a1)), function(x) a2[x, a1roworders[x,]] ))
a2ord
[,1] [,2] [,3]
[1,] 102 101 103
[2,] 102 101 103
[3,] 103 101 102
[4,] 101 102 103
Mark's question about the
Hello
I haven't found better instructions, it's just
not an easy thing to do.
You might also consider joining the r-sig-mac group
and reviewing threads there for additional information.
Rather than try to configure, make and install with
one giant command, I'd suggest breaking the task down
Have you reviewed the old
xlisp machinery that used to
be in R? Check out the RXLisp
library at
http://www.omegahat.org/RXLisp/
and see if this will work. The
PDF
http://www.omegahat.org/RXLisp/examples.pdf
reviews calling xlisp from R.
HTH
Steve McKinney
-Original Message-
From:
1 - 100 of 119 matches
Mail list logo