I have an application in which I want to plot points with different
colours and symbols to distinguish (three) types of point. I thought I
had it figured out, mostly by a trial-and-error approach and by
simple-mindedly following instructions from others, without (sad to say)
really understanding what I was doing.
I have now encountered an example where my syntax screws up, and I can't
see how to fix it in a coherent and generally applicable way. The
problem arises if the data column (factor) specifying the point types is
missing one level. When this happens, the types get out of synch.
I have attached a source-able file "demo.txt" and the two data sets upon
which it depends as "data1.txt" and "data2.txt" (dput() format).
The result from using the data from "data1.txt" (plotObj1) demonstrates
the problem, i.e. the unwanted phenomenon. The result from using the
data from "data2.txt" (plotObj2) is "as it should be".
In the data from data1.txt the "extremes" column, which specifies the
point type, is missing the level "0":
table(Dat1$extremes)
0 intermed 1
0 39 43
As you will be able to see, the plot of these data has the "0" symbol (a
triangle) where there should be an "intermed" symbol (a solid dot), and
the "intermed" symbol where there should be a "1" symbol (an upside-down
triangle).
The data from data2.txt have an extremes column with all three levels
and the corresponding plot is as desired.
I'm pretty sure that the solution to my problem is quite simple, but I
am too stupid to figure it out. Can anyone guide me in the right direction?
Thanks.
cheers,
Rolf Turner
--
Honorary Research Fellow
Department of Statistics
University of Auckland
Phone: +64-9-373-7599 ext. 88276
#
# Script demo.txt
#
library("ggforce")
Dat1 <- dget("data1.txt")
Dat2 <- dget("data2.txt")
shps <- c(24,20,25)
culs <- c("black","blue","black")
tkst <- c("0","intermed","1")
plotObj1 <- ggplot(Dat1) +
geom_point(aes(y = Sfrac, x = Dose,
shape = extremes, colour = extremes)) +
scale_shape_manual("Point type:",values = shps,labels=tkst) +
scale_colour_manual("Point type:",values = culs,labels=tkst) +
xlab("") + ylab("") +
facet_grid_paginate(facets=degC~Duration,nrow=2,ncol=2,page=1)
plotObj2 <- ggplot(Dat2) +
geom_point(aes(y = Sfrac, x = Dose,
shape = extremes, colour = extremes)) +
scale_shape_manual("Point type:",values = shps,labels=tkst) +
scale_colour_manual("Point type:",values = culs,labels=tkst) +
xlab("") + ylab("") +
facet_grid_paginate(facets=degC~Duration,nrow=2,ncol=2,page=1)
structure(list(Dead = c(33, 34, 73, 41, 39, 6, 20, 42, 24, 56,
41, 20, 2, 7, 25, 65, 80, 63, 79, 44, 3, 31, 39, 40, 81, 80,
103, 54, 4, 62, 92, 89, 40, 76, 2, 10, 17, 34, 64, 24, 72, 49,
7, 13, 88, 84, 69, 78, 46, 1, 3, 16, 29, 61, 7, 25, 31, 49, 70,
81, 23, 49, 6, 28, 20, 28, 44, 42, 3, 27, 17, 50, 58, 28, 12,
24, 9, 20, 22, 24, 26, 33), Alive = c(0, 0, 0, 0, 0, 29, 34,
2, 1, 2, 0, 0, 62, 40, 8, 7, 1, 2, 0, 0, 19, 4, 0, 0, 0, 0, 0,
0, 29, 1, 0, 0, 0, 0, 24, 50, 17, 10, 0, 0, 0, 0, 67, 26, 8,
0, 0, 0, 0, 20, 29, 18, 0, 0, 53, 14, 12, 9, 2, 0, 0, 0, 34,
0, 0, 0, 0, 0, 19, 28, 7, 0, 0, 0, 42, 80, 38, 9, 6, 0, 0, 0),
Duration = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("8hour",
"16hour"), class = "factor"), degC = structure(c(1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L), .Label = c("10deg", "20deg"), class = "factor"),
Rep = structure(c(3L, 3L, 3L, 3L, 3L, 7L, 7L, 7L, 7L, 7L,
7L, 7L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L, 5L, 5L, 9L, 9L, 9L, 9L,
9L, 9L, 9L, 9L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L,
8L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 2L, 2L, 2L, 2L,
2L, 2L, 6L, 6L, 6L, 6L, 6L, 6L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 10L), .Label = c("1", "2", "4", "5", "7", "8",
"10", "11", "13", "14", "16", "17"), class = "factor"), Dose = c(30,
40, 45, 50, 55, 0, 15, 20, 25, 30, 35, 40, 0, 10, 15, 20,
25, 30, 35, 40, 0, 10, 15, 18, 21, 24, 27, 37, 0, 10, 18,
21, 24, 37, 0, 8, 12, 15, 18, 22, 20, 24, 0, 6, 10, 14, 18,
26, 30, 0, 4, 8, 13, 18, 0, 4, 8, 10, 12, 14, 16, 18, 0,
6, 9, 13, 15, 18, 0, 3, 6, 9, 13, 18, 0, 2, 4, 6, 8, 10,
12, 14), Trt = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label =
c("16hour10deg",
"16hour20deg", "8hour10deg", "8hour20deg"), class = "factor"),
Sfrac = c(1, 1, 1, 1, 1, 0.171428571428571, 0.37037037037037,
0.954545454545455, 0.96, 0.96551724137931, 1, 1, 0.03125,
0.148936170212766, 0.757575757575758, 0.902777777777778,
0.987654320987654, 0.969230769230769, 1, 1, 0.136363636363636,
0.885714285714286, 1, 1, 1, 1, 1, 1, 0.121212121212121, 0.984126984126984,
1, 1, 1, 1, 0.0769230769230769, 0.166666666666667, 0.5, 0.772727272727273,
1, 1, 1, 1, 0.0945945945945946, 0.333333333333333, 0.916666666666667,
1, 1, 1, 1, 0.0476190476190476, 0.09375, 0.470588235294118,
1, 1, 0.116666666666667, 0.641025641025641, 0.720930232558139,
0.844827586206897, 0.972222222222222, 1, 1, 1, 0.15, 1, 1,
1, 1, 1, 0.136363636363636, 0.490909090909091, 0.708333333333333,
1, 1, 1, 0.222222222222222, 0.230769230769231, 0.191489361702128,
0.689655172413793, 0.785714285714286, 1, 1, 1), extremes = structure(c(3L,
3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 2L, 2L, 2L, 2L,
2L, 2L, 3L, 3L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 3L,
3L, 3L, 3L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 3L,
3L, 3L, 3L, 2L, 2L, 2L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 3L, 3L,
3L, 2L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L,
2L, 2L, 2L, 3L, 3L, 3L), .Label = c("0", "intermed", "1"), class =
"factor")), altReps = structure(list(
newRep = structure(c(6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 18L, 18L, 18L, 18L, 18L, 18L,
18L, 18L, 3L, 3L, 3L, 3L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 15L,
15L, 15L, 15L, 15L, 15L, 15L, 15L, 4L, 4L, 4L, 4L, 4L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 16L, 16L, 16L, 16L, 16L, 16L,
16L, 16L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 7L, 7L, 7L, 7L,
7L, 7L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 11L, 11L, 11L, 11L, 11L, 17L, 17L, 17L, 17L,
17L, 17L, 17L, 17L, 2L, 2L, 2L, 2L, 2L, 2L, 8L, 8L, 8L, 8L,
8L, 8L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L), .Label = c("1",
"2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12",
"13", "14", "15", "16", "17", "18"), class = "factor"), oldRep = c(1,
1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3,
3, 3, 3, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3,
3, 3, 3, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3,
3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3,
3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,
3, 3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,
2, 3, 3, 3, 3, 3, 3, 3, 3), repWithNames = structure(c(6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 12L, 12L, 12L, 12L, 12L, 12L,
12L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 3L, 3L, 3L,
3L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 15L, 15L, 15L, 15L, 15L,
15L, 15L, 15L, 4L, 4L, 4L, 4L, 4L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 7L, 7L, 7L, 7L, 7L, 7L, 13L, 13L,
13L, 13L, 13L, 13L, 13L, 13L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
11L, 11L, 11L, 11L, 11L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
17L, 2L, 2L, 2L, 2L, 2L, 2L, 8L, 8L, 8L, 8L, 8L, 8L, 14L,
14L, 14L, 14L, 14L, 14L, 14L, 14L), .Label = c("1.16.10",
"1.16.20", "1.16.5", "1.8.10", "1.8.20", "1.8.5", "2.16.10",
"2.16.20", "2.16.5", "2.8.10", "2.8.20", "2.8.5", "3.16.10",
"3.16.20", "3.16.5", "3.8.10", "3.8.20", "3.8.5"), class = "factor")),
row.names = c(NA,
124L), class = "data.frame"), data.set.name = "inLogs", trtVars = c("Duration",
"degC"), upn = "Dose", row.names = 43:124, class = c("doseResp",
"data.frame"))
structure(list(Dead = c(4, 4, 67, 39, 73, 68, 32, 39, 2, 3, 16,
66, 99, 26, 29, 5, 5, 20, 113, 87, 56, 30, 121, 7, 2, 21, 52,
0, 15, 29, 21, 21, 48, 54, 2, 19, 10, 53, 75, 88, 48, 66, 33,
34, 73, 41, 39, 6, 20, 42, 24, 56, 41, 20, 2, 7, 25, 65, 80,
63, 79, 44, 3, 31, 39, 40, 81, 80, 103, 54, 4, 62, 92, 89, 40,
76, 2, 10, 17, 34, 64, 24, 72, 49), Alive = c(50, 44, 20, 1,
0, 1, 1, 0, 45, 46, 4, 1, 0, 0, 0, 50, 76, 8, 0, 4, 0, 0, 0,
22, 18, 0, 0, 35, 12, 18, 0, 0, 0, 0, 28, 80, 74, 6, 2, 0, 0,
0, 0, 0, 0, 0, 0, 29, 34, 2, 1, 2, 0, 0, 62, 40, 8, 7, 1, 2,
0, 0, 19, 4, 0, 0, 0, 0, 0, 0, 29, 1, 0, 0, 0, 0, 24, 50, 17,
10, 0, 0, 0, 0), Duration = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("8hour",
"16hour"), class = "factor"), degC = structure(c(1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L), .Label = c("5deg", "10deg"), class = "factor"), Rep = structure(c(4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 12L, 2L, 2L, 2L, 2L, 6L, 6L, 6L,
6L, 6L, 6L, 6L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 3L, 3L,
3L, 3L, 3L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 11L, 11L, 11L, 11L, 11L,
11L, 11L, 11L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L,
5L, 5L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L), .Label = c("1", "3",
"4", "6", "7", "9", "10", "12", "13", "15", "16", "18"), class = "factor"),
Dose = c(0, 15, 30, 45, 55, 65, 75, 85, 0, 10, 30, 40, 50,
60, 80, 0, 10, 20, 30, 40, 50, 60, 80, 0, 10, 35, 50, 0,
10, 20, 35, 40, 45, 50, 0, 10, 20, 25, 30, 35, 45, 40, 30,
40, 45, 50, 55, 0, 15, 20, 25, 30, 35, 40, 0, 10, 15, 20,
25, 30, 35, 40, 0, 10, 15, 18, 21, 24, 27, 37, 0, 10, 18,
21, 24, 37, 0, 8, 12, 15, 18, 22, 20, 24), Trt = structure(c(4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L,
4L, 4L, 4L, 4L, 4L, 4L, 4L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("16hour10deg",
"16hour5deg", "8hour10deg", "8hour5deg"), class = "factor"),
Sfrac = c(0.0740740740740741, 0.0833333333333333, 0.770114942528736,
0.975, 1, 0.985507246376812, 0.96969696969697, 1, 0.0425531914893617,
0.0612244897959184, 0.8, 0.985074626865672, 1, 1, 1, 0.0909090909090909,
0.0617283950617284, 0.714285714285714, 1, 0.956043956043956,
1, 1, 1, 0.241379310344828, 0.1, 1, 1, 0, 0.555555555555556,
0.617021276595745, 1, 1, 1, 1, 0.0666666666666667, 0.191919191919192,
0.119047619047619, 0.898305084745763, 0.974025974025974,
1, 1, 1, 1, 1, 1, 1, 1, 0.171428571428571, 0.37037037037037,
0.954545454545455, 0.96, 0.96551724137931, 1, 1, 0.03125,
0.148936170212766, 0.757575757575758, 0.902777777777778,
0.987654320987654, 0.969230769230769, 1, 1, 0.136363636363636,
0.885714285714286, 1, 1, 1, 1, 1, 1, 0.121212121212121, 0.984126984126984,
1, 1, 1, 1, 0.0769230769230769, 0.166666666666667, 0.5, 0.772727272727273,
1, 1, 1, 1), extremes = structure(c(2L, 2L, 2L, 2L, 3L, 2L,
2L, 3L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 2L, 3L, 2L, 3L,
3L, 3L, 2L, 2L, 3L, 3L, 1L, 2L, 2L, 3L, 3L, 3L, 3L, 2L, 2L,
2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L,
2L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 2L, 2L, 3L, 3L,
3L, 3L, 3L, 3L, 2L, 2L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 3L,
3L, 3L, 3L), .Label = c("0", "intermed", "1"), class = "factor")), altReps
= structure(list(
newRep = structure(c(6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 12L,
12L, 12L, 12L, 12L, 12L, 12L, 18L, 18L, 18L, 18L, 18L, 18L,
18L, 18L, 3L, 3L, 3L, 3L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 15L,
15L, 15L, 15L, 15L, 15L, 15L, 15L, 4L, 4L, 4L, 4L, 4L, 10L,
10L, 10L, 10L, 10L, 10L, 10L, 16L, 16L, 16L, 16L, 16L, 16L,
16L, 16L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 7L, 7L, 7L, 7L,
7L, 7L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 5L, 5L, 5L,
5L, 5L, 5L, 5L, 11L, 11L, 11L, 11L, 11L, 17L, 17L, 17L, 17L,
17L, 17L, 17L, 17L, 2L, 2L, 2L, 2L, 2L, 2L, 8L, 8L, 8L, 8L,
8L, 8L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L), .Label = c("1",
"2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12",
"13", "14", "15", "16", "17", "18"), class = "factor"), oldRep = c(1,
1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3,
3, 3, 3, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3,
3, 3, 3, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3,
3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3,
3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,
3, 3, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,
2, 3, 3, 3, 3, 3, 3, 3, 3), repWithNames = structure(c(6L,
6L, 6L, 6L, 6L, 6L, 6L, 6L, 12L, 12L, 12L, 12L, 12L, 12L,
12L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 3L, 3L, 3L,
3L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 15L, 15L, 15L, 15L, 15L,
15L, 15L, 15L, 4L, 4L, 4L, 4L, 4L, 10L, 10L, 10L, 10L, 10L,
10L, 10L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 7L, 7L, 7L, 7L, 7L, 7L, 13L, 13L,
13L, 13L, 13L, 13L, 13L, 13L, 5L, 5L, 5L, 5L, 5L, 5L, 5L,
11L, 11L, 11L, 11L, 11L, 17L, 17L, 17L, 17L, 17L, 17L, 17L,
17L, 2L, 2L, 2L, 2L, 2L, 2L, 8L, 8L, 8L, 8L, 8L, 8L, 14L,
14L, 14L, 14L, 14L, 14L, 14L, 14L), .Label = c("1.16.10",
"1.16.20", "1.16.5", "1.8.10", "1.8.20", "1.8.5", "2.16.10",
"2.16.20", "2.16.5", "2.8.10", "2.8.20", "2.8.5", "3.16.10",
"3.16.20", "3.16.5", "3.8.10", "3.8.20", "3.8.5"), class = "factor")),
row.names = c(NA,
124L), class = "data.frame"), data.set.name = "inLogs", trtVars = c("Duration",
"degC"), upn = "Dose", row.names = c(NA, 84L), class = c("doseResp",
"data.frame"))
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
