Hi,
I tried to use rmnlIndepMetrop (bayesm package) for my MNL model with 4
choice alternatives, 5 independent variables, 69 observations,
dim(X) [1] 276 5, nu=6. So I run such code:
if(nchar(Sys.getenv(LONG_TEST)) != 0) {R=2000} else {R=10}
set.seed(66)
Le 08.02.2006 04:21, Taka Matzmoto a écrit :
Hi R users
This looks a simple question
Is there any difference between between rnorm(1000,0,1) and running
rnorm(500,0,1) twice in terms of outcome ?
TM
Not here :
R set.seed(1)
R x - rnorm(1000, 0, 1)
R set.seed(1)
R y - rnorm(500, 0, 1)
R z
Dear Colleagues,
I would like to only inspect rules that contain a certain label substring on
the rhs. In this special case the item labels are built like this:
itemtype_itemvalue e.g. Artikelgruppe_E0815 what I want to do is only
show rules where Artikelgruppe is contained in the rhs - has
Romain Francois wrote:
Le 08.02.2006 04:21, Taka Matzmoto a écrit :
Hi R users
This looks a simple question
Is there any difference between between rnorm(1000,0,1) and running
rnorm(500,0,1) twice in terms of outcome ?
TM
Not here :
R set.seed(1)
R x - rnorm(1000, 0, 1)
R
Dear Collegues,
does anybody know how I can place a legend at a sensible position in a star
diagram automatically? The legend should tell the name of the variable (from
colnames) for every color used.
Sincerely, Markus
___
Markus Preisetanz
Consultant
Client
Hi,
How can I install the server version of rpad on a website. I have read
the instructions in http://www.rpad.org/Rpad/ServerNotes.html, but
could not figure out how. Do I have to have something like Apache on
the website? Is there any instructions in addition to those on
Hi list,
I am fitting microarray data (intensity) model using the lme package in R
environment. I have 5 fixed variables in the model. One of the fixed
variables is genes. I am trying to get p-values for different genes. But I
am getting only one p-value for all genes together. I can get a list
Why don't you test it yourself?
E.g.,
set.seed(42)
bob1 - rnorm(1000,0,1)
set.seed(42)
bob2 - rnorm(500,0,1)
bob3 - rnorm(500,0,1)
identical(bob1, c(bob2, bob3))
I won't tell you the answer. :-)
--
Bjørn-Helge Mevik
__
R-help@stat.math.ethz.ch
Hello,
I'm trying to draw a simplified tree of a cluster object
with maptree, but I get an error message:
class(vars.agn)
[1] agnes twins
draw.clust (clip.clust (vars.agn, k=5))
Error in names-.default(`*tmp*`, value = c(merge,
height, order,'names' attribute
Hi R-Experts,
Can anyone tell me how Bloomberg data can be directly downloaded to R?
Is there any package?
Sumanta Basak.
---
This e-mail may contain confidential and/or
Hello list.
Is it possible to use par(srt=45) to rotate text by 45 degrees along the
x-axis of a plot. Using:
code
x_names-c(C57 Nv, C57 Vacc, 129 Nv, 129 Vacc, IFNgR Nv, IFNgR Vacc)
par(srt=45)
mtext(font=2, x_names, side=1, line=1, at=l, cex=1.2)
par(srt=0)
/code
doesn't seem to work in R
Hello!
Could anyone please tell me how should I include R in a text section for
References?
Regards,
Sara Mouro
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
citation()
2006/2/8, Sara Mouro [EMAIL PROTECTED]:
Hello!
Could anyone please tell me how should I include R in a text section for
References?
Regards,
Sara Mouro
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch
Hi Sara.
From the R faq here:
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Citing-R
2.8 Citing R
To cite R in publications, use
@Manual{,
title= {R: A Language and Environment
for Statistical
Computing},
author = {{R Development Core
Dear All,
I have to enter many lines of data in the same object.
I usually use copy-paste to transfer data from an Word file to R.
But, for large lines of data, R gets confused and gives an error message,
i.e. it breaks one line somewhere, and lines get no meaning at all.
Some times I
Sara Mouro wrote:
Dear All,
I have to enter many lines of data in the same object.
I usually use copy-paste to transfer data from an Word file to R.
What is the best way to do that?
Use 'Save As' to save your Word file - or rather just the data section
- as a plain text or
On 2/8/06 6:55 AM, Sara Mouro [EMAIL PROTECTED] wrote:
Dear All,
I have to enter many lines of data in the same object.
I usually use copy-paste to transfer data from an Word file to R.
But, for large lines of data, R gets confused and gives an error message,
i.e. it breaks one
I cleaned up your function a bit but please double check
generate.matrix - function(nr, runs=5){
h - nr/2## half of nr
nc - nr/10 + 1
mat - matrix(0, nr, nc) ## initialize
mat[ ,1] - c( rep(1, h), rnorm(h) ) ## 1st
How does the data look and how are you storing in R (e.g. matrix, list)?
I think this an issue related to Word where it is using either unequal
spaces or different carriage returns. I would not recommend storing
data, especially numerical ones in the form of a matrix, in Word files.
I would
Iain Gallagher wrote:
Hello list.
Is it possible to use par(srt=45) to rotate text by 45 degrees along the
x-axis of a plot. Using:
code
x_names-c(C57 Nv, C57 Vacc, 129 Nv, 129 Vacc, IFNgR Nv, IFNgR
Vacc)
par(srt=45)
mtext(font=2, x_names, side=1, line=1, at=l, cex=1.2)
Please read the posting
1) I think BioConductor mailing list might be better as some of these
could be implemented via LIMMA (I believe)
2) Provide sufficient information and perhaps a simple example.
Regards, Adai
On Wed, 2006-02-08 at 10:42 +0100, Mahdi Osman wrote:
Hi list,
I am
Dear R helper,
I hope that u can help me to sort out my problem
because I sent an E-mail last night to R-list but I
have not receive any help and at the same time I think
this problem is not so hard.
I have used the following functions before
K-10
prime-c(2,3,5,7,11,13,17)
No Excel attachment came through.
Just taking a guess here but there seems to be very little variation the
columns V10 till column V23.
BTW, can you not issue the following call :
mydata[ , 1:7] ~ mydata[ , 8] + mydata[ ,9]
instead of creating y1, y2, ... separately then cbind-ing them ?
On 2/8/2006 4:53 AM, Bjørn-Helge Mevik wrote:
Why don't you test it yourself?
E.g.,
set.seed(42)
bob1 - rnorm(1000,0,1)
set.seed(42)
bob2 - rnorm(500,0,1)
bob3 - rnorm(500,0,1)
identical(bob1, c(bob2, bob3))
I won't tell you the answer. :-)
This isn't really something that can be
Reza,
Yes, you will need Apache and Perl already installed, then follow the
install notes you linked to above. I tried for a long time to get it to
work with Microsoft IIS and could not for the life of me, but when I tried
it with Apache 2.0 it worked pretty easily. I have never set up a web
Paul:
It is a little difficult to understand what you are trying to translate
since you do not show what the model would look like using lme. If you
show lme, then it is easy to translate into lmer syntax.
A few thoughts, first, use lmer in the Matrix package and not in lme4.
Second, see the
Hello Sumanto,
your question might be more appropriately been posted to the R-sig-finance
list:
[EMAIL PROTECTED] (I have cc'ed this mail to this list).
To my knowledge neither a function nor a CRAN-package does exist. However,
on the last useR! conference Dirk Edelbuettel presented a
Hello all,
In a previous posting regarding glmm.admb it is stated that glmm.admb
can handle 2 nested random effects. I can only fit a single random
term at the moment, and wondered if anyone could provide me with some
information on how to specify a model with 2 (nested or
On Wed, 2006-02-08 at 10:33 +, Iain Gallagher wrote:
Hello list.
Is it possible to use par(srt=45) to rotate text by 45 degrees along the
x-axis of a plot. Using:
code
x_names-c(C57 Nv, C57 Vacc, 129 Nv, 129 Vacc, IFNgR Nv, IFNgR
Vacc)
par(srt=45)
mtext(font=2, x_names, side=1,
On 08-Feb-06 Duncan Murdoch wrote:
On 2/8/2006 4:53 AM, Bjørn-Helge Mevik wrote:
Why don't you test it yourself?
E.g.,
set.seed(42)
bob1 - rnorm(1000,0,1)
set.seed(42)
bob2 - rnorm(500,0,1)
bob3 - rnorm(500,0,1)
identical(bob1, c(bob2, bob3))
I won't tell you the answer. :-)
We did get your previous message. See the posting guide about what to do
if you do not get an answer.
See MASS4, chapter 16 (and see the FAQ for what MASS4 is) for fitting
mixtures of normals and testing their fit.
I want to know how to read the data to txt file.
What does that mean? (You
On Wed, 8 Feb 2006, Duncan Murdoch wrote:
On 2/8/2006 4:53 AM, Bj?rn-Helge Mevik wrote:
Why don't you test it yourself?
E.g.,
set.seed(42)
bob1 - rnorm(1000,0,1)
set.seed(42)
bob2 - rnorm(500,0,1)
bob3 - rnorm(500,0,1)
identical(bob1, c(bob2, bob3))
I won't tell you the
Mahdi Osman m_osm at gmx.net writes:
Hi list,
I am fitting microarray data (intensity) model using the lme package in R
environment. I have 5 fixed variables in the model. One of the fixed
variables is genes. I am trying to get p-values for different genes. But I
am getting only one
I estimate some parameters and I want to print them (pretty) on my plot:
# somehow estimated parameters
z-c(1.543523e+00, 1.23453e+00, 3.454000e+00)
x-seq(-1,1,length=100)
plot(x,z[3]*x^2+z[2]*x+z[3],type=l, main=My nice plot of the
estimated function)
I found it immediately after posting :(
substitute is my friend:
text(0.5,5,substitute(f[Sv] ==k*x^2, list(k=zf[1])))
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
On 2/8/2006 8:30 AM, Brian D Ripley wrote:
On Wed, 8 Feb 2006, Duncan Murdoch wrote:
On 2/8/2006 4:53 AM, Bj�rn-Helge Mevik wrote:
Why don't you test it yourself?
E.g.,
set.seed(42)
bob1 - rnorm(1000,0,1)
set.seed(42)
bob2 - rnorm(500,0,1)
bob3 - rnorm(500,0,1)
I am using vioplot to make nice boxplots, which I am outputting to
postscript. I am using the paper=special and width and height to get
the graphs to look nice. I do however have a problem with the size of
the labels and titles. When I set fontsize as an option to the
postscript call too high,
Dirk De Becker wrote:
* Determine the range of the spectrum to be used - For this, I should
be able to calculate the regression coefficients
You can get the regression coefficients from a PLSR/PCR with the
coef() function. See ?coef.mvr However, using the regression
coefficients alone for
I have a dataframe with a column, say x consisting of values, each
value appearing different times, e.g.
x: 1,1,1,1,2,2,4,4,4,9,10,10,10,10,10 ...
and a vector, including e.g.:
y: 2,9,10,...
I need a subset of the dataframe: all rows where x is equal to one of
the values in y. Currently I use a
Only put the expression inside the expression?
plot(x,z[3]*x^2+z[2]*x+z[3],type=l, main=My nice plot)
text(-0.9,5,paste(zs[1],' ',expression(x^3))) # should work
Robert W. Baer, Ph.D.
Associate Professor
Department of Physiology
A. T. Still University of Health
Also see ?bquote
On 2/8/06, Thomas Steiner [EMAIL PROTECTED] wrote:
I found it immediately after posting :(
substitute is my friend:
text(0.5,5,substitute(f[Sv] ==k*x^2, list(k=zf[1])))
__
R-help@stat.math.ethz.ch mailing list
text(-0.9,5,paste(zs[1],' ',expression(x^3))) # should work
this produces a 2x^3 and not 2xsuperscript(3) - just execute it and
you see what I mean. anyway, substitute helps and bquote could shorten
it a bit.
Anywa thanks!
Thomas
__
On 2/8/2006 9:21 AM, Bernhard Baumgartner wrote:
I have a dataframe with a column, say x consisting of values, each
value appearing different times, e.g.
x: 1,1,1,1,2,2,4,4,4,9,10,10,10,10,10 ...
and a vector, including e.g.:
y: 2,9,10,...
I need a subset of the dataframe: all rows where x
Bernhard Baumgartner wrote:
I have a dataframe with a column, say x consisting of values, each
value appearing different times, e.g.
x: 1,1,1,1,2,2,4,4,4,9,10,10,10,10,10 ...
and a vector, including e.g.:
y: 2,9,10,...
I need a subset of the dataframe: all rows where x is equal to one of
Dear Bernhard,
if I understand correctly your question
may be you want something like
df-data.frame(x=sample(1:10, 100, repl=T),
y=sample(1:5, 100, repl=T))
subset(df, x%in%y)
Regards,
Stefano
-Messaggio originale-
Da: [EMAIL PROTECTED]
[mailto:[EMAIL
Sounds like you may need no use match().
On Wed, 2006-02-08 at 15:21 +0100, Bernhard Baumgartner wrote:
I have a dataframe with a column, say x consisting of values, each
value appearing different times, e.g.
x: 1,1,1,1,2,2,4,4,4,9,10,10,10,10,10 ...
and a vector, including e.g.:
y:
Hi
something like
xx-data.frame(x=sample(1:10,100,replace=T))
y-c(2,5,8)
xx[xx$x%in%y,]
HTH
Petr
On 8 Feb 2006 at 15:21, Bernhard Baumgartner wrote:
From: Bernhard Baumgartner [EMAIL PROTECTED]
Organization: Universitaet Regensburg
To:
I installed the MinGW-5.0.0 compilers today (gcc-3.4.2 is now contained in the
'Current' distribution. When the current compilers failed to compile the
referenced Matrix update and not knowing how to check the exact version numbers
for individual files, I applied the Id.exe and f771.exe fixes
Thanks to all,
the %in% function solved my problem!
Bernhard
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Here's one way,
x - data.frame(V=c(1,1,1,1,2,2,4,4,4,9,10,10,10,10,10))
y - data.frame(V=c(2,9,10))
xy - merge(x,y,all=FALSE)
Pay close attention to what happens if you have duplicate values in y, say
y - data.frame(V=c(2,9,10,10))
-Original Message-
From: [EMAIL PROTECTED]
Duncan Murdoch [EMAIL PROTECTED] writes:
This isn't really something that can be proved by a test. Perhaps the
current implementation makes those equal only because 500 is even, or
divisible by 5, or whatever...
I think the intention is that those should be equal, but in a quick
search
On 8 Feb 2006, [EMAIL PROTECTED] wrote:
I would like to only inspect rules that contain a certain label
substring on the rhs. In this special case the item labels are built
like this:
itemtype_itemvalue e.g. Artikelgruppe_E0815 what I want to do
is only show rules where Artikelgruppe is
I am trying to fit a generalised least squares model using gls in the nlme
package.
The model seems to fit very well when I plot the fitted values against the
original
values, and the model parameters have quite narrow confidence intervals
(all are
significant at p5%).
The problem is that the
Please forgive a rather naïve question...
Could someone please give a quick explanation for the differences in conf
intervals achieved via confint.glm (based on profile liklihoods) and the
intervals achieved using the Design library.
For example, the intervals in the following two outputs are
Dear R-Users,
I am searching a way of changing the labels of the scale of a image.plot
function in library(fields). Taking the following example I would like
to add my own labels on the right side of the scale and it would be nice
to decide also where to have the ticks.
In this particular template
On Wed, 8 Feb 2006, Duncan Murdoch wrote:
On 2/8/2006 8:30 AM, Brian D Ripley wrote:
On Wed, 8 Feb 2006, Duncan Murdoch wrote:
On 2/8/2006 4:53 AM, Bjÿÿrn-Helge Mevik wrote:
Why don't you test it yourself?
E.g., set.seed(42) bob1 - rnorm(1000,0,1) set.seed(42) bob2 - rnorm(500,0,1)
Thanks for pointing out the aod package and the beta-binomial logistic
models Renaud.
While I see how betabinom could be applied to some of our other analyses ,
I don't see how it can be used in our habitat selection analysis where
individual locations are coded as 0 or 1 rather than
Hi Sumanta,
1. This messages is much more appropriate for the sig-finance DL
instead. Consider signing up (I read up on Amba, so I am sure you have
good contributions to make in that forum).
2. To my knowledge, there isn't a direct package. However, if you use
Bloomberg's excel plugin, just get
Hi Jarrod ,
I think you are right that neither nested or cross-classified models can be
specified
in glmm.admb() as it stands (except in ad-hoc ways). This is not a restriction
of the underlying
software AD Model Builder, but it is a feature of how I wrote the R-interface.
I will try to figure
Hi,
after looking around I was not able to get info about these two questions:
1. Is there a function to have a jackknifed corrected var/cov estimate (as
described in MacKinnon and White 1985) in a standard OLS regression?
2. Does R possess a LAD (Least Absolute Deviation) regression
Hello all,
Research Systems (www.rsinc.com) have developed and distributes the language
IDL,
and the GIS ENVI, written in IDL.
To my oppinion, R language is superior, compared to IDL, in all aspects.
However, ENVI is the rather convenient and feature rich tool.
Is anyone aware about
R-Helpers,
I am trying to insert a value into a dataframe. This value is a proportion
calculated by counting the number of those individuals with that value and then
inserting the proportion at the end of the dataframe to only those individuals
with the given value. The problem I am
Please read the R-admin manual and as it suggests use the candidate
compilers, now gcc 3.4.5.
You also need the Rtools at the front of your path, so you find find.exe
not FIND.exe.
Finally, it seems to be requiring latex, so you may need to install that
(I am really not sure why though).
Hall Christian
Ich habe Dir hier den Output desselben Models mir den selben Daten mit R
und mit Stata gerechnet angehängt.
Wie Du siehst sind die Unterschiede nicht unbeträchtlich. Da Du bestimmt
schon ähnliche Erfahrungen gemacht hast, wollte ich Dich fragen, ob
diese Unterschiede im Rahmen
On Wed, 8 Feb 2006, Cox, Stephen wrote:
Please forgive a rather naïve question...
Could someone please give a quick explanation for the differences in
conf intervals achieved via confint.glm (based on profile liklihoods)
and the intervals achieved using the Design library.
Well, the Design
Cox, Stephen wrote:
Please forgive a rather naïve question...
Could someone please give a quick explanation for the differences in conf
intervals achieved via confint.glm (based on profile liklihoods) and the
intervals achieved using the Design library.
For example, the intervals in
*
Gary A. Nelson, Ph.D
Massachusetts Division of Marine Fisheries
30 Emerson Avenue
Gloucester, MA 01930
Phone: (978) 282-0308 x114
Fax: (617) 727-3337
Email: [EMAIL PROTECTED]
I have a time series of temperature data recorded every 2 hours and I
would like to generate a spectrogram using the spectrum function that
examines cycles per day. My statement for coding the series using the
ts function is
bevtemp-ts(bevtemp,deltat=0.084)
where 0.084 is 2 hours/24
I have some computers with a massive amount of memory, and I have some
jobs that could use very large matrix sizes. Can R handle matrices of
larger than 2GB? If I were to create a matrix of 1,000,000 x 1,000, it
would use about 8GB. Can R work with an array of that size if I have
compiled R
R-Helpers,
I am trying to insert a value into a dataframe. This value is a proportion
calculated by counting the number of those individuals with that value and then
inserting the proportion at the end of the dataframe to only those individuals
with the given value. The problem I am
Dear R-users,
Sorry for bothering you twice in a day, but I was wondering whether
there is any R-function which can easily plot the single elements of a
Lexis diagram. I though that Lexis.diagram(Epi) could be the case, but
it just plot life times in a frame.
In particular I have been searching
White, Charles E WRAIR-Wash DC [EMAIL PROTECTED] writes:
I installed the MinGW-5.0.0 compilers today (gcc-3.4.2 is now
contained in the 'Current' distribution. When the current compilers
failed to compile the referenced Matrix update and not knowing how
to check the exact version numbers for
Check this out:
https://stat.ethz.ch/pipermail/r-sig-finance/2004q3/51.html
You should consider adapting it with rcom instead of RDCOMClient
Neuro
From: Sumanta Basak [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject: [R] Bloomberg Data Import to R
Date: Wed, 8 Feb 2006 15:51:13
Hi,
I was wondering if anyone could tell me why prcomp() will Invent modes
of variation in a PCA on identical replicates of data? I would have
expected 50 (or whatever number) of identical replicates to return a null
score in such an analysis (or at the least, all variables would share the
same
[EMAIL PROTECTED] wrote:
There are lots of tutorials at
http://cran.r-project.org/other-docs.html.
I have those of course. I was looking for further tutorials.
What do you mean the statistical aspect of R?
I mean using and building statistical models for data analysis
specifying the
RSiteSearch() is your friend.
For 1) Try RSiteSearch('Jackknife',restrict='functions')
For 2) see package quantreg
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P. Box
Dear list,
I've recently came across a problem that I think I've solved and that I wanted
to share with you for two reasons:
- Maybe others come across the same problem.
- Maybe someone has a much simpler solution that wants to share with me ;-)
The problem is as follows: expand.grid() allows
Angelo Secchi wrote:
Hi,
after looking around I was not able to get info about these two questions:
1. Is there a function to have a jackknifed corrected var/cov estimate
(as described in MacKinnon
and White 1985) in a standard OLS regression?
Not sure, but look at package sandwich
Dear members,
this can sound trivial for some people but I don't have experience on
compilation.
I'm trying to install R-2.2.1 on a laptop running Mandriva 2006. I
already installed many of the recommended packages/libraries but it
seems that something is still missing. The problem is that I
Wladimir Eremeev wrote:
Hello all,
Research Systems (www.rsinc.com) have developed and distributes the
language IDL,
and the GIS ENVI, written in IDL.
I find it hard to believe they wrote it all in IDL! I'm guessing its
probably scriptable in IDL, but underneath its written in
HI All,
I have a data frame such as:
test
x y p d
[1,] 1 0 10 21 0
[2,] 2 3 11 12 0
[3,] 3 4 12 23 0
[4,] 3 5 13 24 0
and I want to perfor some operations on the first two coulums,
conditional on the uneqaulity values on the 3rd and 4th columns.
For instance:
j = 3
test[test[,1] ==
Hi,
I am using R 2.2.0 on Mac OS X 10.3.9. But I test the issue also with R
2.2.1 on OS X 10.4 (Tiger)
I have a question regarding the functioning of the history from the
command line.
When I press the 'up-arrow' I call back the last command (everything
ok), let's say that I go back until the
On Feb 8, 2006, at 10:22 AM, Angelo Secchi wrote:
1. Is there a function to have a jackknifed corrected var/cov
estimate (as described in MacKinnon and White 1985) in a standard
OLS regression?
package: sandwich
2. Does R possess a LAD (Least Absolute Deviation) regression
If you don't get a response:
reproducible example, please.
I suspect that you're looking at numerical analysis artifacts, assuming I
have correctly interpreted you.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the
On 2/8/2006 12:50 PM, Laura Quinn wrote:
Hi,
I was wondering if anyone could tell me why prcomp() will Invent modes
of variation in a PCA on identical replicates of data? I would have
expected 50 (or whatever number) of identical replicates to return a null
score in such an analysis (or at
On Wed, 8 Feb 2006, Mike Miller wrote:
I have some computers with a massive amount of memory, and I have some
jobs that could use very large matrix sizes. Can R handle matrices of
larger than 2GB? If I were to create a matrix of 1,000,000 x 1,000, it
would use about 8GB. Can R work with an
Hi
I had similar issue which was resolved by ensuring that I have the lib/include
folder for readline in my LD_LIBRARY_PATH and PATH respectively.
Thanks
Gautam Bhola
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of André Beló
Sent: Wednesday,
Why do I see this error?
library(stats)
require(stats)
[1] TRUE
tHSD - TukeyHSD(aov)
Error in TukeyHSD(aov) : no applicable method for TukeyHSD
In case it helps:
aov
Call:
aov(formula = roi ~ (Cue * Hemisphere) + Error(Subject/(Cue *
Hemisphere)), data = roiDataframe)
Grand Mean:
I have been creating some graphic postscript files using xyplot but now
I find that I need then in either the jpeg or tiff format. I have tried
re-running just the lattice plot so that it is on the screen and then
using the jpeg(filename=myplot. from the bmp, jpec and png
graphics devices
Barry,
I hope you will share your code! We've been using QGIS and I have in
mind just such an interface to R.
We have ENVI in the lab and it is very powerful for image processing. It
would take many many person years to reproduce all of the advanced
algorithms in that package. (It does a lot of
Hi -
I have two response variables 'y1' and 'y2' and a factor 'x'. I would like
to create paired box-whiskers plots for y1~x and y2~x and labeled for the
same x. the b-w plots would be side-by-side in the same panel - almost like
a barchart with two parallel columns for the same x.
the code
On Wed, 8 Feb 2006, Dean Sonneborn wrote:
I have been creating some graphic postscript files using xyplot but now
I find that I need then in either the jpeg or tiff format. I have tried
re-running just the lattice plot so that it is on the screen and then
using the
From: =?iso-8859-1?q?Lu=EDs_Torgo?= [EMAIL PROTECTED]
Date: Wed, 8 Feb 2006 18:08:40 +
Dear list,
I've recently came across a problem that I think I've solved and that I
wanted
to share with you for two reasons:
- Maybe others come across the same problem.
- Maybe someone has a
Hello,
I want to find the parameters mu and sigma that minimize the following
function.
It's important, that mu and sigma are strictly positive.
-
optimiere = function(fmean,smean,d,x,mu,sigma)
{
merk = c()
for (i in 1:length(d))
hello everybody,
i want to plot a count process in the following way, but i don't know
how i can do that. i have data for example x-(0,2,6,2,8,4,.) and
dates y which is a vector of weekly dates for example
(01/01/06, 08/01/06, 15/01/06, 22/01/06, ), now i want to plot the
y's an the
Hi,
I just installed R version 2.2.1, I installed DBI package and worked
fine and then I installed ROracle but I got this error:
[EMAIL PROTECTED] bin]$ R CMD INSTALL
--configure-args='--enable-extralibs' ../../ROracle_0.5-5.tar.gz
* Installing *source* package 'ROracle' ...
creating cache
On 2/8/06, Chris Behr [EMAIL PROTECTED] wrote:
Hi -
I have two response variables 'y1' and 'y2' and a factor 'x'. I would like
to create paired box-whiskers plots for y1~x and y2~x and labeled for the
same x. the b-w plots would be side-by-side in the same panel - almost like
a barchart with
Hi,
I want to maximize a liklihood function with multiple parameters. There
is no closed-form analytical solution to the estimates of the
parameters, and I would like to implement a Newton-Raphson iterrative
approach. Is there a maximization procedure, such as the Newton-Raphson
algorithm,
Hello,
I have a list (mode and class are list) in R that is many elements long and of
the form:
length(list)
[1] 5778
list[1:4]
$ID1
[1] num1
$ID2
[1] num2 num3
$ID3
[1] num4
$ID4
[1] NA
I'd like to convert the $ID2 value to be in one element rather than in two. It
shows up as c(\num2\,
LD == Liz Dem [EMAIL PROTECTED] writes:
LD I have a list (mode and class are list) in R that is many elements long
and of the form:
length(list)
LD [1] 5778
list[1:4]
LD $ID1
LD [1] num1
LD $ID2
LD [1] num2 num3
LD $ID3
LD [1] num4
LD $ID4
LD
1 - 100 of 104 matches
Mail list logo