[R] Help on 'find.BIB' function
Chuck Cleland
ccleland at optonline.net
Wed Apr 25 20:52:18 CEST 2007
Jason Parcon wrote:
> Hello everyone,
>
> I am trying to use the 'find.BIB' function to construct a balanced incomplete block design. When I ran the example given in the help file (find.BIB(10,30,4)), I obtained the following error message:
>
> Error in optBlock(~., withinData = factor(1:trt), blocksize = rep(k, b)) :
> object ".Random.seed" not found
>
> I investigated what the optBlock function is doing but the manual / help files did not give me any information regarding the above error.
>
> I hope somebody can help me regarding this matter.
The following seems to work for me:
> library(crossdes)
Loading required package: AlgDesign
Loading required package: gtools
Loading required package: MASS
> set.seed(671969)
> find.BIB(10,30,4)
[,1] [,2] [,3] [,4]
[1,] 4 5 7 10
[2,] 1 2 3 10
[3,] 1 5 6 10
[4,] 2 8 9 10
[5,] 3 5 6 7
[6,] 3 4 9 10
[7,] 1 5 8 9
[8,] 1 6 7 9
[9,] 1 2 4 7
[10,] 2 6 8 10
[11,] 2 3 5 7
[12,] 1 6 7 9
[13,] 2 6 7 10
[14,] 1 2 3 9
[15,] 2 5 6 8
[16,] 2 4 5 9
[17,] 3 4 6 8
[18,] 1 5 8 10
[19,] 2 4 7 8
[20,] 3 6 9 10
[21,] 1 2 4 6
[22,] 3 7 8 10
[23,] 2 3 5 9
[24,] 1 4 5 10
[25,] 4 6 8 9
[26,] 4 7 9 10
[27,] 1 3 7 8
[28,] 3 4 5 6
[29,] 5 7 8 9
[30,] 1 3 4 8
I get the same error you report if I don't do the set.seed() step.
> sessionInfo()
R version 2.4.1 Patched (2007-03-31 r41127)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
attached base packages:
[1] "stats" "graphics" "grDevices" "utils" "datasets" "methods"
[7] "base"
other attached packages:
crossdes MASS gtools AlgDesign
"1.0-7" "7.2-33" "2.3.1" "1.0-7"
> Best regards,
>
> Jason Parcon
>
>
>
> ---------------------------------
>
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Chuck Cleland, Ph.D.
NDRI, Inc.
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212) 845-4495 (Tu, Th)
tel: (732) 512-0171 (M, W, F)
fax: (917) 438-0894
More information about the R-help
mailing list