similar to: dev.ask

I am trying to create a set of vectors which have similar names. They have a common part and an extra unique part. I am sure I have seen questions on this before but I have been unable to find it in the archive. There are examples of sets of names for columns of a data frame, but not for individual vectors. Suppose I want vectors named myvar.1 myvar.2 etc. I tried: >

I have a problem with expressions. I am trying to create a title where the parameter of interest is displayed as a Greek character. Which parameter is being considered is stored in a character variable. As an example, if I have param <- "alpha" and then do plot(0, 0, main = bquote(Parameter==.(param))) then in the title I get "Parameter = alpha", whereas I want the

I have been struggling trying to write some code to produce all combinations subject to some restrictions. I thought someone might have some bright ideas. I have 11 values which fall into 5 groups. I want all combinations of 2,3, and 4 values where each value must be from a different group. The numbers in the groups are different. Here is a definition of the groups: groups <- list(gp1 =

I have used xlsReadWrite to read data from an Excel spreadsheet. I had a problem with converting times of the day so that I could create POSIXct date-time objects. I was wondering if there was a better solution. Excel stores times of the day as fractions of a day so I wrote a function to convert the fraction to a number of seconds, extract the hours, minutes and seconds and output it in ISO

I am investigating some approaches to reproducible research. I need in the end to produce .html or .doc or .docx. I have used hwriter in the past but have had some problems with verbatim output from R. Tables are also not particularly convenient. I am interested in R2HTML and R2wd in particular, and possibly odfWeave. Does anyone have sample documents using any of these approaches which

I have a client running Microsoft SQL Server. I am interested in ways of accessing data from this server using R. I would welcome any information about how this can be done. I have a reasonable grasp of SQL and have experience with MySQL and RODBC but don't know anything much about Microsoft SQL Server. David Scott -- _________________________________________________________________

I have been playing around with producing tables using xtable and the type = "html" argument when printing. For example, if xtbl is the output of a dataframe which has been run through xtable, using the command: print(xtbl, type = "html", html.table.attributes = "border = '1', align = 'center'") I would be interested to see other examples of

I have a client who uses Stata 11. Can anyone advise me on ways of transferring data from this version of Stata to R? Reading the documentation for read.dta in package foreign it seems it only allows for versions up to Stata 10. I don't know Stata at all but a bit of poking around on their website led me to believe Stata 11 couldn't write in Stata 10 format. However Stata 11 can

I am writing a package for a company for its internal use only. What is an appropriate license statement for the DESCRIPTION file? I would like a statement which reflects the private and proprietary nature of the package, giving copyright to the writer and the company. I also don't want to violate the licensing of R and the packages I am using (RODBC, ggplot2, zoo). David Scott --

I would like to add some vertical lines to a lattice plot of histograms. What I am after is a lattice version of abline(v = 1234). The lattice histogram plot is just: histogram( ~ LTSE | approach, data = arrivals) Can anyone point me in the right direction for this? David Scott _________________________________________________________________ David Scott Department of Statistics The

I am looking for a method of dealing with the modified Bessel function K_\nu(x) for large \nu. The besselK function implementation of this allows for dealing with large values of x by allowing for exponential scaling, but there is no facility for dealing with large \nu. What would work for me would be an lbesselK function in the manner of lgamma which returned the log of K_\nu(x) for large

I am having trouble with creating a POSIXct object. If I create a variable of class Date first out of the date part of my data, I am ok, but if I just paste the date and time parts together and try and create the POSIXct object, I have problems. Here is a toy example created from the actual data which caused the problem. I am using R 2.0.1 on Windows XP. > # Data frame with dates and

Delving into the murky world of dates and times I found this: dates <- c("02/27/92", "02/27/92", "01/14/92", "02/28/92", "02/01/92") > times <- c("23:03:20", "22:29:56", "01:03:30", "18:21:03", "16:56:26") > x <- paste(dates, times) > z <- strptime(x, "%m/%d/%y

I am having trouble with try and tryCatch. I have read the documentation but I just don't get it. Here is what I am trying to do. I am testing a function which has a number of parameters. I want to test it for different values of the parameters and I have some loops, in the middle of which is a test of the function. Sometimes the routine fails and so I have put the bit that might fail

I want to do some recoding of variables: code Age into groups and recode a factor into a smaller number of levels. There are a couple of options for recode functions, in the car package and in memisc, and I think in gmisc. Does anyone have any opinions on the the easiest, most reliable approach for these problems? David Scott _________________________________________________________________

I am looking for a recommendation for an open source competitor to Winzip. I seem to recall Brian Ripley mentioning one in the last year or so, but couldn't find it in the mail archives. (Searching on Ripley there is somehow not terribly useful.) Suggestions? David Scott _________________________________________________________________ David Scott Department of Statistics, Tamaki Campus

Sorry to be completely off topic here, but I know there will be many list members who have an interest in and knowledge of hardware. We have for some time used Sun hardware for our departmental servers, most recently Sun Ultras of various specifications, running Solaris. We share some hardware with Mathematics as well. The applications we use are the usual suspects: R (obviously), S-PLUS, the

I have no experience with incorporating Fortran code and am probably doing something pretty stupid. I want to use the following Fortran subroutine (not written by me) in the file SSFcoef.f subroutine SSFcoef(nmax,nu,A,nrowA,ncolA) implicit double precision(a-h,o-z) implicit integer (i-n) integer l,i,nmax double precision nu,A(0:nmax,0:nmax) A(0,0) =

I have no experience with incorporating Fortran code and am probably doing something pretty stupid. I want to use the following Fortran subroutine (not written by me) in the file SSFcoef.f subroutine SSFcoef(nmax,nu,A,nrowA,ncolA) implicit double precision(a-h,o-z) implicit integer (i-n) integer l,i,nmax double precision nu,A(0:nmax,0:nmax) A(0,0) =

Hello everybody. I have a question that has stumped me and the usual "apply" tricks don't seem to work. I run a course where each student's performance is marked by one or more assessors. I have a data frame containing students' names, assessors' names and their marks, arranged as follows: ID student assessor Q1A Q1B Q1C Q2A Q2B Q3 1 2152833