Displaying 6 results from an estimated 6 matches for "noncent".
2004 Sep 20
2
asypow.noncent: how does it work?
I am trying to do power calculations for the proportional odds model using the asypow library.
The code
noncenta90b10<-asypow.noncent(theta.ha=a9010,info.mat=infomatrixa90b10,constraints=constrt)
returns
Error in max(..., na.rm = na.rm) : invalid "mode" of argument.
the various arguments I've used are:
a9010
[,1]
[1,] -1.7357568
[2,] -0.1928619
specifying the theta.ha array as...
2018 Apr 25
1
Can't Get Lattice Histogram Minor Tick Marks to Work
..., 1.0)
for (i in 1:20) {
maxCount = 100 * counts[i]
for (j in 1:maxCount){
pmat = rbind( pmat, c(i*0.05-0.001,2) )
}
}
# Convert the matrix to a data frame and add column names.
pDat = data.frame(pmat)
colnames(pDat) = c("pValue","nonCent")
# Create a factor to enable titles on the two panels.
nonCent.f = factor(pDat$nonCent,
labels=c("Noncentrality = 0.0", "Noncentrality = 2.0"))
# Use the lattice histogram function to draw the histogram.
histogram( ~ pValue | nonCent.f, data = pDat, nint = 20,
end...
2018 Apr 25
0
Can't Get Lattice Histogram Minor Tick Marks to Work
Per the Posting Guide, why didn't you post the reproducible R code example?
On April 24, 2018 8:22:15 PM PDT, Donald Macnaughton <donmac at matstat.com> wrote:
>I'm drawing a paneled histogram using the lattice package. I've
>succeeded in
>adding minor tick marks to the vertical axis, but I can't get the
>desired
>number of minor tick marks between the major
2018 Apr 25
3
Can't Get Lattice Histogram Minor Tick Marks to Work
I'm drawing a paneled histogram using the lattice package. I've succeeded in
adding minor tick marks to the vertical axis, but I can't get the desired
number of minor tick marks between the major tick marks.
I've attached a self-contained program to illustrate the problem.
Thanks for your help,
Don Macnaughton
Here's my sessionInfo:
R version 3.4.3 (2017-11-30)
Platform:
2003 Aug 01
1
ncp t & Fortran error & power of some tests
Hi everybody, I have three questions to ask us:
a) R incorporates a function for the Non-central T
distribution which unfortunately and, as you know, is
not available in Splus 4.5. In
http://www.stats.ox.ac.uk/pub/Swin I found the Don
MacQueen?s noncent.zip but when I run it in Splus 4.5
the following error message appears: "Error in
.Fortran ("vectnc",: "VECTNC" is not a symbol in the
load table". May be I did not installed it correctly
or (as I suppose) it is incompatible with this version
of Splus. I looked in the...
2006 Mar 08
1
power and sample size for a GLM with Poisson response variable
...I get very different results using asypow package compared
to my other (home grown) approaches.
library(asypow)
pois.mean<-c(0.0065,0.0003)
info.pois <- info.poisson.kgroup(pois.mean, group.size=c(1081,3180))
constraints <- matrix(c(2,1,2), ncol = 3, byrow=T)
poisson.object <- asypow.noncent(pois.mean, info.pois,constraints)
power.pois <- asypow.power(poisson.object, c(1081,3180), 0.05)
print(power.pois)
[1] 0.2438309
asy.pwr() # the function is shown below.
$power
[1] 0.96
sim.pwr() # the function is shown below.
4261000 Poisson random variates simulated
$power
[1] 0.567...