I am trying to replicate ARIMA models in R from SAS. Specifically -
Factored ARMA models, which I can do very easily in SAS (
http://v8doc.sas.com/sashtml/ets/chap7/sect11.htm).
I would like to know if there is a way to do this in any of the R packages.
Thanks in advance,
Indrajit
Hello,
Please have a look at:
http://cran.r-project.org/web/views/TimeSeries.html
Regards,
Pascal
2013/9/20 Indrajit Sengupta indrajitsg2...@gmail.com
I am trying to replicate ARIMA models in R from SAS. Specifically -
Factored ARMA models, which I can do very easily in SAS (
On 19/09/2013 15:28, Prigot, Jonathan wrote:
Sadly, I am limited to the Solaris 10 system. I wish that I could use
Linux, the world uses it.
What does that have to do with this?
The CRAN check farm uses gfortran 4.8 on a Solaris 10 system. g77 is
not the native compiler there, nor is that
Package news (see below for general description of functionality)
depmixS4 version 1.3-0 has been released on CRAN. See the NEWS file
for an overview of all changes. The most important user-visible
changes are:
1) more compact pretty-printing of parameters in print/summary of
(dep)mix objects
Hi
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of David Winsemius
Sent: Friday, September 20, 2013 6:39 AM
To: ivanc010
Cc: r-help@r-project.org
Subject: Re: [R] Need help to find out the name of my columns and rows
in a
Hi,
as far as I know, there is no limitation on data size in regard to foreach. You
should reserve though enough memory for your application on the cluster (via
ulimit -s unlimited and ulimit -v unlimited).
Furthermore I would check the following:
Check if there are two versions of R on the
Hello again,
I have one question on subscripting matrix. Let say I have following matrix:
Mat - matrix(1:9, 3)
colnames(Mat) - c(a, b, a)
Mat
a b a
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
Now I want to fetch data for colnames 'a'.I did following:
Mat[, a]
[1] 1 2 3
However it is not
On 20-09-2013, at 12:52, Christofer Bogaso bogaso.christo...@gmail.com wrote:
Hello again,
I have one question on subscripting matrix. Let say I have following matrix:
Mat - matrix(1:9, 3)
colnames(Mat) - c(a, b, a)
Mat
a b a
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
Now I want to
Hi,
it seems the following works as needed...
Mat[,colnames(Mat)==a]
Olivier.
On Fri, 20 Sep 2013 16:22:37 +0530
Christofer Bogaso bogaso.christo...@gmail.com wrote:
Hello again,
I have one question on subscripting matrix. Let say I have following
matrix:
Mat - matrix(1:9, 3)
Can you point to any function in specific?
On Fri, Sep 20, 2013 at 12:10 PM, Pascal Oettli kri...@ymail.com wrote:
Hello,
Please have a look at:
http://cran.r-project.org/web/views/TimeSeries.html
Regards,
Pascal
2013/9/20 Indrajit Sengupta indrajitsg2...@gmail.com
I am trying to
Hi Elio,
Try this:
library(stringr)
lines1-str_trim(gsub(\t, ,readLines(elio.txt)))
lst1-lapply(split(lines1,cumsum(lines1==)),function(x) x[x!=])
lst2- lapply(lst1[lapply(lst1,length)0],function(x)
as.matrix(read.table(text=x,row.names=1)))
names(lst2)- paste0(m,seq_along(lst2))
dat-
Replying with some code.
input- readLines(textConnection(
6.2 httpd (18)
4.0 httpd (11)
4.0 httpd (11)
3.3 httpd (9)
4.2 httpd (12)
4.2httpd (12)
4.2 httpd (12)
4.2 httpd (12)
4.2 httpd (12)
))
data-input[input!=]
data
[1] 6.2 httpd (18) 4.0
library(stringr)
read.table(text=gsub([()],,str_trim(input[input!=])),sep=,header=FALSE,stringsAsFactors=FALSE)
# V1 V2 V3
#1 6.2 httpd 18
#2 4.0 httpd 11
#3 4.0 httpd 11
#4 3.3 httpd 9
#5 4.2 httpd 12
#6 4.2 httpd 12
#7 4.2 httpd 12
#8 4.2 httpd 12
#9 4.2 httpd 12
A.K.
- Original
Hi,
I guess this is pretty basic.
I have 25 variables in the data file (name: score), i.e. X1,X2,.,X25.
I dont want to use score$X1, score$X2 everytime I use these variables.
Is there a way I can rename all these variables as simply X1,X2,.X25
without writing 25 lines of code, one
You are close--think names, not titles, as in rownames or colnames (no
reason to completely spell out column). Summary already gave you the
column names, so type ?rownames to learn more.
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET:
You haven't said yet, what object your 'data file' is. If you mean a data.frame
I would use colnames(dataName) - c(Col1Name, col2Name, ….)
Best
Simon
On Sep 20, 2013, at 4:10 PM, Preetam Pal lordpree...@gmail.com wrote:
Hi,
I guess this is pretty basic.
I have 25 variables in the data
Dear R help list,
I was just wondering whether there is a way to cluster the documentation files
of data sets in the package documentation index file, so that common prefixes
such as dat... are not necessary.
Best wishes,
Alrik
Dr.
On Fri, Sep 20, 2013 at 10:10 AM, Preetam Pal lordpree...@gmail.com wrote:
I have 25 variables in the data file (name: score), i.e. X1,X2,.,X25.
I dont want to use score$X1, score$X2 everytime I use these variables.
attach(score)
plot(X1, X2) # etc. etc.
-Aaron
[[alternative
On Sep 13, 2013, at 11:21 PM, E Joffe wrote:
Hi David,
First I ordered the levels of each factor in a descending order
based on
frequency.
Then, I used the following code to generate a matrix from the
dataframe with
dummy variables and subsequently run the glmnet (coxnet)
## tranform
or
with(score.plot(X1, X2))
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600FAX:(360) 407-7534
Olympia,
Dear R community,
I am having trouble changing the tick marks on y-axis to every 5 units? I
have the following:
x - c(12, 16, 6, 23, 27, 8, 5, 19, 23, 13, 16, 8)
y - c(29, 29, 23, 34, 38, 24, 22, 34, 36, 27, 33, 27)
plot(x, y, pch=19)
Should I change ylim=c(0,40), and then use axis()?
I kept
Hello everyone,
R beginner, I am confronted with the need to use Rmpf.
In my first scripts I made use of
X=read.table(file.choose(), header=FALSE, sep=,,dec=.)
X=as.matrix(X)
to load into a matrix data from file before matrix use.
How can I do to load the same data in a mpfrMatrix.
Is it
I have data that provide longitude, latitude, and local date and time but
no information about the corresponding time zone. How to identify the time
zone so they can be converted to a common date/time?
Thanks,
Carlisle
[[alternative HTML version deleted]]
Hi everyone!
I have now been using odfWeave() for a while in order to make some word
document reports of my data. I have encountered a problem when I try to
organize my program script into some functions. It seems that if I use the
odfWeave() within a function, then odfWeave() is not able to
please ignore this, I see the error.
greg
hi
probably a silly mistake, but I expected gam to minimise the penalised
deviance.
thanks
greg
set.seed(1)
library(mgcv)
x-runif(100)
lp-exp(-2*x)*sin(8*x)
y-rpois(100,exp(lp))
plot(x,y)
m1-gam(y~s(x),poisson)
hi
probably a silly mistake, but I expected gam to minimise the penalised
deviance.
thanks
greg
set.seed(1)
library(mgcv)
x-runif(100)
lp-exp(-2*x)*sin(8*x)
y-rpois(100,exp(lp))
plot(x,y)
m1-gam(y~s(x),poisson)
points(x,exp(lp),pch=16,col=green3)
points(x,fitted(m1),pch=16,cex=0.5,col=blue)
On 20/09/2013 11:52 AM, C W wrote:
Dear R community,
I am having trouble changing the tick marks on y-axis to every 5 units? I
have the following:
x - c(12, 16, 6, 23, 27, 8, 5, 19, 23, 13, 16, 8)
y - c(29, 29, 23, 34, 38, 24, 22, 34, 36, 27, 33, 27)
plot(x, y, pch=19)
Should I change
I initially posted this question to one of the StackExchange sites, and they
suggested that I repost my problem here.
After using ezANOVA as my primary way of specifying mixed ANOVAs, I've hit a
stumbling block when it come to adding a covariate to the model. I am using an
ANCOVA in order to
Hi,
I am struggling to save my figures as TIFF files (yes, the yournal only accept
TIFF and not any other format).
I can manage to save a simple plot as TIFF at the correct dpi (I need at least
600dpi) and compression works. Only not for my own plot. Everything looks
distorted. I have read
Hi,
May be this helps:
dat1- read.table(text=x1 x2 x3 x4
1 xz ab cd ef
2 ab fz cd ef
3 ab cd dy dx,sep=,header=TRUE,stringsAsFactors=FALSE)
dat1$changes_to_row_above- sapply(seq_len(nrow(dat1)),function(i)
{x1-dat1[,i]%in% dat1[,i-1];if(any(x1))
Hey all,
I've fitted two GAMs to some data using mgcv. The only difference between the
two models is that one includes an additional smooth term (the smooth terms are
s(x), s(y) and s(log(y)), the difference being that one model contains s(y) as
additional term whereas the other one only
If you make no further assumptions then this question is not solvable. For
example we use standard time in our data collection systems even though legal
time here applies daylight savings offset in the summer. In some cases I have
seen data collected from sites in multiple time zones recorded
On Sep 20, 2013, at 11:17 AM, Duncan Murdoch murdoch.dun...@gmail.com wrote:
On 20/09/2013 11:52 AM, C W wrote:
Dear R community,
I am having trouble changing the tick marks on y-axis to every 5 units? I
have the following:
x - c(12, 16, 6, 23, 27, 8, 5, 19, 23, 13, 16, 8)
y - c(29,
Thanks I'm lookin for yur example
-Message d'origine-
De : r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] De
la part de arun
Envoyé : vendredi 20 septembre 2013 18:11
À : R help
Objet : Re: [R] Compare two subsequent rows based on specific values of a
string
Hi,
May
HI,
we are pleased to announce 2 weeks intensive training in:
*SPATIO-TEMPORAL ANALYSIS AND BIG DATA PROCESSING USING FREE AND OPEN
SOURCE SOFTWARE*
16-20 December 2013 20-24 January 2014
University of Twente - Enschede - The Netherlands
4 ECTS
Over the few decades there has been an
On Sep 20, 2013, at 8:02 AM, carlisle thacker wrote:
I have data that provide longitude, latitude, and local date and
time but
no information about the corresponding time zone. How to identify
the time
zone so they can be converted to a common date/time?
Perhaps as an approximation you
Hi,
Are you using Linux? If so, you may use ImageMagick and try your luck
using
convert filename.whatever.format filename.tiff
There are lots of options in ImageMagick. (Read the manual.)
I have done this in the past, and it has worked. However, your exact
situation may be different, so
Thanks Duncan. I have no knowledge of SAS, though many in industry use it.
I hope R would expand its usage to more at the industry level.
Mike
On Fri, Sep 20, 2013 at 12:17 PM, Duncan Murdoch
murdoch.dun...@gmail.comwrote:
On 20/09/2013 11:52 AM, C W wrote:
Dear R community,
I am having
Dear list
I have a matrix composed of islandID as rows and speciesID as columns.
IslandID: Island A, B, C
.O (15 islands in total)
SpeciesID: D0001, D0002, D0003
.D0100 (100 species in total)
The cell of the matrix describes presence (1) or absence (0) of the species
in an island.
Now I
Once you learn to use dput() I am sure someone will be happy to help you.
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
John Kane
Kingston ON Canada
-Original Message-
From: elaine.kuo...@gmail.com
Sent: Sat, 21 Sep 2013 07:14:38 +0800
To:
Depending on what your doing perhaps with() could help?
Or assuming a data.frame or matrix, score[ , 25] will give you Score$X25
If you would supply a bit more information (and code) about what you are doing
we probably can help more.
John Kane
Kingston ON Canada
-Original Message-
Hi,
Try this:
set.seed(248)
lst1- lapply(1:1000,function(i) matrix(
sample(0:1,15*100,replace=TRUE),ncol=100,dimnames=list(paste(Island,LETTERS[1:15]),
paste0(D,sprintf(%04d,1:100)
lst2-lst1[sapply(lst1,function(x) any(colSums(x)==0))]
##The above steps are just to create some matrices
Hello all,
I´m have a very large data frame (more than 5 million lines) as below (dput
example at the end of mail):
Station Antenna TagDateTime Power Events
1 1 2 999 22/07/2013 11:00:2117 1
2 1 2 999 22/07/2013 11:33:4731 1
3 1 2
43 matches
Mail list logo