Dear R-helpers,
I am using plm package to run panel data analysis for different models
with different dependent variables. After running the PLM, I run the
COEFTEST(plm1, vcovBK) to get the panel corrected standard errors. I got
the PCSEs for some of the models but some tests produce error
Just a heads up: I'm pretty sure Josh and the other QS developers have
access to intra-day data, but they can't provide an example using it
in the package because the EULAs of most data providers won't allow
direct redistribution of their data. The examples included all go
online and download data
On 08-04-2012, at 08:28, Navin Goyal wrote:
Dear R users,
I am running a loop with the integrate function. I have pasted the code
below. I am integrating a function from time=0 to the time value in every
row.
I have to perform this integration over thousands of rows with different
Dear R-listers,
I am an MD and clinical epidemiologist developing a measure of comorbidity
severity for patients with liver disease. Having developed my comorbidity score
as the linear predictor from a Cox regression model I want to compare the
discriminative ability of my comorbidity measure
Hi all,
I have this problem with lattice that xyplot() won't draw some of my axis
labels if the type (i.e. the relation argument) of scales is set as free. For
example, in the plot below, I would want it to also show:
1. the labels E1,...E6 below the 10th panel (i.e. 3rd row, 2 col)just as
Sorry, the B data.frame() didn't copy/paste well. Here is a working version
(hopefully)
B-structure(list(yval = c(0.7, 0.61, 0.65, 0.63, 6.08, 0.64, 5.68,6.77,
1.48, 7.71, 0.82, 1.15, 0.54, 1.01, 0.59, 4.84, 0.69, 0.71, 8.7, 0.48, 0.69,
4.81, 1.42, 1.19, 0.84, 4.89, 0.85, 0.67, 7.07, 0.66,
Hi
I use outline mode in ESS and was wondering if it is possible to change the
color of the comments depending on the level.
E.g. # Green
Red
### Blue
Cheers
Renger
Modelworks
Gewerbestrasse 15
3600 Thun - Switzerland
+41 79 818 53 73
Dear R-Users,
I have a newbie-question about producing several plots with five variable
pairs within in one code (loop).
Problem:
I define two objects, erveryone has five vectors (columns /variables):
dimensionen - data.frame(meinspss$attr_diff_gesamt,
meinspss$finanz_diff_gesamt,
On Apr 8, 2012, at 12:56 AM, Berend Hasselman wrote:
On 08-04-2012, at 03:31, nkschmidt wrote:
Hi Maria, you need to name your history file
somefilename.Rapp.history,
example
savehistory(file = history.Rapp.history)
/nkschmidt
You replied to a message almost two years old, which to my
Sorry, I didn't have time to check the speed, indeed.
However - isn't apply the same as a loop, just hidden?
D.
On Fri, Apr 6, 2012 at 6:59 PM, ilai ke...@math.montana.edu wrote:
On Fri, Apr 6, 2012 at 4:02 PM, Dimitri Liakhovitski
dimitri.liakhovit...@gmail.com wrote:
This works great:
Thank you very much, Rui.
This definitely produces the result needed.
Again, I have not checked the speed yet.
input-as.matrix(data.frame(a=c(5,1,3,7),b=c(2,6,4,8)))
f - function(x){
nr - nrow(x)
result - matrix(0, nrow=nr, ncol=ncol(x))
colnames(result) - colnames(x)
thanks for the heads up, I am checking all my providers EULAs to see
which one I can share online... I imagine that as long as it is
delayed 24 hours, there might be no issues... we will see... either
way, worst case, I lost $9 for a registration... best case... I can
share most ETF's with the
Please do your homework first: Read at least the beginning of An
Introduction to R or other R tutorial -- there are many. I say this
because in your post below you do not seem to realize that R is case
sensitive, and you do not seem to know the difference between a data
frame and a column thereof
Preoccupation with speed of execution is typically (certainly not
always) misplaced. Provided you have used sensible basic
vectorization, loops in whatever form work adequately. First get
working code. Then, if necessary, parallelization, byte compilation,
or complex vectorization strategies can
On Sat, Apr 7, 2012 at 8:29 PM, wcheckle wchec...@jhsph.edu wrote:
Thank you David, the bwplot option does what I need:
snip
However, I am interested in also learning how to do it in xyplot as well. I
wasn’t able to follow the last two set of instructions
That was me. Sorry for any
On 02.04.2012 21:35, Jonathan Greenberg wrote:
R-helpers:
I'm curious what support R has for parallel writes to a binary file?
If I want to use snow to have each node write different rows of a
flat binary file (possibly out of sequence), are there any
tricks/issues I should be aware of?
On Apr 8, 2012, at 8:51 AM, slyrs66 wrote:
Dear R-Users,
I have a newbie-question about producing several plots with five
variable
pairs within in one code (loop).
Problem:
I define two objects, erveryone has five vectors (columns /variables):
dimensionen -
On 01.04.2012 20:15, Bazman76 wrote:
Hi there,
When I run the code below I get the error
Error in xy.coords(x, NULL, log = log) :(list) object cannot be coerced to
type 'double'
Any tips how I can resolve this?
The mra output is a list, hence plot.ts cannot handle it.
You probably want to
Export the default method (and otehr methods) from your NAMESPACE? See
Writing R Extensions about the S3method directive.
Note that print.f is the print method for a class f, while f is a
generic function, both fs are unrelated.
Uwe Ligges
On 02.04.2012 12:02, Johannes Radinger wrote:
This sounds neat, and I wish you the best of luck with it. While I'm
sure it'll be great for folks who are already curious enough to dive
into Sweave or knitr, my guess is that coaxing collaborators into using
markdown with SVN or GIT instead of passing around a Word file with
Track Changes is
On Sun, Apr 8, 2012 at 3:42 AM, Kaveh Vakili kaveh.vak...@ulb.ac.be wrote:
Hi all,
I have this problem with lattice that xyplot() won't draw some of my axis
labels if the type (i.e. the relation argument) of scales is set as free. For
example, in the plot below, I would want it to also
Hello -
I want to generate stacked plots with par(mfrow)) function. However,
my axis labels aren't showing.
My script is here:
http://pastebin.com/yXXeMQgb
The plot is here:
http://www.crypticlineage.net/rdisc/strplot.pdf
Thank you for your time.
Vikram
Dear all,
I have been using the graph package until last week, especially the functions
independence.number and largest.independent.vertex.sets. This weekend
I formatted my hard disk and when trying to install packages I
noticied that graph package is no longer on R repository although can
On Sun, Apr 8, 2012 at 11:40 AM, Vikram Chhatre
crypticline...@gmail.com wrote:
Hello -
I want to generate stacked plots with par(mfrow)) function. However,
my axis labels aren't showing.
Your mar (2) are too narrow. You could increase back to the default or
use the lines option in mtext to
Indeed. Changing margins brought back the axis labels. Thank you.
What I would really like to do is stack these plots atop each other so
they share the common x-axis. Is there a way to do it with ggplot2?
Vikram
On Sun, Apr 8, 2012 at 1:33 PM, ilai ke...@math.montana.edu wrote:
On Sun,
First:
sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8LC_COLLATE=en_US.UTF-8
[5] LC_MONETARY=en_US.UTF-8LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=C LC_NAME=C
On Sun, Apr 8, 2012 at 8:26 AM, Dimitri Liakhovitski
dimitri.liakhovit...@gmail.com wrote:
Sorry, I didn't have time to check the speed, indeed.
However - isn't apply the same as a loop, just hidden?
D.
Yes ?apply is a loop but not the same as ?for, see Intro to R. As
Bert Gunter pointed out
On Sun, Apr 8, 2012 at 2:56 PM, Kjetil Halvorsen
kjetilbrinchmannhalvor...@gmail.com wrote:
First:
sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
[3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
Agreed, Rui provided a very elegant vectorized solution - and I am
very thankful to him. Unfortunately (for myself), I am not as
proficient in vectorization - otherwise, I would not have asked the
question.
Why I did not follow on the original hint to use apply? For this reason (quote):
However -
that's it!
Thank you.
--
View this message in context:
http://r.789695.n4.nabble.com/Error-in-xy-coords-x-NULL-log-log-list-object-cannot-be-coerced-to-type-double-tp4524046p4541340.html
Sent from the R help mailing list archive at Nabble.com.
__
It's very simple my boy!! Do you already to play with mar? So..Try to
change the values this object. For example, par(..., mar=c(*5*,2,2,2)).
Bye!
--
View this message in context:
http://r.789695.n4.nabble.com/axis-labels-not-showing-tp4541268p4541354.html
Sent from the R help mailing list
Hi there,
Can someone explain what the difference between spec.pgram and spec.ar is?
I understand that they attempt to do the same thing one using an AR
estimation of the underlying series to estimate teh sensity the other using
the FFT. However when applied to teh same data set they seem to be
Thank you, David.
What I actually would like to do in an more elegant way is the following
boxplot (meinspss$attr_diff_gesamt ~ meinspss$R1_02)
boxplot (meinspss$finanz_diff_gesamt ~ meinspss$R2_02)
boxplot (meinspss$leist_diff_gesamt ~ meinspss$R3_02)
boxplot (meinspss$soz_diff_gesamt ~
Hi, i am writing a function to plot a pdf functions of a Generalized normal
laplace distribution.
The code is this
{
y = x-rho*mu
cf.fn = function(s){
cplex = complex(1,0,1)
temp1 = alpha*beta*exp(-sigma*s^2/2)
temp2 =
Hi all,
I wanted to make sure that I am getting the most from my 8 core computer
running R. I have a distance function like this:
distance - function(pointX, pointY, vecX, vecY)
{
#convert everything to radians
pointX = pointX * pi/180
pointY = pointY * pi/180
vecX = vecX * pi/180
vecY
On Sun, Apr 8, 2012 at 11:57 AM, Duncan Murdoch
murdoch.dun...@gmail.com wrote:
On 12-04-07 4:51 PM, Henrik Bengtsson wrote:
On Sat, Apr 7, 2012 at 1:30 PM, Mark Heckmannmark.heckm...@gmx.de
wrote:
Hello,
using the- assignment operator I do not understand why the following
does not work.
CRAN (and crantastic) updates this week
New packages
* bigdata (0.1)
Maintainer: Han Liu
Author(s): Han Liu, Tuo Zhao
License: GPL-2
http://crantastic.org/packages/bigdata
The big data package is a collection of scalable methods for
large-scale data analysis.
* boss
Michael,
I have not had time to look at this for a while but still wanted to say
thanks for looking into it and sending this solution.
By the way, Jeff mentioned that the version of quantmod on the SVN
(0.3.18) works for this. I tried to figure out how to download that
version, but found the
This can be simplified by using the layering abilities that Felix Andrews
made available in latticeExtra. These are too little known. These pretty
much make it unnecessary to resort to trellis.focus(), at least in such
cases as this. These layering abilities are too little known:
Hi,
I am not sure I follow the scalar part of the suggestion. Could you please
elaborate it a bit more? Each row has a different value of score for each
subject at each time point. Also I simplified the example but there are
other terms that change over each row for each subject.
Does this mean
On Apr 8, 2012, at 7:24 PM, John Maindonald wrote:
This can be simplified by using the layering abilities that Felix
Andrews
made available in latticeExtra. These are too little known. These
pretty
much make it unnecessary to resort to trellis.focus(), at least in
such
cases as
i appreciate all the interest in my question, and thank you Elai for both of
your suggestions, which work very well.
Elai, your last code was particularly simple and helpful to generate the
figure i was looking for.
x11(height=8,width=11)
par(lend=2)
xdat = data.frame(mortality =c(5,
First I'd like to thank Ilaik and Paul for their responses.
I seem to be getting what I need by using LaTeX's tikz package and
including tikz=true in the options for the figure chunks using vcd,
e.g.,
fig=true,tikz=true,width=5,height=5,echo=FALSE,results=hide=
... R commands ...
@
I have
Hello,
I am doing some text mining and text analysis in R. However, I have NO idea
how to install hunspell into R and use the aspell() function correctly,
I keep running into this error:
Error in aspell(f, Rd) : No suitable spell check program found.
Come someone give me a detailed step by
Hi, I would really appreciate all the help I can get. Unfortunately, I am
really new to statistics! I hope you guys don't mind this.
I am trying to find significance levels, beta, R, R squared, adjusted R
squared, standard error and t test.
FILE
Good day fellow R users:
I have routinely received the following message when attempting to
estimate a GMM model for a somewhat square panel (N = 20, T = 9-27,
Obs = 338) using the pgmm function in the plm package:
Error in function (..., deparse.level = 1) :
number of rows of matrices must
On Apr 8, 2012, at 1:28 PM, Guaramy wrote:
Hi, i am writing a function to plot a pdf functions of a Generalized
normal
laplace distribution.
The code is this
{
y = x-rho*mu
cf.fn = function(s){
cplex = complex(1,0,1)
temp1 =
Hi, all:
Can anybody check where is wrong with my code? I tried a lot of times, but
did not find an error. The parameters' estimator is not accurate. It's a
simple model about a multiple regression, with five covariates. rwmetrop is
supposed to give a much more accurate estimand. Thanks a lot.
Dear R Help,
I am attempting to write a function that takes a list of variable groups and a
vector of numbers (e.g., jtsv would indicate one group of variables, jtsv1,
jtsv2, etc, each of which name a variable in a data frame), and returns a
list of variable group sums. For example, given
?get
?rowSums
?rowsum
a group is not a defined data structure in R, btw. You have a much
better chance of getting precise answers if you ask precise questions.
See the Posting Guide.
-- Bert
On Sun, Apr 8, 2012 at 9:01 PM, Frank Tamborello frank.ta...@gmail.com wrote:
Dear R Help,
I am
On 09-04-2012, at 01:26, Navin Goyal wrote:
Hi,
I am not sure I follow the scalar part of the suggestion. Could you please
elaborate it a bit more? Each row has a different value of score for each
subject at each time point.
That last remark has nothing to do with it.
The point is that
51 matches
Mail list logo