Changes between Version 115 and Version 116 of CookBook join BDMAP_CCM v2

Timestamp:

Nov 13, 2010 9:38:54 PM (14 years ago)

Author:

cedric

Comment:

--

CookBook join BDMAP_CCM v2

@@ -350 +350 @@
 ratiosource<-ccm_bd$dist_source_ccm/ccm_bd$dist_source
 ccm_bd$pbdiffsource<-!((diffsource>-20&diffsource< 20)|(ratiosource>0.7&ratiosource<1.3))
-sum(ccm_bd$pbdiffsource,na.rm=TRUE)/nrow(ccm_bd) # 23 % 3281
-sum(ccm_bd$pb_ratio_dist_sea,na.rm=TRUE)/nrow(ccm_bd) # 13 0.00094
-sum(ccm_bd$pbstrahler,na.rm=TRUE)/nrow(ccm_bd) #5% 756
+sum(ccm_bd$pbdiffsource,na.rm=TRUE)/nrow(ccm_bd) # 14 % 2027
+sum(ccm_bd$pb_ratio_dist_sea,na.rm=TRUE)/nrow(ccm_bd) # 12 0.00087
+sum(ccm_bd$pbstrahler,na.rm=TRUE)/nrow(ccm_bd) #3% 457
 #sum(!ccm_bd$pbdiffsource&!ccm_bd$pbstrahler&!ccm_bd$pb_ratio_dist_sea,na.rm=TRUE)/nrow(ccm_bd) # 64 % si on garde strahler
 # finally I don't trust there is a strahler rank pb
-sum(!ccm_bd$pbdiffsource&!ccm_bd$pb_ratio_dist_sea,na.rm=TRUE)/nrow(ccm_bd)  # on garde 67% des données
+sum(!ccm_bd$pbdiffsource&!ccm_bd$pb_ratio_dist_sea,na.rm=TRUE)/nrow(ccm_bd)  # on garde 74% des données
 # histograms
 layout(matrix(c(1,2), 1, 2, byrow = TRUE))
@@ -368 +368 @@
 g<-ggplot(ccm_bd)
 breaks=as.vector(c(1, 2, 5) %o% 10^(-1:3))
-g+geom_point(aes(x=dist_source,y=dist_source_ccm,colour=pbdiffsource,alpha=0.5,size=0.3))+
+g+geom_point(aes(x=dist_source,y=dist_source_ccm,colour=pbdiffsource),alpha=0.5,size=0.3)+
                 geom_abline(slope=1.3,colour="blue",lty=2)+
                 geom_abline(slope=0.7,colour="blue",lty=2)+