Dear Dennis
Thanks for the comments, I read thru the script that was originally
created and gathered that the issue was because of the use of index.cond
in the xyplot function.
I modified the script to drop the index.cond parameter in the function and
now it works to give the subject order as indicated by ordered levels of
the Subj factor.
Following this however, I have another requirement for which i am seeking
the help of R Group again.
I am providing the illsutrative data and the script to run the lattice
plot.
What i require is to get the individual panels in each lattice plot carry
a red or a blue background when the slope of the fitted line is "-ve" or
"+"ve respectively.
Woule be greatful to receive help on how this could be acheived.
library(lattice)
MyData<-structure(list(Subj = structure(c(1L, 2L, 3L, 4L, 1L, 2L, 3L,
4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L,
4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L,
4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L,
4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L), .Label = c("S1", "S2", "S3",
"S4"), class = "factor"), Product = structure(c(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, 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), .Label = c("a",
"b"), class = "factor"), Arm = structure(c(1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("L",
"R"), class = "factor"), Attribute = 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, 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("A",
"B"), class = "factor"), Time = structure(c(1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("T0",
"T1"), class = "factor"), value = c(44, 33, 31, 34, 36, 33, 38,
32, 33, 35, 44, 38, 33, 36, 37, 34, 32, 35, 37, 33, 35, 41, 39,
40, 35, 32, 32, 44, 42, 30, 43, 38, 13, 13, 14, 12, 15, 11, 15,
14, 10, 15, 11, 15, 15, 10, 13, 13, 15, 11, 11, 12, 14, 10, 13,
13, 10, 12, 10, 13, 13, 15, 13, 10)), .Names = c("Subj", "Product",
"Arm", "Attribute", "Time", "value"), row.names = c(NA, 64L), class =
"data.frame", na.action = structure(65:352, .Names = c("65",
"66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76",
"77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87",
"88", "89", "90", "91", "92", "93", "94", "95", "96", "97", "98",
"99", "100", "101", "102", "103", "104", "105", "106", "107",
"108", "109", "110", "111", "112", "113", "114", "115", "116",
"117", "118", "119", "120", "121", "122", "123", "124", "125",
"126", "127", "128", "129", "130", "131", "132", "133", "134",
"135", "136", "137", "138", "139", "140", "141", "142", "143",
"144", "145", "146", "147", "148", "149", "150", "151", "152",
"153", "154", "155", "156", "157", "158", "159", "160", "161",
"162", "163", "164", "165", "166", "167", "168", "169", "170",
"171", "172", "173", "174", "175", "176", "177", "178", "179",
"180", "181", "182", "183", "184", "185", "186", "187", "188",
"189", "190", "191", "192", "193", "194", "195", "196", "197",
"198", "199", "200", "201", "202", "203", "204", "205", "206",
"207", "208", "209", "210", "211", "212", "213", "214", "215",
"216", "217", "218", "219", "220", "221", "222", "223", "224",
"225", "226", "227", "228", "229", "230", "231", "232", "233",
"234", "235", "236", "237", "238", "239", "240", "241", "242",
"243", "244", "245", "246", "247", "248", "249", "250", "251",
"252", "253", "254", "255", "256", "257", "258", "259", "260",
"261", "262", "263", "264", "265", "266", "267", "268", "269",
"270", "271", "272", "273", "274", "275", "276", "277", "278",
"279", "280", "281", "282", "283", "284", "285", "286", "287",
"288", "289", "290", "291", "292", "293", "294", "295", "296",
"297", "298", "299", "300", "301", "302", "303", "304", "305",
"306", "307", "308", "309", "310", "311", "312", "313", "314",
"315", "316", "317", "318", "319", "320", "321", "322", "323",
"324", "325", "326", "327", "328", "329", "330", "331", "332",
"333", "334", "335", "336", "337", "338", "339", "340", "341",
"342", "343", "344", "345", "346", "347", "348", "349", "350",
"351", "352"), class = "omit"))
MyData$Product<-factor(MyData$Product, levels = c("a","b"))
MyData$Time = factor(MyData$Time,levels=c("T0","T1"))
MyData$Subj <- factor(MyData$Subj, levels = c(paste('S', 1:4, sep = '')))
MyData<-MyData
Attributeindata<-levels(MyData$Attribute)
productindata<-levels(MyData$Product)
MyData<-MyData[ order(MyData$Subj) ,]
pdf(file="D:/latticeplotsbysubject.pdf")
for(i in 1:length(Attributeindata))
{
data<-subset(MyData,Attribute==Attributeindata[i])
for(j in 1:length(productindata))
{
data2<-subset(data,Product==productindata[j])
print(xyplot(value ~ Time | Subj, data=data2,layout = c(6,3), aspect =
"xy",as.table = TRUE, type = c("g", "p", "r"),drop.unused.levels=TRUE,ylab
= paste("value of",Attributeindata[i],"for Product",productindata[j]),
xlab = "Time",))
}
}
graphics.off()
Regards
Vijayan Padmanabhan
"What is expressed without proof can be denied without proof" - Euclide.
Dennis Murphy <djmu...@gmail.com>
06/10/2011 02:26 AM
To
Vijayan Padmanabhan <v.padmanab...@itc.in>
cc
r-help@r-project.org
Subject
Re: [R] lattice plot query
Hi:
Try the following - the illustration below is for exactly one set of
panel plots. Adapt as necessary to your situation:
# Reorder the levels of subject - in this case I put them in a new
variable
MyData$subj <- factor(MyData$Subj, levels = c(paste('S', 1:13, sep = '')))
xyplot(value ~ Time | subj, data = MyData, as.table = TRUE,
layout = c(3, 4), aspect = 'xy', type = c('g', 'p', 'r'),
drop.unused.levels=TRUE, xlab = 'Time')
as.table = TRUE starts the panel order from the upper left rather than
the lower left.
HTH,
Dennis
On Thu, Jun 9, 2011 at 4:00 AM, Vijayan Padmanabhan
<v.padmanab...@itc.in> wrote:
> Dear R Group
> I have the following data for which I am trying to create subject wise
> lattice plot for a given attribute and product .
> though the lattice plot is generated, for some reasons that i dont
> understand in each plot the subject panels take a random order, I would
> rather want all the plots to display the subject order in the same order
> as how i have ordered this particular factor level.
>
> Would appreciate if someone would provide me the parameter that i need
to
> pass to the xyplot function of library(lattice) to get this sorted out.
>
>
> ##the data for illustrative purpose..is as follows:
>
> MyData<- structure(list(Subj = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L,
> 3L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
> 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L,
> 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 7L, 7L, 7L, 7L, 7L,
> 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 11L, 11L, 11L, 11L,
> 11L, 11L, 11L, 11L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L,
> 5L, 5L, 5L, 5L, 5L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L,
> 8L, 8L, 8L, 8L, 8L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 3L,
> 3L, 3L, 3L, 3L, 3L, 3L, 3L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 9L,
> 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L,
> 10L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 6L, 6L, 6L, 6L, 6L, 6L,
> 6L, 6L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L,
> 10L, 10L, 10L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L,
> 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L,
> 5L, 5L, 5L, 5L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L,
> 8L, 8L, 8L, 8L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 4L, 4L,
> 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 7L, 7L,
> 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 11L,
> 11L, 11L, 11L, 11L, 11L, 11L, 11L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,
> 3L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 9L, 9L, 9L, 9L, 9L, 9L, 9L,
> 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L), .Label = c("S1",
> "S10", "S11", "S12", "S13", "S3", "S4", "S5", "S6", "S8", "S9"
> ), class = "factor"), Product = 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, 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, 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, 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, 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, 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, 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), .Label = c("a",
> "b"), class = "factor"), Arm = 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, 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, 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, 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, 2L, 2L, 2L, 2L, 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, 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, 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, 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, 2L, 2L, 2L, 2L), .Label = c("L",
> "R"), class = "factor"), Attribute = 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, 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, 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, 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, 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, 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, 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, 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, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label =
c("AttriA",
> "AttriB"), class = "factor"), Time = structure(c(1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L,
> 4L, 5L, 6L, 7L, 8L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L), .Label = c("T0",
> "T1", "T2", "T3", "T4", "T5", "T6", "T7"), class = "factor"),
> value = c(39.5527423278832, 38.2473086972328, 40.2276345888737,
> 37.5267332387073, 40.9688951303111, 39.0569778014684,
> 39.1695786007607,
> 38.9139593338983, 68.8265651458469, 72.0887495716813,
> 67.5611670846489,
> 66.9874884774268, 66.9829433510606, 54.8956638696188,
> 63.0168969989417,
> 65.6194764863213, 66.4129384714848, 66.0961638072511,
> 68.9661166309607,
> 68.7059517735383, 65.7542747314885, 65.5069438367477,
> 64.7683246003315,
> 61.4962919814996, 69.5762947663818, 73.3430582378711,
> 76.8645584330803,
> 68.2802524581975, 77.8816834151128, 74.2317212282283,
> 70.4892680115894,
> 73.5430313565372, 45.0183417657248, 48.1530494831865,
> 46.1558267920891,
> 48.6541405845911, 48.4048804759923, 48.590343912936,
44.7362547527865,
>
> 42.6373136974049, 39.9917820318598, 43.4586441692621,
> 41.6364183563171,
> 42.8906355926352, 45.5495273145838, 41.3577962405852,
41.496338708892,
>
> 40.0151968165421, 45.5648864982887, 46.2956996070121,
> 47.0280057728853,
> 50.1831679877159, 51.061682921112, 45.9831040407059,
44.9747480866787,
>
> 49.6665131229635, 54.8281039296233, 54.1930629171256,
> 54.4341268439816,
> 57.1044051058105, 56.1143459282325, 52.6161245346639,
> 52.7130458240943,
> 55.4358255649899, 47.2220107045472, 47.3389267304398,
> 46.2106468232384,
> 48.9990363616358, 46.8056280265232, 46.5981260727775,
> 46.4927574552312,
> 40.3333491167334, 62.2287823821622, 66.3810749127077,
> 67.2394807083061,
> 62.4429875125554, 69.0376706648671, 60.2812609331127,
61.455201318328,
>
> 59.9863854008917, 62.2997182356239, 68.3165884503308,
> 70.8613200740368,
> 75.5343092931554, 68.0613856491741, 67.8297794393095,
69.723445715177,
>
> 68.2294718401127, 52.8140766072603, 51.4240493093108,
> 49.1567458446157,
> 52.7574433410593, 51.8096008153661, 51.0542236981057,
> 49.1650793229206,
> 46.1961387381546, 39.3501757447082, 40.6370147384917,
> 42.6126353906914,
> 42.8371247270351, 45.6224375518759, 38.1946693028616,
> 42.6841902106385,
> 41.2891110860076, 46.5589402020299, 46.9236938185125,
> 46.7985819317017,
> 52.9563382730808, 52.7678198507564, 42.6784974818896,
> 46.0254807534733,
> 49.3617924018815, 55.0347336430881, 53.0553009508891,
> 59.1754785059011,
> 59.2551361838499, 58.7606127906148, 56.0068865579661,
> 59.7877691476509,
> 59.9894635155904, 45.8292743149289, 44.3975346123602,
> 44.0878763837905,
> 43.0034203212678, 44.7694956749990, 44.5206026642237,
> 44.1129676323369,
> 41.1441655434737, 64.502950186561, 66.2596818647181,
66.9718770730599,
>
> 68.3492812036726, 71.1329167451021, 60.993158766883,
64.0476800619917,
>
> 59.5722548586232, 69.1060115847223, 68.8225246639702,
> 71.3984976174081,
> 80.9630180588849, 73.007714590984, 65.7790674687866,
68.8673355306755,
>
> 66.8813383526656, 40.0423309528637, 44.1216443512537,
> 41.9734395557559,
> 44.8434380300801, 43.7494322571724, 43.3183356186337,
> 45.5420445201619,
> 42.5147260103379, 68.7745002341477, 61.981198669435,
61.7864827933066,
>
> 60.3458594826638, 59.9697661699882, 53.9643718467014,
> 59.1755018103991,
> 59.2276568453173, 62.7915792989106, 71.0620731987689,
> 72.2036228167234,
> 70.9247270575656, 67.1627603989188, 69.8883404495198,
> 63.0049593358823,
> 64.8912785345249, 68.8781924408982, 74.137601748426,
75.2732261184303,
>
> 71.669750963672, 74.4341280041274, 76.2382402210529,
69.7026459277572,
>
> 76.5074726515534, 29.2874363361783, 27.2892266662294,
> 29.5616488492521,
> 29.6852458522043, 30.5496200171305, 29.9966252445121,
> 29.1288540248101,
> 27.4806129668759, 33.5863501946243, 30.4931479044075,
> 29.5862807824957,
> 34.9382017751209, 35.238424232798, 29.4052475587535,
33.0798423075101,
>
> 31.3567180682194, 33.4690209827123, 30.6523035293150,
> 32.6208075805269,
> 33.9407265800815, 31.0650808724465, 36.7534951846668,
> 34.6967011594745,
> 31.9567044598027, 44.5158738295474, 40.4690204233383,
> 38.0597778223743,
> 37.8034937870776, 35.1720350896499, 39.8550085954103,
> 34.9099191361032,
> 42.0149439902365, 26.2151235666232, 23.5679852219729,
> 24.2001786793751,
> 23.6961503280406, 23.7185784815809, 25.3066274885425,
> 24.0292003617619,
> 24.2505451653069, 27.0291768503551, 27.9586095806653,
> 27.1539912074178,
> 27.3800648825246, 27.0086310697033, 29.5006091073482,
> 27.6618986740352,
> 26.6054039275563, 26.7544589952818, 25.7739723955002,
> 24.7493173361881,
> 24.8974539247832, 26.0066344196213, 25.7242948862424,
> 24.2129429825322,
> 26.2810081950791, 32.8501876144134, 27.2087875772192,
> 22.9241752520302,
> 23.9219026236064, 27.8607929641585, 27.7779664335406,
> 25.6974359198254,
> 26.0076906518561, 28.8549852933942, 27.2116802644414,
> 28.5718875220588,
> 30.5721859479267, 29.2724683799154, 29.7833383429998,
> 29.0348479362964,
> 27.1223342602551, 32.5582746876519, 28.4460181623792,
> 30.5181241094846,
> 27.8313626195291, 36.6535660010998, 31.2270378533636,
> 30.7505338399547,
> 28.3032529970344, 32.3568559611009, 34.9553492819238,
> 35.7759456755808,
> 36.0195354653364, 38.5285002594439, 35.7380230921731,
> 37.1484588196843,
> 38.5416308644471, 24.8506692697752, 24.6343520347229,
> 23.6083195613371,
> 24.7607264684248, 22.9854229667107, 23.9036751221927,
> 24.0387963674900,
> 22.7723194255264, 30.3793074740804, 28.8639832951117,
> 28.3030690013896,
> 27.8655679757535, 29.1997887609265, 30.1487803931287,
> 29.2987532171346,
> 28.5871764431944, 24.6444755508312, 24.7330180273872,
> 23.7338446480537,
> 25.8596684061866, 26.0787857361113, 24.8953653206929,
> 23.5625674439808,
> 24.479442599818, 29.6545075345955, 25.6090083380698,
27.1114528433399,
>
> 29.5090563448938, 31.455338578686, 31.5408279032055,
32.0051421579701,
>
> 33.4618122321904, 27.4889692245983, 25.5464291578372,
> 24.5476688346391,
> 25.2081446883923, 26.3672783411296, 25.6265928833026,
> 27.1950139669246,
> 27.130001618485, 28.6347421265670, 26.7649268696797,
28.0626929191564,
>
> 30.4520327581967, 30.4717296368313, 29.8812011276167,
> 30.4064647569007,
> 30.7831770362286, 37.7827415537242, 35.0879233395691,
> 34.7668091883955,
> 37.3427148382295, 38.3775395012458, 36.7485828173724,
38.193636145425,
>
> 34.5786595343394, 28.6171348871096, 27.9470627780942,
> 29.8327633149340,
> 29.6097653173267, 30.0836756702669, 31.0359954571051,
> 31.0662062327447,
> 28.8585430558031, 29.7252130104938, 26.5553131151419,
> 28.0655400327958,
> 28.7575356242266, 28.3515991163835, 26.4610318654439,
> 27.2362207214083,
> 28.0005735091048, 29.1041755824943, 29.5347668690827,
> 32.1870020976934,
> 31.1451531051784, 25.2860128452551, 29.1578475612553,
> 27.0033058437616,
> 24.4368356738120, 39.7697712026109, 35.1675811365275,
> 32.1336514906396,
> 35.5524206935849, 33.5035236183325, 34.2439798004073,
> 33.6508805746601,
> 33.8661875340464)), .Names = c("Subj", "Product", "Arm",
> "Attribute", "Time", "value"), row.names = c(NA, 352L), class =
> "data.frame")
>
>
> ## Here is the code to create lattice plot subject wise, for each
> attribute and product from the study..
> ##This code creates a pdf file output of all plots and saves the file in
> the destination given.
>
> library(lattice)
> Attributeindata<-levels(MyData$Attribute)
> productindata<-levels(MyData$Product)
> pdf(file="D:/Test/latticeplotsbysubject.pdf")
> for(i in 1:length(Attributeindata))
> {
> for(j in 1:length(productindata))
> {
> print(xyplot(value ~ Time | Subj,
> data=subset(MyData,Attribute==Attributeindata[i] & Product=="b"),layout
=
> c(6,3), aspect = "xy",, type = c("g", "p",
> "r"),drop.unused.levels=TRUE,index.cond = function(x,y) coef(lm(y ~
> x))[1],ylab = paste("value of",Attributeindata[i],"for
> Product",productindata[j]),
> xlab = "Time",))
> }
> }
>
>
> ##On examining the pdf file, you wil find that the order of subjects as
> appearing in each page of plot is different.
> ## I would want the order of subjects to appear in the same order as
> ##follows: "S1","S3","S4","S5","S6","S8","S9","S10","S11", "S12","S13"
> ##I have tried reorder function on the factor level "Subj" but it does
not
> give the desired result.
>
> Any help in modifying the code above to meet the requirement would be
> highly appreciated.
>
> Thanks & Regards
> Vijayan Padmanabhan
>
>
> "What is expressed without proof can be denied without proof" - Euclide.
>
> Please visit us at www.itcportal.com
>
******************************************************************************
> This Communication is for the exclusive use of the intended recipient
(s) and shall
> not attach any liability on the originator or ITC Ltd./its
Subsidiaries/its Group
> Companies. If you are the addressee, the contents of this email are
intended for your use only and it shall not be forwarded to any third
party, without first obtaining written authorisation from the originator
or ITC Ltd./its Subsidiaries/its Group Companies. It may contain
information which is confidential and legally privileged and the same
shall not be used or dealt with by any third party in any manner
whatsoever without the specific consent of ITC Ltd./its Subsidiaries/its
Group Companies.
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org 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.
>
Please visit us at www.itcportal.com
******************************************************************************
This Communication is for the exclusive use of the intended recipient (s) and
shall
not attach any liability on the originator or ITC Ltd./its Subsidiaries/its
Group
Companies. If you are the addressee, the contents of this email are intended
for your use only and it shall not be forwarded to any third party, without
first obtaining written authorisation from the originator or ITC Ltd./its
Subsidiaries/its Group Companies. It may contain information which is
confidential and legally privileged and the same shall not be used or dealt
with by any third party in any manner whatsoever without the specific consent
of ITC Ltd./its Subsidiaries/its Group Companies.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
