Demonstrating the Power of F Test with gWidgets | /en/2010/04/demonstrating-the-power-of-f-test-with-gwidgets/

Guest *Tal Galili* @ 2010-04-03 12:56:46 originally posted:

Guest *Max* @ 2010-05-04 22:54:28 originally posted:

Guest *Max* @ 2010-05-04 22:56:53 originally posted:

> library(gWidgetsRGtk2)
Loading required package: gWidgets
> options(guiToolkit = "RGtk2")
> ## if you don't want to install Gtk+, you may use the less fancy tcl/tk
> options(guiToolkit = 'tcltk')
> tbl = glayout(container = gwindow("Power of the F Test"),
+     spacing = 0)
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".glayout", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> tbl[1, 1:4, anchor = c(0, 0), expand = TRUE] = g.f = ggraphics(container = tbl,
+     expand = TRUE, ps = 11)
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".ggraphics", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> tbl[2, 1, anchor = c(1, 0)] = "numerator df"
Error in tbl[2, 1, anchor = c(1, 0)] = "numerator df" : 
  object 'tbl' not found
> tbl[2, 2, anchor = c(0, 0), expand = TRUE] = g.dfn = gslider(from = 1,
+     to = 50, value = 3, container = tbl, handler = function(h,
+         ...) {
+         p.Ftest(dfn = svalue(h$obj))
+     })
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".gslider", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> tbl[2, 3, anchor = c(1, 0)] = "denominator df"
Error in tbl[2, 3, anchor = c(1, 0)] = "denominator df" : 
  object 'tbl' not found
> tbl[2, 4, anchor = c(0, 0), expand = TRUE] = g.dfd = gslider(from = 1,
+     to = 50, value = 20, container = tbl, handler = function(h,
+         ...) {
+         p.Ftest(dfd = svalue(h$obj))
+     })
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".gslider", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> tbl[3, 1, anchor = c(1, 0)] = "delta^2"
Error in tbl[3, 1, anchor = c(1, 0)] = "delta^2" : object 'tbl' not found
> tbl[3, 2, anchor = c(0, 0), expand = TRUE] = g.ncp = gslider(from = 0,
+     to = 100, value = 10, container = tbl, handler = function(h,
+         ...) {
+         p.Ftest(ncp = svalue(h$obj))
+     })
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".gslider", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> tbl[3, 3, anchor = c(1, 0)] = "alpha"
Error in tbl[3, 3, anchor = c(1, 0)] = "alpha" : object 'tbl' not found
> tbl[3, 4, anchor = c(0, 0), expand = TRUE] = g.alpha = gslider(from = 0,
+     to = 1, by = 0.01, value = 0.05, container = tbl, handler = function(h,
+         ...) {
+         p.Ftest(alpha = svalue(h$obj))
+     })
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".gslider", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> tbl[4, 1, anchor = c(1, 0)] = "x range"
Error in tbl[4, 1, anchor = c(1, 0)] = "x range" : object 'tbl' not found
> tbl[4, 2:4, anchor = c(0, 0), expand = TRUE] = g.xr = gslider(from = 1,
+     to = 50, value = 15, container = tbl, handler = function(h,
+         ...) {
+         p.Ftest(xr = svalue(h$obj))
+     })
Loading required package: gWidgetstcltk
Error in function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ".gslider", for signature "guiWidgetsToolkittcltk"
In addition: Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE,  :
  there is no package called 'gWidgetstcltk'
> ## draw the graph
> p.Ftest = function(dfn = svalue(g.dfn), dfd = svalue(g.dfd),
+     ncp = svalue(g.ncp), alpha = svalue(g.alpha), xr = svalue(g.xr)) {
+     x = seq(0.001, xr, length.out = 300)
+     yc = df(x, dfn, dfd)
+     yn = df(x, dfn, dfd, ncp = ncp)
+     par(mar = c(4.5, 4, 1, 0.05))
+     plot(x, yc, type = "n", ylab = "Density", ylim = c(0, max(yc,
+         yn)))
+     xq = qf(1 - alpha, dfn, dfd)
+     polygon(c(xq, x[x >= xq], xr), c(0, yn[x > xq], 0), col = "gray",
+         border = NA)
+     lines(x, yc, lty = 1)
+     lines(x, yn, lty = 2)
+     legend("topright", c(as.expression(substitute(F[list(df1,
+         df2)] ～ " density", list(df1 = dfn, df2 = dfd))), as.expression(substitute(F[list(df1,
+         df2)](ncp) ～ " density", list(df1 = dfn, df2 = dfd, ncp = ncp))),
+         as.expression(substitute("Power = " ~ p, list(p = round(1 -
+             pf(xq, dfn, dfd, ncp = ncp), 4))))), lty = c(1:2,
+         NA), fill = c(NA, NA, "gray"), border = NA, bty = "n")
+ }
> p.Ftest()
Error: object 'g.xr' not found
Error in svalue(g.xr) : 
  error in evaluating the argument 'obj' in selecting a method for function 'svalue'

Guest *j verzani* @ 2010-05-10 04:40:23 originally posted:

require(traitr)                         # also uses cairoDevice, gWidgetsRGtk2, gWidgets
options(guiToolkit = "RGtk2")           # not tcltk -- using ggraphics


drawGraph <- function(.) {
  vals = xq], xr), c(0, yn[x > xq], 0), col = "gray",
            border = NA)
    lines(x, yc, lty = 1)
    lines(x, yn, lty = 2)
    legend("topright", c(as.expression(substitute(F[list(df1,
        df2)] ～ " density", list(df1 = dfn, df2 = dfd))), as.expression(substitute(F[list(df1,
            df2)](ncp) ～ " density", list(df1 = dfn, df2 = dfd, ncp = ncp))),
                         as.expression(substitute("Power = " ~ p, list(p = round(1 -
                             pf(xq, dfn, dfd, ncp = ncp), 4))))), lty = c(1:2,
                                                                    NA), fill = c(NA, NA, "gray"), border = NA, bty = "n")
    return(1 - pf(xq, dfn, dfd, ncp = ncp))
  })
}

dlg <- aDialog(items=list(
                 graph=graphicDeviceItem(),
                 dfn=rangeItem(value=3, from=1, to=50, label="numerator d.f."),
                 dfd=rangeItem(value=20, from=1, to=50, label="denominator d.f."),
                 ncp=rangeItem(value=10, from=0, to=100, label="delta^2"),
                 alpha=rangeItem(value=0.05, from=0, to=1, by=0.01, label="alpha"),
                 xr = rangeItem(value=15, from=1, to=50, label="x range")
                 ),
               title="Power of the F Test",
               model_value_changed=drawGraph,
               Ok_handler=drawGraph,
               help_string="Change the value to see the power of the F test")

dlg$make_gui()                          # can define a view to adjust layout of widgets
drawGraph(dlg)                          # initial graphic