Re: [R] Svar: Re: help with simple R-question

From: John Kane <jrkrideau_at_yahoo.ca>
Date: Wed, 06 Jun 2007 07:46:56 -0400 (EDT)


I think that you have the same situation as before though I have never used ranef(). The boar ids are acting as the row names and are not really part of the data.frame. It just looks like that when R prints the data.frame.

Try

boars <- rownames(far1)
far1 <- cbind( boars, far1)

The results may look a bit funny with two printed columns looking the same but one will be a column in the data.frame and the other will be the row names column. I hope :)

>
> The left column is boar id number, and the right is
> the random effect
> estimate. I need the numbers in the left column when
> i merge far1
> together with other data.frames based on the id
> numbers. When i use
> ranef the output is the class list and R only sees
> the intercepts, but i
> need a data.frame with 'boar id' and 'niveau', two
> columns not just
> one...
>
> fx
> > far1 <- ranef(resultat1)[1]
> > far1
> $farid
> (Intercept)
> 1437987 -3.775851e-03
> 10000100 -3.220044e-03
> 10007995 1.848914e-02
> 10008295 -4.583903e-03
> 10008792 -9.518371e-03
> 10033401 -7.538132e-03
> 10048900 1.540309e-02
> > far1 <- as.data.frame(far1)
> > far1
> X.Intercept.
> 1437987 -3.775851e-03
> 10000100 -3.220044e-03
> 10007995 1.848914e-02
> 10008295 -4.583903e-03
> 10008792 -9.518371e-03
> 10033401 -7.538132e-03
> 10048900 1.540309e-02
> Thanks again
>
> Rina
>
>
> >>> John Kane <jrkrideau_at_yahoo.ca> 06/06/07 11:38
> >>>
>
> --- Rina Miehs <rmi_at_danishmeat.dk> wrote:
>
> > thanks, that works great!!
> >
> > just have another thing...i the same area
> > What if the class is list instead of array, how
> can
> > you name the first unrecognized column?
>
> I am not sure that I understand the question. You
> don't really have an unrecognised column in the
> dataframe but an array of rownames ( I don't know
> how
> they are stored).
>
> I think you can do the same thing as you did for the
> data.frame but as I say, I am not sure I understand
> the question. Would you post a little example?
>
>
> >
> > Rina
> >
> > >>> John Kane <jrkrideau_at_yahoo.ca> 06/05/07 3:17
> >>>
> >
> > --- Rina Miehs <rmi_at_danishmeat.dk> wrote:
> >
> > > hello
> > >
> > > what do i write for R to recognize both columns?
> > >
> > > In the R-script downunder you can see that i use
> > > tapply to get my
> > > information out of my data, and then i need to
> use
> > > it as a dataframe
> > > with both columns! It is like R is using the
> first
> > > column as an
> > > observationnumber or something, how can i change
> > > that??
> >
> > It is using the names of the variables as
> rownames.
> >
> > try
> > n.ant <- names(antall)
> > antal1 <- data.frame(n.antal1, antal1)
> >
> >
> > >
> > > > antal1 <-tapply(l1$omlob1,l1$farid,length)
> > > > antal1
> > > 1437987 10000100 10007995 10008295 10008792
> > > 10010203 10018703
> > > 10033401
> > > 2 3 3 2
> 3
> >
> > > 1 1
> > > 2
> > > 10048900 10050492 10055897 10076495
> 10081892
> > > 10094801 10100692
> > > 10101395
> > > 3 1 3 3
> 6
> >
> > > 2 5
> > > 20
> > > 10101495 10101595 10104692 10113592
> 10113697
> > > 10114297 10120797
> > > 10120897
> > > 1 5 4 2
> 6
> >
> > > 11 1
> > > 4
> > > 10121697 10121897 10121997 10133592
> 10142892
> > > 10142995 10146495
> > > 10150497
> > > 16 3 6 1
> 1
> >
> > > 6 4
> > > 4
> > > 10150692 10157092 10157292 10164792
> 10170892
> > > 10171795 10171895
> > > 10172300
> > > 5 2 4 4
> 4
> >
> > > 4 4
> > > 1
> > > 10175195 10187802 10192499 10192897
> 10198295
> > > 10200493 10201693
> > > 10211593
> > > 1 2 2 3
> 5
> >
> > > 1 3
> > > 5
> > > > antal1 <- data.frame(antal1)
> > > > antal1
> > > antal1
> > > 1437987 2
> > > 10000100 3
> > > 10007995 3
> > > 10008295 2
> > > 10008792 3
> > > 10010203 NA
> > > 10018703 NA
> > > 10033401 2
> > > 10048900 3
> > > 10050492 1
> > > 10055897 3
> > > 10076495 3
> > > 10081892 6
> > > 10094801 2
> > > 10100692 5
> > >
> > > Thanks
> > > Rina
> > >
> > > [[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 ( http://www.r/ ) ( http://www.r/
> > )-project.org/posting-guide.html
> > > and provide commented, minimal, self-contained,
> > > reproducible code.
> > >
> >
> >
> >
> > Get news delivered with the All new Yahoo!
> > Mail. Enjoy RSS feeds
> > right on your Mail page. Start today at
> > http://mrd.mail.yahoo.com/try_beta?.intl=ca
> >
> >
>
>
>
> Be smarter than spam. See how smart SpamGuard
> is at giving junk

> http://mrd.mail.yahoo.com/try_beta?.intl=ca
>
>
>

=== message truncated ===



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. Received on Wed 06 Jun 2007 - 12:10:29 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 06 Jun 2007 - 13:31:43 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.