A million thanks Peter, subsetting fixed it.
George

On Mon, Oct 12, 2009 at 11:51 PM, Peter Ehlers <ehlers@ucalgary.ca> wrote:

> Hi George,
>
> Your problem is not with xyplot, but with the NA occurrences in
> your data. Try adding
>
>  subset = {!is.na(MSE)},
>
> to your xyplot call, or (better), subset the data before
> calling xyplot.
>
>  -Peter Ehlers
>
>
> George Kalema wrote:
>
>> Hi Peter (and anyone else willing to help me out),
>> Many thanks for your help. Having used your code plus a few other
>> modifications, I only get the points plotted but without the two lines. I
>> just cannot figure out what the problem is.
>>
>> My code is as follows:
>>
>> library(lattice)
>> datos2 <- subset(datos, samplesize != 10 & parm != "Theta0")
>> unq <- sort(unique(datos2$samplesize))
>> datos2$fsamplesize <- factor(datos2$samplesize, labels = paste("Sample
>> size =", unq))
>> datos2$parm <- factor(datos2$parm, levels = c("Intercept", "time",
>> "trt", "time*trt"))
>> tp1.sim <- xyplot(MSE ~ ntimes | fsamplesize + parm, group = group, data
>> = datos2,
>>    type = "b", lty = 1:2, pch = 1:2,
>>    scales = list(x = list(at = c(2, 4, 8, 16)), alternating = 1),
>>    as.table = TRUE, key = list(text = list(c("GNA", "PNA")), points =
>> list(pch = 1:2))
>> )
>> plot(tp1.sim)
>>
>> I have attached my real dataset (called datos) as well.
>>
>> Kind appreciations to your efforts.
>>
>> George
>>
>>
>> On Wed, Oct 7, 2009 at 9:20 AM, Peter Ehlers <ehlers@ucalgary.ca> wrote:
>>
>>  see below
>>>
>>> George Kalema wrote:
>>>
>>>  Dear R users,
>>>> I am trying to have an xyplot of a data set which has the following
>>>> variables:
>>>> case (n=10,20,30)
>>>> parameter (parm=a,b)
>>>> group (grp=g1,g2)
>>>> y (y values)
>>>> x (x=2,4,8)
>>>>
>>>> My plot should be parameter by case such that I have 2 rows (each row=
>>>> each
>>>> parameter) and 3 columns (each column=each case). My R-code is as
>>>> follows
>>>> but I am not able to get what I want to:
>>>>
>>>> tp1.sim <- xyplot(y~ x | case + parm , group=group, data = data, lty =
>>>> 1:4
>>>> ,
>>>> pch = 1:4)
>>>> print(tp1.sim)
>>>>
>>>> How can I have two lines (for g1 and g2) in each plot (each box)?
>>>>
>>>>  include the type="b" argument
>>>
>>>  How do I label the x-axis with only values 2, 4, 8?
>>> include the scales= argument or make x a factor
>>>
>>>  How do I label each column with the corresponding case number?
>>> make 'case' a factor
>>>
>>> The following should do what you want:
>>>
>>> xyplot(y ~ x | factor(case) + parm, group=group, data=data,
>>>   type='b', lty=1:2, pch=1:2,
>>>   scales=list(x=list(at=c(2,4,8)))
>>> )
>>>
>>> I don't understand why you want 4 line types/point chars.
>>>
>>>  -Peter Ehlers
>>>
>>>
>>>  My hypothetical data set is as follows:
>>>>
>>>> parm x case y group
>>>> a 2 10 0.03 g1
>>>> b 2 10 0.02 g1
>>>> a 4 10 0.03 g1
>>>> b 4 10 0.02 g1
>>>> a 8 10 0.03 g1
>>>> b 8 10 0.02 g1
>>>> a 2 20 0.03 g1
>>>> b 2 20 0.02 g1
>>>> a 4 20 0.03 g1
>>>> b 4 20 0.02 g1
>>>> a 8 20 0.03 g1
>>>> b 8 20 0.02 g1
>>>> a 2 30 0.03 g1
>>>> b 2 30 0.02 g1
>>>> a 4 30 0.03 g1
>>>> b 4 30 0.02 g1
>>>> a 8 30 0.03 g1
>>>> b 8 30 0.02 g1
>>>> a 2 10 0.13 g2
>>>> b 2 10 0.12 g2
>>>> a 4 10 0.13 g2
>>>> b 4 10 0.12 g2
>>>> a 8 10 0.13 g2
>>>> b 8 10 0.12 g2
>>>> a 2 20 0.13 g2
>>>> b 2 20 0.12 g2
>>>> a 4 20 0.13 g2
>>>> b 4 20 0.12 g2
>>>> a 8 20 0.13 g2
>>>> b 8 20 0.12 g2
>>>> a 2 30 0.13 g2
>>>> b 2 30 0.12 g2
>>>> a 4 30 0.13 g2
>>>> b 4 30 0.12 g2
>>>> a 8 30 0.13 g2
>>>> b 8 30 0.12 g2
>>>>
>>>> Many thanks in advance for your response.
>>>>
>>>> George
>>>>
>>>>       [[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.
>>>>
>>>>
>>>>
>>>>
>>
>>


-- 
........................................
George Williams KALEMA,

Schapenstraat 37/282,
3000 Leuven,
Belgium.

Cell: +32 495 33 13 02
........................................

	[[alternative HTML version deleted]]


