Displaying 5 results from an estimated 5 matches for "adjp".
Did you mean:
adj
2008 Dec 14
1
how to convert factors to numbers
...of the data as
"factors" and I do not know how to treat the data as numbers:
> spy<-read.table("c:\\StockData\\SPY.txt")
> attach(spy)
> names(spy)<-c("QDate","OpenP","HighP","LowP","CloseP","Vol","AdjP")
> spy[1,]
QDate OpenP HighP LowP CloseP Vol AdjP
1 12/14/2006, 141.86, 143.24, 141.84, 143.12, 64755200, 138.34
> ChangeFromOpen<-spy[1,5]-spy[1,2]
Warning message:
In Ops.factor(spy[1, 5], spy[1, 2]) : - not meaningful for factors
>
As you can see, I c...
1999 Nov 25
0
desperate!
...---------------
Code for Bonferroni, Holm and Hochberg post hoc tests:
p.adjust.holm <- function(p, n=length(p)) {
##n <- length(p)
r <- rank(p)
index <- order(p)
qi <- p*(n+1-r)
for (i in 2:length(p)) {
qi[index[i]] <- max(qi[index[i+1]])
}
list(adjp=pmin(qi, 1), p=p, method="Holm")
}
p.adjust.hochberg <- function(p) {
n <- length(p)
r <- rank(p)
index <- order(p)
qi <- p*(n+1-r)
for (i in (n-1):1) {
qi[index[i]] <- min(qi[index[i]], qi[index[i+1]])
}
list(adjp=qi, p=p, method="...
1998 Jul 14
1
Are post-hoc tests being developed for R?
Hi-
Is anyone working on multiple comparisons of means or post-hoc tests
(ie: Tukey, Bonferroni) for R? I saw in the winter 98' archives of the R
mailing lists that these tests had not been implemented yet, I was just
wondering if I could look foward to having them. ;-)
I also looked through the contributed packages and didn't see anything
that offered such tests.
I guess I could check
1998 Jul 16
1
R-beta: Re: Post-hoc tests
Matt,
Here's a Bonferroni-corrected multiple one-sample t-test that I wrote
some years ago. It took a while to get it into R, as na.omit doesn't
seem to handle vectors and I had to write a quick kludge (na.remove).
Another more general point was that I discovered that the help page for
t.test gives the name "parameters" for the degrees of freedom, as in S.
However, the name
2010 Mar 03
1
Help with multtest (rawp2adjp)
Hello R experts,
I am trying to analyze this dataset and am stuck on this problem for quite
some time now.
I am using mt.rawp2adjp.
the output that came out was a matrix with two colums since I had asked it
to calculate the adjusted p values using one method.
so it has the two columns as: rawp BH
I combined these using cbind with my actual dataframe.
checked using head all was fine.
thereafter I am trying to extract the rows...