Hi R users,
Currently I want to fit a FIGARCH model to a dataset. The only package that
allow for it that I could find is fGarch. However it seems that the FIGARCH
model class fitting of this package has been moved to Oxmetrics. I tried to
install the old versions of it using 'tar.gz' files from
Original post: On 13/10/2013 18:53, Byron Dom wrote:
Due to convention a script l - $$\ell$$ (LaTeX \ell) is used to
represent a certain quantity in something I'm working on. I'm
unable to figure out how to use it in R. It's not included in the
list on ?plotmath.
Can anyone tell me how
Could you please show me what the code should be then?
--
View this message in context:
http://r.789695.n4.nabble.com/Function-not-working-as-I-d-like-tp4678878p4678898.html
Sent from the R help mailing list archive at Nabble.com.
__
Thanks for the help guys. Sorry, I will not post further RStudio questions
here. Yihui, thanks for the tip. Looks like these commands would do:
grDevices::pdf.options(useDingbats = FALSE); require(knitr);
opts_knit$set(concordance = TRUE); knit('filename.rnw',
encoding='UTF-8')
pdflatex
Dear All:
I have
v = c(xy, xy^2, z)
dotv = sum(v*v)
I hope, I have the following result:
dotv = x^2*y^2 + x^2*y^4 + z^2 (still in variable x, y and z, not
numeric result)
How to make simple script, so that I have: dotv = x^2*y^2 + x^2*y^4 + z^2
Thanks a lot for your help.
Best regards,
Yes, that is pretty much it. Setting the options useDingbats or
concordance is optional. You may or may not really need them. What is
essential is the knit() function.
Regards,
Yihui
--
Yihui Xie xieyi...@gmail.com
Web: http://yihui.name
Department of Statistics, Iowa State University
2215
Hi,
With the 'merge' function using different tolerance level, I've come up with a
solution.
For e.g. consider r1 and r2 as two raster layers, attributes are as follows:
r1 - raster(xmx=-150, ymn=60, ncols=30, nrows=20)
r1[]-1:ncell(r1)
r1
class : RasterLayer
dimensions : 20, 30, 600
Hi list,
In June 2013, using R.2.15.1 (i386) on Windows 7, I calculated a set of ecdf
which I stored in lists of records having a size.ecdf field, with following
structure :
[skip previous structure...]
$ size.ecdf :function (v)
..- attr(*, class)= chr [1:3] ecdf stepfun function
..-
Hi again, this is a follow-up of my previous post :
I've noticed that
attr(first.list[[1]]$size.ecdf,call) gives ecdf(test.moms[, m.pxs]) and
that there is no test.moms field in my records (a local variable named
test.moms was used inside a function, but it was stored under another name in
On 24/10/2013 07:54, Olivier Eterradossi wrote:
Hi list,
In June 2013, using R.2.15.1 (i386) on Windows 7, I calculated a set of ecdf which I
stored in lists of records having a size.ecdf field, with following structure
:
Note that R 2.15.1 was not current even then.
[skip previous
We would like to announce the following stats course;
Course: Introduction to MCMC, Linear mixed effects models and GLMM with R
When: 20-24 January, 2014
Where: Bangor University, UK
Info: http://www.highstat.com/statscourse.htm
Flyer: http://www.highstat.com/Courses/Flyer2014_02Bangor.pdf
Kind
this works for me:
plot(1,main=\u2113)
cheers
Am 24.10.2013 01:39, schrieb Byron Dom:
Original post: On 13/10/2013 18:53, Byron Dom wrote:
Due to convention a script l - $$\ell$$ (LaTeX \ell) is used to
represent a certain quantity in something I'm working on. I'm
unable to figure
Thanks Jim, maybe now I can start learning.
Next I want to use my choice of n digits or width instead of round.
--
View this message in context:
http://r.789695.n4.nabble.com/labeling-abscissa-using-a-function-of-the-plotted-scale-tp4678075p4678946.html
Sent from the R help mailing list
Dear, Patrick.
I'm using a workaround to work with RWeka and multicore package. The most
important point is not loading the RWeka package and use the namespace in a
encapsulated function. Take a look on this code:
# C4.5 classifier. Return the prediction for a test dataset.
cl.c45 =
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of wayne.w.jo...@shell.com
Sent: 22 October 2013 18:52
To: R-help@r-project.org
Subject: [R] nls model definition help
Hi fellow R users,
I'm trying to fit a model using nls
Hi Arun,
That seemed to do the trick - thanks!!
Jonathan
On Wed, Oct 23, 2013 at 11:12 PM, arun smartpink...@yahoo.com wrote:
HI,
Better would be:
res1 - dcast(df,gene~case,value.var=issue,paste,collapse=,,fill=0)
str(res1)
#'data.frame':2 obs. of 4 variables:
# $ gene : chr
Dear all,
I was wondering whether (or better: how) I can use R to read recursively a
directory to get all the sub-folders and files located in the root folder and
put it into a tree like structure where the leaves are files and intermediate
nodes are the directories? The idea is that I'd like
Hi all,
I am trying to plot a raster object (I can explain why but the point is that it
would be a raster objeçt)..
I have selected a small code to show you exactly the problem
require(raster)
test-matrix(data=runif(1),nrow=100)
m-raster(test)
plot(m,axes=FALSE)
A wild guess -- take a look at the CRAN phylohenetics task view, as
that sounds like the sort of thing that might have tree generation and
manipulation functions.
... but you may do better with some non-R tool out there.
(Hopefully, you'll get a better response, though).
Cheers,
Bert
On Thu,
Thanks Jim, maybe now I can start learning.
Here is a run of my trying to learn:
xvalue-c(5.2,1.3,9.7,2.8,8.1,4.7,6.6,7.4)
yvalue-c(9,3,4,7,2,5,3,6)
plot(xvalue,yvalue)
axis(1,at=NULL,labels=1/xvalue,digits=5)
Error in axis(1, at = NULL, labels = 1/xvalue, digits = 5) :
'labels' is supplied
One R package that might be of interest would be 'diagram':
http://cran.r-project.org/web/packages/diagram/
I would also agree with Bert here and would point you in the direction of
PSTricks, which can handle these sorts of complex figures. It would of course
require learning LaTeX, but that
perhaps as a somewhat tamer guess than Bert's, it's probably not
very hard to write a function that pieces together a dendrogram
(S3 class provided by the stats package), using the list.dirs and
list.files functions.
This could then be plotted via the plot function, and you could
even designate
This is by no means an explanation, but as a possible approach to consider:
plot(0:1, 0:1, xaxt='n', asp=1)
axis(1,at=c(0,1),labels=c(a,b))
plot(m,add=TRUE)
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 10/24/13
I'm writing a function that needs an optional variable. If the variable is
given, then a subset must be made using that variable, but if the variable
is not given the subset must be done using all the values for that variable
(ignoring that variable).
Something like this:
*band - function(file,
Dear all,
Running the following code produces an error in dummy.coef. I guess it might be
a bug..
ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt - c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group - gl(2, 10, 20, labels = c(Ctl,Trt))
weight - c(ctl, trt)
x -
Dear R-users,
I am trying to analyze data with ordinal independent variables and
binary response . Could some one suggest me an appropriate statistic
for this kind of analysis.
Thanks,
Regards,
Kishir
__
R-help@r-project.org mailing list
Hi professors
I'm finding the best package in R for manage a big database about vessels. I
need to divided the travel of each boat into tracks. I have GPS points.
I would join these points in a line on the map. I make that with the function
lines, but I can't divided by vessel or by track,
I
Berend Thanks. I am aware of the missing kb's in the download. However,
my problem is that it only happens with this package. I download other
packages manually and via command-line install.packages() and
Do not have any problems with completing full downloads of the entire file.
Hi All
I've gotten some awesome help getting a formular that finds the intersection of
two vectors. This works brilliantly, but I can't figure out how to make it run
over another factor. A simple example looks likes this:
df - data.frame(
id = factor(rep(c(supply, demand), each =
On 24-10-2013, at 16:15, David Mora dmv...@gmail.com wrote:
I'm writing a function that needs an optional variable. If the variable is
given, then a subset must be made using that variable, but if the variable
is not given the subset must be done using all the values for that variable
Kishor,
It sounds like you will need to use logistic regression. Although the following
URL might help you better understand logistic regression,
http://ww2.coastal.edu/kingw/statistics/R-tutorials/logistic.html it would
probably be most helpful if you would see help from a local statistician.
It is difficult to advise you when you don't provide reproducible data.
For example, you could include the first 6 lines of your data using
dput(head(mydata))
If you have GPS points stored as easting, northing, and if you have a
variable that identifies the vessel, and if your data are
On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
Hi R users,
Currently I want to fit a FIGARCH model to a dataset. The only package that
allow for it that I could find is fGarch. However it seems that the FIGARCH
model class fitting of this package has been moved to Oxmetrics. I tried to
install
Dear John,
Thanks for the information. It is quite helpful.
Regards,
Kishor
On Thu, Oct 24, 2013 at 9:24 PM, John Sorkin
jsor...@grecc.umaryland.edu wrote:
Kishor,
It sounds like you will need to use logistic regression. Although the
following URL might help you better understand logistic
Using your csv file, I get no error running the following code.
SPIV2 - read.csv(C:/Temp/SPI V2.csv)
heatmap(as.matrix(SPIV2), na.rm=TRUE)
Perhaps you modified SPIV2 in some way prior to running heatmap()?
Jean
On Tue, Oct 22, 2013 at 10:29 AM, David Romero
Sorry if this subject has been already dealt here.
Which are some common tests for nonparametric k-way ANOVA?
I have read about Kruskal-Wallis test as a kind of nonparametric one-way
ANOVA, but I have not found anything about a general-setting (I mean k-way)
nonparametric ANOVA.
Can you
Have you tried to download from another CRAN mirror, e.g.
http://cran.r-project.org/web/packages/XML/
More mirrors at http://cran.r-project.org/mirrors.html
As already others said, it's very unlikely that this is not an issue
on your end.
/Henrik
On Thu, Oct 24, 2013 at 4:54 AM, Steven
On Oct 24, 2013, at 11:38 AM, David Winsemius dwinsem...@comcast.net wrote:
On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
Hi R users,
Currently I want to fit a FIGARCH model to a dataset. The only package that
allow for it that I could find is fGarch. However it seems that the FIGARCH
Rick,
This uses a different approach than what you propose, but it gets the job
done and perhaps you will find it helpful.
Jean
library(jpeg)
# download image from internet for use in example
# http://www.tootsie.com/wallpaper/wp_concord2_1344x1008.jpg
img -
On Oct 24, 2013, at 18:59 , Vicent Giner-Bosch wrote:
Sorry if this subject has been already dealt here.
Which are some common tests for nonparametric k-way ANOVA?
I have read about Kruskal-Wallis test as a kind of nonparametric one-way
ANOVA, but I have not found anything about a
On 24/10/2013 18:25, Marc Schwartz wrote:
On Oct 24, 2013, at 11:38 AM, David Winsemius dwinsem...@comcast.net wrote:
On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
Hi R users,
Currently I want to fit a FIGARCH model to a dataset. The only package that
allow for it that I could find is
Thanks. That works for me too.
I had worked that out based on Brian's response to my original post.
If you read thru my response to him in detail, you'll see:
Here are a couple of examples using it that worked:
plot(1:10,xlab=\u2113)
plot(1:10,xlab=\u2113(\u2113 + 1))
On Oct 24, 2013, at 1:11 PM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote:
On 24/10/2013 18:25, Marc Schwartz wrote:
On Oct 24, 2013, at 11:38 AM, David Winsemius dwinsem...@comcast.net wrote:
On Oct 23, 2013, at 7:53 PM, Long Vo wrote:
Hi R users,
Currently I want to fit a FIGARCH
Obviously row 11 is incomplete.
It may contain nothing more than a space or an EOL.
Using the fill=TRUE argument to 'scan' or 'read.table' is your friend
here.
--
View this message in context:
http://r.789695.n4.nabble.com/Error-in-scan-tp4678976p4678982.html
Sent from the R help mailing
Dear R users,
Can I add a text to a plot without x axis?
I saw a function that enables placing the text with mouse
text(location, text to place, pos, ...)
locator(1)
but it did not work out.
--
Ahmed M. Attia
Research Assistant
Dept. Of SoilCrop Sciences
Texas AM University
Hi,
Try:
fun1 - function(dat){
mat1 - combn(colnames(dat1),2)
res - sapply(seq_len(ncol(mat1)),function(i) {x1- dat[,mat1[,i]];
wilcox.test(x1[,1],x1[,2])$p.value})
names(res) - apply(mat1,2,paste,collapse=_)
res
}
set.seed(432)
dat1 -
David,
You should cc r-help on all correspondence so others can follow the thread.
The error appears to be a result of the missing values yielded by the
dist() function on your matrix.
See https://stat.ethz.ch/pipermail/r-help/2009-July/204443.html.
Jean
On Thu, Oct 24, 2013 at 1:25 PM, David
Dear R users,
How can I customize x axis in a boxplot?
xaxt=n to suppress x axis, but I would like to customize x axis for
different groups.
Thanks
--
Ahmed M. Attia
Research Assistant
Dept. Of SoilCrop Sciences
Texas AM University
ahmed.at...@ag.tamu.edu
Cell phone: 001-979-248-5215
Depends on what the original poster really wants with a nonparametric
k-way ANOVA. If what they are really looking for is an approach that has
fewer distributional assumptions, then it is possible to perform
permutation tests for linear models parameterized for k-way anova that
would eliminated
You should provide an example of what you are trying to do. Try
this
plot(0, 1, xlab=, xaxt=n)
text(locator(1), This is my text, xpd=TRUE)
The xpd= parameter controls the clipping region. The default
xpd=FALSE clips text to the plot region. You were not clear, but
if you were trying to plot
On Oct 24, 2013, at 7:00 AM, Hurr wrote:
Thanks Jim, maybe now I can start learning.
Here is a run of my trying to learn:
xvalue-c(5.2,1.3,9.7,2.8,8.1,4.7,6.6,7.4)
yvalue-c(9,3,4,7,2,5,3,6)
plot(xvalue,yvalue)
axis(1,at=NULL,labels=1/xvalue,digits=5)
Error in axis(1, at = NULL, labels =
?axis
-
David L Carlson
Associate Professor of Anthropology
Texas AM University
College Station, TX 77840-4352
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Ahmed Attia
Sent: Thursday, October
Hello,
There's a bug in your function, it should be 'dat', not 'dat1'. In the
line marked, below.
fun1 - function(dat){
mat1 - combn(colnames(dat),2) # Here, 'dat' not 'dat1'
res - sapply(seq_len(ncol(mat1)),function(i) {x1- dat[,mat1[,i]];
wilcox.test(x1[,1],x1[,2])$p.value})
Hi,
You may try:
?list.files()
nm1 - list.files(pattern=.txt)
res - lapply(nm1,function(x) {
ln1 - readLines(x)
indx1 - grep(DATE PROCESSED,ln1)
indx2 - grep([A-Z],ln1)
There are several packages for track data on CRAN and they don't all
force the data to have specific variables beyond position,
date-time, track ID. With trip you can do this easily enough, but
depending on your needs adehabitatLT or move might be better. I wrote
trip so I know that you can do
A good choice might be to store your data in objects of class
SpatialLinesDataFrame
Then you can assign any attributes you want to each line, such as your
vessel, type, destination, or status.
This will make it relatively easy to, for example, plot subsets of tracks
specified by different
Hi,
You may try
?rollmean from library(zoo)
?SMA from library(TTR)
dat1 - structure(...
library(zoo)
res1 - rollmean(dat1[-c(1:5),2],8)
library(TTR)
res2 - SMA(dat1[-c(1:5),2],8)
all.equal(res1,res2[!is.na(res2)])
#[1] TRUE
A.K.
Hi,
How can I calculate running 8-hour averages of
I agree that this seems to be a bug, but OTOH the help for dummy.coef()
*does* say:
The method used has some limitations, and will give incomplete results
for terms such as |poly(x, 2)|.
(Note that it says it will give incomplete results not throw an error!)
Looking into the code of
On 10/25/2013 01:00 AM, Hurr wrote:
Thanks Jim, maybe now I can start learning.
Here is a run of my trying to learn:
xvalue-c(5.2,1.3,9.7,2.8,8.1,4.7,6.6,7.4)
yvalue-c(9,3,4,7,2,5,3,6)
plot(xvalue,yvalue)
axis(1,at=NULL,labels=1/xvalue,digits=5)
Error in axis(1, at = NULL, labels = 1/xvalue,
Thanks David, did I use incorrect english to call labels = 1/xvalue a
formula or function?
Seems to me an obvious first try would be to place the labels where
they would have been placed without the formula.
I haven't yet, but I will study from your suggestions.
--
View this message in
On Oct 24, 2013, at 1:51 AM, Alexandra Kuznetsova wrote:
Dear all,
Running the following code produces an error in dummy.coef. I guess it might
be a bug..
ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt - c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group -
On Oct 24, 2013, at 3:56 PM, Hurr wrote:
Thanks David, did I use incorrect english to call e a
formula or function?
Incorrect R terminology. A function in R is a specific sort of object and a
formula is a different sort of object. R terminology for the use of labels =
1/xvalue inside the
You could store your first approach in a function and lapply it to
your by_hour variable:
df - data.frame(
hour = factor(rep(1:5,4)),
id = factor(rep(c(supply, demand), each = 10)),
price= c(5,7,9,11,13,15,17,19,21,23,
20,18,16,14,12,10,8,6,4,2 ),
quantity =
Dear all,
Running the following code produces an error in dummy.coef. I guess it might be
a bug..
ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt - c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group - gl(2, 10, 20, labels = c(Ctl,Trt))
weight - c(ctl, trt)
x -
Hi Jean,
Problem result: the 1st column value must be complete for heatmap to run
correctly.
This solved the same problem I had with vegdist.
Thank you for your help.
Regards
David
De: Adams, Jean [mailto:jvad...@usgs.gov]
Enviado el: 24 de octubre de 2013 14:12
Para: David
Sorry for the html code:
I'm writing a function that needs an optional variable. If the variable is
given, then a subset must be made using that variable, but if the variable
is not given the subset must be done using all the values for that variable
(ignoring that variable).
Something like
Hi,
I'm trying to install R on CentOS 6.4.
Following some instructions online, I've done this:
rpm -Uvh
http://download.fedoraproject.org/pub/epel/6/i386/epel-release-6-8.noarch.rpm
yum install R
But yum fails, with this (full output below):
Error: Package: R-core-2.10.0-2.el5.x86_64
I'm trying to place labels so as to avoid covering points (and other labels) in
an ordination plot. I've been trying to use orditkplot() for this purpose, but
get an error message, even when I try to replicate the example:
require(vegan)
data(varespec)
ord - cca(varespec)
orditkplot(ord,
Hi,
Check out this function:-
pairwise.wilcox.test {package=stats}.
example(pairwise.wilcox.test)
On Fri, Oct 25, 2013 at 2:15 AM, Rui Barradas ruipbarra...@sapo.pt wrote:
Hello,
There's a bug in your function, it should be 'dat', not 'dat1'. In the line
marked, below.
fun1 -
It looks much better than mine.
with p value adjustment:
p.adjust(fun1(dat1), method = holm, n = 153)
#
dat1$id - 1:10
library(reshape2)
dat2 - melt(dat1,id.var=id)
with(dat2,pairwise.wilcox.test(value,variable))
with(dat2,pairwise.wilcox.test(value,variable,p.adj=none))
A.K.
On
Hello,
Do you have headers installed (generally tcl-devel and tk-devel)?
Regards,
Pascal
On 25 October 2013 10:33, Michael Stauffer mgsta...@gmail.com wrote:
Hi,
I'm trying to install R on CentOS 6.4.
Following some instructions online, I've done this:
rpm -Uvh
71 matches
Mail list logo