similar to: Reading in data

Hello, from a software I have the following output in xml (see below): It is a series of matrices, for each age one. I have 3 categories (might vary in the application), hence, 3x3 matrices where each element gives the probability of transition from i to j. I would like read this data into R (preferably in a list, where each list element is one of the age specific matrices) and - after altering

Hi! I have a dataframe like this: dat=data.frame(Age=c(rep(30,8),rep(40,8),rep(50,8)),Period=rep(seq(2005,2008,1),3),Rate=c(seq(1,8,1),seq(9,16,1),seq(17,24,1)),Sex=rep(c(rep(0,4),rep(1,4)),3))attach(dat)dat    Age Period Rate Sex1   30   2005    1   02   30   2006    2   03   30   2007    3   04   30   2008    4   05   30   2005    5   16   30   2006    6   17   30   2007    7   18   30   2008  

Jim Lemon <bitwrit <at> ozemail.com.au> writes: > I have been thinking about how to reduce the number of basic questions that > elicit the ...ahem... robust debate that has occurred about how to answer The traffic on r-help could be reduced by creating a second list where more elementary questions are asked. There may be other ways to partition the universe of questions

dear£º in the example£¨nomogram£©£¬I don't understand the meanings of the program which have been marked by red line.And how to compile the program(L <- .4*(sex=='male') + .045*(age-50) + (log(cholesterol - 10)-5.2)*(-2*(sex=='female') + 2*(sex=='male'))). n <- 1000 # define sample size set.seed(17) # so can reproduce the results age <- rnorm(n, 50, 10)

Thank you for reading. I am trying to get sphericity values, and I understood I need to use mlm, but how do I implement a nested within subject design in mlm? I already read the R newsletter, fox chapter appendix, EZanova, and whatever I could find online. My original ANOVA anova(aov(resp ~ sucrose*citral, random =~1 | subject, data = p12bl, subset = exps==1)) Or anova(aov(resp ~

Dieter Menne pointed out that the (small) xml attachment didn't make it. Here is an in-line version (see end of message). Let's hope it works this time. I'm struggling with interpreting XML files created by ADODB as data.frames and I'm looking for advice. Note: This xlm contains a result set which comes from a rectangular data array. I've been trying to play with

Hi I have a dataset (.txt file and .dat file) in which the length of one record is 144. There is no header in the .txt or .dat file itself. When I read this file using the read.table command, and want to drop some coloumns by setting the argument colClasses to "NULL" for the columns that I want to drop, it does not work, because the entire record of the length 144 is being treated as

I am using rw0641. In my continuing quest to understand repeated measures analysis, I again return to lme. I exported the Potthoff and Roy data Orthodont.dat from S-PLUS 4.5 to avoid capture errors and ran the examples in the R help. I imported the data.frame with data <- read.table("Orthodont.dat",header=T) attach(data) and created the objects Orthodont.fit1 <-

Hi, I have recently been attempting to find the LD50 from two predicted fits (For male and females) in a Generalised linear model which models the effect of both sex + logdose (and sex*logdose interaction) on proportion survival (formula = y ~ ldose * sex, family = "binomial", data = dat (y is the survival data)). I can obtain the LD50 for females using the dose.p() command in the MASS

To whom this may concern, I would be very grateful if someone could give me some advice on where I am going wrong with a logistic regression I am trying to run. I am trying to run a logistic regression on an aggregated data set and have input the command: logistic.mod<-glm(x~Frequency+Location+Sex+Age.Group,family=binomial(link="logit"),data=earsag1.dat) where x is the count of my

Dear list, I have following plot, where I have set the color (red and green) and lines (lty=2:3) in the panel.groups but can't not figure out how change the lines and color of the legend in the "key" to the same lines and color as in the panel.groups. bwplot(means ~ age | scales , dat, panel = "panel.superpose", groups = sex,scales = list(x = list(rot =

I would like to read the following hierarchical data set. There is a family record followed by one or more personal records. If col. 7 is "1" it is a family record. If it is "2" it is a personal record. The family record is formatted as follows: col. 1-5 family id col. 7 "1" col. 9 dwelling type code The personal record is formatted as follows: col.

Hi everyone, I have a certain number of samples and I want to visualize the groups those samples belong to. For example, suppose to have three variables, age, sex, and smoker/nonsmoker, and three samples, S1, S2, S3. S1 is 35, male, nonsmoker S2 is 24, female, nonsmoker S3 is 24, female, smoker at the end I have the following data frame: S1 S2 S3 age 35 24 30 sex M F F smk N N S What I

I''m using" superredcloth (1.160) Using ''_.'' to indicate table headers doesn''t work: $ irb irb(main):001:0> require ''superredcloth'' => true irb(main):005:0> SuperRedCloth.new("|_. name |_. age |_. sex |\n| joan | 24 | f |\n| archie | 29 | m |\n| bella | 45 | f |\n").to_html =>

Hi, can anybody tell me how to contact the maintainer of the "concord" package? The address given in the help file is not valid anymore. ,---- | > help(package=concord) | | Information f?r Paket 'concord' | | Description: | | Package: concord | Version: 1.4-2 | Date: 2005-05-10 | Title: Concordance and reliability | Author: Jim Lemon

It is frustrating to see the labels I want in the dimensions of a list but not be able to extract those labels into titles for plots generated from component objects. If someone could set me straight, I would appreciate it. For your amusement, I have provided an example of the Byzantine code I am currently using to avoid loops: # Simulate ANOVA type test data sex<-c(rep(1,8),rep(0,8))

Hello, I have written code to plot an xyplot as follows: library(lattice) xyplot(len~ageJan1|as.factor(cohort),groups=sex,as.table=T,strip=strip.c ustom(bg='white',fg='white'),data=dat, xlab="Age (January 1st)",ylab="Length (cm)",main="Linear models for male and female cod, by cohort",type='p',

How do you label Xtick and Ytick marks with Currency symbols: $2000 instead of 2000? I would like to add dollar symbols to tickmarks on boxplots, histograms and back-to-back histograms. My Examples (requiring the lattice and Hmisc packages): data(case0102, package="Sleuth2") str(case0102) boxplot(Salary~Sex, case0102) histogram(~ Salary | Sex, data=case0102) require(Hmisc) #

Dear R users, my dataframe looks like this head(dat) Id sex byear age 1 300 m 2003 50 2 300 m 2003 36 3 402 f 2003 29 4 402 f 2003 21 5 402 f 2003 64 6 150 m 2005 43 ... ...(where Id is just the Identification number of Individual, sex (male or female), byear (=birthyear)) now, I 'd like to add a column, where each Individual gets an automated number starting

Hello all, I am working on a longitudinal study of children in the UK and trying the PAN package for imputation of missing data, since it fulfils the critical criteria of taking into account individual subject trend over time as well as population trend over time. In order to validate the procedure I have started by deleting some known values ?we have 6 annual measures of height on 300 children