Re: [R] signal handling

On Fri, 12 Aug 2005, Paul Gilbert wrote: > Omar Lakkis wrote: > > > Is ther a signal handling model in R? similar to Perl's > > %SIG hash. I want to do fast clean up in my R code before > > exit when a kill signal is issued. > > I'm not sure about perl's signals, but Unix signals can be > passed w

Re: [R] boot error: Error in statistic(data, original, ...) : unused argument(s) ( ...)

Please look at the examples and the book which package boot supports (see its DESCRIPTION file). From ?boot statistic: A function which when applied to data returns a vector containing the statistic(s) of interest. When 'sim="parametric"', the first argument to 'statistic'

Re: [R] g77 and R

> [EMAIL PROTECTED] wrote: >> Jerry, >> >> Here are the two R packages. >> >> Thanks. >> >> Rick > All the fortran codes compile OK so the problem is invoking the compiler > you > want. One approach which is a work around is to rename the gfortran > executable > to somthing else like gfortran-save

[R] boot error: Error in statistic(data, original, ...) : unused argument(s) ( ...)

I have added on additional line to the code below, the line the defines the matrix passed to AdjForBase2. I am having difficulty setting up a boot strap. My code is listed below. I am getting the following error message: Error in statistic(data, original, ...) : unused argument(s) ( ...) I neither

[R] boot error: Error in statistic(data, original, ...) : unused argument(s) ( ...)

I am having difficulty setting up a boot strap. My code is listed below. I am getting the following error message: Error in statistic(data, original, ...) : unused argument(s) ( ...) I neither understand the error I have made that leads to the error message, nor do I know how to fix my error. Any

Re: [R] g77 and R

> [EMAIL PROTECTED] wrote: >> Jerry, >> >> Here are the two R packages. >> >> Thanks. >> >> Rick > All the fortran codes compile OK so the problem is invoking the compiler > you > want. One approach which is a work around is to rename the gfortran > executable > to somthing else like gfortran-save

[R] Broken tkrplot on Mac OS X

Dear all-- I have been trying to get the tkrplot package to work in order to use the bioconductor package genArise. I am trying to build it on a Mac running OS 10.4.1 with R 2.1.1. Following a stanadard install from R here the error I get when I try to load tkrplot > library(tkrplot) Loading re

Re: [R] question

What are RGN$G and RGN$R? Consider the following: > .Machine$double.xmax [1] 1.797693e+308 > log(.Machine$double.xmax, base=2) [1] 1024 > 2^(log(.Machine$double.xmax, base=2)-1) [1] 8.988466e+307 > 2^log(.Machine$double.xmax, base=2) [1] Inf If you want somethin

Re: [R] Problems with numeric variable containing ? values

I'm very sorry, but I do not understand your question. If you are still waiting for an answer, may I suggest you PLEASE do read the posting guide! "http://www.R-project.org/posting-guide.html";. This may help you formulate your question so it will be easier for someone else to under

[R] R/S-Plus/SAS yield different results for Kendall-tau and Spearman nonparametric regression

Colleagues, I ran some nonparametric regressions in R (run in RedHat Linux), then a colleague repeated the analyses in SAS. When we obtained different results, I tested S-Plus (same Linux box). And, got yet different results. I replicated the results with a small dataset: DATA: 37.5 23

Re: [R] R-devel and bad magic number

Sean Davis wrote: > Is anyone else having issues like this (r-devel)? The files look fine in an > editor > > >>load('stuff1.R') > > Error: bad restore file magic number (file may be corrupted) -- no data > loaded > >>load('/Users/sdavis/Desktop/Downloads/biomaRt/R/biomaRt.R') > > Error: b

Re: [R] Manually Calculating Odds from POLR Model

Dear Tate, If I understand correctly what you're asking, the formulas are on p. 21 of the paper at . But why do you want to do this when you can get the fitted probabilities from predict()? I hope this helps. John

Re: [R] quotes

col.names <- strsplit("one, two, three", ", ")[[1]] is an easy way to do this. But you can do paste('"', gsub(', ', '","', "one, two, three"), '"', sep="") R accepts either single or double quotes and it can make life easier to use both. On Fri, 12 Aug 2005, Jean-Pierre Gattuso wrote: > Hi:

Re: [R] multiple responses classification/regression problem

polymars() in polspline pakage can do the work. HTH. On 8/12/05, Weiwei Shi <[EMAIL PROTECTED]> wrote: > Hi, there: > I am wondering if anyone knows about this? the response variable is a > vector. I knew mvpart might be able to do this. anyone would like to > share some examples? > > of course,

Re: [R] need help

Hi, there: here is some part from my previous email: [,1] [,2] [,3] [,4] [,5] [1,] 34.216166 96.928587 330.125990 330.183222 330.201215 [2,] 2.819183 8.134491 8.275841 8.525256 8.828448 [3,] 2.819183 7.541680 7.550333 8.374636 8.690998 [4,] 4.67255

Re: [R] Dating Objects

Barry Rowlingson wrote: > Rolf Turner wrote: > > >> I don't follow what you aren't following It seems to me >> to be an eminently reasonable thing to want to do. Answers >> questions such as ``Was x modified since y was modified?'' > > > I wanted to have some sort of

Re: [R] need help

> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > On Behalf Of Weiwei Shi > Sent: Friday, August 12, 2005 2:05 PM > To: r-help > Subject: [R] need help > > Hi, there: > I think i need to re-phrase my question since last time I did not get > any reply but i think th

Re: [R] signal handling

I'm not sure about perl's signals, but Unix signals can be passed with q("yes/no", status=whatever) See ?q. This is pretty useful for passing signal to make, for example. Paul Omar Lakkis wrote: > Is ther a signal handling model in R? similar to Perl's %SIG hash. > I want to do fast clean

[R] multiple responses classification/regression problem

Hi, there: I am wondering if anyone knows about this? the response variable is a vector. I knew mvpart might be able to do this. anyone would like to share some examples? of course, nnet can do that too, but what else? Thanks, you guys have a good weekend! Weiwei -- Weiwei Shi, Ph.D "Did you

[R] quotes

Hi: I have been struggling with gsub to no avail and am seeking help from the list. I want to make a vector of the string "one, two, three" and have tried to: - replace each comma by "," - add " at the beginning of the string - add " at the end of the string to issue the comm

[R] need help

Hi, there: I think i need to re-phrase my question since last time I did not get any reply but i think the question is not that hard, probably i did not make the question clear: I want to find cases like 35, 90, 330, 330, 335 from the rest which look like 3, 3, 3, 3.2, 3.3 4, 4.4, 4.5, 4.6, 4.7 .

[R] converting a t statistic to r2

HI I wonder if anyone can help. I have a longitudinal sample of 100 subjects: 200 data points were acquired starting at different ages and at irregular intervals (subjects have different numbers of repeated data points, so some have only one data point). I have been examining the relationship o

[R] R-devel and bad magic number

Is anyone else having issues like this (r-devel)? The files look fine in an editor > load('stuff1.R') Error: bad restore file magic number (file may be corrupted) -- no data loaded > load('/Users/sdavis/Desktop/Downloads/biomaRt/R/biomaRt.R') Error: bad restore file magic number (file may be

[R] Manually Calculating Odds from POLR Model

Hello, I am using polr(...) to generate a model. The summary shows the coefficients and the intercepts. For example: coefficient for x1 = c1 coefficient for x2 = c2 intercept A|B = i1 intercept B|C = i2 I can then run predict(..., type="p") with the model and see the odds for

Re: [R] as.formula and lme ( Fixed effects: Error in as.vector(x, "list") : cannot coerce to vector)

You can use substitute() to get the actual formula in your call. demo.lm2 <- eval(substitute(lme(fma, data=demo, random=~1|subj), list(fma=fma))) works for me in your example. On Fri, 12 Aug 2005, Liu Song wrote: > This is a continuing issue with the one on the list a long time ago (I > couldn

Re: [R] coercing created variables into a new data frame using na.omit()

On Fri, 12 Aug 2005, Kjetil Brinchmann Halvorsen wrote: > Prof Brian Ripley wrote: > >> I don't know if you can read your message, but I find it exceedingly >> difficult and there seem to be several typos. Please use the space and >> return keys ... and only send a message once. >> >> You prob

[R] as.formula and lme ( Fixed effects: Error in as.vector(x, "list") : cannot coerce to vector)

This is a continuing issue with the one on the list a long time ago (I couldn't find a solution to it from the web): -- > Using a formula converted with as.formula with lme leads > to an error message. Same works ok with lm, a

Re: [R] HPD credible sets

There are 2 very basic functions in the package TeachingDemos. They only work on univariate unimodal (single interval) situations. There is also a function in the boa package, but I don't know anything about that one. Greg Snow, Ph.D. Statistical Data Center, LDS Hospital Intermountain Health Ca

Re: [R] coercing created variables into a new data frame using na.omit()

Prof Brian Ripley wrote: >I don't know if you can read your message, but I find it exceedingly >difficult and there seem to be several typos. Please use the space and >return keys ... and only send a message once. > >You problem is perhaps that you are not looking at the data frame, but at >th

Re: [R] Issues with tcltk for Tiger Mac OSX

Dear Bernard, The Rcmdr web site has some instructions prepared by Rob Goedman for Mac users at that might prove useful. (I'm not a Mac user myself so can't offer specific advice.) I hope this helps, John --

Re: [R] How to insert a certain model in SVM regarding to fixed kernels

> David, Please correct me if I am wrong but I think svm partially works > with dyn although I don't remember what the specific limitations were. Yes, the fitted values / residuals can be extracted from the trained model. The 'newdata' argument of predict() is not functional yet for time series.

Re: [R] Dating Objects

Rolf Turner wrote: > I don't follow what you aren't following It seems to me > to be an eminently reasonable thing to want to do. Answers > questions such as ``Was x modified since y was modified?'' I wanted to have some sort of 'make'-like functionality in R. Suppose

[R] Bug in order with na.last=FALSE?

I scanned the archives a bit and found discussion of a bug in order when na.last=NA, but nothing about when na.last=FALSE. Assign a (character) vector x as follows: x <- c("2005-08-12 12:38:35","2005-08-12 12:38:46",NA, "2005-08-12 12:38:27",NA) Then executing o <

Re: [R] Dating Objects

Paul Roebuck wrote: On Fri, 12 Aug 2005, Jason Skelton wrote: > > I'm trying to carryout the equivalent of ls -l in R > > to give some date/time label to each of my objects > > > > If the case is that there is no equivalent, is it > > possible to list all objects in an environment that > > share

Re: [R] Dating Objects

Hi Jason, Perhaps you get do what you want to do using attributes? Does this get close to what you need? > ## Create some objects > xx <- "adsfasdf" > zz <- 3027 > yy <- "abdsf" > > ## Assign a time attribute to them. > attr(zz, "time") <- Sys.time() > attr(yy, "time") <- Sys.time() > > ## To se

Re: [R] Concerning reading of SAS-files

bogdan romocea wrote: > The first one is an index, not a data set. Anyway, just use SAS to > export the data sets in text format (CSV, tab-delimited etc). You can > then easily read those in R. (By the way, the help for read.xport says > that 'The file must be in SAS XPORT format.' Is .sas7bdat an

Re: [R] evaluating string variables

Given my own language inadequacies (I speak only one), I must, of necessity, refrain from criticizing the multilingual. However, in this one instance my sensibilities have been overwhelmed: "symmetric" , please :-) . Incidentally, you may find that file.path() is also useful. Cheers, Bert --

Re: [R] evaluating string variables

Thanks a lot!!! using m<-matrix(scan(paste("C:/R/",i,".mat",sep=""), n=5*5),5,5, byrow=TRUE) worked very well!!! Leonardo __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://

Re: [R] evaluating string variables

Le 12.08.2005 16:44, Leonardo Sepulveda Durán a écrit : >Hello!!! > >I have a folder (C:/R/) with matrix files, named by number, i.e. >0.mat, 1.mat,...,1250.mat. In this case, they are 5x5 simetric >matrices. I would like to compute a property for each matrix and put >calculated values into a data

Re: [R] evaluating string variables

look at ?paste(), e.g., paste("C:/R/", i, "mat", sep = "") I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/16/336899 Fax: +32/16/337015 Web: h

Re: [R] evaluating string variables

use paste() to construct the file name e.g. fn <- paste(i,".mat",sep="") Sean On 12/08/05, Leonardo Sepulveda Durán <[EMAIL PROTECTED]> wrote: > Hello!!! > > I have a folder (C:/R/) with matrix files, named by number, i.e. > 0.mat, 1.mat,...,1250.mat. In this case, they are 5x5 simetric > matri

[R] evaluating string variables

Hello!!! I have a folder (C:/R/) with matrix files, named by number, i.e. 0.mat, 1.mat,...,1250.mat. In this case, they are 5x5 simetric matrices. I would like to compute a property for each matrix and put calculated values into a data frame for posterior ploting and printing. Below there is an ex

Re: [R] Dating Objects

On Fri, 12 Aug 2005, Jason Skelton wrote: > I'm trying to carryout the equivalent of ls -l in R > to give some date/time label to each of my objects > > If the case is that there is no equivalent, is it > possible to list all objects in an environment that > share a common component so that the co

Re: [R] Creating new columns inside a loop

Thanks, everyone! I knew it had to be simple, just didn't know the magic word to search on. I'm off! On 8/11/05, Berton Gunter <[EMAIL PROTECTED]> wrote: > > ?assign > > -- Bert Gunter > Genentech Non-Clinical Statistics > South San Francisco, CA > > "The business of the statistician is to cat

[R] Ecartis command results: -- Binary/unsupported file stripped by Ecartis --

>> The original message was received at Fri, 12 Aug 2005 09:18:43 -0500 Unknown command. >> from 175.107.114.173 Unknown command. >> - The following addresses had permanent fatal errors - Unknown command. >> [EMAIL PROTECTED] Unknown command. >> - Transcript of the session follows

[R] General expression of a unitary matrix

Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, s

Re: [R] Groups in histograms?

On 8/9/05, Fredrik Karlsson <[EMAIL PROTECTED]> wrote: > Dear list, > > Sorry for answering my own post, but I have had partial sucess in > this. With the panel function below, I get separate histograms in each > panel using the group argument. > > histogram(~vot | agem, nint=50,data=work,groups=

[R] Issues with tcltk for Tiger Mac OSX

Dear R-helpers, I have installed the latest version of R 2.1.1 in the Tiger OSX. However, when I load up R and use the following command: > library(tcltk) I encounter the following error: Loading Tcl/Tk interface ... Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to l

[R] HowTo derive a correct likelihood-ratio chi-squared statistic from lrm() with a rsc() ?

Dear R helpers, >From the lrm( ) model used for binary logistic regression, we used the L.R. model value (or the G2 value, likelihood-ratio chi-squared statistic) to evaluate the goodness-of-fit of the models. The model with the lowest G2 value consequently, has the best performance and the hi

Re: [R] Finding regions of overlap

On 8/12/05 9:30 AM, "Liaw, Andy" <[EMAIL PROTECTED]> wrote: > Am I missing something? It seems quite easy to me: > > Suppose the intervals are (l1, u1), (l2, u2), ..., (ln, un). The > intersection, if non-empty, would be (max(l1, ..., ln), min(u1, ..., un)), > and the union (min(l1, ..., ln), m

Re: [R] Finding regions of overlap

Am I missing something? It seems quite easy to me: Suppose the intervals are (l1, u1), (l2, u2), ..., (ln, un). The intersection, if non-empty, would be (max(l1, ..., ln), min(u1, ..., un)), and the union (min(l1, ..., ln), max(u1, ..., un)). Andy > From: Sean Davis > > I have a number of reg

Re: [R] Problem with lme4

On 8/12/05, Peter Ho <[EMAIL PROTECTED]> wrote: > Hi, > > I cannot seem to get lme4 to work. I have installed the lme4 and Matrix > package with apt-get. and both can be found in /usr/lib/R/site-library. > When I tried an example for lmer, R could not find the function lmer(), Try using install.

Re: [R] Concerning reading of SAS-files

The first one is an index, not a data set. Anyway, just use SAS to export the data sets in text format (CSV, tab-delimited etc). You can then easily read those in R. (By the way, the help for read.xport says that 'The file must be in SAS XPORT format.' Is .sas7bdat an XPORT file? Hint: no.) > ---

[R] Finding regions of overlap

I have a number of regions consisting of (start,end) pairs. I would like to find regions of overlap among various sets of these regions. In other words, I would like to find unions/intersections of arbitrary sets of regions. Does this smell of anything anyone has already done? Thanks, Sean ___

Re: [R] coercing created variables into a new data frame using na.omit()

I don't know if you can read your message, but I find it exceedingly difficult and there seem to be several typos. Please use the space and return keys ... and only send a message once. You problem is perhaps that you are not looking at the data frame, but at the variable in the workspace. at

[R] coercing created variables into a new data frame using na.omit()

Hi, I am an R newbie and one thing I am having trouble with binding variables that I have created within one data frame into a new data frame when using na.omit(). To illustrate this problem I will give the example I am working on and the approah I have been using:- data.frame1<-filepath att

Re: [R] error tracing

Steve Adams yahoo.com> writes: > Hi, I am running some Cox modeling on large number of > variables (thousands) using apply. For some reasons, > some of the variables have problems in Cox regression, > so the the run was stopped automatically. Is there a > way to keep the running for all the var

[R] coercing created variables into a new data frame using na.omit()

Hi, I am an R newbie and one thing I am having trouble with binding variables that I have created within one data frame into a new data frame when using na.omit(). To illustrate this problem I will give the example I am working on and the approah I have been using:- data.frame1<-filepath att

Re: [R] chisq warning

Stephen Choularton bymouth.com> writes: > I am running chisq as below and getting a warning. Can anyone tell me > the significance or the warning? > > > chisq.test(c(10 ,4 ,2 ,6 ,5 ,3 ,4 ,4 ,6 ,3 ,2 ,2 ,2 ,4 ,7 ,10 ,0 ,6 > ,19 ,3 ,2 ,7 ,2 ,2 ,2 ,1 ,32 ,2 ,3 ,10 ,1 ,3 ,9 ,4 ,10 ,2 ,2 ,4 ,5 ,7 ,6

[R] Dating Objects

Hi I know this subject has been mentioned before but from the mail archives I'm under the impression that this is not possible ? I'm trying to carryout the equivalent of ls -l in R to give some date/time label to each of my objects If the case is that there is no equivalent, is it possible to l

Re: [R] How to insert a certain model in SVM regarding to fixed kernels

On 8/12/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > David, Please correct me if I am wrong but I think svm partially works > with dyn although I don't remember what the specific limitations were. > Its possible that what works already is enough for Amir. For example, > > library(e1071) > l

Re: [R] How to insert a certain model in SVM regarding to fixed kernels

David, Please correct me if I am wrong but I think svm partially works with dyn although I don't remember what the specific limitations were. Its possible that what works already is enough for Amir. For example, library(e1071) library(dyn) set.seed(1) y <- ts(rnorm(100)) y.svm <- dyn$svm(y ~ lag(y

Re: [R] Exponential, Weibull and log-logistic distributions in glm()

>> Dear R-users! >> >> I would like to fit exponential, Weibull and log-logistic via glm() like >> functions. Does anyone know a way to do this? Bellow is a bit longer >> description of my problem. > > I think you want to use survreg(). It will still work when there is no > censoring. > > > Addin

[R] Problem with lme4

Hi, I cannot seem to get lme4 to work. I have installed the lme4 and Matrix package with apt-get. and both can be found in /usr/lib/R/site-library. When I tried an example for lmer, R could not find the function lmer(), > library(lme4) Attaching package: 'lme4' The following object(s

Re: [R] R_X11 module cannot be loaded

On 8/12/05 6:50 AM, "Pat Meyer" <[EMAIL PROTECTED]> wrote: > Hi, > > I encountered a problem while working in R on our Server (4 Itanium 2 > processors). We have redhat for IA64 installed, and R version R-2.0.1 > (2004-11-15). > > Here is my problem. > > When I run R in a console, the png comma

Re: [R] Digest reading is tedious

Hi Adrian, Here's what I used to use for a list I used to subscribe to (it is the first example in the man page for formail, it worked for me and I never went further playing with formail): :0: * ! * [EMAIL PROTECTED] | formail +1 -ds >>new/expert Here is a rough summary of what it does: :0:

[R] R_X11 module cannot be loaded

Hi, I encountered a problem while working in R on our Server (4 Itanium 2 processors). We have redhat for IA64 installed, and R version R-2.0.1 (2004-11-15). Here is my problem. When I run R in a console, the png command works flawlessly and it outputs a png file. Nevertheless, when I run R

Re: [R] Digest reading is tedious

On Tuesday 09 August 2005 23:47, Martin Maechler wrote: > > "Trevor" == Trevor Hastie <[EMAIL PROTECTED]> > > on Tue, 9 Aug 2005 10:27:32 -0700 writes: > > Trevor> [...snip...] > > But that has been an option in mailman, the software behind our > mailing lists --- for ages --- > >

Re: [R] How to insert a certain model in SVM regarding to fixed kernels

Amir, > > Suppose that we want to regress for example a certain autoregressive > model using SVM. We have our data and also some fixed kernels in > libSVM behinde e1071 in front. The question: Where can we insert our > certain autoregressive model ? During creating data frame ? Yes, I think.

[R] handling warning messages

Hi, I have query regarding R & Rserve. In Rserve, there is a way to capture Errors by RSrvException class, but is there any way to capture warning messages? I have found that there is "warnings()" command in R, which lists the last warning message, but I am not able to get the warning message

[R] R Help

Hi, I have query regarding R & Rserve. In Rserve, there is a way to capture Errors by RSrvException class, but is there any way to capture warning messages? I have found that there is "warnings()" command in R, which lists the last warning message, but I am not able to get the warning messa

Re: [R] RGUI crash when opening script in XP Home enviroment

>Prof Brian Ripley schrieb: > > > >>However, as I don't know what the `open file button' is, I cannot >>reproduce it. Is this MDI or SDI mode, >> >> found it : MDI Mode Pager Style multiple windows Knut Krueger __ R-help@stat.math.ethz.ch ma

[R] help on cross hedge optimal hedge variance ratio

Hi everyone I am trying to estimate the optimal hedge variance ratio for cross hedging two commodities. the price levels are used (compared to price change and % price change) and used the OLS with dummy variable for estimating the co-efficients. the equation looks like this Y = B + B1*D1 + B2*X

[R] Concerning reading of SAS-files

Hi! I'm trying to start a credibility estimation study with a coule of data sets that are created for SAS. The data sets are saved as ".sas7bndx" and ".sas7bdat". I've tried reading them to R with the function 'read.xport' but this returns the error message 'Error in lookup.xport(file) : unabl