X = U D V' ## D are the singular values of X
X'X = V D^2 V' ## D^2 are the eigenvalues of X'X
V is the same in both factorizations.
2010/9/22 OndÅej Mikula onmik...@gmail.com
Dear R-helpers,
could anybody explain me briefly what is the difference between
eigenvectors returned by 'eigen'
See the FAQ 7.22 Why do lattice/trellis graphics not work?
The most likely reason is that you forgot to tell R to display the graph.
Lattice functions such as xyplot() create a graph object, but do not display
it (the same is true of
*ggplot2*http://cran.r-project.org/package=ggplot2graphics, and
tmp - data.frame(matrix(rnorm(30), 10, 3,
dimnames=list(letters[1:10],
c(company, person, salary
tmp
company person salary
a -1.04590176 -0.7841855 1.07150503
b -1.06643101 0.6545647 0.43920454
c
Jake,
You can easily use glht for that.
See ?MMC in the HH package for examples.
You may need
install.packages(HH) ## if you don't already have HH.
Specifically, look at the examples showing the use of focus.lmat
On Thu, Sep 16, 2010 at 11:43 AM, Jake Kami jakejk...@gmail.com wrote:
dear
Is this a Windows system behind a firewall?
Can you reach the CRAN mirrors directly from Internet Explorer?
If yes to both, then you need to use
setInternet2(TRUE)
before doing the package installations from the RConsole menu.
See the Windows FAQ and the ?setInternet2 for more details.
On
If you are on windows, and running manually, then
just draw any plot to create the graphics window and then drag the window to the
other monitor.
Rich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
On Windows, use the 8.3 name for file and directory names that have
embedded blanks.
You can discover the 8.3 name directly from R.
system(paste(Sys.getenv(COMSPEC),/c, dir/x c:\\pr*))
Volume in drive C has no label.
Volume Serial Number is B869-26CD
Directory of c:\
02/12/2010 03:52 PM
Also, see R for Windows FAQ 2.16 R can't find my file, but I know it is there!
On Tue, Sep 7, 2010 at 5:01 PM, Oliver Soong osoon...@gmail.com wrote:
Can somebody confirm some unexpected behavior? This is under Windows,
with R 2.11.0 and 2.11.1.
__
Zdeněk,
The way to isolate the problem is to start
emacs with
emacs -q --no-site-files
then manually load ess-site.el
and start R.
We hope that runs normally.
Then add additional lines from your site-start.el
and your .emacs until you find the problem line.
Followup should be on the ESS
Since you are already on emacs, try using ESS. I am pretty sure we do that.
http://ess.r-project.org
Please follow up on the ESS mailing ess-h...@stat.math.ethz.ch
Rich
On Mon, Aug 30, 2010 at 12:48 PM, Barry Rowlingson
b.rowling...@lancaster.ac.uk wrote:
On Mon, Aug 30, 2010 at 5:29 PM, Erik
See
?table
for details.
tmp - ' yx
+ 1 1 2008
+ 2 1 2008
+ 3 0 2008
+ 4 0 2009
+ 5 1 2009'
data - read.table(textConnection(tmp), header=TRUE)
data
yx
1 1 2008
2 1 2008
3 0 2008
4 0 2009
5 1 2009
table(data$x, data$y)
0 1
2008 1 2
2009 1
Kay,
doe this do what you want?
dotplot(y1+y2 ~ facs$Treatment|facs$Sites,
outer=TRUE,
scales = list(x = list(rot = 90, tck=c(1,0))),
ylab=c(y1, y2),
xlab=c(Site 1, Site 2),
strip=FALSE)
On Wed, Aug 25, 2010 at 11:04 AM, Kay Cichini
The multiple y axes are protecting you in this situation.
z - cbind(rnorm(100,c(1,10),1), rnorm(100,c(20,30),1))
dotplot(z[,1]+z[,2] ~ facs$Treatment|facs$Sites,
outer=TRUE,
scales = list(
y = list(
relation=free)),
ylab=c(y1, y2),
Reduce(`+`, k)
On Mon, Aug 23, 2010 at 12:18 PM, Thaler, Thorn, LAUSANNE, Applied
Mathematics thorn.tha...@rdls.nestle.com wrote:
Hi all,
In R it is possible to sum tables:
(a - table(rep(1:3, sample(10,3
1 2 3
2 5 7
a+a
1 2 3
4 10 14
Now suppose that I have a list of
sorry, cancel that. I will try again.
On Mon, Aug 23, 2010 at 12:35 PM, RICHARD M. HEIBERGER r...@temple.eduwrote:
Reduce(`+`, k)
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
Please see this email
https://stat.ethz.ch/pipermail/r-help/2009-March/190541.html
On Sat, Aug 21, 2010 at 11:42 PM, FishLover goldenpuppy_p...@hotmail.comwrote:
Hello
I'm trying to make boxplots
However I'm having issues because my x axis is distance downstream on a
river. The
The question isn't completely clear. I am guessing you want something
like Figure 1.7 or Figure 7.18 in Paul Murrell's book.
library(party)
example(ctree)
Rich
On Sat, Aug 21, 2010 at 3:40 PM, Bernard Leemon bernie.lee...@gmail.comwrote:
I want to make a graph where each element plotted is
cbind(A=x, B=y)
On Sat, Aug 21, 2010 at 6:53 PM, r.ookie r.oo...@live.com wrote:
Is there a way to rename the columns to something like A and B in the cbind
function?
x - rnorm(n = 10, mean = 0, sd = 1)
y - rnorm(n = 10, mean = 0, sd = 1)
cbind(x,y)
Bert,
we can save a lot of time by using paste and then only one call to eval and
parse.
x2 - c(1,2:5, 3:6, 4,8,5:7, 10)
system.time(for (i in 1:100) unlist(lapply(parse(text=x2),eval)))
user system elapsed
0.060.000.03
system.time(for (i in 1:100)
tmp - VariablePARPlot1Plot2
Plot3Plot4
+ ParasiteA3114
+ ParasiteB1235
+ ParasiteC2113
+ ParasiteD2
Yes,
x[order(order(z)),]
Two uses of order are needed, as shown.
Rich
On Sun, Aug 15, 2010 at 11:47 PM, Leon Yee yee.l...@gmail.com wrote:
Dear all,
I have a data frame with several columns, and I have the specified
order of a given column. How can I re-order the rows of my data
Your colleagues was careless with back quotes. See
?'`'
for a discussion of non-syntactic names. This next transcript shows how to
use the
back quote to create and then remove an object with a non-symantic name.
?which
?'`'
`females[\dp\]` - 123
ls(pat=fe)
[1] females[\dp\]
As far as I can tell, the largest number R can take has 308 digits
(1+E308).
This is not a correct statement. The magnitude of the largest 64-bit double
precision floating
point number in R is approximately 1E308. The internal storage in R (and
in most computers today)
contains only the
Please look at RExcel, available through the RExcelInstaller package from
CRAN or
included with a full R download from rcom.univie.ac.at (go to the download
page and take
the most recent RAndFriends).
You can use the supplementary book R through Excel by Erich Neuwirth
result - data.frame(samp=samp, merged=ifelse(is.na(student2), student1,
student2))
result
samp merged
1 A 25.50
2 B 29.40
3 C 22.00
4 D 34.30
5 E 13.60
6 F 10.50
7 G 19.70
8 H 20.22
9 I 17.00
10J 12.90
Rich
[[alternative HTML
From the S+8.1 documentation graphics.pdf
*
uin *information numeric inches per *usr *unit *c(0.73, 0.05)
*
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Michael,
This will get you started. What you are doing with the seven rows isn't
clear
from your description. I made the dates into Date objects. I called your
data
mydata as data is potentially ambiguous.
Rich
mydata - read.table(header=TRUE, textConnection(
ds c1 c2
1
if (FALSE) {
codeline 1
codeline k
}
On Fri, Jul 16, 2010 at 7:59 PM, Joshua Wiley jwiley.ps...@gmail.comwrote:
Hello,
Is there an way to easy comment of sections of code? I was thinking
something along the lines of
\dontrun{
codeline 1
codeline k
}
but that could be
a - rbind(1:3, 4:6, 7:9, 10:12, 13:15)
a1 - a[1:3,]
a2 - a[4:5,]
a1
a2
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Please get a copy of
R for SAS and SPSS Users
*by*
*Muenchen*, Robert A.
http://www.springer.com/statistics/computanional+statistics/book/978-0-387-09417-5
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
You missed FAQ 7.22
7.22 Why do lattice/trellis graphics not work?
The most likely reason is that you forgot to tell R to display the graph.
Lattice functions such as xyplot() create a graph object, but do not display
it (the same is true of *ggplot2* graphics, and Trellis graphics in S-Plus).
Thorn,
All the formulas you listed make sense. Remember that the notation x:y is
an indexing scheme
for dummy variables. It is not a statement about degrees of freedom. The
dummy variables
specified may be linearly dependent on earlier dummy variables in the system
x/y is correctly expanded
Use the timer in tcl. See this email from Phillipe Grosjean with details.
https://stat.ethz.ch/pipermail/r-help/2009-July/203151.html
Rich
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
You need three more lines of code
abline(h=0)
xx - seq(1.96, 4, length=25)
polygon(x=c(xx[1],xx,xx[25]), y=c(0,dnorm(xx),0), col='red')
Please see the normal.and.t.dist function in the HH package for more
detailed control of the graph.
## install.packages(HH) ## if you don't have it yet.
On Thu, Jul 8, 2010 at 3:14 PM, Addi Wei addi...@gmail.com wrote:
Hopefully simple question: What is the best way to name, and treat factor
columns for data that has lots of columns?
This is my column list:
id pID50 D.1 D.2 D.3 D.4 D.5 , etc. all the way to D.185
It would be much better
blah
Error: object 'blah' not found
blah - c(1, 2, 3)
Error in c(1, 2, 3) : 'file' must be a character string or connection
It used to work, and 6 months later I'm back on R, and now it will not
work.
Most likely, you inadvertently masked the base function c with your own
function c
Assuming you are on Windows, you can probably do this within R, without
involving the IT department.
?setInternet2
for the details and discussion. Enter the line
setInternet2(use = TRUE)
before doing the install.packages()
On Wed, Jul 7, 2010 at 7:55 PM, Katie Ashton
Hadley,
The S language modeling language was designed with Wilkinson and
Rogers in mind. The notation was changed from their paper to
retain consistency with the parsing rules for ordinary algebra in
S. I think of : as an indicator of an indexing system into the
dummy variables. It is not an
Michael,
Look at the ancova function in the HH package.
## install.packages(HH) ## if not there yet.
library(HH)
library(nlme) ## for the Orthodont data
ancova(distance ~ age*Sex, data=Orthodont)
ancova(distance ~ age+Sex, data=Orthodont)
ancova(distance ~ age, groups=Sex, data=Orthodont)
You have only one Age value in each Lot. Hence the Age is aliased with one
of the
5 degrees of freedom between lots. You can see it in the picture
xyplot(Conc ~ Age|Lot, data=blots, pch=16)
There is no opportunity for a slope within each of the groups.
Rich
[[alternative HTML version
RExcel on Windows is designed for this task.
See rcom.univie.ac.at for examples including a video.
While you can download the RExcelInstaller package from CRAN, we recommend
installing R with RExcel and other packages it needs with the RAndFriends
installer available
on the download page at
You can normally get through the firewall by using
the internet2 option.
Use ??internet for the exact function name. I am not at my computer now so I
can't check for you.
Rich
__
R-help@r-project.org mailing list
Uwe,
Very nice, thank you.
Rich
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide
Is there a cleaner way of combining two expressions.
This example works and gives what I want
plot(1:10)
aa - expression(alpha==.05)
bb - expression(beta ==.80)
aabb - expression(alpha==.05 ~ , ~ beta ==.80)
text(5, 10, aa)
text(5, 9, bb)
text(5, 8, aabb)
text(5,1,
Seb,
Thanks. That doesn't solve the problem of combining two expressions.
My aa and bb are expressions constructed somewhere else and passed to the
current function which wants to use them together. Your solution moves the
construction
of aa and bb into the function and is equivalent to my
The / is used for nesting and is defined by
A/B == A + (B %in% A)
thus
(a+b)/c == (a+b) + c %in% (a+b) == a + b + a:b:c
Rich
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
ess.r-project.org
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Tukey's pairwise comparisons test is included in the multicomp package.
Also, see the
MMC function in the HH package, available from CSAN.
Rich
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
I will guess that your use of cat without a line end leaves ESS waiting for
a prompt.
cat(b, \n)
would be the way to avoid that.
The way to unfreeze ESS (and emacs more generally) is Ctrl-g
Please put followup ESS questions on the ESS list
ess-h...@stat.math.ethz.ch
Rich
Karl,
The definition and interpretation of contrasts is part of any intermediate
design of
experiments text.
Contrasts for interactions say that the effect of moving
from level 1 of A to level 2 of A depends on the level of B.
I will use notation YAB to indicate the levels of A and B.
For
I think you are asking for M-x ess-transcript-clean-buffer
Save your output file in a file named something.rt which will open
in ess-transcript-mode. It opens read-only and you will need to C-x C-q
to make it writable. Then run
M-x ess-transcript-clean-buffer
On Wed, Jun 2, 2010 at 2:35 PM,
This is a guess since you didn't say where it comes from.
It looks like a data.frame where some rows were deleted.
Also, note that the column names in your example are sorted alphabetically
instead of numerically. This example shows one way to get that behavior.
Rich
tmp -
tmp - matrix(1:24, 6, 4, dimnames=list(letters[1:6], LETTERS[1:4]))
tmp
A B C D
a 1 7 13 19
b 2 8 14 20
c 3 9 15 21
d 4 10 16 22
e 5 11 17 23
f 6 12 18 24
tmp[c(a, c, d, f), ]
A B C D
a 1 7 13 19
c 3 9 15 21
d 4 10 16 22
f 6 12 18 24
See
?`[`
for discussion.
#Error that is returned from above code
..
'C:\Program' is not recognized as an internal or external command,
operable program or
This statement means that the Windows CMD processor does not realize that
'C:\Program Files' is a single word. Use the 8.3 version of the filename
'C:\progra~1'
Please look at the maiz example, the last example in ?MMC in the HH package.
Follow the example all the way to the end. It illustrates a problem and
then the resolution
of the problem.
install.packages(HH) ## if you don't have it yet.
library(HH)
?MMC
Rich
[[alternative HTML version
M.Ribeiro,
You can install on Windows from the tar.gz using the RTools.
You need to set the PATH to find RTools and it looks like you skipped that
step.
Full details are in the R Extensions manual.
Here are my notes, from several years ago as you can see by the R-2.8.0
paths
and the Rtools29.exe
Use lattice.
require(lattice)
?lattice
?xyplot
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
Usually you will want to look at simple effects of the factors when there is
interaction. Please look at the WoodEnergy demos in the HH package.
These examples use glht.
install.packages(HH) ## if you don't currently have HH
library(HH)
demo(MMC.WoodEnergy-aov)
demo(MMC.WoodEnergy)
Rich
Download and install
install.packages(RcmdrPlugin.HH)
library(RcmdrPlugin.HH)
Then there are two options.
The Rcmdr menu item
Statistics Means One-way ANOVA...
has a checkbox for pairwise comparison of means
It uses glht in the multcomp package.
The second option, which I prefer, is to use
format(1.4, nsmall=2)
Rich
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide
Knut,
With the assumption that you are asking about Fisher's LSD for all pairwise
contrasts in ANOVA,
then it is available, along with many superior tests, in the multcomp
package in R.
See ?glht.
Since you mention Excel, please look at RExcel, a package which seamlessly
merges
R and Excel. Any
Here is the full repair for the latex functions in Hmisc to make pdflatex
work in Windows.
This version is still slightly awkward. I hope that Charles and Frank will
smooth it out
and put it in their next release.
I added two new options() and revised show.dvi so it will use them.
Rich
On Thu, May 13, 2010 at 12:46 PM, Erik Iverson er...@ccbr.umn.edu wrote:
Looks good.
Thank you.
I came up with something on my own which was basically to re-write
print.latex. I call pdflatex a few times because with the longtable
package, sometimes things don't line up right until you
At least I thought it did, since I believe show.dvi is an exported
function, can't you just overwrite it in your global environment and
anything that calls it will use your copy? I'll have to look into how I
got your code working on my system, because I know I didn't use that trick!
I
Your example showed up in landscape for me on the pdf file.
I needed to rotate the image in the Adobe Reader.
I like the appearance better with
show.dvi(dvi(x.tex - latex(report,landscape=TRUE), width=8.5,
height=11))
where I changed the paper size in the dvi() command.
Rich
I didn't see that in the Adobe help file. There was some reference in the
Adobe online forum
http://forums.adobe.com/
but the first few hits were queries, not responses to the query. If you
pursue it, please
report success back to the list.
Rich
On Thu, May 13, 2010 at 4:08 PM, Felipe
Tao Shi,
Here is the repair for your query. Use this version of show.dvi temporarily
(until Hmisc includes it in its official version of show.dvi). This
function and
option setting will let the command
latex(x)
display the dvi file in a YAP window on Windows.
This version is not to be used
You can use the bwplot panel in the HH package.
You might need to install HH first, by uncommenting the line below.
girafe.txt - textConnection(
day nmgml
1 29 2.72655
2 29 2.48820
3 15 2.85945
4 15 2.58915
5 15 2.88345
6 15 2.66675
7 47 3.29125
8 15 2.44085
9 29 2.43065
Shi, Tao wrote:
so, I always get the error window (attached) poped out.
The image was stripped by the mailer. Please type the text of the error
into the body of the email.
I tried to include the following into my PATH:
C:\Program Files\MiKTeX 2.7\tex\latex
Where? In the Windows
I recommend the ancova function in the HH package.
install.packages(HH)
library(HH)
example(ancova)
?ancova
For your example,
tmp - textConnection(
day A B
0 10.010.0
7 9.0 9.1
14 8.0 8.2
21 7.0 7.3
28 6.0 6.4
35 5.0 5.5
42
I recommend that you use RExcel to get your files directly from the
xls or xlsx file into R. RExcel correctly transfers dates from Excel to R.
rcom.univie.ac.at for full information.
You can begin with the RExcelInstaller package from CRAN.
Rich
[[alternative HTML version deleted]]
It is exactly the same
tmp - matrix(1:24,6,4)
tmp[4,] - NA
tmp
apply(tmp, 2, sum, na.rm=TRUE)
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
tmp - 'R is a free software environment for statistical computing and
graphics'
strsplit(tmp, )
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
z[yz] - y[yz]
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal,
attr(terms(~B*A), term.labels)
[1] B A B:A
attr(terms(~A/C), term.labels)
[1] A A:C
attr(terms(~B*A/C), term.labels)
[1] B A B:A B:A:C
attr(terms(~(B*A)/C), term.labels)
[1] B A B:A B:A:C
attr(terms(~B*(A/C)), term.labels)
[1] B A A:C B:A B:A:C
The
Santosh,
continuing with your example, I recommend several functions in the HH
package.
## install.packages(HH) ## needed once, if you don't already have it
require(HH)
tmp - data.frame(y=rnorm(100),
category=rep(factor(letters[1:5]), each=20),
Please see the FAQ on R for Windows2.16 where the distinction between the
standard file
separator / and the usual Windows file separator \ is discussed.
On Fri, Apr 23, 2010 at 12:37 PM, weix1 weix1_2...@hotmail.com wrote:
I have string as follows:v:\work\gene
however, i can not write this
Serdal,
You still don't want explicit dummy variables. You have instead indicated
that you need two
factors, one for industryGroup, and one for industry nested within
industryGroup.
The model you suggested
lnQ~lnC+lnM+lnL+lnE+eco_inno+inno+(sum)ind_3d
will be easier to read and interpret as
The spurious condition (ds==2) is actually interpreted as the action taken
by the else statement.
Then the braced plot statement is seen as the next statement after the
conclusion of the
if () {} else ()
statement.
[[alternative HTML version deleted]]
On Wed, Apr 7, 2010 at 9:25 PM, Eric Scott ersco...@illinois.edu wrote:
Thank you for your reply. The WoodEnergy example helped a lot. I
understand now that it is inappropriate to make all pairwise comparisons
with an interaction present and better to make comparisons between levels of
one
Covelli,
From
?rep
Non-integer values of times will be truncated towards zero. If times is a
computed quantity it is prudent to add a small fuzz.
Thus, the example is a very interesting special case of FAQ 7.31,
sprintf(%17.17a,66)
[1] 0x1.08000p+6
sprintf(%17.17a,pr)
[1]
Michael and others,
Here is my complete ancova example
http://astro.ocis.temple.edu/~rmh/HH/hotdog.pdf
This example, especially in Figure 6, places them in a context of a
Cartesian
product of models with the intercept having two levels and slope having
three levels.
It is based on the ancova
## Steve,
## please use the ancova function in the HH package.
install.packages(HH)
library(HH)
## windows.options(record=TRUE)
windows.options(record=TRUE)
# hypothetical data
beak.lgth -
c(2.3,4.2,2.7,3.4,4.2,4.8,1.9,2.2,1.7,2.5,15,16.5,14.7,9.6,8.5,9.1,
Joaquin,
p -levelplot(frame, colorkey=TRUE)
print(update(p, legend=NULL), position=c(0/4,0,1/4,1), more=TRUE)
print(update(p, legend=NULL), position=c(1/4,0,2/4,1), more=TRUE)
print(update(p, legend=NULL), position=c(2/4,0,3/4,1), more=TRUE)
draw.colorkey(p$legend$right$args$key, draw=TRUE,
## this is much easier to read
riffle3 - function(a, b) {
mlab - min(length(a), length(b))
seqmlab - seq(length=mlab)
c(rbind(a[seqmlab], b[seqmlab]), a[-seqmlab], b[-seqmlab])
}
riffle3((1:10),(50:55))
## [1] 1 50 2 51 3 52 4 53 5 54 6 55 7 8 9 10
riffle3((50:55),(1:10))
##
Márcio,
Think matrix!
Do you really want b to be 100 copies of the same numbers?
You are asking for a strange crossproduct with the main diagonal zeroed out.
a2 - crossprod(a)
a2[cbind(1:1000, 1:1000)] - 0
all.equal(a, a2)
Rich
[[alternative HTML version deleted]]
On Thu, Mar 25, 2010 at 10:28 PM, Joaquin Rapela rap...@usc.edu wrote:
Your reply was very useful Rich. Thanks!
I would like to set the height of the viewport in draw.colorkey to match
the
height of the levelplots, as it is done when I add colorkeys to all the
levelplots. The height of the
xyplot(1:10 ~ 1:10, ylab=list(expression(atop(PM[2.5] ~ Concentration ~
frac(mu*g,m^3),(random text in brackets, crt=90)
See ?plotmath for atop and many more functions
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing
PL1 used semi-colons, therefore so does SAS (Yes, SAS was originally written
in PL1).
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Please look at the Error term in ?aov
On Wed, Mar 17, 2010 at 4:25 AM, Galanidis Alexandros a...@env.aegean.grwrote:
Greetings to all,
This is my model: aov.fit-aov(Y~A+B+C+D+E+A:C+A:E)
In summary(aov.fit) all F values are comptuted by eg MS(A)/MS(Residuals).
This is not correct (or what I
Please look at the latex function in the Hmisc package.
The default display is very good. And there are many optional
arguments that give you very fine control over the appearance of
the table.
Rich
[[alternative HTML version deleted]]
__
Why use a csv dataset as an intermediary? Use RExcel and get
the dataset directly from your Excel source.
See http://rcom.univie.ac.at for full details.
You can download the RExcelInstaller package from CRAN with
install.packages(RExcelInstaller)
library(RExcelInstaller)
installRExcel()
Rich
barchart(yield ~ year, data = barley,
stack=TRUE,
col=black,
ylim=c(0,85),
scales=list(x=list(labels=c(expression(atop(March, 1932)),
expression(atop(April, 1932)
)
This gets most of what you want.
[[alternative HTML version deleted]]
Please see the maiz example in ?MMC in the HH package.
maiz is the last example in the help file. Keep going all the way to the
end of
the help file. See also the
demo(MMC.WoodEnergy-aov, HH)
These examples show how to use glht in the presence of interactions and
covariates.
Rich
In addtition to the example I mentioned previously,
demo(MMC.WoodEnergy-aov, HH)
Please also see
demo(MMC.WoodEnergy, HH)
In this example, since anova(energy.aov.4),
shows that the Wood factor and Stove:Wood interaction are significant,
all possible pairwise comparisons of the 12 Stove:Wood
I have several responses, not necessarily consistent with each other.
1. An appendix should be a user manual with screen shots.
2. The archival version of the package should be Yourpackage_x.y-z.tar.gz
stored
at your university's archival web location. The actual usable version of
your package
## It looks like something like this will meet your needs
## Rich
tmp - data.frame(trat=letters[1:3], y=rnorm(3), x=1:3)
xyplot(y ~ x | trat, data=tmp)
strip.levelnames - c(expression(H[2]*O), expression(x^2+y^2),
expression(fraction(x^3, 1+x^3)))
xyplot(y ~ x | trat, data=tmp,
## Sorry, typo in my last email. frac is the correct spelling, not
fraction
## RIch
tmp - data.frame(trat=letters[1:3], y=rnorm(3), x=1:3)
xyplot(y ~ x | trat, data=tmp)
strip.levelnames - c(expression(H[2]*O), expression(x^2+y^2),
expression(frac(x^3, 1+x^3)))
xyplot(y ~ x | trat, data=tmp,
It looks like you need something like this. The data you attached didn't
make
it through the email system.
## you might need
install.packages(HH) ## do this once
library(HH) ## do this every time you start a new R session
interaction2wt(y ~ Block + CO2 + temp + moisture)
## the rest is
On Sat, Mar 13, 2010 at 11:51 PM, Megh megh700...@yahoo.com wrote:
Hi, I was trying to run your code in my VBA editor, however could not
succeed. The execution stumbled in the line Call Rinterface.StartRServer
itself. I have RCOM package installed into my R environment. Do I need to
install
601 - 700 of 843 matches
Mail list logo
