stepHeat {astrochron} | R Documentation |

The stepHeat function will evaluate data from stepwise heating experiments, producing an Ar/Ar age spectrum, a weighted mean age with uncertainty, and other helpful statistics/plots (with interactive graphics for data culling). The function includes the option to generate results using the approach of IsoPlot 3.70 (Ludwig, 2008) or ArArCALC (Koppers, 2002).

```
stepHeat(dat,unc=1,lambda=5.463e-10,J=NULL,Jsd=NULL,CI=2,cull=-1,del=NULL,output=F,
idPts=T,size=NULL,unit=1,setAr=95,color="black",genplot=T,verbose=T)
```

`dat` |
dat must be a data frame with seven columns, as follows: (1) %Ar39 released, (2) date, (3) date uncertainty (one or two sigma), (4) K/Ca, (5) %Ar40*, (6) F, and (7) F uncertainty (one or two sigma). NOTE: F is the ratio Ar40*/Ar39K (see Koppers, 2002). |

`unc` |
What is the uncertainty on your input dates? (1) one sigma, or (2) two sigma. DEFAULT is one sigma. This also applies to the F uncertainty, and the J-value uncertainty (if specified) |

`lambda` |
Total decay constant of K40, in units of 1/year. The default value is 5.463e-10/year (Min et al., 2000). |

`J` |
Neutron fluence parameter |

`Jsd` |
Uncertainty for J-value (neutron fluence parameter; as one or two sigma) |

`CI` |
Which convention would you like to use for the 95% confidence intervals? (1) ISOPLOT (Ludwig, 2008), (2) ArArCALC (Koppers, 2002) |

`cull` |
Would you like select dates with a graphical interface? (0=no, 1=select points to retain, -1=select points to remove) |

`del` |
A vector of indices indicating dates to remove from weighted mean calculation. If specified, this takes precedence over cull. |

`output` |
Return weighted mean results as new data frame? (T or F) |

`idPts` |
Identify datum number on each point? (T or F) |

`size` |
Multiplicative factor to increase or decrease size of symbols and fonts. The default is 1.4 |

`unit` |
The time unit for your results. (1) = Ma, (2) = Ka |

`setAr` |
Set the %Ar40* level to be illustrated on the plot. The default is 95%. |

`color` |
Color to use for symbols. Default is black. |

`genplot` |
Generate summary plots? (T or F) |

`verbose` |
Verbose output? (T or F) |

This function performs weighted mean age calculations for step-heating data, including estimation of age uncertainties, mean square weighted deviation, and probability of fit.

The following plots are produced:

(1) %Ar40* versus %Ar39 released

(2) K/Ca versus %Ar39 released

(3) Ar/Ar age spectrum, with 2 sigma uncertainties for each step, and weighted mean with 95% confidence interval (in red)

If the J-value and its uncertainty are input, stepHeat will calculate and include the uncertainty associated with J. The uncertainty is calculated and propagated following equation 18 of Koppers (2002).

A NOTE regarding confidence intervals: There are two conventions that can be used to calculate the confidence intervals, selected with the option 'CI':

(1) ISOPLOT convention (Ludwig, 2008). When the probability of fit is >= 0.15, the confidence interval is based on 1.96*sigma. When the probability of fit is < 0.15, the confidence interval is based on t*sigma*sqrt(MSWD).

(2) ArArCALC convention (Koppers, 2002). When MSWD <=1, the confidence interval is based on 1.96*sigma. When MSWD > 1, the confidence interval is based on 1.96*sigma*sqrt(MSWD).

ADDITIONAL ADVICE: Use the function readMatrix to load your data in R (rather than the function read).

A.A.P. Koppers, 2002,
*ArArCALC- software for 40Ar/39Ar age calculations*: Computers & Geosciences, v. 28, p. 605-619.

K.R. Ludwig, 2008,
*User's Manual for Isoplot 3.70: A Geochronological Toolkit for Microsoft Excel*: Berkeley Geochronology Center Special Publication No. 4, Berkeley, 77 p.

I. McDougall and T.M. Harrison, 1991,
*Geochronology and Thermochronology by the 40Ar/39Ar Method*: Oxford University Press, New York, 269 pp.

K. Min, R. Mundil, P.R. Renne, and K. Ludwig, 2000,
*A test for systematic errors in 40Ar/39Ar geochronology through comparison with U/Pb analysis of a 1.1-Ga rhyolite*: Geochimica et Cosmochimica Acta, v. 64, p. 73-98.

I. Wendt and C. Carl, 1991,
*The statistical distribution of the mean squared weighted deviation*: Chemical Geology, v. 86, p. 275-285.

```
## Not run:
# Check to see if this is an interactive R session, for compliance with CRAN standards.
# YOU CAN SKIP THE FOLLOWING LINE IF YOU ARE USING AN INTERACTIVE SESSION.
if(interactive()) {
# Sample MT-09-09 incremental heating Ar/Ar data from Sageman et al. (2014).
perAr39 <- c(4.96,27.58,19.68,39.9,6.25,1.02,0.42,0.19)
age <- c(90.08,89.77,89.92,89.95,89.89,89.55,87.71,86.13)
sd <- c(0.18,0.11,0.08,0.06,0.14,0.64,1.5,3.22)
KCa <- c(113,138,101,195,307,27,17,24)
perAr40 <- c(93.42,99.42,99.64,99.79,99.61,97.99,94.64,90.35)
Fval <- c(2.148234,2.140643,2.144197,2.145006,2.143627,2.135163,2.090196,2.051682)
Fsd <- c(0.00439,0.00270,0.00192,0.00149,0.00331,0.01557,0.03664,0.07846)
ex <- data.frame(cbind(perAr39,age,sd,KCa,perAr40,Fval,Fsd))
stepHeat(ex)
# plot without points identified
stepHeat(ex,size=0,idPts=FALSE,cull=0)
}
## End(Not run)
```

[Package *astrochron* version 1.2 Index]