Hi,
embed() seemed well-suited, but I couldn't figure out an elegant way to use it
embed(c(A,A), 4)[1:4, 4:1]
HTH,
baptiste
On 6 January 2011 22:34, ADias wrote:
>
> Hi
>
> Suppose we have an object with strings:
>
> A<-c("a","b","c","d")
>
> Now I do:
>
> B<-matrix(A,4,4, byrow=F)
>
> and I
Hi:
This is an abridged version of the reply I sent privately to the OP.
Generate an artificial data frame
# function to randomly generate one of the Q* columns with length 1000
mysamp <- function() sample(c(-1, 0, 1, NA), 1000, prob = c(0.35, 0.2, 0.4,
0.05), replace = TRUE)
# use above fu
You could re-do part of your code to run with mpfr-class variables, and use
this function:
# mpfr choose(n,k)
rmpfac<-function(n,k,prec=50)
factorial(mpfr(n,prec))/factorial(mpfr(k,prec))/factorial(mpfr(n-k,prec))
Converting into and out of mpfr may not be worth it, but you wi
Dear all,
how can I perform a post hoc analysis for ANOVA with repeated measures (in
presence of a balanced design)?
I am not able to find a good example over internet in R...is there among you
someone so kind to give
me an hint with a R example please?
For example, the aov result of my analysis
Dear everyone:
I am doing a research on several stock markets. And I need to
construct an Bartlett HAC covariance matrix estimator for Sigma(Cov(Y0,Yj)),
j is from 0 to T. Can you tell me how to do it.
Your Sincerely!
Nigel Gregory
01/09/11
[[alternati
My suggestion is that you take the weekend lull in readership to now
read the Posting Guide which it appears you nave not yet done. Had you
done so you should have encountered the advice about how to
successfully attach data files. (None came through.) You would also
have read the request t
Make sure your data is a matrix. There are many examples of expression heatmaps
available on the bioconductor list. After checking out these examples, I would
post to the bioconductor list if you are still having problems. Also
consider a small
example to get you a working heatpmap. You have to ins
Hello, thanks for all your replies, it was a helpful lesson for me
(and hopefully for my colleagues, too). Bests, Sören
On 11-01-07 11:23, David Winsemius wrote:
Date: Fri, 7 Jan 2011 11:23:04 -0500
From: David Winsemius
To: sovo0...@gmail.com
Cc: r-help@r-project.org
Subject: Re: [R] Differe
Please read the maiz example in HH
install.packages("HH")
library(HH)
?MMC
Read all the way to the end of the maiz example.
On Sun, Jan 9, 2011 at 7:10 AM, Frodo Jedi wrote:
> Dear all,
> how can I perform a post hoc analysis for ANOVA with repeated measures (in
> presence of a balanced design
Good day R-listers,
My question is more a statistical question than an R related question,
so please bear with me
i'm currently applying the chow test of poolability
in fact i'm working with panel N=17 T=5 , and my model looks like
this : Yit= a0+B1X1+B2X2+B3X3+B4X4+eit
My question is the foll
Hello everyone,
I have a simple histogram of gasoline prices going back a few years that
I want to insert a line graph of consumer price index (cpi) over the
histogram. I have looked through the "Lattice" book by Deepayan Sarkar
but don't see anything there. How might this be done? An example
Have 40,000 rows of buy/sell trade data and am trying to add up the buys for
each second, the code works but it is very slow. Any suggestions how to
improve the sapply function ?
secEP = endpoints(xSym$Direction, "secs") # vector of last second on an XTS
timeseries object with multiple entries
Ive been trying to help someone with using R scripts with audio on Windows
Vista/7. The problem is that all the packages Ive tried fail to work
properly when being used on Windows 7.
I couldnt get "audio" to work at all. It kept coming up with the message
that wmm output could not happen. Worked
This is a generalize question: basically, say you are typing a long line of
command in R console, and then you realize you forgot to add something in
the beginning, is there a way to get to the beginning of the line without
pressing the "left" key on the keyboard and waiting for the cursor to get
Hi,
I have been trying to get the label under the total column - i.e. a mean value
of columns 2 to 6 - in a barplot I generate with this script:
t1 <- tapply(A, B, sum)
t1[8] <- mean(t1[2:6])
t1 <- as.table(t1)
barplot(t1, ylim=c(0,3000))
mtext("Var1", side = 1, line = 3)
mtext("Var2", side = 2,
try "home" or better yet use a text editor and cut/paste your command; let you
type a lot of shorter lines that are more readable and easy to change.
Sent from my iPad
On Jan 9, 2011, at 16:00, casshyr wrote:
>
> This is a generalize question: basically, say you are typing a long line of
> co
On Jan 9, 2011, at 8:13 PM, Jim Burke wrote:
Hello everyone,
I have a simple histogram of gasoline prices going back a few years
that I want to insert a line graph of consumer price index (cpi)
over the histogram. I have looked through the "Lattice" book by
Deepayan Sarkar but don't see
On Jan 9, 2011, at 4:00 PM, casshyr wrote:
This is a generalize question: basically, say you are typing a long
line of
command in R console, and then you realize you forgot to add
something in
the beginning, is there a way to get to the beginning of the line
without
pressing the "left" k
On 01/09/2011 04:00 PM, casshyr wrote:
This is a generalize question: basically, say you are typing a long line of
command in R console, and then you realize you forgot to add something in
the beginning, is there a way to get to the beginning of the line without
pressing the "left" key on the key
Try looking at the barplot function and notice that is should be
returning values for the mid-points. You should use those instead of
the at=values
On Jan 9, 2011, at 12:03 PM, Luca Meyer wrote:
Hi,
I have been trying to get the label under the total column - i.e. a
mean value of column
split the data by truncating the time to a second, then process each group.
this will save the subsetting you are doing. also merge the data with direction
and size in the same frame. it looks like you can subset by "buy" to begin
with.
Sent from my iPad
On Jan 9, 2011, at 19:10, rivercode w
Hi Jim,
Some example data would help us. I typically think of a histogram as
the frequency of values falling within a certain range (determined by
bins). Since they are univariate plots, I'm not sure how you are
planning on adding a line graph to that. If you just want bars of the
average gasol
Just to add to the silly solutions, here's how I would have done it...
mu <- 40
sdev <- 10
days <- 100:120 # range to explore
p <- 0.8
days[ match(TRUE, qnorm(0.2, mu*days, sqrt(sdev * sdev * days)) >= 4000) ]
Michael
On 9 January 2011 08:48, Bert Gunter wrote:
> If I understand what you have
On Sun, Jan 9, 2011 at 8:13 PM, Jim Burke wrote:
> Hello everyone,
>
> I have a simple histogram of gasoline prices going back a few years that I
> want to insert a line graph of consumer price index (cpi) over the
> histogram. I have looked through the "Lattice" book by Deepayan Sarkar but
> don
Dear list members,
I extensively searched in the previous threads of this mailing list an example
easy to understand for me,
and able to fit my problem, but I didn´t succed to find a solution for which I
feel certain.
At the moment I am stuck with this solution for my unbalanced design, and I
Thanks Josh, Gabor, and David,
I appreciate your suggestions and the time you took to think about this.
This was all most helpful. Gabor I will look at the zoo package soon.
Sounds interesting.
Below is what worked for me from Josh to overlay a line graph on a
histogram.
obj1 <- xyplot( r
Dear Dennis,
Thank you so very much for your quick reply. What an introduction to
R-help!! Especially I appreciated the time you put to explain the code
privately.
After a few hick-ups I got it working on my data as well.
Regards,
- Kari
Quoting Dennis Murphy :
Hi:
This is an abridged v
Dear all,
Thanks for your input. It works fine now. All it took was for me to tidy up
the workspace. Simple solution, that I should have considered earlier.
Regards,
Andreas
--
View this message in context:
http://r.789695.n4.nabble.com/Problems-with-glht-function-for-lme-object-tp3179128p32
Altogether I got five more or less silly solutions (not my judgment!), some of
them further discussed in private mail, for a problem where my expectation was
to get a simple one-liner back: "Check ?clt" or so...
Fortunately, with all of them I seem to arrive at a result that is consistent
with
Thomas Jensen-6 wrote:
>
> ... data set ...
> The data set has in total 27 countries for the years 1999 to 2008, but
> with unbalanced panels.
>
> I want to be able to estimate a model and do forecasting for each
> country in the data set. I have been looking into the YourCast package
> f
I have started exploring potential of R in applying IRT to a dataset. I have
a data of about 20k students who took a Maths test, a diagnostic in nature.
I find that I don't get stable solution while using tpm() even after passing
an argument "start.val = RANDOM". What should be done in this case t
I am trying to automate scoring done by applying 3 parameter model of IRT to
the response data. I call R from PHP to do this using exec(). Since I face
convergence issues while using R, it is important that the program captures
the error messages like "the solution is not stable" given by R while
r
Hi,
I have a fairly simple linear regression using the lm function. There
are about 100 variables and 30,000 rows of data. It runs fine and
produces a decent looking R2 value. I'm interested in performing a
stepwise variable selection to see if things can be cleaned up a bit.
Calling the step
33 matches
Mail list logo
