Hi,
I am trying to print several dataset in several different sheets of the
same xlsx file.
I sometimes run the codes more than once, and I need to erase all things in
the target file, and I do not use "append" in the first function call.
However I want to write it in a loop. Is there a way to
Is it possible? I am expecting the result to be the second row of the data
frame ...
d<-data.frame(a=1:3, b=3:5,c=9:11)
> d
a b c
1 1 3 9
2 2 4 10
3 3 5 11
> d[identical(d[c("b","c")],c(4,10)),]
[1] a b c
<0 rows> (or 0-length row.names)
[[alternative HTML version deleted]]
Hi,
Could someone suggest a way to rename a data frame column? For example,
I want to rename the column beta, but I can't do it this way
> d <- data.frame(alpha=1:3, beta=4:6, gamma=7:9)
> mm<-"beta"
> rename(d, c(mm="two", "gamma"="three"))
The following `from` values were not present in
Thanks.
Could we print the row/column names, "alpha1" and "alpha2" to the csv file?
2016-04-30 17:06 GMT-07:00 Jim Lemon <drjimle...@gmail.com>:
> Hi jpm miao,
> I think you can get what you want like this:
>
> alpha1<-sample(LETTERS[1:3],50,TRUE)
> a
table(alphatab)
> sink()
> ?delim.xtab
2016-04-26 16:14 GMT-07:00 Jim Lemon <drjimle...@gmail.com>:
> Hi jpm miao,
> You can get CSV files that can be imported into Excel like this:
>
> library(prettyR)
> sink("excel_table1.csv")
> delim.table(table(df
Hi,
How could we print the frequency table (produced by "table") to an Excel
file?
Is there an easy way to do so? Thanks,
Miao
> df <- data.frame(x = 1:3, y = 3:1, z = letters[1:3])
> table(df[,c("y","z")])
z
y a b c
1 0 0 1
2 0 1 0
3 1 0 0
> test<-table(df[,c("y","z")])
>
Hi,
I have a data frame with two variables x, y, both of which take values
in the set {1,2,3}. I'd like to count the frequency by the command "table",
but exclude the value "1" in variable x, but keep "1" in variable y. Is it
possible? When I use "exclude", value 1 in both x and y are
Hi,
I tried to read a (small) xlsx file by "readWorksheetFromFile" function
of "XLConnect" package and "read.xlsx" function in "xlsx" package, but I
got this error message:
Error: OutOfMemoryError (Java): Java heap space
I tried to follow the solution on the web
Hi,
I am trying to add a vertical arrow (from top to bottom or from bottom
to up) to a time series plot using ggplot2 and xts. It seems that the
vertical line command "geom_vline" does not work for this purpose (Correct
me if I am wrong). I try the command "geom_segment" as follows, but I got
Hi,
I made a few graphs by ggplot. The following codes produce a pdf file
with graphs, sometimes in landscape orientation, sometimes in portrait
orientation. I am using both Mac and Windows PC.
Question: how can I control the orientation of the pdf file? I try to
add a line
Hi,
I made a few graphs by ggplot. The following codes produce a pdf file
with graphs in landscape orientation on my Windows PC, while they produce a
pdf file with the same graphs, but in portrait orientation:
*p2 <- lapply(1:(2*n), function(.x) xyz_outl[.x][[1]]) #a sequence of
graphs made
Hi,
I try to read a large xlsx file from by the function
"readWorksheetFromFile" in "XLConnect" package. Sometimes it works, but
sometimes it gives an error message. When it does not work and I try to run
the program line by line (in RStudio), then it would sometimes work, but
sometimes it
Hi,
The function "ggplot" does plot and return a plot. For example, we can
write:
y = ggplot(.. .) Then y is a plot.
How can we modify the multiplot function so that it can also return a
plot? Multiplot.R is here:
Hi,
I am outputting a large number of graphs (made by ggplot2) to a Word
file. The function "addPlot" in the package "rtf" does help output the
graphs sequentially in one column; can it be modified so that it outputs
graphs in more than one columns?
(I try R2wd package, but it does not
My code is:
if(type=="none")
type2<-"nc"
if(type=="drift")
type2<-"c"
if(type=="trend")
type2<-"ct"
I am wondering if there's a concise way to write a mapping from type to
type2, especially when the number of categories is high. Thanks!
[[alternative HTML version deleted]]
Hi,
I understand that there're many great graphic packages in R (e.g.,
ggplot2) . Nevertheless, my office uses Excel extensively. Is there any
package in R that produces Excel graphs by R codes? Thanks!
[[alternative HTML version deleted]]
Hi,
Given that I frequently use A, B, C, D, E packages.
In each of the 5 packages, I intend to modify one of the functions.
Two ways:
1st: Modify A, B, C, D, E to A', B', C', D', E'
Then I will use library(A') . library(E') instead of A, B, C, D, E.
or Just keep the names A,
Hi,
The package r2wd is good at writing MS Word document from R.
Can Chinese be written into MS Word by this package or any other method
from R? Are Chinese fonts (e.g., Kai 楷書) available in r2wd?
Thanks!
[[alternative HTML version deleted]]
Hi,
I try to build a package myself on RStudio. The command
package.skeleton is successfully run, but the command build is not. In
my case, the username and the package name are both abcd.
Error message:
R CMD build abcd
Error: unexpected symbol in R CMD
Could someone help me to
Hi,
If the relation between y and x is piecewise linear, the package
segmented can be used to model it. It works pretty well.
My situation is a little different: y and x are both I(1). Is there an
R-package dealing with this situation? If not, is there an econometrics
paper on this issue?
Hi,
I run a statistical test function in the package fUnitRoots that
returns a S4 object but I am wondering how to extract the p-value, one of
the output elements.
The document of the function urdfTest:
.
All tests return an object of class fHTEST with the following slots:
@call
.
)
2015-03-13 1:23 GMT+08:00 David Winsemius dwinsem...@comcast.net:
On Mar 12, 2015, at 1:04 AM, jpm miao wrote:
Hi,
I run a statistical test function in the package fUnitRoots that
returns a S4 object but I am wondering how to extract the p-value, one of
the output elements
Sorry. Let me modify the question: Does there exist any unit root test
function (with trend or intercept) where the p-value can be extracted? The
function adf.test in tseries package does return the p-value, but there's
no choice of trend or intercept. Thanks.
2015-03-13 10:49 GMT+08:00 jpm miao
Hi ,
I have a zoo object with several rows fx, r3m, etc. I would like to
create in the zoo object the lag of one column. One thing I can think of is
to create the lag first and then merge it into the original zoo object,
which is tedious. Can I add an NA column to the zoo object first and then
Hi,
My question is about NAs in the function gmm, but I believe that the
same issues occur in the case of lm.
I try to estimate a model by gmm function (GMM, generalized method of
moments). Each of the variables has 94 rows, but the resulting fitted
model has only 89 rows. Then the
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey
-Oorspronkelijk bericht-
Van: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
Namens jpm miao
Hi,
I could not find a nice lag operator on zoo object. Perhaps there is,
but I just couldn't find it. Basically I want the operator to return the
lagged zoo object (with one or more variables ) with the original date. For
example, if I write lag(x, -3), then I got the lagged series, but the
Hi,
I am plotting time series by ggplot2, but I believe that my question
applies to other plotting tool as well.
I want to make my x-axis the quarterly scale, e.g:
2000Q1 2000Q2.
However, scale_x_date and date_format(%m/%d) support all time formats
BUT QUARTERs
Hi,
I am working on zoo (time series) objects.
Is there any way to do a time series regression with a lag period?
E.g., Y(t) = b1*X1(t)+b2*X(t-1)+b3*X2(t)
Is dynlm the default one to use? Anything else
Thanks!
[[alternative HTML version deleted]]
Hi,
Is there any way to find the growth rate of a quarterly/monthly time
series? For example, I have a quarterly CPI level series, and I wonder how
to find the CPI growth rate (inflation rate). Is there easier way to
convert to another class and convert back?
Thanks!
dput(dtz[1:16,CPI])
I tried to convert a monthly dataset to quarterly by averaging
dt_mz-zoo(dt_m, dt_m$date_m)
dt_mz2q-aggregate(dt_mz, as.yearqtr, FUN=mean)
However, the resulting zoo object are just NAs
How can I get a quarterly average dataset?
Thanks!
dput(dt_mz)
structure(c(Jan 1981, Feb 1981, Mar 1981,
Hi,
I've not used R for about one year and don't know well about the updates
on the time series-related package.
My primary job is to do economic/financial time series data analysis -
annual, monthly, daily, etc. I usually read data by the package
XLConnect, which can read xls or xlsx files
Hi,
timeDate package create a date vector like this:
Dates - c(1989-09-28,2001-01-15,2004-08-30,1990-02-09)
I have a date whose size is large. Could this package read the dates
from xls or txt files? Could we convert the read vector (e.g., I usually
use XLConnect to read xls files) to
Thanks Pascal! It certainly helps!
2014-09-29 17:10 GMT+08:00 Pascal Oettli kri...@ymail.com:
Hi Miao,
You certainly will find useful answers here :
http://cran.r-project.org/web/views/TimeSeries.html
Regards,
Pascal Oettli
On Mon, Sep 29, 2014 at 6:05 PM, jpm miao miao...@gmail.com
#. .OO#. rocks...1k
---
Sent from my phone. Please excuse my brevity.
On September 29, 2014 2:47:01 AM PDT, jpm miao miao...@gmail.com wrote:
Hi,
timeDate package create a date vector like this:
Dates - c
Hi,
Is there any optimization package in R recommended for general purposes
(constrained, unconstrained, etc)?
Has anyone used the function in NAGFWrapper package
http://www.nag.co.uk/numeric/R/r-package
Are the optimization functions therein available for free?
Thanks,
Miao
Hi,
I would like to fit my data with a 4th order polynomial. Now I have only
5 data point, I should have a polynomial that exactly pass the five point
Then I would like to compute the fitted or predict value with a
relatively large x dataset. How can I do it?
BTW, I thought the model
Hi,
I find a few Pearson correlation test functions like
fBasics::correlationTest or stats::cor.test
which give the p-value of the test result. Is there a function that
calculate the cutoff correlation value for a specific p-value , e.g., p =
0.05?
I have a plot for the cross
Hi,
I am building a bivariate SVAR model
y_1t=c_1+Ã_1 (1,1) y_(1,t-1)+Ã_1 (1,2) y_(2,t-1)+Ã_2 (1,1) y_(1,t-2)+Ã_2
(1,2) y_(2,t-2)+å_1t
b y_1t+ y_2t=c_2+Ã_1 (2,1) y_(1,t-1)+Ã_1 (2,2) y_(2,t-1)+Ã_2 (2,1)
y_(1,t-2)+Ã_2 (1,2) y_(2,t-2)+å_2t
Now y1 is relatively exogenous in that y1
Hi,
I attempt to conduct the unit root test with the two functions
unitrootTest(x, lags = 1, type = c(nc, c, ct), title = NULL,
description = NULL)
adfTest(x, lags = 1, type = c(nc, c, ct), title = NULL,
description = NULL)
According to the document,
lagsis the maximum number of
Hi,
I am wondering if C++ programs could be compiled in RStudio. I search on
the web, and I find many yes's and no's. It looks like Rcpp can do it (or
partially?). I just wonder if it can replace the C++ IDE, e.g., Eclipse or
Visual Studio since RStudio is much easier to use.
One says yes:
Is B a reserve word in vars package?
I tried to run an SVAR-AB model by SVAR function and find the IRF in vars
package. The problem is that when B matrix is named by B, an error
message occurs. However, if the same matrix is named by Bm, then things
run smoothly. What's wrong? Is B a reserve
Hi,
I have a matrix whose columns are named as in and out. Then I coerce
it to be a data.frame. However the system seems to forbid me from using the
name in, but I am not aware of it until I call it by the dollar sign $.
Is there something R should remind me but it does not?
Is there any
Just wonder why R does not remind me so when I use it as a dimname...
2013/7/11 Pascal Oettli kri...@ymail.com
Hello,
in is a reserved word.
?Reserved
Hope this clarifies,
Pascal
2013/7/11 jpm miao miao...@gmail.com
Hi,
I have a matrix whose columns are named as in and out
, at 7:18, jpm miao miao...@gmail.com wrote:
Hi,
I try to build a toy package by running the following codes in an R
program
require(stats)
f - function(x,y) x+y
g - function(x,y) x-y
d - data.frame(a=1, b=2)
e - rnorm(1000)
package.skeleton(list=c(f,g,d,e), name=test1pkg,
path
,
Did you run R CMD build test1pkg in a command prompt window?
Regards,
Pascal
On 20/06/13 17:01, jpm miao wrote:
Hi,
Following your advice, I install R tools from
http://cran.r-project.org/bin/**windows/Rtools/http://cran.r-project.org/bin/windows/Rtools/
Rtools215.exe http://cran.r
/6/20 Berend Hasselman b...@xs4all.nl
On 20-06-2013, at 10:01, jpm miao miao...@gmail.com wrote:
Hi,
Following your advice, I install R tools from
http://cran.r-project.org/bin/windows/Rtools/
Rtools215.exe
http://cran.r-project.org/bin/windows/Rtools/Rtools215.exe
Files
arch.test in vars package with multivariate.only = FALSE
The function
arch.test(x, lags.single = 16, lags.multi = 5, multivariate.only = FALSE)
produces the univariate and multivariate parts of the ARCH test
By_t = C+A_1 y_{t-1}+
.. A_n y_{t-n}+\epsilon (homoskedastic VAR)
And then
Hi,
I try to build a toy package by running the following codes in an R
program
require(stats)
f - function(x,y) x+y
g - function(x,y) x-y
d - data.frame(a=1, b=2)
e - rnorm(1000)
package.skeleton(list=c(f,g,d,e), name=test1pkg,
path=D:/R/pkgtest)
Then the program runs smoothly
Creating
Hi,
I am trying to build a package with package.skeleton function.
I already have the function quadprod2.R in the current folder. After
running the program,
library(frontier)
source(quadprod2.R)
package.skeleton(name=sfa_ext)
package.skeleton(name=sfa_ext)
Creating directories ...
Hi,
I have a structure, which is the result of a function
How can I access the elements in the gradient?
dput(test1)
structure(-1.17782911684913, gradient = structure(c(-0.0571065371783791,
-0.144708170683529), .Dim = 1:2, .Dimnames = list(NULL, c(x1,
x2
test1[[1]]
[1] -1.177829
test1
Hi,
From time to time I need to do the aggregation. To illustrate, I present
a toy example as below. In this example, the task is to aggregate x and y
by z with the function mean.
Could I call the aggregation function with x_test, where
x_test=c(x,y)? Thanks
Miao
Hi,
I have a few graph objects created by some graphic package (say, ggplot2,
which I use frequently). Because of the existent relation between the
graphs, I'd like to index them in two dimensions as p[1,1], p[1,2], p[2,1],
p[2,2] for convenience.
To my knowledge, the only data type capable
Hi,
I have a sequence whose 1st, 3rd, 4th, 6th are non-NAs. How could I let R
return 1,3,4,6, the indices?
I know only how to find the non-NA elements. Thanks,
Miao
test-c(2,NA,6,8,NA,12)
test[is.na(test)==FALSE]
[1] 2 6 8 12
[[alternative HTML version deleted]]
Hi,
How could I see the source code of functions in an R package?
If we type ?function_name , we will see documentations of the
function_name.
If we type function_name, is what returns just the source code? Could we
just save it in an .R file and modify as we want? However, it seems that
Hi,
As an example, how can I get the data such that field a of ab, ab[a],
equals 3? I expect the answer to be the union of 2 and 4, as
Thanks,
a-c(1,3,4,3,5,6,5)
b-c(2,4,6,7,3,1,2)
ab-data.frame(a,b)
ab
a b
1 1 2
2 3 4
3 4 6
4 3 7
5 5 3
6 6 1
7 5 2
ab[a==3]
Error in
and the median of the data with the
assumption that all data points inside each interval are evenly
distributed. Are there such functions in R?
70-80480-90590-1008100-1107110-1203
2013/5/8 David Winsemius dwinsem...@comcast.net
On May 7, 2013, at 8:40 PM, jpm miao wrote:
Is there a function
Thanks, but why does f[f=4] yield 16 instead of 4?
f
[1] -24 -8 16 -32 64 -128
f[f0]
[1] -2 -8 -32 -128
f[f0]
[1] 4 16 64
f[f=4]
[1] 16
2013/5/8 Jorge I Velez jorgeivanve...@gmail.com
f [ f 0 ]
On Wed, May 8, 2013 at 11:54 AM, jpm miao miao...@gmail.com wrote:
Hi
Sorry. I just got it
f[f==4]
[1] 4
2013/5/8 jpm miao miao...@gmail.com
Thanks, but why does f[f=4] yield 16 instead of 4?
f
[1] -24 -8 16 -32 64 -128
f[f0]
[1] -2 -8 -32 -128
f[f0]
[1] 4 16 64
f[f=4]
[1] 16
2013/5/8 Jorge I Velez jorgeivanve
Hi,
From time to time I merge two dataframes with possibly a common field.
Then the common field is no longer present,but what are present fieldname.x
and fieldname.y. How can I fix the problem so that I can still call by the
orignal fieldname? If you don't understand my problem, please see
Hi,
I have a vector f with some negative columns. I remember that there is
an easy expression that can find out negative items. Can someone tell me
how I can do it?
It seems to be
f[i such that f[i]0 ...]
Thanks,
Miao
d-1:7
f-(-2)^d
f
[1] -24 -8 16 -32 64 -128
Is there a function in R that calculate the mean and median for a grouped
data?
For example, a survey shows the oil price outlook in the future. How can I
calculate the mean/median?
(Of course, I understand that the groups below 80 and above 110 must be
defined more specifically)
below 80 4
Hi,
Below is the output from an R package. The first column (4, 5, 6, 7,
which is unnamed) is the company name (code), while the second column
efficiency is the performance of each company, which is the primary output
result. How can I access the first column, rowname, 4, 5, 6, 7, 8, 9, 11?
)
# convert all NA fields to zeroes if you prefer
temp[ is.na( temp ) ] - 0
On Wed, May 1, 2013 at 11:55 PM, jpm miao miao...@gmail.com wrote:
Hi,
Attached are two datasheet to be read.
My raw data 130502temp.xlsx contains numbers with ' symbols, and they
can't be read as numbers
Hi,
I would like to read several datasets and would like to create a set
(list? sequence?) of many empty dataframes. How could this be done? How
could I declare a set (list? sequence?) of many empty matrices?
Thanks,
Miao
[[alternative HTML version deleted]]
I have a program, when I write
if(num!=NA)
it yields an error message.
However, if I write
if(is.na(num)==FALSE)
it works.
Why doesn't the first statement work?
Thanks,
Miao
[[alternative HTML version deleted]]
__
R-help@r-project.org
/crp4700/coelli_Intro_effic.pdf
According to my friend working on panel data, the case where efficiency
1 is seldom/never encountered, but the book says the opposite. Thanks,
Miao
2013/4/25 Arne Henningsen arne.henning...@gmail.com
Dear Miao
On 25 April 2013 03:26, jpm miao miao...@gmail.com
Hi,
Attached are two datasheet to be read.
My raw data 130502temp.xlsx contains numbers with ' symbols, and they
can't be read as numbers. Even if I copy and paste as numbers to form a new
file 130502temp_number1.xlsx, they could not be read smoothly.
1. How can I read the datasheet as
Hi,
I am using Extrafont package to install more fonts for my graphs. My
primary graphic tool is ggplot2. I seem to have problem installing the
package, but could not pinpoint where it is.
I try to follow the instruction here:
https://github.com/wch/extrafont
I guess (but am not sure)
Hi,
I'd like to print a string vertically. For example, I would like to print
abcd as a\nb\nc\nd
Is there a function in R such that
Input: abcd
Output: a\nb\nc\nd?
Thanks,
Miao
[[alternative HTML version deleted]]
__
Hi,
I am trying to find out the scale efficiency and optimal scale of banks
by stochastic frontier analysis given the panel data of bank. I am free to
choose any model of stochastic frontier analysis.
The only approach I know to work with R is to estimate a translog
production function by
Hi,
I have two questions on the function readNamedRegionFromFile in
XLConnect pacakge.
1.
In the documentation,
# multiregion xlsx file from demoFiles subfolder of package XLConnect
demoExcelFile - system.file(demoFiles/multiregion.xlsx,
package = XLConnect)
I have many xls grade report sheets with the same format
XXX High School Grade Report
Confidential Yes
Math English Science
John 90 85 90
Mary 75 88 93
Since the reports are prepared on a regular basis, I have many reports with
identical format. Without the title XXX High
HI,
I have a dataframe with two variable A, B. I transform the two variable
and name them as C, D and save it in a dataframe dfcd. However, I wonder
why can't I call them by dfcd$C and dfcd$D?
Thanks,
Miao
A=c(1,2,3)
B=c(4,6,7)
dfab-data.frame(A,B)
C=dfab[A]*2
D=dfab[B]*3
Hi,
Is it possible to create Dual y axes with ggplot2 package? If so, is
there any sample codes that reaches (or partially) this goal?
The author of the package, Dr Hadley Wickham, seem to discourage the
dual axes, but some remedies were provided by other users (If I am wrong,
please
of avg_cost_2012 ~ asset_2012, but do
not see the line for y1 against x2. y1 and x1 are numeric vectors of the
same size. How can I fix it?
Thnaks
Miao
2013/4/10 David Winsemius dwinsem...@comcast.net
On Apr 9, 2013, at 8:21 PM, jpm miao wrote:
Thank you very much.
Could it be done in Lattice package
, 4365158.32240167, 4466835.92150963,
4570881.89614875, 4677351.41287198, 4786300.92322638, 4897788.19368447,
5011872.33627272), .Dim = c(201L, 2L), .Dimnames = list(NULL,
c(y1, x2)))
2013/4/11 David Winsemius dwinsem...@comcast.net
On Apr 10, 2013, at 12:44 AM, jpm miao wrote:
Hi David,
Many
Hi David,
Thanks. How can I make the graph in log scale?
Thanks
Miao
2013/4/11 David Winsemius dwinsem...@comcast.net
On Apr 10, 2013, at 6:11 PM, jpm miao wrote:
Thanks for your comments, David
avg_cost_2012 and asset_2012 are numeric vectors of length 39
while x2 and y1
Hi,
How can I plot two curves with distinct x and y vectors? I would like to
join one of them by regular lines and plot the other just by points (no
lines). Can it be done in regular R graphic tools, say plot function?
Can it be done in Lattice package, say xyplot function?
Thanks,
Miao
use the following command.
plot(x1,y1,type=l)
points(x2,y2)
Hope it helps.
On Tue, Apr 9, 2013 at 8:24 PM, jpm miao miao...@gmail.com wrote:
Hi,
How can I plot two curves with distinct x and y vectors? I would like
to
join one of them by regular lines and plot the other just by points
Taking the Grunfeld data, which is built-in in R, for example,
(1)How can I construct a dataset (or dataframe) that consists of the data
of all firms in 1951?
(2)How can I calculate the average capital in each form over the period
1951-1954?
What I can imagine is to categorize the data by firm,
Hi,
I conduct a panel data estimation and obtain estimators for two of the
coefficients beta1 and beta2. R tells me the mean and covariance of the
distribution of (beta1, beta2). Now I would like to find the distribution
of the quotient beta1/beta2, and one way to do it is to simulate via the
Hi,
I am graphing with the following command in the Lattice and LatticeExtra
package
xyplot(xts,lty=c(1,2),col=c(blue,red),type=c(l,g),par.settings =
list(layout.heights = list(panel = c(2, 2))), aspect=xy,xlab=,ylab=%,
key=key1,screen=list(a,a,b,b,c,c,d,d), layout=c(2,2),
Dear All,
I attempted to fit a DVEC and a BEKK multivariate GARCH model, but am
wondering which package to use.
1. I tried to use rmgarch package in R, but I couldn't find the
subroutines for DVEC and BEKK.
2. I tried to find rmgarch package of R, which is not located on the
official R
Hi ,
Are there any multivariate GARCH package (e.g., BEKK) in R?
I do see a few projects constructing multivariate GARCH package, but
wonder if there is anyone that is ready for use.
Thanks,
Miao
[[alternative HTML version deleted]]
Hi everyone,
Do we have GARCH-in-mean VAR in R?
Is there any GARCH package in R?
Thanks,
Miao
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hello,
I have a table (in a txt file) which look like this:
Monday 12 78 89
Tuesday 34 44 67
Wednesday 78 98 2
Thursday 34 55 4
Then the table repeats Monday , Tuesday, ... followed by several numbers
My goal is to read values after the table. My problem is a little more
complicated,
Hi,
I am using xylot function in Lattice package. I find it excellently
written, but I don't know how to find resources on this function. One thing
I would like to change is the line type. To my knowledge, type =l means a
regular line, type g means broken lines, but I can't find a complete
Hello,
I work with time series data. From time to time I run programs to
produce results that are in time series form (e.g., quarterly or monthly
data). After a few days I might need to access part of the results and to
run another program. Is there any function or package (like dataframe or
Hello,
Is there a function in R by which one can run AR model with Bayesian
Information Criteria (BIC)? To my knowledge, functions ar and ar.ols could
select the order only by AIC.
Thanks,
Miao
[[alternative HTML version deleted]]
__
population, visit the Age UK
Knowledge Hub
http://www.ageuk.org.uk/professional-resources-home/knowledge-hub-evidence-statistics/
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of jpm miao
Sent: 27 September 2012 09:43
To: r-help
Hi,
I am trying to fit an AR model, maximum order =4, order selection
criterion is aic. I wonder why these two give different results:
m1-ar.ols(x, aic=TRUE, method=ols, order.max=4)
m1-auto.arima(x,d=0, D=0, max.p=4, max.P=0, max.q=0, max.Q=0, ic=aic)
Could they both use the
To dear Dr Sarkar and anyone that knows about Lattice package,
I make 4 graphs by Lattice package. Each of the graphs has two time
series. All the series are plotted in plain lines by default, and I would
like one series to be in plain line and the other to be in dotted line in
each graph. How
Hello,
I would like to write a program that compute the principal components of
a set of data and then
1. Run the dependent variable against the principal components (lagged
value)
2. Do prediction
, following Stock and Watson (1999) Forecasting Inflation. All data
are time series.
Hi,
I am reading several .csv file of time series. Some values of the time
series are unavailable. Some of the unavailable series are marked NA, and
some are just blank. The NA value can be read, but the blank cannot be
read. Then I change the change the blank to NA, but they can't be
properly
Hi,
I find the problem. Some of the numbers like 1357 are presented as
1,357. After deleting the comma, everything is ok. Thank you for your
attention!
Miao
2012/8/22 jpm miao miao...@gmail.com
Hi,
I am reading several .csv file of time series. Some values of the time
series
Hi ,
To my knowledge, there're two functions that can do principal component
analysis, princomp and prcomp.
I don't really know the difference; the only thing I know is that when
the sample size number of variable, only prcomp will work. Could someone
tell me the difference or where I can
Hi,
To my knowledge, Ryacas can do symbolic computation in R, like
Mathematica or Maple.
I wonder if it (or any other R package) can solve symbolically a system
of equations? I have a system of four equations in four variables
(non-linear, but not very hard to solve) and wonder if Ryacas
Hi,
Part of my program is to calculate the number of time series in a zoo
object. It works well if it has more than one time series, but it fails if
it has only one. How can I access the number of column (i.e. the number of
time series) when I have only one column? Why is the number of an
names of the zoo
objects? In case of more than one time series, the function colnames
works, but not for the univariate time series.
Thanks,
Miao
2012/8/9 Achim Zeileis achim.zeil...@uibk.ac.at
On Thu, 9 Aug 2012, jpm miao wrote:
Hi,
Part of my program is to calculate the number
1 - 100 of 132 matches
Mail list logo