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 diasan...@gmail.com 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
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
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
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
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
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
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 dwinsem...@comcast.net
To: sovo0...@gmail.com
Cc: r-help@r-project.org
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 frodo.j...@yahoo.com wrote:
Dear all,
how can I perform a post hoc analysis for ANOVA with repeated measures (in
presence of a
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
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, line
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 cass...@hotmail.com wrote:
This is a generalize question: basically, say you are typing a
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
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
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
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
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
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 gunter.ber...@gene.com wrote:
If I
On Sun, Jan 9, 2011 at 8:13 PM, Jim Burke j.bu...@earthlink.net 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
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(
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 djmu...@gmail.com:
Hi:
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:
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
from
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 to
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
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
