# Compute tracheidograms

### Description

Function to compute tracheidograms from microscopic measurements.

### Usage

1 2 3 4 5 6 |

### Arguments

`traq` |
microscopic light profile. A vector with the light measurements or a |

`val50` |
"y-value" at wich measurements should be made. |

`mw` |
integer. Width of the rolling window to smooth the original data. |

`plotit` |
logical. should the original data and the computed tracheidogram be drawn? |

`xlim` |
extent of the x-axis. A vector of length 2. |

`ylim` |
extent of the y-axis. A vector of length 2. |

`x` |
an object of class tgram, resulting from |

`colores` |
a vector of length 2, with the colors to draw the lumen and wall measurements, respectively. |

`leyenda` |
a vector of length 2 with the legend to appear in the plot. By default leyenda =c("lumen","double wall"). |

`lwd` |
width of the lines in the legend. |

`add` |
logical. If TRUE, add to a current plot. |

`traq.0` |
logical. If TRUE, draw the original measurements. |

`bg.legend` |
background color for the legend. |

`...` |
additional graphical parameters. |

### Details

The purpouse of this fucntion is obtaining cell anatomical data from microscopic light measurements (see DeSoto et al. for details of data adquisition). The microscopic lectures are first smoothed with a rolling window (using function `rollmean`

of package `zoo`

, using the selected width `mw`

. Then, the smoothed curve is "cut" at the treshold value `val50`

and the distances among the intersection points are computed. This provides an ordered sequence of lumen diameters (LD) and double wall thikness (DWT) measurements. From this sequence some other anatomical measurements are computed. Radial cell wall thickness is computed as * CWT[t] = 1/2 * (DWT[t]/2 + DWT[t+1]/2) *. Tracheid diameter is computed as * TD[t] = DWT[t]/2 + LD[t] + DWT[t+1]/2*.

### Value

tgram returns an object of class *tgram*, bassically a list with

`traq ` |
original data. |

`traq0 ` |
if |

`cut.points ` |
two column matrix with the coordinates of the intersection of |

`what ` |
vector indicating if the computed distances are of lumen (1) or double wall (2) . |

`distances ` |
ordered sequence of the computed distances (both of lumen and double wall). |

`LD ` |
ordered sequence of lumen diameters. |

`DWT ` |
ordered sequence of double wall thickness. |

`mw ` |
width of the rolling window employed to smooth the data. |

`CWT ` |
ordered sequence of radial cell wall thikness'. |

`TD ` |
ordered sequence of tracheid diameters. |

`LD_CWT_ratio ` |
ordered sequence of LD/CWT ratio. |

### Author(s)

Marcelino de la Cruz marcelino.delacruz@upm.es and Lucia DeSoto

### References

DeSoto, L., De la Cruz, M. & Fonti, P. 2010. Intra-annual pattern of tracheid size in the Mediterranean Juniperus thurifera as indicator for seasonal water stress. *Submited*

### Examples

1 2 3 4 5 6 7 8 9 | ```
## Not run:
data(traq.profile)
plot(tgram(traq.profile, mw=10),leyenda=c("lumen","double wall"),
xlab="distance pixel", ylab="grey value",ylim=c(0,250),
bg.legend="white")
## End(Not run)
``` |