-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathACFL-NoInteractions.Rnw
97 lines (88 loc) · 2.69 KB
/
ACFL-NoInteractions.Rnw
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
\documentclass{article}
\usepackage{pdflscape}
\usepackage{longtable}
\title{Follow-up analysis ACFL nest success along trails}
\author{Max A Henschell}
\begin{document}
% \SweaveOpts{concordance=TRUE}
% % \SweaveOpts{concordance=TRUE}
% \SweaveOpts{width=10, height = 8}
\maketitle
<<include=FALSE, echo=FALSE, results='hide'>>=
source("ACFL_NoInteractionModels_MJP.R")
source("ACFL_NoInteractionModels_MAH.R")
lapply(packages.list, require, character.only=T)
@
%\section*{Variable Correlations}
%\begin{landscape}
\begin{figure}[!htbp]
\centering
<<echo=FALSE, results='hide'>>=
####################################
#
#Correlation for MJP
#
####################################
corrplot.mixed(cor(ACFL.MJP.nl[,c(12:26)]))
@
\caption{Correlation of environmental variables for MJP ACFL nests}
\end{figure}
\begin{center}
\end{center}
%\end{landscape}
%\begin{landscape}
\begin{figure}[!htbp]
\centering
<<echo=FALSE, results='hide'>>=
####################################
#
#Correlation MAH
#
####################################
corrplot.mixed(cor(ACFL.MAH.nl[,c(12:26)]))
@
\caption{Correlation of environmental variables for MAH ACFL nests}
\end{figure}
\begin{center}
\end{center}
%\end{landscape}
% variables:\\
% Forest1km: proportion forest cover within 1 km of nest -> all studies\\
% Road1km: km of roads within 1 km of nest -> all studies\\
% Homesteads1km : number of homesteads within 1 km of nest -> all studies\\
% RoadDist: distance to closest road -> all studies\\
% W2: trail width at 1.3 m from ground -> all studies\\
% GPH: groups per hour -> MAH,MJP\\
% TPH: trees pre hectare around nest -> MAH\\
% BA\_ha: basal area per hectare around nest -> MAH,MJP\\
% TYPE: trail type: constructed, old road -> all studies\\
% PLOT: trail name -> all studies\\
% HEIGHT: nest height -> all studies\\
% DIST: distance from nest to trail -> all studies\\
% ORIENT: Nest orientation from trail: 90 = opposite side of tree from trail\\
% OrdDate: Ordinal date from Jan 1 - 140 -> all studies\\
% STUDY: study\\
% YEAR: year of nest\\
% NEST: nest ID\\
%\section*{Analysis of data from AL, MJP, MAH}
\begin{figure}[!htbp]
<<echo=FALSE, results='hide'>>=
imageplot.bma(ACFL.MJP.BMA, color = c("blue", "red", "white"), order = "mds")
@
\caption{Variable selection for BMA from MJP data}
\end{figure}
%\section*{Analysis of data from MAH}
\begin{figure}[!htbp]
<<echo=FALSE, results='hide'>>=
imageplot.bma(ACFL.MAH.BMA, color = c("blue", "red", "white"), order = "mds")
@
\caption{Variable selection for BMA from MAH data}
\end{figure}
\begin{landscape}
<<echo=FALSE, results= 'asis'>>=
xtable(MJP.BMA.df[,c(1:3)], caption = "BMA MJP data")
xtable(MAH.BMA.df[,c(1:3)], caption = "BMA MAH data")
#print(MJP.xtab)
@
\end{landscape}
\end{document}