similar to: THANK YOU: Updating R version

I can't seem to find just what I'm looking for in R help, Everitt and Hothorn HSAUR, Murrell's book, or the R graphics gallery at http://addictedtor.free.fr/graphiques/. Probably not looking efficiently, but anyway, If my data look like this: > head(data) cat startyear studentid 1 other 2001 12 2 UHS 2001 17 3 Lourdes 2001 10 4

R 2.5.1 on WinXP I'm trying to create new variables in a dataframe called jaw, as powers of jaw$age up to the sixth power, and name them jaw$age.(--the digit corresponding to the power--) Obviously none of these work (silly for me to try, I suppose): for (i in 2:6) { jaw$age.'i' <- jaw$age^i } for (i in 2:6) { jaw$age.i <- jaw$age^i } for (i in 2:6) {

Suppose I have a Master.Rnw file that looks something like this: \documentclass[12pt]{mypaper} \usepackage[margin=1in]{geometry} \usepackage{setspace} \usepackage{url} \usepackage{indentfirst} \usepackage{fancyhdr} \usepackage{Sweave} \pagestyle{fancy} \lhead{sonographic rectal diameter and ADHD} \rhead{ } \usepackage{abbrevs} %\usepackage{natbib} %\usepackage{apacite}

I start with a dataframe called xrays. It contains scores on films from each of two radiologists. It is in "long" format. I used the reshape package to melt a data frame and then cast it into "wide" format, one line for each patient (identified by redlognumb) with scores from both radiologists for a given patient on the same line. I named the result of the casting xrays.data.

Is there a way to automate fitting and assessing loglinear models for several nominal variables . . . something akin to step or drop1 or add1 for linear or logistic regression? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu "If you want to build a ship, don't drum

I'm trying to make a factor with levels for day of week, appear as Monday-Sunday, rather than alphabetically. I'm using reorder.factor from gdata package, but obviously not using it properly. What am I doing wrong? Thanks. >data <- read.table("DataEarly2008.txt", sep=" ", header=T, as.is=1:2) >str(data) 'data.frame': 7953 obs. of 9 variables: $

Running R 2.5.1 on WinXP with Tinn-R, and using Sweave via MikTex. I am having trouble with the code below: ------------------------------------------------------- <<lowerCI95>>= lowerCI95713 <- t.test(data713, alternative="greater") lowerCI95713.bound <- lowerCI95713$conf.int[1] @ Thus a lower 95\% confidence bound on the mean octane concentration is

I've got a 20 slide Powerpoint set up against the expected day someone says "Mitch, could you talk for a short while to the residents on how to use statistics?" It looks good at this end, but I'd love to get at least one more pair of eyes on it. Any volunteers? -- Due to the recession, requests for instant gratification will be deferred until arrears in scheduled

Dear all: Just to make sure I am on the same page as everyone else. I just got a new computer and had to reinstall MikTex and of course R and Tinn-R. In order to be able to generate PDF's with Sweave I had to install notepad++ to tell Sweave where Latex.exe and yap.exe where located under MikTex path. So,I have two questions: 1- Do I need to have a Miktex editor in order to use Sweave? I am

Dear all, I like the way the Rcmdr package computes reliability. E.g reliability(cov(d[,c("q1", "q2", "q3", "q4", "q5", "q6")], use="complete.obs")) will not only give me the alpha score, but also for each variable, alpha.score if deleted. However - when writing scripts it's very tiresome to load the whole Rcmdr GUI just

I am wondering if there is a better alternative than Excel for data storage that does not require database knowledge (I will eventually have to learn this, but it is not on my immediate todo list). I need something that is not limited to 256 columns... I don't need any of the built in functions in excel just a spreadsheet like program with cells that hold data in a data.frame format for a

I can't seem to use the point-and-click identify() function properly. I'm running R 2.5.1 (I know, I need to get around to upgrading) under Win XP. The problem is, when I click on a point on the graph, I get an error, "no point within 0.25 inches." But in some areas, I can click where there is no visible point anywhere close, and an identify() label will pop up. The troublesome

Hello, I'm looking for a textbook that can explain some of the math behind the intro-to-intermediate stuff like ANOVA, multiple regression, non- parametric tests, etc. A little background: I took an intro stats course last year and would like to further my education. Being as that was the highest (and only) stats class the local community college offers, it looks like I'm on

Every day I get a csv file containing the names of the 64 schools in our county, the number of students sent home ill, and the number of students absent (plus lots of other variables). The file is cumulative since fall of 2009. It is in "long" format: one line per school per day. Each line is also supposed to contain the total number of students enrolled in the school. That number

A quick question. I am trying to understand how I could move the installed packages in my R 2.3 version to the newly installed R 2.5 version, without having to install all the packages again. I copied the files under the old library subdirectory to the new library subdirectory. But still the newer version is not recognizing the packages that were copied over. Thanks. [[alternative HTML version

I'm trying to learn Sweave. So far things are going well with the chunks of code identified by << >>= But I'm having trouble with the in-line text use of \Sexpr. Here is a short example .Rnw file: \documentclass[12pt]{article} \usepackage[margin=1.25in]{geometry} \usepackage{graphicx} \usepackage{Sweave} \begin{document} <<test>>= pnorm(1) @ \Sexpr{2+6}

suppose I have some logical vector x <- as.logical(c(0,0,0,1,0,0,1,1,0)) x How would I make the words TRUE appear on the screen in a different color from the words FALSE? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu "If you want to build a ship, don't drum

Dear all, Please, please I really need your help, since I am not able to solve the problem myself: It seems that until now most Windows users of my Bioconductor package "xps" use Windows XP, but recently a user tried to use it on Windows Vista where it crashed R immediately. Since I assume that some users may migrate to the upcoming Windows 7 where the problem will be the same, I

Running R 2.13.1 on Windows XP. I would like to get week of the year (1-52) for each date. library(chron) dts <- dates(c("02/27/92", "02/27/92", "01/14/92","02/28/92", "02/01/92")) dts dts.chron <- as.chron(dts) dts.chron class(dts.chron) # all of these component extractions work: months(dts.chron) weekdays(dts.chron) years(dts.chron)

Hello I use the function .libPaths() > .libPaths() [1] "/usr/local/lib/R/library" > .libPaths(c("/usr/local/lib/myRlib","/usr/local/lib/library")) > .libPaths() [1] "/usr/local/lib/R/library" but it simply does not pick up the new path. Can anybody help? (The problem originates from the inability to install packages into