Matthew,
I would try first looking at the Rcdk package. From it sound like what you
want to do they already have some very efficient code to that. If you don't
want all of CDK then you could try Rdisop from Bioconductor. The package
will take a formula and give you it's mass. It won't do a
R help,
I would like to change the labels of my tick marks on the x-axis, and I am
having difficulty understanding how. Basically, my current axis labels
are: 0, 50, 100, 150, 200, and I would like: 1962Q1, 1974Q2, 186Q4, 1999Q2,
2011Q4.
With all respect to generality, could you please simply
I have modified a C function in the `phmm` package. Now I want to rebuild the
package, but the following error occurred. How can I fix the error?
$ cd ~/Desktop/phmmw
$ R CMD build phmm
* checking for file ‘phmm/DESCRIPTION’ ... OK
* preparing ‘phmm’:
* checking DESCRIPTION
On 01/30/2012 11:09 PM, Roopashree Shrivastava wrote:
Dear all,
I am new to plotting Taylor Diagram using plotrix package within R, hence
this post. I have written a script which plots Taylor Diagram with one
reference and 7 model values. However the font size, line width and line
type are not
thanks don
I have here enough to study for a while
thank you for your help
max
- Original Message -
From: MacQueen, Don macque...@llnl.gov
To: Massimo Bressan mbres...@arpa.veneto.it; r-help@r-project.org
Sent: Monday, January 30, 2012 4:47 PM
Subject: Re: [R] merge multiple data
You are using Mac OS X without telling us! See the appropriate FAQ at
http://cran.r-project.org/bin/macosx/RMacOSX-FAQ.html
And please post Mac-specific questions on R-sig-mac .
On 31/01/2012 06:21, moli wrote:
I have modified a C function in the `phmm` package. Now I want to rebuild the
On Mon, 30-Jan-2012 at 08:03PM -0800, Ryan Murphy wrote:
| R help,
|
| I would like to change the labels of my tick marks on the x-axis, and I am
| having difficulty understanding how. Basically, my current axis labels
| are: 0, 50, 100, 150, 200, and I would like: 1962Q1, 1974Q2, 186Q4,
On 01/30/2012 11:56 PM, gianni lavaredo wrote:
Hey Thanks Jim,
is It possible to plot in the same plot mean, standard error AND
standard deviation. As in figure 3 of the paper I attached
Hi gianni,
Here is a rather rough function that will do more or less what you want.
You can choose any
On 01/31/2012 08:54 AM, rkevinbur...@charter.net wrote:
Sorry. I am not sure how to post a link. Basically the legend looks like:
* * * * * * * * * * *
-4 -3 -2 -1 0 1 2 3 4
Where ' * ' are colored boxes that are right next to each other. Kind of
like a gradient.
Hi Kevin,
The color.legend
Hi Andrew
You could plot each section of your plot as a mini-line, something like this:
x = seq(1, 40)
y = seq(1, 40)^2
g = rep(c(A, B, A, B), each = 10)
newdf - data.frame(x = x, y = y, g = g)
newpanel - function(x, y, cols, ...) {
panel.xyplot(x, y, type = n, ...)
for(i in
On 01/30/2012 08:52 PM, Josiane NJIWA wrote:
Hello all,
I am very new to R and i am facing two problems. First i didn't succeed
changing the konsole language in english even after trying the line command set
language='en'.
I would like to plot ROC curves. I have a serie of 10 threshold
upgrade to the current mgcv version (see last item in changeLog for mgcv
1.7-12)
On 01/30/2012 10:04 PM, Casey Olives wrote:
Hello,
I'm trying to use the bam() function in the R mgcv package for a large set of
grouped binary data. However, I have found that this function does not take
data
On 01/31/2012 09:56 AM, rantree wrote:
Hi there ! I've got 7 Samples that may contain 29 differents kinds of
material every line below corresponds to a differente kind of material,
example: Sample 1, 2,3,5 and 6 has material 1 So I want to know how can
I do a Venn Diagram with the data
dear r-helpers,
i created an object named A, which looks like this:
A - c(4.4,1.9,4.1)
now i needed to get numbers instead of characters and for this i used the
function:
as.integer(A)
which resulted in:
[1] 4 1 4
My question is, why the numbers are rounded or more impotantly if
there is
as.numeric(A)
[1] 4.4 1.9 4.1
Integers are just integers...
On Tuesday 31 January 2012 12:21:48 Marion Wenty wrote:
dear r-helpers,
i created an object named A, which looks like this:
A - c(4.4,1.9,4.1)
now i needed to get numbers instead of characters and for this i used the
They are rounded to integers because you used as.integer(). Perhaps you meant
to use as.double()
Michael
On Jan 31, 2012, at 6:21 AM, Marion Wenty marion.we...@gmail.com wrote:
dear r-helpers,
i created an object named A, which looks like this:
A - c(4.4,1.9,4.1)
now i needed to get
Thank you very much for your help!
I was under the impression that numeric vectors and integer vectors are the
same!
Reading about integers just now I did not understand, what integer vectors
are. However, for my problem, I didn't need them, and used as.numeric
respectively as.double().
Thanks
The reason that it's not implemented is because of computational cost. Some
users had done it on their own using the same idea. It's just that it takes
too much memory for even moderately sized data. It can be done much more
efficiently in MART because computational shortcuts were used.
I'm not exactly sure if this is a problem with indexing by name; i.e., is the
following behavior by design? The problem is that names or dimnames that are
empty seem to be treated differently, and one can't index by them:
R junk = 1:3
R names(junk) = c(a, b, )
R junk
a b
1 2 3
R junk[]
NA
On Jan 30, 2012, at 8:07 PM, Eugene Gallagher wrote:
How do you label Xtick and Ytick marks with Currency symbols: $2000
instead of 2000? I would like to add dollar symbols to tickmarks on
boxplots, histograms and back-to-back histograms.
My Examples (requiring the lattice and Hmisc
arima is one of the largest R functions and not an easy one to
explain, but you can see its source simply by typing arima with no
parentheses at the prompt. To get the C functions that are called,
look here: https://svn.r-project.org/R/trunk/src/library/stats/src/
You might also check the
Hi All,
I want generate data using R that follows the shape of graphs (A and B) in
the attached file. Can anybody suggest me what function fits for each
graph?
Your help is highly appreciated in advance
Val
graphs.pdf
Description: Adobe PDF document
Hi Val
Look at the help file for
?curve
To get the plot.
You also need to decide what relevant function will fit the graphs you need.
Without knowing what your purpose is, I can not help you on that.
Contact
Details:---
I have a function of two arguments, (m, n) that returns a list with the
same structure as the function makefoo below.
When m=1, it returns a list of vectors, each of length n.
When m1, it returns a list of square matrices, each of size m x m.
makefoo - function(m, n) {
if (m==1) {
A
Hello Terry,
thank you for your help. I've tried your suggestion now and worked out the
following code. What I forgot to mention in my last mail is that I have
interval-censored data (I'm testing with discrete increasing fluid-volumes). I
would be thankful if you reviewed my solution:
myfun -
Hi Michael,
I must admit to reading this quickly, so I may have missed the point.
match.arg() takes character arguments. I think that you should be able to
use identical() to test whether the argument given is det or tr, as in
d - det
identical(d, det)
[1] TRUE
tr - function(x) sum(diag(x))
On Tue, Jan 31, 2012 at 10:03:39AM -0500, Val wrote:
Hi All,
I want generate data using R that follows the shape of graphs (A and B) in
the attached file. Can anybody suggest me what function fits for each
graph?
Hi.
The graphs leave open several questions. Are the distributions
Hi
I have a function
y= x^2
min =10
max=20
increment=0.1
I want to find the value of y for the value of x between min and max by step
increment.
how to get the values using apply function instead of for loop
-
Thanks in Advance
Arun
--
View this message in context:
Dear R users,
I am struggling to fit expo-linear equation to my data using nls
function. I am always getting error message as i highlighted below in yellow
color:
Theexpo-linear equation which i am interested to fit my data:
response_variable = (c/r)*log(1+exp(r*(Day-tt))), where
Dear All,
I'm new to R and my problem is this:
I'm reading a table which has list objects such as this:
col1 col2 col3 col4
1 354,5
2 2,3 6 5
NumericTable - as.numeric(Table)
This command doesn't work as it can't
Hi,
I have encountered this error when attempting a One-way Repeated-measure ANOVA
with my data.
I have read the Anova in car: SSPE apparently deficient rank thread
by I'm not sure the within-subject interaction has more degrees of freedom
than subjects in my case.
I have prepared the
I think I sorted it myself (sorry).
Whether this is the best way I dont know but it seemed to work. I changed
the table to a matrix using the as.matrix command. Then used the as.numeric
command on the matrix to change the lists to NA. Then I had to change the
numeric matrix back to a matrix of
Hi,
You can do that even without apply():
v.v-seq(10,20, by=0.1)
y-v.v^2
But if you want it with apply...
apply(as.matrix(v.v), 1, function(x) x^2)
Regards,
Carlos Ortega
www.qualityexcellence.es
2012/1/31 arunkumar akpbond...@gmail.com
Hi
I have a function
y= x^2
min =10
There must be a trick because I am not seeing a legend:
library(maps)
library(plotrix)
colors - grey(0:10/10)
mo - map('county', projection='polyconic')
leg.txt - c(-1, 0, 1)
color.legend(mo$range[1]-1,mo$range[2]-3,mo$range[2]+1,mo$range[2]-2,
+
Hi there,
This is a problem I've run into and do not know how to avoid. It
happens when I make an assignment using the dimension names as the
subscript of the array. The end result is a dimenensionless array
(i.e. a vector) which I don't want. See:
out - array(0, 5, list(1:5))
dim(out)
[1] 5
-Original Message-
From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk]
Sent: Tuesday, January 31, 2012 3:14 PM
To: Alexandris Theodoros
Cc: c...@r-project.org; Vlachos Nikolaos; Chatzivasiloglou Ioannis; Konrat
Aggeliki-Elfride
Subject: Re: FW: R project
This is the address for
Hello!
I have a dataset with monthly temperatures for 4 different years and I would
like to test if they are significantly different between the years. As I
collected the data for different sites I wondered if there were some
possibility to perform the calculation at once or if I have to repeat
Petr,
Thank you very much for help.
Graph A,
Are the distributions restricted to the shown intervals?
Not necessarily.
Based on your suggested R-code,
x - seq(0, 3, length=1001)
y1 - dnorm(x, mean=0.5, sd=1)
y2 - dnorm(x, mean=2.5, sd=1)
plot(x, (y1+1.05*y2)/2.05, type=l)
The following
Sorry this is the commands I am using:
colors - grey(0:10/10));
leg.txt - c(-1, 0, 1)
mo - map('county', projection='polyconic'));
color.legend(mo$range[1]-1,mo$range[2]+1,mo$range[3]-3,mo$range[3]-1, +
tmp - matrix(1:12, 3, 4, dimnames=list(letters[1:3], LETTERS[1:4]))
tmp
A B C D
a 1 4 7 10
b 2 5 8 11
c 3 6 9 12
tmp[a,, drop=FALSE]
A B C D
a 1 4 7 10
tmp[,A,drop=FALSE]
A
a 1
b 2
c 3
you need the correct number of commas.
For your example
dim(out[1, drop=FALSE])
[1] 1
out -
On Tue, Jan 31, 2012 at 12:40:35PM -0500, Val wrote:
[...]
What I want is,
1- let the plot star from 0.2 in Y-axis rather than the minimum value, Then
goes up to 0.23 then stay flat. A slow drop when it reaches to 0.25 on
X-axis. Finally, when it reaches at the coordinate of (0.21,0.3)
Hi petr,
Can the required density be understood as a piecewise
linear function going through 4 or 5 given points?
That is my problem. The function should be nonlinear. However, we can break
it down to the first 3 or 4 points could be linear and then nonlinear
function. On the later points can we
I found the problem. Operator error :(
Thank you.
-Original Message-
From: Kevin Burton [mailto:rkevinbur...@charter.net]
Sent: Tuesday, January 31, 2012 11:53 AM
To: 'Jim Lemon'
Cc: 'r-help@r-project.org'
Subject: RE: [R] Different type of legend?
Sorry this is the commands I am
Suppose I have two data frames A and B
A has three variables and B also has three variables.
I would like to merge these two database but the requirement to merge
is that the value of the second column in database A is less than the
value of the second column in database B. Is there a R code to do
can someone tell me if there is an easier way to do this in R - create a
design matrix? thanks.
design - model.matrix(~
-1+factor(c(1,1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5)))
# Creates an appropriate design matrix
can someone tell me
Which part? This might help:
rep(1:5, each=9)
On Tue, Jan 31, 2012 at 1:22 PM, Daniel Negusse
daniel.negu...@my.mcphs.edu wrote:
can someone tell me if there is an easier way to do this in R - create a
design matrix? thanks.
— design - model.matrix(~
On Tue, Jan 31, 2012 at 01:59:13PM -0500, Val wrote:
Hi petr,
Can the required density be understood as a piecewise
linear function going through 4 or 5 given points?
That is my problem. The function should be nonlinear. However, we can break
it down to the first 3 or 4 points could be
Hi Petr,
Thank you very much for the help,
Consider also to specify directly the inverse distribution
function using an increasing piecewise polynomial. Generating
numbers from the distribution is then immediate and computing
the graph of the density may be obtained using the formula
for the
Hello Everyone,
I am familiar with the use of 2SLS to mimic SEM. I even used this approach once
to estimate a model with a latent interaction term.
Does anyone know how to extend this approach to nested data? I have cancer
patients with measures of cancer symptoms, functional impairment, and
Hi,
Thanks for your comments, the issue remains unsolved though. I insist
the dropping of dimensions only appears to occur in assignments
out - array(0, 5, list(1:5))
out
1 2 3 4 5
0 0 0 0 0
dim(out)
[1] 5
dim(out[as.character(2:3)])
[1] 2
There is no dropping here, but then
I am trying to get the DSC proceedings from 2003.
The computer is not responding since last night (18 hours ago).
Is this transient? Have the DSC files been placed elsewhere?
Rich
http://www.ci.tuwien.ac.at/Conferences/DSC-2003/Proceedings
[[alternative HTML version deleted]]
GIYF - http://www.r-project.org/conferences/DSC-2003/Proceedings/
-Roy
On Jan 31, 2012, at 1:18 PM, Richard M. Heiberger wrote:
I am trying to get the DSC proceedings from 2003.
The computer is not responding since last night (18 hours ago).
Is this transient? Have the DSC files been placed
Roy,
Thank you, that is 80% of the way. Once you get to that page, then all the
abstracts
and papers are still linked back to www.ci.tuwien.ac.at
Fortunately, putting the name of the paper at the end of the link you gave
me
does construct a valid reference.
This now flips from a query to a bug
Hi there,
The command for merging is:
merge
Can you please clarify by which variable you wish to merge by?
Can you maybe use ?dput function in order to give a self contained R code
example of your data - and the output you wish to get?
Best,
Tal
Contact
On Tue, Jan 31, 2012 at 2:24 PM, Tal Galili tal.gal...@gmail.com wrote:
Hi there,
The command for merging is:
merge
A fortunes candidate?
-- Bert
Can you please clarify by which variable you wish to merge by?
Can you maybe use ?dput function in order to give a self contained R code
There are two packages which may be useful: venneuler and VennDiagram.
One of them allows up to four sets; the other is pretty much
unlimited. I forget which is which :-( .
Oops: there's another, evenn. I haven't used that one.
--
Sent from my Cray XK6
Pendeo-navem mei anguillae plena
An alternative is to use the ffDesMatrix(..) procedure in the BHH2
library (from CRAN mirrors) - it's drawn from _Statistics for
Experimenters_ by Box, Hunter, and Hunter (hence the name). It will
general a variety of design matrices - best for factorial and fractional
factorial. It has a
Ram
These errors often show up with starting values that are poor and wander
into 'illegal' territory (divide by zero, exp(large number),
log(negative or zero), etc. Consider putting constraints on the
parameters (?nls will guide you) and working harder to get starting
values by plotting the
Thanks for the reply. I am aware of those packages, as well as of a few
others (epiR, epinet, epibasix, epicalc). Unfortunately, they don't do this.
I'll try to get in touch with the authors of the Stata package, who knows.
Regards,
D.C.
-Message d'origine-
De : MacQueen, Don
I am having trouble generating a graph.
I want to know the % of respondents who answered that they strongly
agree or agree the America owes R's ethnic group a better chance
(BTTRCHNC) and I want to organize it by racial group (RACESHRT).
BTTRCHNC is organized ordinally from 1 through 5 with
Hi,
I have downloaded a bunch of bz2 files. I wonder if R will be able to
unzip them in a batch mode or even one at a time?
I was looking at the unz function. But it didn't work well. Say I have
a bz2 file in H:/Temp/65502805_532.pair.bz2. Anyone has any
suggestion?
thanks
John
Well i really need to do a Venn diagram :(
--
View this message in context:
http://r.789695.n4.nabble.com/Venn-Diagram-help-tp4342882p4345655.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
All,
I am attempting to select all of the contiguous elements of a matrix that meet
some criterion. I.e., values that would be contained within an irregular area
defined by a contour applied around point of interest. So, if I have a matrix
x as follows:
x - matrix(rnorm(25), nrow=5, ncol=5,
I may be wrong, but I don't think unz() handles bz2 files - only zip files.
See bunzip2() of the R.utils package (which utilizes bzfile connections).
/Henrik
On Tue, Jan 31, 2012 at 2:34 PM, ql16717 ql16...@gmail.com wrote:
Hi,
I have downloaded a bunch of bz2 files. I wonder if R will be
Try the bzfile() function instead of unz(). You'll find that you're
taken to the same help page for both of those functions.
?bzfile
On Tue, Jan 31, 2012 at 2:34 PM, ql16717 ql16...@gmail.com wrote:
Hi,
I have downloaded a bunch of bz2 files. I wonder if R will be able to
unzip them in a
Without knowing more about what you are trying to accomplish, and what
you have tried so far (code) it's difficult to say. But I'll venture
this: are you sure you need a graph? Seems like a table might suffice.
BT - sample(1:5, 50, replace=TRUE)
RA - sample(1:5, 50, replace=TRUE)
table(BT)
Prof Brian Ripley wrote:
'library' in R has a different meaning: I've altered the subject to be
more accurate 'libR'.
This is what R CMD SHLIB is for: it does all this for you in a portable way.
But if you want to DIY, you can use R CMD config to find out the
appropriate linker incantation.
Hi,
I am looking for a way to plot bar on a map instead of the standard points.
I have been using ggplot2 and maps libraries.
The points are added with the function geom_point. I know that there is a
function
geom_bar but I can't figure out how to use it.
Thank you for your help,
Simon
Can optimisation (simplex etc) be done through R?
__
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,
Dear R helpers,
I need to estimate a probit model with box constraints placed on several of
the model parameters. I have the following two questions:
1) How are the standard errors calclulated in glm
(family=binomial(link=probit)? I ran a typical probit model using the
glm probit link and the
On 02/01/2012 08:13 AM, Rebecca Lisi wrote:
I am having trouble generating a graph.
I want to know the % of respondents who answered that they strongly
agree or agree the America owes R's ethnic group a better chance
(BTTRCHNC) and I want to organize it by racial group (RACESHRT).
BTTRCHNC is
Dear R users,
I am new to Latex and I am using the R package xtable to generate tables.
I want to produce a table that is very long. in the landscape format, but
I would need to rescale the table so that it fits in the page. xtable
enables me to have the landscape format, but I cannot rescale it,
72 matches
Mail list logo
