The R Project SVN R

Rev

Rev 27542 | Rev 30463 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 27542 Rev 27717
Line 1... Line 1...
1 
 
 1 
 
2 
R : Copyright 2003, The R Foundation for Statistical Computing
 2 
R : Copyright 2003, The R Foundation for Statistical Computing
3 
Version 1.9.0 Under development (unstable) (2003-12-12), ISBN 3-900051-00-3
 3 
Version 1.9.0 Under development (unstable) (2003-12-20), ISBN 3-900051-00-3
4 
 
 4 
 
5 
R is free software and comes with ABSOLUTELY NO WARRANTY.
 5 
R is free software and comes with ABSOLUTELY NO WARRANTY.
6 
You are welcome to redistribute it under certain conditions.
 6 
You are welcome to redistribute it under certain conditions.
7 
Type 'license()' or 'licence()' for distribution details.
 7 
Type 'license()' or 'licence()' for distribution details.
8 
 
 8 
 
Line 54... Line 54...
54 
 
 54 
 
55 
> cat("\nNumber of base objects:\t\t", length(ls.base),
55 
> cat("\nNumber of base objects:\t\t", length(ls.base),
56 
    "\nNumber of builtin functions:\t", sum(base.is.f), "\n\t starting with 'is.' :\t ",
56 
    "\nNumber of builtin functions:\t", sum(base.is.f), "\n\t starting with 'is.' :\t ",
57 
    length(is.bi <- bi[substring(bi, 1, 3) == "is."]), "\n")
 57 
    length(is.bi <- bi[substring(bi, 1, 3) == "is."]), "\n")
58 
 
 58 
 
59 
Number of base objects:		 991
59 
Number of base objects:		 992
60 
Number of builtin functions:	 964
60 
Number of builtin functions:	 965
61 
	 starting with 'is.' :	  42
61 
	 starting with 'is.' :	  42
62 
 
 62 
 
63 
> is.primitive <- function(obj) is.function(obj) &&
63 
> is.primitive <- function(obj) is.function(obj) &&
64 
    is.null(args(obj))
 64 
    is.null(args(obj))
65 
 
 65 
 
Line 965... Line 965...
965 
 
 965 
 
966 
 
 966 
 
967 
	demo(glm.vr)
 967 
	demo(glm.vr)
968 
	---- ~~~~~~
 968 
	---- ~~~~~~
969 
 
 969 
 
- 
 
 970 
> require(stats)
- 
 
 971 
[1] TRUE
- 
 
 972 
 
970 
> Fr <- c(68, 42, 42, 30, 37, 52, 24, 43, 66, 50, 33,
973 
> Fr <- c(68, 42, 42, 30, 37, 52, 24, 43, 66, 50, 33,
971 
    23, 47, 55, 23, 47, 63, 53, 29, 27, 57, 49, 19, 29)
 974 
    23, 47, 55, 23, 47, 63, 53, 29, 27, 57, 49, 19, 29)
972 
 
 975 
 
973 
> Temp <- gl(2, 2, 24, labels = c("Low", "High"))
 976 
> Temp <- gl(2, 2, 24, labels = c("Low", "High"))
974 
 
 977 
 
Line 1023... Line 1026...
1023 
 
 1026 
 
1024 
> detg.mod <- glm(terms(Fr ~ M.user * Temp * Soft +
1027 
> detg.mod <- glm(terms(Fr ~ M.user * Temp * Soft +
1025 
    Brand * M.user * Temp, keep.order = TRUE), family = poisson,
1028 
    Brand * M.user * Temp, keep.order = TRUE), family = poisson,
1026 
    data = detg)
 1029 
    data = detg)
1027 
 
 1030 
 
1028 
> summary(detg.mod, correlation = FALSE)
 1031 
> summary(detg.mod)
1029 
 
 1032 
 
1030 
Call:
 1033 
Call:
1031 
glm(formula = terms(Fr ~ M.user * Temp * Soft + Brand * M.user *
1034 
glm(formula = terms(Fr ~ M.user * Temp * Soft + Brand * M.user *
1032 
    Temp, keep.order = TRUE), family = poisson, data = detg)
 1035 
    Temp, keep.order = TRUE), family = poisson, data = detg)
1033 
 
 1036 
 
Line 1228... Line 1231...
1228 
sm     3    3    3    3    4    4    4    4    2     2     2
 1231 
sm     3    3    3    3    4    4    4    4    2     2     2
1229 
 
 1232 
 
1230 
> par(op)
 1233 
> par(op)
1231 
>
1234 
>
1232 
> cat("Time elapsed: ", proc.time() - .ptime, "\n")
 1235 
> cat("Time elapsed: ", proc.time() - .ptime, "\n")
1233 
Time elapsed:  3.1 0.12 3.22 0 0
1236 
Time elapsed:  3.07 0.13 3.2 0 0
1234 
>
1237 
>