similar to: Using variable names in for loops - Generating plots semi-automatically from a series of variables Partly solved

Hi, I have a valid SPSS .sav file (which I can open happily in SPSS v11 on Windows XP). Opening it in R2.41 on Linux we get this message :- > HSE3023 <- read.spss("HSE.sav") Error in read.spss("HSE.sav") : error reading system-file header In addition: Warning message: HSE.sav: Variable X234 indicates variable label of invalid length 256 Now variable X234 has indeed a

Dear colleagues, I have a problem which has bitten me occasionally. I often need to prepare graphs for many variables in a data set, but seldom for all. or for any large number of sequential or sequentially named variables. Often I need several graphs for different subsets of the dataset for a given variable. I run into similar problems with other needs besides graphing. What I would like to

Hi, I think I'm missing something very obvious, but I am missing it, so I would be very grateful for help. I'm using rpart to analyse data on skull base morphology, essentially predicting sex from one or several skull base measurements. The sex of the people whose skulls are being studied is known, and lives as a factor (M,F) in the data. I want to get back predictions of gender, and

Hi, As part of the analysis of a longitudinal study, I have a list, called wg, on R 1.3.1 for windows created as follows - wg <- by(xx, xx$mrn, function(xx) lm(xx$maternal.wt ~ xx$gestation, data=xx)) xx is a data.frame, and so wg contains the results of 293 linear regressions, as lm objects. Now what I want is to extract the 293 intercepts and regression coefficients as a vector, and

The Statistics & Decision Sciences community at Johnson & Johnson is searching for individuals to help in the mission of discovering, developing, manufacturing, and delivering beneficial medicines for patients. The group resides within the Janssen R&D organization in Johnson & Johnson's pharmaceutical segment. We need individuals with skills and motivation to consistently and

Hi, I'm analysing some anthropometric data on fifty odd skull bases. We know the gender of each skull, and we are trying to develop a predictor to identify the sex of unknown skulls. Rpart with cross-validation produces two models - one of which predicts gender for Males well, and Females poorly, and the other does the opposite (Females well, and Males poorly). In both cases the error

Position: Statistician in a systems biology team focused on metabolic disorders such as diabetes and obesity. The position is part of a statistics group supporting Pfizer global research and development. Role and Responsibilities: The statistician will be an integral member of a systems biology team which develops and/or uses computational and statistical approaches to manage and derive

Dear colleagues, This may be a question with a really obvious answer, but I can't find it. I have access to a large file with real medical record identifiers (mixed strings of characters and numbers) in it. These represent medical events for many thousands of people. It's important to be able to link events for the same people. It's much more important that the real record numbers

Dear colleagues, I would be very grateful for your help with the following. I have banged my head off this question several times in the past, and repeatedly over the last week. I have looked in the usual places and found no obvious solution. I fear that this just means I didn't recognize it, but I'd be very grateful for your help. I am scoring 8000 psychometric tests - the SCQ, if

Hi, I would be grateful for help with a problem which is irritating me. I am quite sure that I am doing something stupid, but I can't see what it is. I am running R 1.7 on Windows 2000. The graphics device is the PC screen. The graphics from the nlme demonstration in Bates an Pinheiro's manual work just as advertised. The CO2 data and the Orthodont data dsiplay

Dear Sirs. I am working with the R packages Hmisc and PredictABEL to make NRI estimates from my Cox models with and without a specific biomarker. According to Pencina et al (Statistics in Medicine 2010, DOI: 0.1002/sim.4085 ), a continuous/non-categorical NRI (NRI>0) is to be used when there are no obvious reason to categorize risk, such as the risk of future cardiovascular events in patients

Edwards Lifesciences (NYSE: EW) is a global leader in products and technologies to treat advanced cardiovascular disease, the global leader in acute hemodynamic monitoring and the number-one heart valve company in the world. Headquartered in Irvine, California, Edwards has more than 5,000 employees worldwide, selling medical technologies in more than 100 countries. Edwards Lifesciences leverages

Statistical Programmer (SAS or Splus/R) About the Company: Edwards Lifesciences (NYSE: EW) is a global leader in products and technologies to treat advanced cardiovascular disease, the global leader in acute hemodynamic monitoring and the number-one heart valve company in the world. Headquartered in Irvine, California, Edwards has more than 5,000 employees'' worldwide, selling medical

Thanks Andris, Michael and Petr for your prompt and kind feedbacks. I will try generating my own biplot from low-level graph commands... I hope it will work. Best regards, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) PO Box 616 6200 MD Maastricht E-mail: m.manca at path.unimaas.nl Office

About the Company: Edwards Lifesciences (NYSE: EW) is a global leader in products and technologies to treat advanced cardiovascular disease, the global leader in acute hemodynamic monitoring and the number-one heart valve company in the world. Headquartered in Irvine, California, Edwards has more than 5,000 employees worldwide, selling medical technologies in more than 100 countries. Edwards

Two Postdoctoral fellows (fulltime) Leiden University, Leiden, The Netherlands Job description ________________________________________ TOP Institute Pharma (TI Pharma) has granted our proposal to set up a mechanism-based PK-PD modelling platform. This platform focuses on the transfer of knowledge from academia to the pharmaceutical industry and is a collaborative effort of

Post-Doctoral Fellowship Modelling of Cardiovascular Safety (full time) Leiden University, Leiden, The Netherlands <javascript:doredirect();> Job description TOP Institute Pharma (TI Pharma) has granted our proposal to set up a mechanism-based PKPD modelling platform. This platform focuses on the transfer of knowledge from academia to the pharmaceutical industry and is a

Hello, I have a longitudinal data with about 30 subjects. I used xyplot() to plot the longitudinal data. One problem is that xyplot() recycles the color of auto.key so that every 7th subject has the same color (symbol if setps() was used). Is there a way so that every subject will have a unique color or symbol? Thanks Osman -- Osman O. Al-Radi, MD, MSc, FRCSC Fellow, Cardiovascular Surgery

Dear R-help, I would like to add text to each of four panels in a plot generated by xyplot in lattice library. A sample code is given below, the plot generated has the first label repeated in all panels! How can I get the labels to be different in each panel? library(lattice) x <- rnorm(400) y <- rnorm(400) a <- gl(4, 100) xyplot(y~x|a, panel=function(...){

Hello, I would like to create a plot composed of stacked boxes (squares or rectangles), where the size of the box would represent the frequency of observations based on a categorical variable (group), the color would represent the proportion of success (binary) within that group (outcome) on a predetermined color scale. Ideally the boxes can be stacked from the bottom left to the top right based