similar to: Fortran vs R

Hello,   I am currently working on a modeling project using Fortran to run large repetitive loops (many DO loops). As part of this process I would like to use a model fit in R and currently stored as an R object. This is a rather complex model, a GAMM, and I am curious whether there is a way to call this model from Fortran. I am not sure "call" is correct terminology, but I

Hello! Thanks for reading this request for assistance. I have a question regarding creating a histogram-like figure from data that are not currently in the correct format for the "hist" command. Specifically, my data have been processed and are in a matrix with columns containing the variables of interest and separate columns containing the number of times this variable was

Hello, I'm trying to install a package that requires a Fortran compiler (Hmisc) using R CMD INSTALL. I downloaded the package source onto my Desktop, unzipped it, and then typed: R CMD INSTALL /Users/brianbeckage/Desktop/Hmisc * Installing *source* package 'Hmisc' ... ** libs g77 -fno-common -g -O2 -c cidxcn.f -o cidxcn.o g77 -fno-common -g -O2 -c cidxcp.f -o cidxcp.o g77

Greetings! I have one quick question: How do you do a Tukey test on an ANCOVA? Thanks for any tips! -Paul

Hello, I have a question regarding how to subset/select parts of a data frame (matrix) in order to plot data associated only with this subset. Specifically I have a large data frame in which one column contains ID values (dates), and other columns contain data I would like to plot (temperature, light, etc.). I would like to break up this large matrix so as to plot data associated

Hi, I am doing a very informal presentation for my office about R capabilities to deal with and analyze spatial data, display data and maps, and connections with GIS. I've used in my presentation info from the CRAN, the spatial Task view, and the more striking graphics examples from http://addictedtor.free.fr/graphiques/thumbs.php and NCEAS

Dear R-ers: Recognizing that it is better to remain silent and have people think you stupid than to speak and remove any doubt, I fear that I must speak nonetheless. I am having the Devil's own time trying to make a simple R package for my own use. Yes, I have perused the FAQs and searched the archives. Yes, I have read "Writing R Extensions" and eventually divined what Rcmd

Hello, I am writing regarding log transformation of data in a single matrix column, and subsequent use of these data in a glm model fit. I have a data matrix in which I am using the log function to transform the values. This transformation results in -Inf values in some places, though. I then receive an error when this matrix is used in the glm function, and would like to know this

Greetings R-ians: I know what doesn?t work but I don?t know why, nor how to remedy things. I have a character string containing "." which I want to replace with " " gsub(".", " ", file.label) replaces the every character with a blank. However gsub(".xls", " ", file.label) replaces ".xls" with a blank as expected. It appears

Greetings, R-ians: I would like to know which version or R was used to create a given package. I think I remember seeing that topic discussed recently but cannot find it among my notes. Can anyone tell me how to determine which version of R created a package? Thanks. Charles Annis, P.E. <mailto:Charles.Annis@StatisticalEngineering.com> Charles.Annis@StatisticalEngineering.com

R-ians: After some effort I coerced my code to do what I want but my syntax is a kludge. Suggestions on more elegant syntax? par <- NIM.results$par do.call("Draw.NIM.POD.curve", list(par[1], par[2], par[3], par[4], par[5], a.hat.decision, .... et cetera ... It seems that I should be able to avoid defining the variable "par" and then specifying each of its elements,

Greetings R-ians: Searching the Searchable Mail Archives I discovered that ctrl L will clear the Rgui screen, which is what I'd like to do from a print (or some similar) statement. Is there a mechanism to use the ctrl L clear-screen sequence in a script, or print statement? Thanks for your counsel. Charles Annis, P.E. Charles.Annis at StatisticalEngineering.com phone: 561-352-9699 eFax:?

Greetings, R-ians: I'm using R 2.3.1 on WindowsXP. I need to find the name of a file at the end of a sting that contains the path + file, with the problematic "\\" as separators. The string looks something like this: "C:\\Documents and Settings\\myName\\My Documents\\R Projects\\Project1\\file.name.csv" What I want is "file.name.csv" Currently I use the

Greetings, R-ians: I'm using R 2.3.1 on WindowsXP. I need to find the name of a file at the end of a sting that contains the path + file, with the problematic "\\" as separators. The string looks something like this: "C:\\Documents and Settings\\myName\\My Documents\\R Projects\\Project1\\file.name.csv" What I want is "file.name.csv" Currently I use the

Dear R Helpers: I am trying to annotate a plot. The following code snippet works, but it is kind of a kludge since it adds the partial derivative symbols after creating the plotmath frac(). Is there a more elegant way to write a partial derivative? dev.off() plot(NA, xlim=c(-3,3), ylim=c(0,1.6), xlab="", ylab="", tck=-0.015) text(1.6, 1, expression(paste("slope =

Dear R-ians: I'm running R2.9.2 on a 6 year old Windows XP DELL with 2 Gig RAM and a 3MHz Pentium 4 chip. I've written a package using the User Menus Under Windows commands (winMenuAdd, etc). It works very well. I have 6 test cases and running any one of them requires many selections form the menus and some keyboard entry, and it takes me hours to exercise them all to see

Greetings R-wizards: For historical reasons I have filenames with the character "?" and have successfully used "\u00e2" in its place, with the hoped-for result on all my on-screen plots. However since R2.7.0 I have trouble with savePlot() when the file name includes that character as it does in this example: savePlot(paste("diagnostic ? vs a ", file.label,

Greetings, Amigos: I have been trying without success to convert a character string, > repeated.measures.columns [1] "3,6,10" into c(3,6,10) for subsequent use. as.numeric(repeated.measures.columns) doesn't work (likely because of the commas) [1] NA Warning message: NAs introduced by coercion I've tried many things including strsplit(repeated.measures.columns, split =

Greetings: I'm running R2.6.2 on a WinXP DELL box with 2 gig RAM. I have created a new glm link function to be used with family = binomial. The function works (although any suggested improvements would be welcome), logit.FC <- function(POD.floor = 0, POD.ceiling =1) { if (POD.floor < 0 | POD.floor > 1) stop ("POD.floor must be between zero and one.") if

I?ve been using this routine for several years. I?m sorry, I don?t remember where I got it. It works as it should, viz. it blanks the R console. But it requires package rcom and now that requires rscproxy. cls <- function () { require(rcom) wsh <- comCreateObject("Wscript.Shell") comInvoke(wsh, "SendKeys", "\f") invisible(wsh) } > cls()