pvalue1

How to download and try this example

List input data set pvalue1.dat

Usage: pl -gif pvalue1.htm


#proc page
  pagesize: 18 10

#proc annotate
  location: 1 10
  textdetails: align=left size=18
  text: Comparison of SEA/GnJ  Bmp5 genotype   se/se  vs.  +/se    (means +/- 1SEM)

#proc legendentry
  sampletype: color
  details: purple
  label: se/se

#proc legendentry
  sampletype: color
  details: green
  label: +/se

#proc legend
  location: 1.2 9.7
  format: across


#set slide = 0.1


#proc getdata
  file: pvalue1.dat


#proc categories
  axis: x
  categories:   f
		m

#set xleft = 1
#set bot = 7
#loop

  #proc processdata
    action: breaks
    fields: 5 
  #endproc 

  #if @NRECORDS = 0
    #break
  #endif

  #if @BREAKFIELD1 like "lower_jaw"
    #set bot = 4
    #set xleft = 1
  #endif
  #if @BREAKFIELD1 like "skull_width"
    #set bot = 1
    #set xleft = 1
  #endif

  #set xright = $arith( @xleft+1.5 )
  #set top = $arith( @bot+2 )

  #proc areadef
    title: @BREAKFIELD1
    titledetails: size=12
    rectangle: @xleft @bot @xright @top
    xrange: categories
    yautorange: datafield=6
    xaxis.stubs: usecategories
    xaxis.stubdetails: size=12 adjust=0,0.08
    yaxis.stubs: inc
    yaxis.stubdetails: size=6 adjust=0.1,0

  #proc categories
    axis: x
    slideamount: -@slide
 
  #proc scatterplot
    select: @@4 = (se/se)
    xfield: 2
    yfield: 6
    symbol: shape=circle style=filled color=gray(0.7) radius=0.02
    cluster: yes

  #proc categories
    axis: x
    slideamount: @slide
 
  #proc scatterplot
    select: @@4 = (+/se)
    xfield: 2
    yfield: 6
    symbol: shape=circle style=filled color=gray(0.7) radius=0.02
    cluster: yes
    
  #proc processdata
    action: summary
    fields: 2 4
    valfield: 6

  #proc categories
    axis: x
    slideamount: -@slide
 
  #proc bars
    select: @@id2 = (se/se)
    locfield: id1
    lenfield: mean
    errbarfield: se
    tails: 0.05
    thinbarline: color=purple
    // color: purple

  #proc scatterplot
    select: @@id2 = (se/se)
    xfield: id1
    yfield: mean
    symbol: shape=pixcircle radius=0.03 color=purple style=filled

  #proc categories
    axis: x
    slideamount: @slide
 
  #proc bars
    select: @@id2 = (+/se)
    locfield: id1
    lenfield: mean
    errbarfield: se
    tails: 0.05
    thinbarline: color=green
    // color: green

  #proc scatterplot
    select: @@id2 = (+/se)
    xfield: id1
    yfield: mean
    symbol: shape=pixcircle radius=0.03 color=green style=filled

  #proc processdata
    action: join
    field: id1
    // should be able to use @@id2 in these select statements but can't due to bug....
    leftselect: @@2 = (se/se)
    rightselect: @@2 = (+/se)
    fieldnames: Lid1 Lid2 Lmean Lsd Lse Ln_obs Lmin Lmax Lsum Rid1 Rid2 Rmean Rsd Rse Rn_obs Rmin Rmax Rsum 

  #proc categories
    axis: x
    slideamount: 0
 
  // print p values
  #proc pvalue
    xfield: Lid1
    statfields1: Lmean Lsd Ln_obs
    statfields2: Rmean Rsd Rn_obs
    textdetails: size=8 adjust=0,-0.4

    
  #set xleft = $arith( @xleft+2 )

  #proc usedata
   original: yes

#endloop


#exit