IF you provide us with the reproducible code that is requested in the
posting guide, you may get some valuable help. I gues sthe is from ggplot2?
Uwe Ligges
Felipe Carrillo wrote:
Hi: See the attached pdf graphic. The legend overlaps a little bit with the
strip. Is there a way to move it
Example:
boxplot(d1$b ~ d1$a, pars = list(boxwex = 0.15), at=c(1:3 - 0.1), xaxt=n)
boxplot(d2$b ~ d2$a, add=T, pars = list(boxwex = 0.15), at=c(1:3 + 0.1),
xaxt=n)
axis(1, at=1:3)
Uwe
Malcolm Ryan wrote:
Is there any way to get a boxplot of several data sets beside one
another on the same
Thank you again for all the R help folks who responded. I again appreciate all
the help and insight and will investigate the options suggested.
I guess I still doing a little head scratching at how the division occurred:
It looks like the default hist(...) behavior is doing the following:
The zoo package has lattice methods and does not use require
and it puts lattice in Imports in the DESCRIPTION file. If a user
wants to the use the zoo lattice functions then they must issue
a library(lattice) call. As a result users do not have to load
lattice but if they want to use it they
It looks to me like you need to add x to the arguments in your
call to constrOptim, something like the following:
constrOptim(beta_i, myFunction, NULL, ui, ci, mu = 1e-04, control = list(),
method = Nelder-Mead,outer.iterations = 100, outer.eps = 1e-05, x=x)
The first x in x=x
Or even:
library(chron)
ch - chron(HouseDates)
years(ch)
[1] 1990 1991 1992 1993 1994 1995 1996
Levels: 1990 1991 1992 1993 1994 1995 1996
# or
with(month.day.year(ch), year)
[1] 1990 1991 1992 1993 1994 1995 1996
On Thu, Jun 4, 2009 at 4:02 AM, ONKELINX, Thierry
edw...@web.de wrote:
Hi,
Unfortunately, they can't handle NA. Any suggestion? Some row for Ip
don't have ip address. This cause an error/ wrong result.
A quick fix could be to substitute ... or 0.0.0.0 for the NA
entries. (Use something like
ipch - as.character(df$ip)
ipch[is.na(df$ip)] -
Does anyone know which package include the computation of Cochrans Q
statistic in R?
jlfmssm
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Dear Malcom,
Another option is to merge both dataset to one big dataset and then plot
the big dataset.
d1 - data.frame(a = c(rep(1:3, each = 3)), b = c(1:9), d = A)
d2 - data.frame(a = c(rep(1:3, each = 3)), b = c(9:1), d = B)
Dataset - rbind(d1, d2)
library(ggplot2)
Dataset$a -
Suggest you be more specific. The solution I posted does handle missing values
sorting them to the beginning:
Lines - id rank color status ip
+ 138 29746 yellow no 162.131.58.26
+ 138 29746 red yes 162.131.58.16
+ 138 29746 blue yes 162.131.58.10
+ 138 29746 red no 162.131.58.17
+ 138 29746
Hello,
Can anyone help me with the following:
if one enters a function name in the R console then usually one sees the code
of that function. But there are functions that one cannot see. For example I
want to see the code of print.htest or t.test.default. These functions are
non-visible. Is it
-- begin included message --
I am trying to set up a data set for a survival analysis with time-varying
covariates. The data is already in a long format, but does not have a variable
to signify the stopping point for the interval. The variable DaysEnrolled is
the
variable I would like to use
which package have Hedges g statistic calculation in R?
[[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
jlfmssm wrote:
Does anyone know which package include the computation of Cochran’s Q
statistic in R?
install.packages(outliers)
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
On Jun 4, 2009, at 7:48 AM, jlfmssm wrote:
Does anyone know which package include the computation of Cochrans Q
statistic in R?
jlfmssm
See this thread:
https://stat.ethz.ch/pipermail/r-help/2006-September/113139.html
HTH,
Marc Schwartz
__
Yes. Uwe Ligges' R Help Desk article in this Rnews should answer your
questions:
www.r-project.org/doc/Rnews/Rnews_2006-4.pdf
On Jun 4, 2009, at 9:08 AM, Karin Martijn wrote:
Hello,
Can anyone help me with the following:
if one enters a function name in the R console then usually one
Dear,
I am using cwt function from Rwave package to perform Morlet wavelet analysis.
d1-c(1.31673209591515, -0.171333455797346, -1.67618546079420,
-0.931604651010362, -0.954614183224057, -1.19821794547758, 0.516096363353144,
-0.0432811977638559, 0.737764943619919, 0.438046629673177,
jlfmssm schrieb:
Does anyone know which package include the computation of Cochran?s Q
statistic in R?
jlfmssm
[[alternative HTML version deleted]]
I might be wrong but the packages rmeta and meta report a Q-statistic.
For an example, see
library(meta)
meta1
metagen(meta1$TE,
On Thu, Jun 4, 2009 at 12:36 PM, Patrick Burns pbu...@pburns.seanet.com wrote:
How about:
order(x, runif(length(x)))
Thanks - that is really elegant.
Rainer
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of The R Inferno and A Guide for the
Hello r-help,
I'm fitting a model with lm() and using the orthogonal polynomials
from poly() as my basis:
dat - read.csv(ConsolidatedData.csv, header=TRUE)
attach(dat)
nrows - 1925
Rad - poly(Radius, 2)
ntheta - 14
Theta - poly(T.Angle..deg., ntheta)
nbeta - 4
Beta - poly(B.Beta..deg., nbeta)
I want to make a log-log plot with a regression line, but I can't figure
out what I'm doing wrong. What I'm trying is:
plot(mass,area, log=xy, pch=as.numeric(food))
abline(lm(mass~area))
or
plot(mass,area, log=xy, pch=as.numeric(food))
islands$logmass - log(mass)
islands$logarea - log(area)
Hello,
I need to know how to import ARIMA coefficients. I already determined the
coefficients of the model with other software, but now i need to do the
forecast in R.
For Example: I have a time series named x
and i have fitted an ARIMA(1,0,1) (with other software)
AR coef =
Hi UseRs,
I recently started working with the RGtk2 library. The documentation
is comprehensive and I've tried learning from the examples in most Gtk
tutorials (which have C code). This is a little problematic, but has
helped quite a bit in getting me started.
I would like to create a GUI for
I use gWidgetsRGtk2, which provides widgets such as notebook etc. You
can start with gWidgetsRGtk2, and if you want to tailor the GUI
according to your needs, you can modified them with RGtk2 package.
Ronggui
2009/6/4 Harsh singhal...@gmail.com:
Hi UseRs,
I recently started working with the
alamoboy wrote:
Newbie here. Many apologies in advance for using the incorrect lingo.
I'm new to statistics and VERY new to R.
I'm attempting to group or bin data together in order to analyze them
as a combined group rather than as discrete set. I'll provide a simple
example of the
The following code is giving me problems. I want to export densities
of a distribution to a csv file. At the bottom of the code I use the
hist function to generate the densities. But hist is returning values
greater than 1. I don't understand, why. Any help you can supply is
greatly
dear list,
i'm trying to install the package RCurl into a linux system running SUSE
Linux Enterprise Server 10 SP2 but i get compilation errors which i
guess should be due to some missing additional software. i've been
searching about the R-help archives without success so i hope somebody
can
On Thu, Jun 4, 2009 at 12:45 PM, Rainer M Krug r.m.k...@gmail.com wrote:
On Thu, Jun 4, 2009 at 12:36 PM, Patrick Burns pbu...@pburns.seanet.com
wrote:
How about:
order(x, runif(length(x)))
Thanks - that is really elegant.
One thing:
it is saver to use sample(length(x)) instead of
Dear R users,
I'm using the rpart package to build classification trees. I'm interested in
implementing the ordered Twoing as a splitting criterion.
Does anyone have experience with this task ?
Thank you for your help
Paolo
Paolo Radaelli
Dipartimento di Metodi Quantitativi per le Scienze
Hi,
If I understand your problem correctly, you didn't need to send us
your entire code.
A simple reproducible example that showed just your problem would have been
sufficient.
Take a look at this, and reread the help for hist().
testdata - runif(1000)
testdata.hist - hist(testdata)
Dear list,
I have problems involving passing a subset= argument through
I'm trying to augment the set of defined analyses for mice (homonymous
package) with a call to lme. This package create multiple imputations of
missing data in a mids object, each completed data set may be obtained
alamoboy wrote:
Newbie here. Many apologies in advance for using the incorrect lingo.
I'm new to statistics and VERY new to R.
I'm attempting to group or bin data together in order to analyze them
as a combined group rather than as discrete set. I'll provide a simple
example of the
Hi Steve,
Thanks for your response and also for the useful information about the
latest version warning.
On the earlier version that I was using (2.6.2 Win), I was expecting
at least an error or warning in response to submitting an obviously
corner case condition that should result in NaN, but
If you want to play R golf then:
sapply(HouseDatesList, '[[', 1)
[1] 1990 1991 1992 1993 1994 1995 1996
Does the same thing in fewer keystrokes, but Ronggui's solution is more
readable.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
Dear Helpers,
I need to fit a gamma function on a distribution. I want to use the Method
of the Least Squares for minimizing the sum of squared residuals (SSE). I
don't know how to do this. I guess I need to calculate the best fit
parameter values and then somehow comparing my empirical
Hi,
I tried fitting
loglinear model using the glm(catspec). The data used is FHtab. . An
independence model was fitted. Here summary() and fitmacro( ) give different
values for AIC.
I understand that fitmacro( ) takes the likelilhood ratio L2(deviance)
to calculate AIC and uses the
Hi there,
I trying to solve this problem for the whole day not going anywhere,
so I really hope maybe somebody can help
me in this community...
I've got an object coefficient2 which I want to plot in differerent
ways, with colors and labels added to the points,
but somehow there seems to be a
On Thu, 4 Jun 2009 15:21:52 +0100 Daniel Mail d20...@live.com.pt
wrote:
DM I need to know how to import ARIMA coefficients. I already
DM determined the coefficients of the model with other software, but
DM now i need to do the forecast in R.
So why then don't you fit the model then in R
Hi,
I'm a beginner in the field, I have to perform the ridge regression with
lm.ridge for many datasets, and I wanted to do it in an automatic way.
In which way I can automatically choose lambda ?
As said, right now I'm using lm.ridge MASS function, which I found quite simple
and fast, and
(copy of earlier reply which was not addressed to r-help. Also added
informative Subject:)
yes .. is not responsive to the question of how you propose to
display or examine such a mathematical object. My signature was
perhaps a lame (certainly an ineffective) effort at getting you to
No attachment came through the mailserver. Did you follow the
directions in the posting guide regarding acceptable types of
attachments? I also do not see any code that would let us reproduce an
input process.
One method that might work is to offer the results of
dput(coefficients) as
I still think that fortune(181) applies here. Someday you (or another user
that you give your function to) will run this, then realize that you/they had
an A, B, or c variable that has just been overwritten that you/they wanted to
keep. (also 'c' is one of the variable names recommended
Your output example below looks the same as the input. But I think the 'aperm'
function may be what you are looking for, read its help page and run the
example to see if that will work for you.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
sorry, the attachment can be found here for download:
http://www.wzw.tum.de/waldinventur/fileadmin/coefficient2.RData
katharina wrote:
Hi there,
I trying to solve this problem for the whole day not going anywhere,
so I really hope maybe somebody can help
me in this community...
I've
Also see this post at
https://stat.ethz.ch/pipermail/r-help/2008-July/169149.html
The same idea is discussed in a SAS proceedings paper (but it is NOT
specific to SAS)
Hello R users
http://cran.r-project.org/doc/manuals/R-intro.html
An Introduction to R
6.2 Constructing and modifying lists
After a short but successful struggle with nested
associative arrays (using named lists), I think the
following documentation change might be beneficial.
The
Could you provide a reproducible example even with fake data would be
fine or dput() yours.
On Thu, Jun 4, 2009 at 10:11 AM, Par Leijonhufvud
p...@hunter-gatherer.org wrote:
I want to make a log-log plot with a regression line, but I can't figure
out what I'm doing wrong. What I'm trying is:
Try adding a version that uses sort with the partial argument, that should be
faster than regular sort (for long enough test vectors) and possibly faster
than the max solutions when finding more than just the largest 2.
Also, for your max solutions, what happens when the 2 largest values are
Was not able to find NCStats on CRAN or BioConductor repositories. You
are wasting our (or mine at any rate) time by giving us dead ends.
You have been bitten by the factor snake in the grass. Color should
be text rather than a factor if you want the desired results.
See if this is
stephen sefick ssef...@gmail.com [2009.06.04] wrote:
Could you provide a reproducible example even with fake data would be
fine or dput() yours.
Sorry, I don't understand what you mean. Im my post was the first 8 lines of
my data, imported into R with
islands - read.table(islands.csv, sep=,,
Sorry for the problem with NCStats, I used it from this page not realizing
that this is not 'official':
http://www.rforge.net/NCStats/files/
Thank you for helping me with the color issue, I finally know now what has
happened with the
colors of the points my plots.
David Winsemius wrote:
Hi,
Is there an easy way to catch errors, in order to arrange for
r-scripts to exit gracefully?
I'm thinking of something along the lines of using is.na with an
if/else statement, but for errors.
Thanks!
__
R-help@r-project.org mailing list
Suppose I have a vector of strings. I'd like to convert this to a vector of
numbers if possible. How can I tell if it is possible?
as.numeric() will issue a warning if it fails. Do I need to trap this
warning? If so, how?
In other words, my end goal is a function that takes a vector of strings
On Thu, Jun 4, 2009 at 6:56 PM, Brigid Mooney bkmoo...@gmail.com wrote:
Hi,
Is there an easy way to catch errors, in order to arrange for
r-scripts to exit gracefully?
I'm thinking of something along the lines of using is.na with an
if/else statement, but for errors.
Try 'try'.
Barry
On Jun 4, 2009, at 12:56 PM, Brigid Mooney wrote:
Hi,
Is there an easy way to catch errors, in order to arrange for
r-scripts to exit gracefully?
I'm thinking of something along the lines of using is.na with an
if/else statement, but for errors.
Thanks!
See ?try and ?tryCatch
HTH,
Marc
I can't copy and paste your example right out of you email and into an
R session.
?dput
and then see if you can copy and paste it into an R session and make
it work. That way it is easier for everyone and you have a better
chance of getting helpful responses.
HTH
Stephen Sefick
On Thu, Jun 4,
On Jun 4, 2009, at 1:01 PM, Steve Jaffe wrote:
Suppose I have a vector of strings. I'd like to convert this to a
vector of
numbers if possible. How can I tell if it is possible?
as.numeric() will issue a warning if it fails. Do I need to trap this
warning? If so, how?
In other words, my
Roslina Zakaria wrote:
Hi Ravi,
I did ask you some question regarding newton method sometime ago.. Now I
have fixed the problem and I also wrote 2 looping code (ff1 and ff2) to
evaluate the modified Bessel function of the first kind and call them in
the newton code. But I dont't
Greg Snow wrote:
Try adding a version that uses sort with the partial argument, that should be
faster than regular sort (for long enough test vectors) and possibly faster
than the max solutions when finding more than just the largest 2.
I find the documentation for the partial argument in
I think the problem is that plot's log axes are to the base 10
so the lm() call needs to use log10, not log. E.g.,
x-101:200
y-sqrt(x)+1+runif(100)
plot(x,y,log=xy)
abline(col=red, lm(log(y)~log(x))) # nothing plotted
abline(col=blue, lm(log10(y)~log10(x))) # passes through points
Suppose I have an n-dimensional array and a logical vector as long as the
first dimension. I can extract an n-dimensional subarray with
a[ i, , , , .. ,] where there are n-1 commas (ie empty indices)
Is there an alternative notation that would better lend itself to more
generic use, e.g. to
Dear all,
I feel like I've been reinventing the wheel with this code (implementing
type = 'b' for Grid graphics),
http://econum.umh.ac.be/rwiki/doku.php?id=tips:graphics-grid:linesandpointsgrob
Has anyone here attempted this with success before? I found suggestions
of overlapping large
All,
There is a simple solution to this problem using R's matrix algebra
commands. I describe in detail at
http://www.stanford.edu/~messing/Affiliation%20Data.html
http://www.stanford.edu/~messing/Affiliation%20Data.html .
Best wishes,
-Solomon
--
View this message in context:
Thank you Peter,
This solved the problem.
edw...@web.de wrote:
Hi,
Unfortunately, they can't handle NA. Any suggestion? Some row for Ip
don't have ip address. This cause an error/ wrong result.
A quick fix could be to substitute ... or 0.0.0.0 for the NA
entries. (Use something
I don't want to plot or visualize the plot. I just want to clculate the
density at certain prespecified points.For that,knowin the funcional form is
enough.I just want a function that returns me back the estimate of density
at some point
On Thu, Jun 4, 2009 at 10:04 PM, David Winsemius
To address the creation of a 5-D density object:
http://finzi.psych.upenn.edu/R/library/ks/html/kda.kde.html
http://finzi.psych.upenn.edu/R/library/ks/html/kde.html
http://finzi.psych.upenn.edu/R/library/locfit/html/locfit.raw.html
--
David Winsemius
On Jun 4, 2009, at 12:34 PM, David
On Mon, Jun 1, 2009 at 2:18 PM, stephen sefick ssef...@gmail.com wrote:
library(ggplot2)
melt.updn - (structure(list(date = structure(c(11808, 11869, 11961, 11992,
12084, 12173, 12265, 12418, 12600, 12631, 12753, 12996, 13057,
13149, 11808, 11869, 11961, 11992, 12084, 12173, 12265, 12418,
thank you.
This is just what I was looking for.
On Fri, Jun 5, 2009 at 3:16 AM, David Winsemius dwinsem...@comcast.netwrote:
To address the creation of a 5-D density object:
http://finzi.psych.upenn.edu/R/library/ks/html/kda.kde.html
Hi,
I recently upgraded to R 2.9.0 on linux x86. After doing so, I switched to the
RPostgreSQL package for interfacing with a postgresql database. I am using
postgresql 8.3.7.
A query that works from the postgresql terminal is causing a segfault when
executed from R.
My sessionInfo, the
On 4 June 2009 at 16:17, Dylan Beaudette wrote:
| Hi,
|
| I recently upgraded to R 2.9.0 on linux x86. After doing so, I switched to
the
| RPostgreSQL package for interfacing with a postgresql database. I am using
| postgresql 8.3.7.
|
| A query that works from the postgresql terminal is
On Thursday 04 June 2009, Dirk Eddelbuettel wrote:
On 4 June 2009 at 16:17, Dylan Beaudette wrote:
| Hi,
|
| I recently upgraded to R 2.9.0 on linux x86. After doing so, I switched
| to the RPostgreSQL package for interfacing with a postgresql database. I
| am using postgresql 8.3.7.
|
|
HI,
I am Angshuman a postdoc in Buck Institute, Novato, CA. I am using random
forest in R. I have a problem. I have a training file and a test file. I need
to generate model file to classify a set of data of the test file. I need to
know the command for that. Please let me know.
Thank you.
Hi all,
I've been playing with ROracle (0.5-9) for a few days
and I can't wrap my mind around this one.
Here's a sample of my R (2.4.0) session.
my.df-data.frame(prd_id=c(123,456),vol_factor=c(.123,.456))
my.df
prd_id vol_factor
1123 0.123
2456 0.456
library(ROracle)
Loading
Thanks... This will be very usefull.
Warren Young wrote:
Jorge Cornejo wrote:
Hi, I am looking any way to communicant with Arduino
(http://www.arduino.cc/) Wiring (http://www.wiring.org.co) boards and
read data generate with sensor on these.
On Linux you can simply open the file
I have a question about antialiasing when R generates bitmaps. (This follows
a thread on the ggplot2 mailing list.)
I mostly use R on Linux, although I sometimes use it in Mac and Windows as
well. On Linux, I've found that plotting shapes 15-18 via cairo results in
bad-looking output. The points
On Thu, Jun 4, 2009 at 7:53 AM, Harsh singhal...@gmail.com wrote:
Hi UseRs,
I recently started working with the RGtk2 library. The documentation
is comprehensive and I've tried learning from the examples in most Gtk
tutorials (which have C code). This is a little problematic, but has
helped
I just checked every project in the Web interfaces section of the FAQ
(http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-Web-Interfaces ).
There are many dead links in this section. Of the links that work,
many projects themselves contain links that don't work, appear to have
not been maintained in
Hello all,
I had both curl and curl-devel
(both 7.15) installed on my x86_64/CentOS machine. However, I still got problem
when I do
R CMD INSTALL RCurl
Â
Error message below,
since library and/or include path is missing. I tried
R CMD INSTALL RCurl
Dear All,
Slightly off -non technical topic ( but hey it is Friday)
Following last week's interview with REvolution Computing which makes
enterprise versions of R, here is another interview with the rapidly
growing company Blue Reference CEOPaul van Eikeren at
Hi all,
I am new to R-project my problem is I tried to get subset from two different
tables its giving error
but if i m tring for geting results from one table its working
actually i have to take values from two tables with applying different
conditions on two tables like
kk- is an object of
Is it really necessary to further advertise this company which already
spams R-help subscribers?
Hadley
On Thu, Jun 4, 2009 at 10:41 PM, Ajay ohri ohri2...@gmail.com wrote:
Dear All,
Slightly off -non technical topic ( but hey it is Friday)
Following last week's interview with REvolution
It seems you need a later curl -- 7.15 is rather old (and I think
incomplete as a version number) and 7.19.5 is current.
It is helpful to tell us what version of packages (here RCurl) you are
trying to use, as well as your precise OS.
On Thu, 4 Jun 2009, heyi xiao wrote:
Hello all,
82 matches
Mail list logo