Displaying 2 results from an estimated 2 matches for "callrate".
Did you mean:
calldate
2007 Nov 26
3
determine the number of digits printed into graph with "text(...)"
...ot;text(...)" - see code
below: "fraction" is my number. I wish to print only two digits of the
number and struggle to do this... options(digit=3) does not work...
Thank you for your suggestions! Cheers, Georg.
**********************************
Georg Ehret
Johns Hopkins
Baltimore
callrates<-read.table("RSG",col.names=c("variant","callrate"))
hist(callrates[,2],breaks=100)
fraction<-nrow(callrates[callrates$callrate>=0.9,])/nrow(callrates)
hist(callrates[,2],breaks=100)
text(0.4,500,"proportion of callrates >= 0.9:")
text(0.4,450,fr...
2009 Mar 23
2
error statement: missing value where TRUE/FALSE needed
...at the error message
means?
The specific line in the code is:
if (r<vQ[i]){ vS[i]=0 } else vS[i]=1
The code is available at:
http://www.michael-curran.com/gibbs.html
and the two data sets are available at:
http://www.geocities.jp/atsmatsumoto/ci.txt
and
http://www.geocities.jp/atsmatsumoto/callrate.txt
Note: I saved the first data set as ci.txt and the second as boj.txt and so
these are the file names that the code loads. If anyone uses the code and
the data and manages to get it to work, I would gratefully appreciate if
they could tell me what adjustments they made to the code in order t...