1
(1/.001132502831257078142695356738391845979614949037372593431483578708946772366930917327293318233295583238958097395243488108720271800679501698754246885617214043035107587768969422423556058890147225368063420158550396375990939977349943374858437146092865232163080407701019252548131370328425821064552661381653454133635334088335220838052095130237825594563986409966024915062287655719139297848244620611551528878822197055492638731596828992072480181200453001132502831257078142695356738391845979614949037372593431483578708946772366930917327293318233295583238958097395243488108720271800679501698754246885617214043035107587768969422423556058890147225368063420158550396375990939977349943374858437146092865232163080407701019252548131370328425821064552661381653454133635334088335220838052095130237825594563986409966024915062287655719139297848244620611551528878822197055492638731596828992072480181200453001132502831257078)

1
(1/.1132502831257078142695356738391845979614949037372593431483578708946772366930917327293318233295583238958097395243488108720271800679501698754246885617214043035107587768969422423556058890147225368063420158550396375990939977349943374858437146092865232163080407701019252548131370328425821064552661381653454133635334088335220838052095130237825594563986409966024915062287655719139297848244620611551528878822197055492638731596828992072480181200453001132502831257078142695356738391845979614949037372593431483578708946772366930917327293318233295583238958097395243488108720271800679501698754246885617214043035107587768969422423556058890147225368063420158550396375990939977349943374858437146092865232163080407701019252548131370328425821064552661381653454133635334088335220838052095130237825594563986409966024915062287655719139297848244620611551528878822197055492638731596828992072480181200453001132502831257078)

1
2
3
4
5
6
7
8
9
song<-function(n=99) {
  a<-" bottle"
  b<-" of beer on the wall, "
  c<-" of beer. "
  d<-"Take one down and pass it around: "
  #
  ## set "beer counter" vector
  counter<-c(as.character(n:1),"No more")
  # 
 ...

1
mean(1,2,3,4)

1
 

1
2
3
4
5
6
7
8
9
# Goal: To make a panel of pictures.
 
par(mfrow=c(3,2))                       # 3 rows, 2 columns.
 
# Now the next 6 pictures will be placed on these 6 regions. :-)
 
# Let me take some pains on the 1st
plot(density(runif(100)), lwd=2)
text(x=0, y=0.2, "100 uniforms")        # Showing you how to place text at will
 ...

1
2
3
4
5
6
7
rnd<-function(mp,sd)
  ((1:(sd[3]^mp*sd[4]^mp)*sd[1])%%sd[3]^mp)*sd[4]^mp+
  ((1:(         sd[4]^mp)*sd[2])%%sd[4]^mp)
n<-ceiling(log10(m<-3333));pn<-c(3L,7L,5L,2L)
r<-rnd(n,pn)
(r<-r[r<m]) # ; diff(r)#
all(table(r)==1)

1
2
3
4
5
6
7
8
9
nloadRegister(function (){Quickling.init("452359;0", 10, {"prefetch_logging":false});});
onloadRegister(function (){window.intl_locale_rewrites = {"meta":{"\/_B\/":"^(.*[.,!?\\s]|)","\/_E\/":"([.,!?\\s].*|)$"},"patterns":{"\/\u0001(.*)('|')s\u0001(?:'|')s(.*)\/":"\u0001$1$2s\u0001$3","\/_\u0001([^\u0001]*)\u0001\/e":"mb_strtolower(\"\u0001$1\u0001\")","\/_\u0001([^\u0001]*)\u0001\/":"javascript"}};});
onloadRegister(function (){DocumentTitle.set("Facebook");});
onloadRegister(function (){JSCC.get('j4e8ac50018dc329714131771').init($("contentArea"), $('rightCol'), null, $('toolbarContainer'));});
onloadRegister(function (){JSCC.get('j4e8ac50018dc329714131773').getData().onLoad(false, true);});
onloadRegister(function (){new CounterDisplay("messages_unread", $("jewelInnerInboxCount"), $("jewelOuterInboxCount"), "fbMessagesItemList", "jewelUnread", true);});
onloadRegister(function (){window.presenceRequests = new RequestsJewel().init(JSCC.get('j4e8ac50018dc329714131774'), "[fb]requests", false)});
onloadRegister(function (){JSCC.get('j4e8ac50018dc329714131774').init("requests", $("fbRequestsJewel"));});
onloadRegister(function (){new MessagesJewel().init(JSCC.get('j4e8ac50018dc329714131775'), "[fb]messages")});
 ...

1
2
3
4
5
6
7
8
9
#H0: p1 >= p2 v.s. H1: p1 < p2
#Assume sum(X) ~ Bin(N, p1), and sum(Y) ~ Bin(M, p2)
#Let p1 = sum(x)/N, p2 = sum(y)/M
N = 20000
M = 40000
p1 = 0.995
p2 = 0.997
 
B = 10000 #number of bootstrap times
 ...

1
2
3
4
5
6
7
8
9
#H0: p1 <= p2 v.s. H1: p1 > p2
#Assume sum(X) ~ Bin(N, p1), and sum(Y) ~ Bin(M, p2)
#Let p1 = sum(x)/N, p2 = sum(y)/M
N = 20000
M = 40000
p1 = 0.995
p2 = 0.997
 
B = 10000 #number of bootstrap times
 ...

1
2
3
4
5
6
7
8
9
#H0: p1 <= p2 v.s. H1: p1 > p2
#Assume sum(X) ~ Bin(N, p1), and sum(Y) ~ Bin(M, p2)
#Let p1 = sum(x)/N, p2 = sum(y)/M
N = 20000
M = 40000
p1 = 0.995
p2 = 0.997
 
B = 10000 #number of bootstrap times
 ...

1
2
3
4
5
6
7
8
9
N1 <- 20000 # 實際的樣本數量
B <- 1000 # bootstrap 的次數
p1 <- 0.995 # Under H0 的良率
 
success_times1_0.975 <- rep() #良品的個數
 
for (i in 1:B)
{
  success_times1_0.975[i] <- quantile(rbinom(B,N1,p1)/N1, 0.975)
 ...

1
2
3
4
5
6
7
8
9
 
N <- 40000 # 實際的樣本數量
B <- 1000 # bootstrap 的次數
p <- 0.995 # Under H0 的良率
 
sucess_times <- rep() #良品的個數
 
for (i in 1:B)
{
 ...

1
2
3
4
5
6
7
8
9
N <- 40000 # 實際的樣本數量
B <- 1000 # bootstrap 的次數
p <- 0.995 # Under H0 的良率
 
sucess_times <- rep() #良品的個數
 
for (i in 1:B)
{
  sucess_times [i] <- quantile(rbinom(B,N,p)/N, 0.975)
 ...

1
2
3
4
5
6
7
8
9
N <- 40000 # 實際的樣本數量
B <- 1000 # bootstrap 的次數
p <- 0.995 # Under H0 的良率
 
sucess_times <- rep() #良品的個數
 
for (i in 1:B)
{
  X[i] <- quantile(rbinom(B,N,p)/N, 0.975)
 ...

1
2
3
4
5
6
7
8
9
N = 20000#sample size
X = rbinom(1,N,0.997)
x = rbinom(1000, N, X/N)
 
M = 40000
Y = rbinom(1, M, 0.995)
y = rbinom(1000, M, Y/M)
 
z = x/N-y/M

1
Print"12234456"

1
2
3
4
5
6
7
8
9
song<-function(n=99) {
  a<-" bottle"
  b<-" of beer on the wall, "
  c<-" of beer. "
  d<-"Take one down and pass it around: "
  #
  ## set "beer counter" vector
  counter<-c(as.character(n:1),"No more")
  # 
 ...

1
2
3
4
5
6
7
8
9
song<-function(n=99) {
  a<-" bottle"
  b<-" of beer on the wall, "
  c<-" of beer. "
  d<-"Take one down and pass it around: "
  #
  ## set "beer counter" vector
  counter<-c(as.character(n:1),"No more")
  # 
 ...

1
2
3
4
5
6
7
8
9
song<-function(n=99) {
  a<-" bottle"
  b<-" of beer on the wall, "
  c<-" of beer. "
  d<-"Take one down and pass it around: "
  #
  ## set "beer counter" vector
  counter<-c(as.character(n:1),"No more")
  # 
 ...

1
 

1
2
3
4
5
 
 
x1<-c(1,2,3)
y1<-c(1,2,10)
plot(x1,y1)

1
2
3
4
5
6
7
8
 
 
aA<-matrix(scan(),byrow=T,ncol=3)
1 2 3
4 5 6
7 8 10
 
 
 ...

1
2
3
4
5
6
7
8
 
 
aA<-matrix(scan(),byrow=T,ncol=3)
1 2 3
4 5 6
7 8 10
 
 
 ...

1
2
3
4
5
6
7
8
 
 
aA<-matrix(scan(),byrow=T,ncols=3)
1 2 3
4 5 6
7 8 10
 
 
 ...

1
2
3
4
5
6
7
8
 
 
aA<-matrix(scan(),byrow=T,col=3)
1 2 3
4 5 6
7 8 10
 
 
 ...

1
2
3
4
5
6
7
8
 
 
aA<-matrix(scan(),byrow=T,col3)
1 2 3
4 5 6
7 8 10
 
 
 ...

1
2
3
4
5
 
 
c(1,2,3)
 
 

1
2
3
4
5
 
 
c(1,2,3)
 
print c

1
2
3
4
5
6
7
8
9
song<-function(n=99) {
  a<-" bottle"
  b<-" of beer on the wall, "
  c<-" of beer. "
  d<-"Take one down and pass it around: "
  #
  ## set "beer counter" vector
  counter<-c(as.character(n:1),"No more")
  # 
 ...