Let us be clear that TIFF is not a format, it is a collection of formats
and for example my camera's RAW images are TIFFs, as are Adobe's DNG
files. (See e.g. http://en.wikipedia.org/wiki/TIFF.) So the journals
concerned are misusing the term. It is very easy to convert either JPEGs
or PNGs
col.region changes both colors in plot and colorkey.
try:
x - seq(pi/4, 5 * pi, length = 100)
y - seq(pi/4, 5 * pi, length = 100)
levelplot(z~x*y) #default
levelplot(z~x*y,col.regions=rainbow(24)) #custom color
On 1/11/07, Bram Kuijper [EMAIL PROTECTED] wrote:
Hi all,
I try to make a
Thanks in advance for any suggestions:
I am using the 'multhist' function in the 'plotrix' package to
display histograms of some variable 'x' given some value of another
variable 'z' -- for example, separate histograms for some variable
according to sex (or another dichotomous) variable. Thus
Petr Pikal wrote:
Hi
On 11 Jan 2007 at 23:11, Farrel Buchinsky wrote:
Not that I know of.
When I get back to the office I will check it out.
I certainly do not recall having created one.
I wonder if a library I am working with maybe created a mean function.
How do I figure out if I
Dear R-users,
Following is part of my data, where slide has 36
levels and block 48 levels. I have done boxplot for
each slide on the same graph. There are outliers for
each slide and I tried to use indentify functtion to
identify outliers in such a way that when I click on
an outlier or point,
My simpleminded understanding of simple regression is that when
plotting regression lines for x on y and y on x in the same plot, the
lines should cross each other at the respective means. But, given the
R function below, abline (lm(y~x)) works fine, but abline (lm(x~y))
does not. Why?
Try this version of your function and then think about it
tst - function () {
attach (attitude)
x - rating
y - learning
detach (attitude)
plot (x, y)
abline(v=mean(x))
abline(h=mean(y))
abline (lm(y~x))
cc - coef(lm(x ~ y))
abline (-cc[1]/cc[2], 1/cc[2])
}
My simpleminded understanding of
because given data is a part of your data, I cannot examine,
however, try:
##out.block-identify(tb_ncs$y,tb_ncs$Slide)
out.block-identify(tb_ncs$Slide,tb_ncs$y)
On 1/11/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Dear R-users,
Following is part of my data, where slide has 36
levels and
On Fri, 12 Jan 2007, Tom Backer Johnsen wrote:
My simpleminded understanding of simple regression is that when
plotting regression lines for x on y and y on x in the same plot, the
lines should cross each other at the respective means. But, given the
R function below, abline (lm(y~x))
On Fri, 12 Jan 2007, Tom Backer Johnsen wrote:
My simpleminded understanding of simple regression is that when
plotting regression lines for x on y and y on x in the same plot, the
lines should cross each other at the respective means. But, given the
R function below, abline (lm(y~x)) works
Tom Backer Johnsen wrote:
My simpleminded understanding of simple regression is that when
plotting regression lines for x on y and y on x in the same plot, the
lines should cross each other at the respective means. But, given the
R function below, abline (lm(y~x)) works fine, but abline
How do I draw non-square matrices with image() and get the axes right?
Try 1:
a - matrix(rnorm(100),20,5)
image(1:20,1:5,a,asp=1,xlab=label here)
# No good because the axes don't touch the image
Try 2:
image(1:20,1:5,a,asp=1,axes=F,xlab=label here)
axis(side=1,pos=0)
# No good because the
Dear R users,
I am trying to include C code into R via the .C interface. I have read
that arguments passed to a C function have to be correctly DEreferenced.
This is something that can be easily done for numbers (integers or
float) by adding
a * before the reference like for instance:
Prof Brian Ripley wrote:
Where did you tell it 'x' was the abscissa and 'y' the ordinate?
(Nowhere: R is lacking a mind_read() function!)
Please stop complaining about missing features. Patches will be considered.
Oh, it's you, Brian. Never mind then. You'll get to it, I'm sure.
;-)
--
This should do the trick:
mind_reader - function() {
ll - letters[round(runif(6, 1, 26))]
ff - ll[1]
for (ix in 2:length(ll)) {
ff - paste(ff, ll[ix], sep = )
}
if (exists(ff)) {
cat(The function that you were
Robin Hankin wrote:
How do I draw non-square matrices with image() and get the axes right?
Try 2:
image(1:20,1:5,a,asp=1,axes=F,xlab=label here)
axis(side=1,pos=0)
# No good because the x axis label is floating far from the x axis.
Its only no good if your plot device isnt a similar
ken knoblauch wrote:
This should do the trick:
mind_reader - function() {
ll - letters[round(runif(6, 1, 26))]
I see my paraNormal distribution package hasn't found its way to CRAN yet:
http://tolstoy.newcastle.edu.au/R/help/05/04/1701.html
Barry
Barry Rowlingson wrote:
ken knoblauch wrote:
This should do the trick:
mind_reader - function() {
ll - letters[round(runif(6, 1, 26))]
I see my paraNormal distribution package hasn't found its way to CRAN yet:
http://tolstoy.newcastle.edu.au/R/help/05/04/1701.html
LOL! Nice!
On 1/12/2007 5:56 AM, Barry Rowlingson wrote:
ken knoblauch wrote:
This should do the trick:
mind_reader - function() {
ll - letters[round(runif(6, 1, 26))]
I see my paraNormal distribution package hasn't found its way to CRAN yet:
Stephane Cruveiller said the following on 1/12/2007 4:15 AM:
Dear R users,
I am trying to include C code into R via the .C interface. I have read
that arguments passed to a C function have to be correctly DEreferenced.
This is something that can be easily done for numbers (integers or
Hello!
I am looking for a function which computes the maximum likelihood
estimator of the autocorrelation function for a gaussian time series.
Does a such function already exist in R?
The estimator by default in R, acf(), uses the method of moments.
Thanks a lot,
Alain
--
Alain Guillet
Dear talepande,
Thanks for your suggestion,
I have already tried to use it, but the identify
function gave me only the observation number
everytime I clicked on any point.What I want is
instead of obervation numbers it would be block and/or
slide numbers.
Any other idea how I can make it
Hello
I want to repeatedly extract coefficients and standard errors from a GLM and
write them into a file (1row=all coefficients of model A, 2 row=all
coefficients of model B, etc.). I can extract coefficients but not standard
errors. furthermore I fail to write extracted values line by line
You will need to give us a reference, as the acf is not a parameter in a
model in your description and MLEs apply to model parameters.
Just possibly ar.mle is what you are looking for, perhaps plus ARMAacf?
On Fri, 12 Jan 2007, Alain Guillet wrote:
Hello!
I am looking for a function which
Hello List.
please find uploaded to CRAN a new package, Brobdingnag.
This package does two things:
(1) allows computation of very large numbers using a logarithmic
representation.
(2) provides a Hello, World example of S4 methods in use: there are
two classes of object
(brob and glub)
Indermaur Lukas said the following on 1/12/2007 7:55 AM:
Hello
I want to repeatedly extract coefficients and standard errors from a GLM and
write them into a file (1row=all coefficients of model A, 2 row=all
coefficients of model B, etc.). I can extract coefficients but not standard
Fortune?
On 1/12/07, Peter Dalgaard [EMAIL PROTECTED] wrote:
Prof Brian Ripley wrote:
Where did you tell it 'x' was the abscissa and 'y' the ordinate?
(Nowhere: R is lacking a mind_read() function!)
Please stop complaining about missing features. Patches will be considered.
Oh, it's
On Friday 15 December 2006 22:38, BRENDAN KLICK wrote:
For summary(GAM) in the mgcv package smooth the degrees of freedom for
the F value for test of smooth terms are the rank of covariance matrix
of \hat{beta} and the residuals df. I've noticed that in a lot of GAMs
I've fit the rank of the
Have you used Tinn-R and what landmines await the inexperienced?
I could not understand why a script that used to work stopped working.
Look at these two scenarios
I opened an excel spreadsheet and copied several cells to the clipboard
Then Scenario 1
Executed from Tinn-R
Thanks to all the responders. Here are some replies to the comments:
1) Concerning the term TIFF format.
It may be that the journals are misusing the term TIFF, but it would
also appear that wikipedia is as well. The first sentence in the wiki
link sent below states:
Tagged Image File FORMAT
Farrel Buchinsky wrote:
Have you used Tinn-R and what landmines await the inexperienced?
Depending on which button you press in Tinn-R, the clipboard is used to
transfer the commands to R, so sometimes you can't rely on the previous
contents of the clipboard while using Tinn-R. If you use
I have had good results post processing xfig output
with fig2dev:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/60735.html
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/60762.html
On 1/12/07, Inman, Brant A. M.D. [EMAIL PROTECTED] wrote:
Thanks to all the responders. Here are some
Inman, Brant A. M.D. wrote:
Thanks to all the responders. Here are some replies to the comments:
1) Concerning the term TIFF format.
It may be that the journals are misusing the term TIFF, but it would
also appear that wikipedia is as well. The first sentence in the wiki
link sent below
Dear R-users
I'm considering a repeated measures experiment where two
within-subject factors A (2 levels) and B (3 levels) have been
measured for each of 14 subjects, S. I wish to test the effect of
factor A. I know that a variance component model with random effects
S, S:A, S:B and S:A:B can be
Dear friends,
I want to do a unequal probability sampling, that is, Probability
Proportionate to size, Is it right for the following programs?
Say my original dataset is:
ID Population
1 100
2 200
3 300
IF the population is large ,then the corresponding ID has the large
Prof. Brian Ripley,
You are right, my question was not clear.
In fact, I want to estimate the k first components of the acf, i.e. I
want to estimate the k parameters (c(0),c(1),...c(k-1)), where c is the
autocorrelation function, by a maximum likelihood estimator.
Alain
Prof Brian Ripley a
Hello R-Users!
Does anyone know of a package to draw/analyze silicon wafer maps?
Here are some examples
http://www.java2s.com/Code/Java/Chart/JFreeChartWaferMapChartDemo.htm
http://dp.pdf.com/site/products/wafermap/binmap.html
I've google'd and searched CRAN with no luck.
It
On Fri, 12 Jan 2007, Alain Guillet wrote:
Prof. Brian Ripley,
You are right, my question was not clear.
In fact, I want to estimate the k first components of the acf, i.e. I
want to estimate the k parameters (c(0),c(1),...c(k-1)), where c is the
autocorrelation function, by a maximum
Thank you Peter Dalgaard.
When I open a DOS box and type gswin32c, I do indeed get an error message
saying that it can't find the program. I edited the Windows system
environmental variable Path and the user environmental variable PATH
(wasn't sure which to edit), to contain the follwing
Dear Kim,
as far as I understandyour problem correct the specification of the model in
lme is:
lme( fixed=y ~ A*B, random=~1|S)
Thilo
On Friday 12 January 2007 15:54, Kim Mouridsen wrote:
Dear R-users
I'm considering a repeated measures experiment where two
within-subject factors A (2
Inman, Brant A. M.D. wrote:
Thank you Peter Dalgaard.
When I open a DOS box and type gswin32c, I do indeed get an error message
saying that it can't find the program. I edited the Windows system
environmental variable Path and the user environmental variable PATH
(wasn't sure which to
How about dev.off() after your plot! It is easy as you predicted :-)
Janusz.
Inman, Brant A. M.D. wrote:
Thank you Peter Dalgaard.
When I open a DOS box and type gswin32c, I do indeed get an error message
saying that it can't find the program. I edited the Windows system
environmental
Peter Dalgaard wrote:
Inman, Brant A. M.D. wrote:
Thank you Peter Dalgaard.
When I open a DOS box and type gswin32c, I do indeed get an error message
saying that it can't find the program. I edited the Windows system
environmental variable Path and the user environmental variable PATH
On Fri, 12 Jan 2007, Inman, Brant A. M.D. wrote:
Thanks to all the responders. Here are some replies to the comments:
[...]
3)Lack of gratitude by R users.
It is interesting to note upon reviewing the R-help files that many
queries (perhaps even the majority?) do, in fact, convey
At 8:07 AM + 1/12/07, Prof Brian Ripley wrote:
Let us be clear that TIFF is not a format, it is a collection of formats
and for example my camera's RAW images are TIFFs, as are Adobe's DNG
files. (See e.g. http://en.wikipedia.org/wiki/TIFF.) So the journals
concerned are misusing the term.
Dear Thilo
Thanks for your suggestion. I guess the model you are fitting here has
only a single random effect term, namely subject. If the effect of A
depends on S, one needs to include an additional random effects term
for the S:A interaction.
With lme I can get output for the effect of A which
Respected Sir
I am very new user of R. I am trying to use RFA package for Regional
frequency analysis of rainfall data. Sir I am not understanding the example
given in refference manual
data(region)
reganalysis(region)
Sir in this example how can I make my own data(region) and what are the
rows
The only button I pressed was the send line button in version 1.19.1.5. I
changed my command to
prelim-read.delim(source(clipboard))
and still got the same problem.
I do not know how to use the source versions of the buttons. Can you please
tell me more?
Thanks
Farrel
Martin Becker [EMAIL
Dear all,
I have a problem with plotting.
I have an irregular set of data, where latitude is listed from bigger to
smaller values. So, I cannot plot with image.plot.
On the other hand, quilt.plot solves the problem, but it can't take the
definition of color. It gives the feedback:
Error in
Several R Helpers pointed out that I forgot to include the dev.off()
statement in my code. This solved my problem with one caviat: the
output file address cannot have any spaces in it (as pointed out by
Chuck Cleland). For instance:
# This file location works great
Zitat von Farrel Buchinsky [EMAIL PROTECTED]:
The only button I pressed was the send line button in version 1.19.1.5. I
changed my command to
prelim-read.delim(source(clipboard))
The problem is not using the wrong R-command, but using a (Tinn-R)
button that leads to a replacement of the
Prof Brian Ripley sent the following at 12/01/2007 08:07:
Let us be clear that TIFF is not a format, it is a collection of formats
... snip ...
[I don't recall anyone ever writing to thank us for the
PNG or JPEG or bitmap drivers, and lack of appreciation does play a part.]
... rest snipped
On Thu, 11 Jan 2007, David Forrest wrote:
Thanks. The xaxs|yaxs='i' works well for the base graphics. Is there an
additional parameter in play for lattice graphics? The closest I could
gotten is the below which still leaves a bit of a margin:
xy-data.frame(x=c(0,1,1,0,0),y=c(0,1,0,0,1))
On 1/12/07, Martin Becker [EMAIL PROTECTED] wrote:
Depending on which button you press in Tinn-R, the clipboard is used to
transfer the commands to R, so sometimes you can't rely on the previous
contents of the clipboard while using Tinn-R. If you use the
(source)-versions of the buttons, the
Instead of discussing this odd behaviour of TINN-R, I would prefer a discussion
on importing data through the clipboard. In my opinion it isn't a good a idea
to import data with the clipboard. I know that it's a quick and dirty way to
get your data fast into R.
But I see two major drawbacks.
Thierry:
Instead of discussing this odd behaviour of TINN-R, I would prefer a
discussion on importing data through the clipboard. In my opinion it isn't a
good a idea to import data with the clipboard. I know that it's a quick and
dirty way to get your data fast into R.
But I see two major
Hello everybody,
I'm handling a matrix dataset composed by a number of variables much higher
than the objects (900 vs 100) and performing a prcomp (centered and scaled) PCA
on it. What I get is a Loadings (rotation) matrix limited by my lower number of
objects and thus 900x100 instead of
Dear All
1)
I'm fitting spatial CAR models
using R2Winbugs and although everything seems to go reasonably well (or I
think so)
the next message appears from WINBUGS 1.4 window:
gen.inits()
Command #Bugs: gen.inits cannot be executed (is greyed out)
The question is if this message means that
On Thu, 11 Jan 2007, David Forrest wrote:
Thanks. The xaxs|yaxs='i' works well for the base graphics. Is there an
additional parameter in play for lattice graphics? The closest I could
gotten is the below which still leaves a bit of a margin:
xy-data.frame(x=c(0,1,1,0,0),y=c(0,1,0,0,1))
I would say rather that for binary data (binomial data with n=1) it
is not possible to detect overdispersion from examination of the
Pearson chi-square or the deviance. Overdispersion may be, and
often is, nevertheless present. I am arguing that overdispersion is
properly regarded as a
ONLINE VENDOR NEUTRAL INTRO TO DATA MINING FOR ABSOLUTE BEGINNERS
(no charge)
A non-technical data mining introduction for absolute beginners
January 26, 2007, 10AM - 11AM PST
Future Sessions (May 23, June 14, Sept 7)
To register: contact [EMAIL PROTECTED]
This one-hour webinar is a perfect
Zitat von BBands [EMAIL PROTECTED]:
On 1/12/07, Martin Becker [EMAIL PROTECTED] wrote:
Depending on which button you press in Tinn-R, the clipboard is used to
transfer the commands to R, so sometimes you can't rely on the previous
contents of the clipboard while using Tinn-R. If you use the
John Maindonald wrote:
I would say rather that for binary data (binomial data with n=1) it
is not possible to detect overdispersion from examination of the
Pearson chi-square or the deviance. Overdispersion may be, and
often is, nevertheless present. I am arguing that overdispersion
I just do not understand this. You are going to be amazed by what I tell
you.
Firstly thank you for orienting me to the Tinn-R source buttons.
I selected the first two lines in my script that were as follows.
prelim-read.delim(clipboard)
str(prelim)
Then I went over to Excel and copied the
There is some information in Programming with Data by John Chambers. He
has reference to a site (
http://cm.bell-labs.com/cm/ms/departments/sia/project/icmanuf/index.html)
that you might want to check out.
On 1/12/07, Walker, Sam [EMAIL PROTECTED] wrote:
Hello R-Users!
Does anyone know of
In fact, I need it in the general case, not only for an ARMA process.
Unfortunately, I have no reference to give so I will code it. Sorry for
the trouble.
Alain
Prof Brian Ripley a écrit :
On Fri, 12 Jan 2007, Alain Guillet wrote:
Prof. Brian Ripley,
You are right, my question was not
Dunja Drvar wrote:
Dear all,
I have a problem with plotting.
I have an irregular set of data, where latitude is listed from bigger to
smaller values. So, I cannot plot with image.plot.
On the other hand, quilt.plot solves the problem, but it can't take the
definition of color. It gives the
This is a general question to people who've installed R on a UNIX
sparc-sun-solaris platform:
Have you had any issues related to maintaining R on this platform, e.g.,
installations that didn't work, instances of R crashing and possibly
requiring a new installation, etc?
I'm especially interested
Hi all,
Sorry about the simple question, but I have searched the web with
prompt , input etc. and never got the answer .
thanks a lot
tong
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
you can get them from return value, try:
out.id-identify(tb_ncs$Slide,tb_ncs$y)
out.block-tb_ncs[out.id,]$Block
On 1/12/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Dear talepande,
Thanks for your suggestion,
I have already tried to use it, but the identify
function gave me only the
try:
readLines(n=1)-str
On 1/13/07, Tong Wang [EMAIL PROTECTED] wrote:
Hi all,
Sorry about the simple question, but I have searched the web with
prompt , input etc. and never got the answer .
thanks a lot
tong
__
71 matches
Mail list logo