similar to: How to select specific rows from a data frame based on values

Dear all, I am estimating a mixed-model in Ubuntu Raring (13.04¸ amd64), with the code: fm0 <- lme(rt ~ run + group * stim * cond, random=list( subj=pdSymm(~ 1 + run), subj=pdSymm(~ 0 + stim)), data=mydat1) When I check the approximate variance-covariance matrix, I get: > fm0$apVar [1] "Non-positive definite

I've been running a glmer logit on a very large data set (600k obs). Running on a 10% subset works correctly, but for the complete data set, R completes apparently without error, but does not display the results. Given these jobs take about 200 hours, it's very hard to make progress by trial and error. I append the code and the sample and complete output. As is apparent, I upgraded R

Hi, I have sort of a newbie question. I've seriously put a lot of effort into how to handle simple replicates in a repeated ANOVA design, but haven't had much luck. I really liked reading "Notes on the use of R for psychology experiments and questionnaires", by Jonathan Baron and Yuelin Li ( http://www.psych.upenn.edu/~baron/rpsych/rpsych.html ) but still didn't run across

Hi I compouted a multiple linear regression with repeated measures on one explanatory variable: BOLD peak (blood oxygenation) as dependent variable, and as independent variables I have: -age.group (binaray:young(0)/old(1)) -and task-difficulty measured by means of the reaction-time 'rt'. For 'rt' I have repeated measurements, since each subject did 12 different tasks. -> so

Hello: I have a tab delim file with 100 rows and 100 columns. I have numerical values in this table. What I want is to create an image color map with color gradation. my values range from -5 to 0. max value is 0. to acheive this, I need to convert my data.frame into matrix. I tried following : mydf <- read.table('mytable',sep='\t',header=T) mydf1 <- mydf[,2:100]

Hi All: I have written an R search plug-in based on Prof Jonathan Baron's R search site. This Firefox/Netscape plug-in lets you search R sites and R discussion group archives directly from your Mozilla (Seamonkey) or Firefox or Netscape browser. Here is the link to the installer for the plug-in: http://arinbasu.bravehost.com/installer.html If you want the source files or have other comments

Hello! # I have a data frame: mydf<-data.frame(c1=rep(NA,5),c2=rep(NA,5),c3=rep(NA,5)) # I have an index whose length is always the same as nrow(mydf): myindex<-c(1,2,3,2,1) # I need c1 to have 1s in rows 1 and 5 (based on the information in myindex) # I need c2 to have 1s in rows 2 and 4 (also based on myindex) # I need c3 to have 1 in row 3 # In other words, I am trying to achieve this

Hi All: Is a distribution of R version 2.3.1 (the current version) available for the Ubuntu Linux? If I want to build a .deb package for R-2.3.1 for Ubuntu, are there specific requirements that I should be aware of? Is a tutorial available? I searched the web, and the repositories, but oculd not locate an appropriate file. TIA, Arin Basu

Hi All, I am scheduled to teach a graduate course on research methods in health sciences at a university. While drafting the course proposal, I decided to include a brief introduction to R, primarily with an objective to enable the students to do data analysis using R. It is expected that enrolled students of this course have all at least a formal first level introduction to quantitative methods

Hi All: This follows my earlier post on webized slides on lecture presentation on introducing R. I learned that in Mozilla (Firefox) browsers, the slides did not show up. Sorry for the no show. As a reluctant windows user, I kind of carelessly clicked through Powerpoint to convert the presentation file to its html form, unwittingly leading to the mess. See if it got corrected now (I do not have

Hi All: I came across the following problem while working with a dataset, and wondered if there could be a solution I sought here. My dataset consists of information on 402 individuals with the followng five variables (age,sex, status = a binary variable with levels "case" or "control", mma, dma). During data check, I found that in the raw data, the data entry

Hi All: I tried to update the package mgcv from my current version 1.1-8 to mgcv version 1.2. Using the command update.packages("mgcv") to update mgcv failed in my computer and the returning message was: "Error in old.packages(lib.loc=lib.loc, contriburl=contriburl, method=method, :no installed.packages for (?invalid) lib.loc=mgcv" When I typed installed.packages(), I got

Hi All: A couple of weeks back, I asked a question on the list that I was invited to provide an introductory lecture on R to a group of academicians in Kolkata. I thank all of you who had generously guided me in providing me web links and words to the wise. Time to give back. I did the presentation on introduction to R and uploaded the presentation files at the following site:

Dear All, Folks at the West Bengal University of Technology has set up a mirror to distribute R and associated packages. Here is the URL: http://mirror.wbut.ac.in/CRAN This will be helpful for R users in South Asia and the adjoining countries. I have yet to see the site featured in the CRAN list of mirrors. What needs to be done? For those of you who are impatient and would like to give it a

Hi All: I was looking for, but could not locate in the packages, or in the R archive searches if there exists an R implementation of Cuzick's test of trend. The test is described as follows: An extension of the Wilcoxon rank-sum test is developed to handle the situation in which a variable is measured for individuals in three or more (ordered) groups and a non-parametric test for trend

Hi Group: >From an R session, I wanted to update packages, and issued the command: update.packages() R started to update the package "foreign" and proceeded. After downloading, while it attempted to install the packages, it came up with the following error message. I could not understand the error message. Does it indicate that the problem lies with the gcc? How do I install the

Hi All: Greetings, and best wishes from the festive times here at Kolkata, India -- the time of Durga Puja celebrations. I seek your advice as I try plotting lines for my data. The problem: I have created a dataframe that looks like this (name: myFrame): lowest second third fourth highest significance INAS 0.107 0.115 0.123 0.115 0.166 0.000 MMA 0.091 0.107 0.115

Dear List: This is in search of advices/opinions for a lecture presentation on R. I have this presentation scheduled on Saturday (16-October-2004); the topic is "introduction to R". The audience is a group of researchers and scholars who have either never used R/are beginning to use R, but otherwise have used a package called limdeb, and SPSS. There is also a request to cover Rgnumeric,

Hi All: Recently, I got a query how R compared cost-wise (cost of purchase, maintenance, and training manuals) with other proprietary statistical software programs (or programming environments). The request, verbatim: <--- begin quote ----> >> dear arin, >> >> can you please find out the price of the various >> prop. statistical >> packages that are in

Hi All: I am interested in tools available in R that will enable me to run analysis of case crossover studies. I searched the web (http://www.google.com), and found that a tool called "funfit" was available for Splus. Searching the R archives, I also found mention of a "funfit" library. However, when I searched in CRAN, I could not find it. Is "funfit" now part of