[R] Help with debugging
Phillip Good
pigood at verizon.net
Thu Apr 22 23:52:44 CEST 2004
Once the IF triggers in the following program, it just won't turn off:
#Get p-value from permutation distribution
N = 64
cnt = 0
for (i in 1:N){
pdata = sample (data)
statp=F1(size,pdata,gmean,samps)
if (stat0 <= statp ){
cnt=cnt+1
print(i)
print (statb)
}
}
cnt/N
#print statements are soley for debugging purposes
#Here are the function code and test data needed to run the program:
F1=function(size,data, gmean, samps){
stat=0
start=0
end=0
for (i in 1:samps){
end=end+size[i]
term=mean(data[start:end])
stat= stat +abs(term-gmean)
start=end+1
}
list(stat=stat)
}
samps = 4
size = c(4,4,3,5)
data = rnorm(16, 2*rbinom(16,1, .4))
data = data +c(2,2,2,2, 0,0,0,0, 0,0,0, 0,0,0,0,0)
gmean = mean(data)
stat0=F1(size,data,gmean,samps)
More information about the R-help
mailing list