Wow! [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2013-01-09 17:10 (4552 d 17:56 ago) – Posting: # 9805
Views: 8,347

Dear Detlew!

❝ Do you believe the subject line?


First: no.

❝ If not, have a look on the following code :cool::


Second: yes. Unbelievable.

power.sim2x2(CV=0.3, n=12, GMR=1.25, details=T)

1e+06 sims. Time elapsed (sec):

   user  system elapsed

   0.50    0.05    0.55

[1] 0.033578


Amazing. With roundit=T I got 0.033605. Same value if I switch lines 65–66 to R’s round(). With my original code I got 0.033236 (unrounded) and 0.033786 (R’s rounding), but run-time was ~30 minutes.
Concerning the question at the top of your source code: I think so but I’m not sure as well. BTW, couldn’t figure out how to modify the code for imbalanced sequences.

Just discovered that I’m not the first simulator. See Chapter 4.2.1 of Chow & Liu (any edition).

Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,688 registered users;
33 visitors (0 registered, 33 guests [including 15 identified bots]).
Forum time: 12:06 CEST (Europe/Vienna)

Complex, statistically improbable things are by their nature
more difficult to explain than
simple, statistically probable things.    Richard Dawkins

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5