For an answer to your question 2. look at help(try)
and help(tryCatch).
--- Yong Wang [EMAIL PROTECTED] wrote:
Dear Rlist
I am runing a for loop on a large dataset to do
exploring
investigation. Code embedded in the loop include the
lm routine.
Unfortunately, for some specification of
Hello,
I have some doubts on TukeyHSD application.
I want to investigate the effects of depth, latitude and month variation on
the length of a fish. These are orthogonal and observational data.
For this, I have made an aov model (L~month+lat+prof+month*lat), after
applying drop1 and step
This is a repost of
https://stat.ethz.ch/pipermail/r-help/2007-September/141727.html
Please do study the posting guide to see why you did not get an answer and
what to do when you do not.
It is nothing to do with TukeyHSD, as the differences are there in the
means.
It is clear that your
Dear forum,
I'm trying to install RMySQL but I'm having problems in loading it, here is
the message that I'm getting in R
install.packages(RMySQL, dependencies=TRUE)
trying URL
'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.5/RMySQL_0.6-0.zip'
Content type 'application/zip' length
Hello,
I create a matrix:
best - matrix(0, ncol=2, nrow=num.selected,
dimnames=list(the.best$.Name, c(Probability(%), Inside)))
best[,1] - as.numeric(the.best$Total*100)
best[,2] - ifelse(the.best$Weight==0, No, Yes)
What I want is the second column of mode numeric, but it is of mode
Dear List,
is there a way to send the contents of a variable (e.g. a list) directly to a
laser printer. I am
aware that I could save the variable to a ascii-file and then print the file,
but I would prefer to
print from R directly.
I have been looking for a way to do so for the better part
Hi Jeff,
That looks like a nice initiative. However, if you are interested in
getting contributions from the community, it might be good to spell
out how others might use the content of the site. Currently you have
copyright r-cookbook.com, but maybe you could consider a creative
commons
S Bina wrote:
Hi,
I cannot find the function kurtosis. Is it sth additional I am meant to
download? I use the MacOS X version of R.
Many thanks
Samira
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On 2007-September-28 , at 16:57 , Frank E Harrell Jr wrote:
jiho wrote:
On 2007-September-28 , at 15:18 , Paul Smith wrote:
On 9/28/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
I know how to export graphics as pdf files and then how to include
them in LaTeX documents. However, I do not
Hi Philippe,
I don't want to be competing in any way with any of the fantastic
official R resources that exist; I only want to supplement them.
Although this list is probably not the best place to discuss the
proper use of wiki's for documentation/learning, I'll make a few
comments on
Hi,
I don't know if it will help you but to retrieve the slope of lm(y~x), I use
:
coefficients(lm(y~x))[2]
([1] for the intercept)
I should not tell this but It took me a long time to find this obvious
thing.
Have a nice week-end,
Ptit Bleu.
On 2007-September-28 , at 15:18 , Paul Smith wrote:
On 9/28/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
I know how to export graphics as pdf files and then how to include
them in LaTeX documents. However, I do not know how to do in
order to
have the text of the graphics written with
On 9/28/07, Manuel Ramon [EMAIL PROTECTED] wrote:
Hi to everyone,
I am starting to work with a model that is not familiar to me.
The model would be like that:
y = . - a*max(Q(i) - Q(0), 0) + .
where Q(i) is the accumulated effect of a variable at time i and Q(0) a
threshold
R Users,
Emine Bayman sent this out earlier and we do not think it went through.
Appologies if it did.
We want to fit GLMM with lmer with binomial data and a one-way random
effects model (overall mean is a fixed effect and there are random
effects for each binomial).
We are using the
Dear All,
Can R draw plots of functions on a Cartesian coordinate system with
axes like the ones shown at
http://en.wikipedia.org/wiki/Image:Cartesian-coordinate-system-with-circle.svg
?
I have already searched the R web-site, but found nothing.
Thanks in advance,
Paul
I am trying to use ROracle (v 0.5-8) with R on Windows (v 2.5.1)
I get an error message when I run the following code
library(DBI)
library(ROracle)
m - dbDriver(Oracle)
con - dbConnect(m, dsn='', uid='', pwd='')
Note masked real values with ###
The
Hello,
I'm starting with data analysis with R, and I'd like to know if there is in
any package a function to have price elasticity of an item.
For example I have information about sale quantity and price, and I want to
know if the series is sensitive to the price.
Sale Quantity Price
34
Adam Wilson-4 wrote:
I am running R 2.5.1, RMySQL 0.6 , and DBI 0.2-3 on Windows XP
Like others, I am having trouble with NA/Null value conversions between R
and a MySQL database via DBI, but I could not find my exact problem in the
archives. Most of the time NA values in R get
On 9/28/07, Frank E Harrell Jr [EMAIL PROTECTED] wrote:
I know how to export graphics as pdf files and then how to include
them in LaTeX documents. However, I do not know how to do in
order to
have the text of the graphics written with the font selected for the
LaTeX document. Is that
According to the docs package e1071 has kurtosis.
__
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,
On Fri, 28 Sep 2007, Paul Smith wrote:
Dear All,
I know how to export graphics as pdf files and then how to include
them in LaTeX documents. However, I do not know how to do in order to
have the text of the graphics written with the font selected for the
LaTeX document. Is that possible?
On 9/27/07, marcg [EMAIL PROTECTED] wrote:
Thanks a lot
This already looks nice and I already checked the ?wireframe, but with no
examples and as newcommer its hard to find out a correct code.
If we set drape=F in the example:
g - expand.grid(x = 1:10, y = 5:15, gr = 1:2)
g$z -
On Fri, 28 Sep 2007, Sergey Goriatchev wrote:
Hello,
For my functions I want to create output similar in appearance to that
of what you get when you print a summary of lm model:
Residuals:
Min1Q Median3Q Max
-0.209209 -0.043133 0.001793 0.044105
Not sure how you want to handle the NAs, but you could try the
following:
#start
MalVar29_37 - read.table(textConnection(V1 V2 V3 V4 V5 V6 V7 V8 V9
0 0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0 0
0 0 0 0 0 1 0 0 0
NA NA NA NA NA NA NA NA NA
0 1 0 0 0 1 0 0 0), header=TRUE)
Hi Everyone,
I would like to add errorbars to a lattice plot. I already have a function that
adds error bars to a plot (originally written by George Gilchrist) which takes
as arguments the coordinates of the point and then the size of the error bar. I
would like to integrate it in my lattice
Your problem is that a matrix can't be a mix of data types, but a data
frame can. Try the following code.
newBest = as.data.frame(best)
Now you can convert the column of numbers in newBest to numbers from
the factors they start as. Check the help file for data.frame(...)
for more details of
On 9/28/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
Can R plot graphs like the one at
http://www.mathwords.com/f/f_assets/floor_graph.gif
with the balls at the discontinuity points?
You can use segments() to draw the segments and symbols() to draw the
balls. For example,
Hello Jeff,
Good initiative,... but why not to put this in the official R Wiki
(http://wiki.r-project.org)? There is a section named 'tips' dedicated
to such little recipes
(http://wiki.r-project.org/rwiki/doku.php?id=tips:tips). It should be
better to centralize all these little tips, don't
Hi,
I have a data set like this:
MutantRepTime OD
02H02100.029
02H02200.029
02H02300.023
02H02180.655
02H02280.615
02H02380.557
02H021121.776
02H02212 1.859
02H023121.668
02H02
I have a large number of textfiles, and one matix in every texfile. I now
want to find an easy way to select those textfiles that contain at least one
value over a certain limit, and collect those matrices in a vector or
something alike. How do I do that?
--
View this message in context:
On 9/27/07, Moisan Yves [EMAIL PROTECTED] wrote:
My experience is that cairoDevice is a lot slower than Cairo,
especially
on Windows, and about equally flaky. If you see how many layers are
involved with Cairo on Windows you will not be surprised.
Actually, I ended up trying to load
Yes, R can do that. Well, actually YOU can do that using R.
But it is hard to believe that you looked very hard before writing. Did you
look at these R functions?
?plot
?line
?points
?arrows
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
Sorry!
paste(x*100,%)
Marc Schwartz [EMAIL PROTECTED] 27/09/2007 15:56:05
On Thu, 2007-09-27 at 14:36 +0100, stat stat wrote:
I am wondering if there is any procedure to write a particular value
in Percentage format, still maintaining it's numeric character. for
example I want to write
Thanks a lot for both solutions of my problem.
I tried it immediately and I understood how they are working.
The next problem for me is now to deal with the NAs. I thought
perhaps it is possible to exclude the variable from the row
comparison if in one of the rows is an NA?
Furthermore it
Samuel Okoye wrote:
Hello,
if we suppose that
times - c(2006-05-14, 2006-06-12, 2006-06-12, 2006-05-14,
2006-05-14, 2006-06-12)
value - c(2,3,1,4,3,1)
then with
plot(times, value)
we have two boxplots in one graph for 2006-05-14 and 2006-06-12
vittorio wrote:
Compilation of MCMCpack under freebsd 6.2 i386 fails because of the following
cryptic error:
* Installing *source* package 'MCMCpack' ...
checking for C++ compiler default output file name... a.out
checking whether the C++ compiler works... yes
checking whether we are
jiho wrote:
On 2007-September-28 , at 16:57 , Frank E Harrell Jr wrote:
jiho wrote:
On 2007-September-28 , at 15:18 , Paul Smith wrote:
On 9/28/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
I know how to export graphics as pdf files and then how to include
them in LaTeX documents.
Sorry sending this again - is anyone familiar with multiple transparent
wireframe plots?
I already checked the ?wireframe, but with no examples and as newcommer its
hard to find out a correct code.
If we set drape=F in the example:
g - expand.grid(x = 1:10, y = 5:15, gr = 1:2)
g$z -
Paul Smith schrieb:
Dear All,
I know how to export graphics as pdf files and then how to include
them in LaTeX documents. However, I do not know how to do in order to
have the text of the graphics written with the font selected for the
LaTeX document. Is that possible?
Thanks in advance,
Is there any way to password-protect script files (either within R or
otherwise)?
The question seemed to me to be more about password protection against
modification, rather than encryption.
I'd have thought this was something a decent OS could take care of. It seems as
daft to try to get a
Dear all,
I would like to know if there is any R package that has the Barnard´s exact
test. I also like to have any opinions about which test is more powerful:
Fisher´s exact test or Barnard´s exact test.
Best regards
João Fadista
Ph.d. student
UNIVERSITY OF AARHUS
Hello Vittorio,
I am running FreeBSD 7.0-CURRENT (i386) from yesterday and I am able to
compile MCMCpack_0.9-1.tar.gz without any problem on R-2.6.0 RC
(2007-09-27 r43002).
Native FreeBSD devel/pcre-7.3 and lang/gcc42 packages are installed.
Please look if your library versions are ok.
On 9/28/07, Federico Calboli [EMAIL PROTECTED] wrote:
Hi Everyone,
I would like to add errorbars to a lattice plot. I already have a function
that
adds error bars to a plot (originally written by George Gilchrist) which
takes
as arguments the coordinates of the point and then the size of
You do need the RMySQL client DLLs in your PATH. Also make sure that your
MySQL version is not one listed as problematic at
http://www.stats.ox.ac.uk/pub/RWin/ReadMe
On Fri, 28 Sep 2007, andreiabb wrote:
Dear forum,
I'm trying to install RMySQL but I'm having problems in loading it, here
On 9/28/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
Can R plot graphs like the one at
http://www.mathwords.com/f/f_assets/floor_graph.gif
with the balls at the discontinuity points?
You can use segments() to draw the segments and symbols() to draw the
balls. For example,
On 9/28/2007 11:20 AM, Abu Naser wrote:
Hi All user,
I have been using R-2.5.1. dose orientlib support this version? I would like
to try. it.
I have been wondering how to install the library.
As far as I know it is fine. You install it as you would install any
other package, but that
I think he wants the axes crossing at 0,0 not on the outer edges like
the default.
You can put the axes in the plot (though it tends to distract rather
than help in many cases) by:
axis(1, pos=0)
axis(2, pos=0)
You will need to draw the arrowheads yourself. There are options (under
?par) for
The source code for print.summary.lm and summary.lm is in the file
/src/library/stats/R/lm.R of the R source code, which you can download
from CRAN if you don't have it already. The file lm.R is also at
https://svn.r-project.org/R/trunk/src/library/stats/R/lm.R
On Fri, 28 Sep 2007, Sergey
Hi,
I would like to extract the fitted values from a model using LMER but
only for the fix portion of the model and not for the fix and random
portion (e.g it is the procedure outpm or outp in SAS). I am aware of
the procedure fitted() but I not sure it give the fitted values both for
the
On 9/28/07, Greg Snow [EMAIL PROTECTED] wrote:
I think he wants the axes crossing at 0,0 not on the outer edges like
the default.
You can put the axes in the plot (though it tends to distract rather
than help in many cases) by:
axis(1, pos=0)
axis(2, pos=0)
You will need to draw the
Sorry for my silly questions. I'm a beginner with R and most statistics
concepts.
I carried out a simple linear regression where the dependent variable is
explained through a combination of powers
of cos(independent variable).
I can see R returned a good R^2 factor ( 0.99) but I have a hard time
[1]voyages et vacances
Destination vacances
voyages et vacances
[2]Merveilles de l'Egypte
voyages et vacances
[3]La Sicile mystérieuse
voyages et vacances
[4]Athènes et ses joyaux
vacances et séjour chateau de
Hi,
I have been trying to run the code below. In the event of non-convergence,
the statistic in the boot function returns NA as a value. To obtain a BCA
confidence interval, I use boot.ci but, if NA recorded as a value for the
statistic for one the replicates the following error appears:
Look at the las argument in ?par for the easiest solution (can be passed
to axis).
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of
On 9/28/07, hadley wickham [EMAIL PROTECTED] wrote:
Yes there is harm. But to make bold lines, easy to read titles is fine.
See the spar function in
http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see
the setps, ps.slide, and setpdf functions in the Hmisc package.
On 9/28/07, Greg Snow [EMAIL PROTECTED] wrote:
Look at the las argument in ?par for the easiest solution (can be passed
to axis).
Thanks a lot, Greg! That is it!
Paul
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Paul Smith
Sent: Friday,
Does this give you what you want?
fit - lm( Petal.Width ~ Petal.Length, data=iris)
tmp1 - resid(fit)
tmp2 - pnorm( tmp1, 0, summary(fit)$sigma )
par(mfrow=c(2,1))
qqnorm(tmp1)
qqline(tmp1)
plot( ppoints(length(tmp1)), sort(tmp2), xlab='Theoretical Percentiles',
ylab='Sample Percentiles')
I find them helpful for displaying the differences in coupled biomarkers
between therapies with different therapeutic targets. You can quickly
see/show which marker is affected quicker and then the other compensating in
response and how they can be reversed depending on target. Having each
On 9/28/07, Anouk Simard [EMAIL PROTECTED] wrote:
Hi,
I would like to extract the fitted values from a model using LMER but
only for the fix portion of the model and not for the fix and random
portion (e.g it is the procedure outpm or outp in SAS). I am aware of
the procedure fitted() but I
Hi,
I am trying to create a plot of a simple PL graph for an option I'm
pricing using Rmetrics. I'm not trying to do anything fancy. I just
want to figure out if it's possible to put the x-axis at 0 rather
than at the bottom the the entire plot. It seems like something that
would be
Use xaxt='n' in the original plot to suppress the default axis (or
axes=FALSE to supress both), then axis(1, pos=0) to draw the axis at 0.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original
On Sep 28, 2007, at 3:22 PM, Greg Snow wrote:
Use xaxt='n' in the original plot to suppress the default axis (or
axes=FALSE to supress both), then axis(1, pos=0) to draw the axis
at 0.
That was what I was looking for. Thank you. I just need to do a bit
of refining now.
James
Is this easier?
x.index - duplicated(x.sample)==FALSE
cbind(x.sample[x.index],y[x.index])
- Matt
On 9/28/07, Marc Schwartz [EMAIL PROTECTED] wrote:
On Fri, 2007-09-28 at 17:48 -0400, Brian Perron wrote:
Hello all,
An elementary question that I am sure can be easily cracked by an R
hadley wickham wrote:
Yes there is harm. But to make bold lines, easy to read titles is fine.
See the spar function in
http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see
the setps, ps.slide, and setpdf functions in the Hmisc package.
I was interested to see that you
Here is yet another approach using aggregate(), which internally,
basically does what my first solution did:
aggregate(z[, 2], list(z[, 1]), [, 1)
Group.1 x
1 1 -1.2006469
2 2 -0.1614918
3 3 -0.5717729
4 4 -0.2398887
5 5 1.1690564
See ?aggregate
Note
On Fri, Sep 28, 2007 at 12:13:46PM -0500, hadley wickham wrote:
Yes there is harm. But to make bold lines, easy to read titles is fine.
See the spar function in
http://biostat.mc.vanderbilt.edu/SgraphicsHints for a starter. Also see
the setps, ps.slide, and setpdf functions in the
Dear List
In a package I want to import the mApply function from the Hmisc package, and I
would like to import only that function.
1) If I write Depends: Hmisc in the DESCRIPTION file I get the whole Hmisc
package, so that is not the way to go ahead.
2) According to Writing R extensions,
Hi Soren,
What I do in cases like this is just copy the function and place it
into my script at the top (or write it into its own source file and
call it from script). Best,
Matt
On 9/28/07, Søren Højsgaard [EMAIL PROTECTED] wrote:
Dear List
In a package I want to import the mApply function
On 9/28/07, Søren Højsgaard [EMAIL PROTECTED] wrote:
Dear List
In a package I want to import the mApply function from the Hmisc package, and
I would like to import only that function.
1) If I write Depends: Hmisc in the DESCRIPTION file I get the whole Hmisc
package, so that is not the
Hi WizaRds,
I'm experiencing a problem connecting to an Oracle 10g database via RODBC
(I'm getting this on Microsoft XP).
The same SQL queries via PL/SQL Developer work just fine, but when I pump the
query through sqlQuery in RODBC then I get a data frame back with 0 rows.
I cut the query down
Hadley,
This was something I forgot to change before making the site public.
Thanks for pointing out the creative commons license--it's definitely
what I was thinking: share-alike content.
I'd be glad to hear any other suggestions you might have,
Jeff.
On Sep 28, 2007, at 9:15 AM, hadley
On Fri, 2007-09-28 at 21:58 -0500, Derek Ogle wrote:
Hello all,
I am attempting to use the seq() function to define breaks for the
cut() function but am finding performance that I do not understand.
The root of my problem appears to be contained in the following simple
example.
brks -
Hi To r-help
Emergency report. Check DMXC!
Price up 21% in 30 minutes!
5 day price: ~$0.50
__
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
73 matches
Mail list logo