Nothing
PANDabove <-
function(behavior,phaseX,v1,v2){
t1<-table(phaseX)
tmaxA<-t1[names(t1)==v1]
startA<-match(v1,phaseX)
endA<-tmaxA+startA-1
A<-behavior[startA:endA]
minA=(min(A,na.rm=T))-.1
vx=which.min(A)
tmaxB<-t1[names(t1)==v2]
startB<-match(v2,phaseX)
endB<-tmaxB+startB-1
#tsxB<-behavior[startB:endB]
B=(behavior[startB:endB])
cdcl<-c(A,NA,B)
y<-na.omit(cdcl)
total=length(y)
iv=(1:total)
end<-which(is.na(cdcl))
iv<-insert(iv,NA,end)
maxy=which.max(cdcl)
max<-cdcl[maxy]+1
numx<-sum(!is.na(cdcl))+3
par(mfrow=c(3,3))
maxy=which.max(behavior)
max<-behavior[maxy]+1
numx<-sum(!is.na(behavior))+3
#graphics.off()
layout(rbind(1,2), heights=c(4,1))
plot(iv,cdcl, ylim=c(0,max),lwd=2,type="o",col="red",bty="l", xlab="time", ylab="behavior", main="PAND" )
#yo<-readline("Y ordinate for your reference line " )
yo=max(A)-.1
abline(h=yo,col="gray",lwd=3)
#***********************above baseline
#nbaboveline<-B<as.numeric(yo)
#nump=sum(nbaboveline)
#nx= (length(B)+length(A))-nump
# p=nx/(length(B)+length(A))
#maxy=which.max(cdcl)
#max<-cdcl[maxy]+1
#numx<-sum(!is.na(cdcl))+3
#par(mfrow=c(3,3))
#maxy=which.max(behavior)
#max<-behavior[maxy]+1
#numx<-sum(!is.na(behavior))+3
#PAND<-c("PAND Above = ",round(p,2),"Minimum overlapping to remove",nump)
#numphalf = nump/2
# a=numphalf
#b=length(B)-numphalf
#c=length(A)-numphalf
#d=numphalf
#rpandphi=phi(c(a,b,c,d))
#print(PAND)
pand=PAND(A_data = A, B_data = B,improvement="increase")
print(pand)
writeLines("-------------------------------------------")
writeLines(".90 or above = very effective" )
writeLines(".70 to .89 = moderate effectiveness")
writeLines(".50 to .69 = debatable effectiveness")
writeLines(" below .50 = not effective")
writeLines("--------------------------------------------")
#phir<-c("Robust Phi = ",abs(rpandphi))
#print(phir)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.