Re: [R] RGui crashes on wle call

On Mon, 13 Jun 2005, Chris Bergstresser wrote: Hi all -- I'm seeing the following commands reliably produce a crash in RGui, version 2.0.1, for both my home and office machine: rm(list = ls(all = TRUE)); load(dataset.R); library(wle); data.wle = wle.lm(abortion ~ year * lib.con +

[R] protection stack overflow??

Hi dear Rers, I am using SSOAP package to access SOAP service at NCBI. I followed the example code in SSOAP but failed. z - .SOAP(http://www.ncbi.nlm.nih.gov/entrez/eutils/soap/soap_adapter.cgi;, method=run_eInfo, db=pubmed, action = I(einfo)) Error: protect(): protection stack overflow

Re: [R] superscript in figures - basic question

Peter == Peter Alspach [EMAIL PROTECTED] on Tue, 14 Jun 2005 14:11:47 +1200 writes: Peter Ben Peter Others have pointed out plotmath. However, for some Peter superscripts (including 2) it may be easier to use Peter the appropriate escape sequence (at in Windows):

Re: [R] superscript in figures - basic question

Chuck Cleland [EMAIL PROTECTED] writes: See ?plotmath. plot(1:10, 1:10, ylab=expression(BA (m^{2}/ha))) Ick... The parser will think that BA is a function and m^2/ha its argument. This probably has consequences for the spacing. Try expression(BA * (m^2/ha)) -- O__ Peter Dalgaard

[R] New Family object for GLM models...

Dear R-Users, I wish to create a new family object based on the Binomial family. The only difference will be with the link function. Thus instead if using the 'logit(u)' link function, i plan to use '-log(i-u)'. So far, i have tried to write the function following that of the Binomial and

Re: [R] protection stack overflow??

On Tue, 14 Jun 2005, Hu Chen wrote: Hi dear Rers, I am using SSOAP package to access SOAP service at NCBI. I followed the example code in SSOAP but failed. z - .SOAP(http://www.ncbi.nlm.nih.gov/entrez/eutils/soap/soap_adapter.cgi;, method=run_eInfo, db=pubmed, action = I(einfo)) Error:

[R] why does the unsubscribe not work

Hi I have tried to unsubscribe from this list now three times, I desperately need to unsubscribe from this list from this address because the list is choking this mailbox. Below are the various suggested ways of unsubscribing To subscribe or unsubscribe via the World Wide Web, visit

[R] ordinary polynomial coefficients from orthogonal polynomials?

How can ordinary polynomial coefficients be calculated from an orthogonal polynomial fit? I'm trying to do something like find a,b,c,d from lm(billions ~ a+b*decade+c*decade^2+d*decade^3) but that gives: Error in eval(expr, envir, enclos) : Object a not found decade - c(1950, 1960, 1970,

Re: [R] why does the unsubscribe not work

[EMAIL PROTECTED] writes: Hi I have tried to unsubscribe from this list now three times, I desperately need to unsubscribe from this list from this address because the list is choking this mailbox. Below are the various suggested ways of unsubscribing To subscribe or unsubscribe via

[R] bs() function of the splines package

Laurent 14 juin 12:02 afficher les options De : Laurent [EMAIL PROTECTED] - Rechercher les messages de cet auteur Date : Tue, 14 Jun 2005 03:02:37 -0700 Local : Mar 14 juin 2005 12:02 Objet : bs() function of the splines package Répondre | Répondre à l'auteur | Transférer |

Re: [R] Learning, if and else

You can also use switch() instead of ifelse(), which makes the code a bit easier to read. The downside to this is that switch does not take vectorised input and thus you need a loop. For example # data dbh - c(30,29,15,14,30,29) form - factor(c(tree, tree, liana, liana, palm, palm)) df -

Re: [R] ordinary polynomial coefficients from orthogonal polynomials?

On Tue, 14 Jun 2005, James Salsman wrote: How can ordinary polynomial coefficients be calculated from an orthogonal polynomial fit? Why would you want to do that? predict() is perfectly happy with an orthogonal polynomial fit and the `ordinary polynomial coefficients' are rather badly

Re: [R] superscript in figures - basic question

On Tue, 14 Jun 2005, Martin Maechler wrote: Peter == Peter Alspach [EMAIL PROTECTED] on Tue, 14 Jun 2005 14:11:47 +1200 writes: Peter Ben Peter Others have pointed out plotmath. However, for some Peter superscripts (including 2) it may be easier to use Peter the

[R] update.packages() - gregmisc

Dear all, I have a problem to update package gregmisc. After I update, update.packages(ask='graphics') trying URL 'http://cran.at.r-project.org/bin/windows/contrib/2.1/gregmisc_2.0.8.zip' Content type 'application/zip' length 2465 bytes opened URL downloaded 2465 bytes package 'gregmisc'

Re: [R] update.packages() - gregmisc

On 6/14/05, Muhammad Subianto [EMAIL PROTECTED] wrote: Dear all, I have a problem to update package gregmisc. After I update, update.packages(ask='graphics') trying URL 'http://cran.at.r-project.org/bin/windows/contrib/2.1/gregmisc_2.0.8.zip' Content type 'application/zip' length 2465

[R] bs() function of the splines package with intercept=FALSE

Hello, I'm implementing a function using non uniform B-Splines. Looking at the code of the bs() function, I realized that if the intercept was set to FALSE, the behavior of the function was the following (df is the number of degrees of freedom that I believe can be interpreted as the number of

Re: [R] update.packages() - gregmisc

Thanks. I do like this, remove.packages(gregmisc, .libPaths()[1]) remove.packages(gtools, .libPaths()[1]) install.packages(gregmisc, .libPaths()[1]) update.packages() update.packages() install.packages(gtools, .libPaths()[1]) update.packages() update.packages()

Re: [R] why does the unsubscribe not work

PD Peter Dalgaard [EMAIL PROTECTED] on 14 Jun 2005 12:06:53 +0200 writes: PD [EMAIL PROTECTED] writes: Hi I have tried to unsubscribe from this list now three times, I desperately need to unsubscribe from this list from this address because the list is choking this

Re: [R] why does the unsubscribe not work

You can test this suggestion by asking for a password reminder from the https site. If you do get an email that says what your password is, then you are still subscribed. I actually asked the website to send me a password reminder to an account that is not subscribed. Although it says that A

[R] load ing and saving R objects

Does anyone know a way to do the following: Save a large number of R objects to a file (like load() does) but then read back only a small named subset of them . As far as I can see, load() reads back everything. The context is: I have an application which will generate a large number of large

[R] Puzzled in utilising summary.lm() to obtain Var(x)

I have a program which is doing a few thousand runs of lm(). Suppose it is a simple model y = a + bx1 + cx2 + e I have the R object d where d - summary(lm(y ~ x1 + x2)) I would like to obtain Var(x2) out of d. How might I do it? I can, of course, always do sd(x2). But it would be much

[R] RES: Learning, if and else

Thank you very much! All methods are very useful. Paulo M. Brando Instituto de Pesquisa Ambiental da Amazonia (IPAM) Santarem, PA, Brasil. Av. Rui Barbosa, 136. Fone: + 55 93 522 55 38 www.ipam.org.br E-mail: [EMAIL PROTECTED] -Mensagem

Re: [R] load ing and saving R objects

Richard Mott wrote: Does anyone know a way to do the following: Save a large number of R objects to a file (like load() does) but then read back only a small named subset of them . As far as I can see, load() reads back everything. Save them to individual files when you generate them?

Re: [R] load ing and saving R objects

I would suggest saving each object to an individual file with some sort of systematic file name. That way, you can implement a rudimentary key-value database and load only the objects you want. You might be interested in the 'serialize()' and 'unserialize()' functions for this purpose. If

Re: [R] ordinary polynomial coefficients from orthogonal polynomials?

Prof Brian Ripley wrote: On Tue, 14 Jun 2005, James Salsman wrote: How can ordinary polynomial coefficients be calculated from an orthogonal polynomial fit? Why would you want to do that? predict() is perfectly happy with an orthogonal polynomial fit and the `ordinary polynomial

Re: [R] ordinary polynomial coefficients from orthogonal polynomials?

I think this is covered in the FAQ: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Models -roger James Salsman wrote: How can ordinary polynomial coefficients be calculated from an orthogonal polynomial fit? I'm trying to do something like find a,b,c,d from lm(billions ~

[R] Manipulating dates

Hello, Given a vector of characters, or factors, denoting the date in the following way: 28/03/2000, is there a method of 1) Computing the earliest of these dates; 2) Using this as a base, then converting all the other dates into merely the number of days after this minimum date Many

Re: [R] Manipulating dates

Use POSIX. To convert: my.dates - strptime(your.characters, format='%d/%m/%Y') once you have that, you can use 'min' to find the minimum. 'difftime' will give you the differences. Jim __ James HoltmanWhat is the problem you

Re: [R] Manipulating dates

try this: x - c(28/03/2000, 27/08/2001, 29/05/2002, 15/12/2003) y - as.Date(x, %d/%m/%Y) y - min(y) 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

[R] some suggestion

it seems R has no function to sort the data.frame according to some variable(s),though we can do these by order() and indexing.but why not make sort() the a generic function,making it can sort vector and other type of objects? maybe this is a silly suggestion,but i think it is quite usefull.

Re: [R] ordinary polynomial coefficients from orthogonal polynomials?

On Tue, 14 Jun 2005, Frank E Harrell Jr wrote: Prof Brian Ripley wrote: On Tue, 14 Jun 2005, James Salsman wrote: How can ordinary polynomial coefficients be calculated from an orthogonal polynomial fit? Why would you want to do that? predict() is perfectly happy with an orthogonal

Re: [R] some suggestion

take a look at this function by Kevin Wright RSiteSearch(sort.data.frame) 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:

[R] problems with the fitted values of the function gls

Dear helpers, I am a beginner user of R and probably I am not using properly the function gls() of the library nlme. I estimated a model with AR(1) residuals but apparently the fitted and the predicted values arent constructed considering the autocorrelation of the residuals. It seems that

Re: [R] load ing and saving R objects

On Tue, 14 Jun 2005, Barry Rowlingson wrote: Richard Mott wrote: Does anyone know a way to do the following: Save a large number of R objects to a file (like load() does) but then read back only a small named subset of them . As far as I can see, load() reads back everything. Save them

Re: [R] New Family object for GLM models...

On Tue, 14 Jun 2005, Abatih E. wrote: Dear R-Users, I wish to create a new family object based on the Binomial family. The only difference will be with the link function. Thus instead if using the 'logit(u)' link function, i plan to use '-log(i-u)'. So far, i have tried to write the

[R] c(recursive=TRUE)

Hi R users, I am currently using c(...,recursive=TRUE) to handle list-structured objects. This allows to represent something like: l1 = list(level1=1,level2=list(sub1=1,sub2=2)) as: (l1names = names(c(l1,recursive=TRUE))) [1] level1 level2.sub1 level2.sub2 Then, one can use: (l1names

Re: [R] use of gam

On Mon, 13 Jun 2005, Kerry Bush wrote: Suppose I fit the following model: library(gam) fit - gam(y~x1+x2+s(x3),family=binomial) and then I use fitf$coef x1x2 s(x3) 4.1947460 2.7967200 0.0788252 are the coefficients for x1 and x2 the estimated coefficients?

[R] plotting confidence intervals for polynomial models? (was Re: ordinary polynomial coefficients from orthogonal polynomials?)

What I really want now are the 95% confidence intervals that I mistakenly thought that linesHyperb.lm() would provide: pm - lm(billions ~ I(decade) + I(decade^2) + I(decade^3)) library(sfsmisc) linesHyperb.lm(pm) Error in names-.default(`*tmp*`, value = c(I(decade), I(decade^2), :

Re: [R] load ing and saving R objects

Thanks everyone for help on this. It looks like the solution is a zillion files. Richard -- Richard Mott | Wellcome Trust Centre tel 01865 287588 | for Human Genetics fax 01865 287697 | Roosevelt Drive, Oxford OX3 7BN

[R] t.test using RSPerl

Hi, I've just started using R and RSPerl. I have some code as follows: R::initR(--no-save); R::call(t.test, ([EMAIL PROTECTED], [EMAIL PROTECTED])); where @array1 and @array2 are both 1-dimensional arrays in Perl having 54675 elements each. On execution the output is as follows: Calling R

Re: [R] t.test using RSPerl

This has nothing to do with RSPerl, instead it has to do what kind of object you obtain and how these are print():ed. Typing the name of an object, say, 'res', at R prompt and pressing ENTER; res is equivalent as typing print(res) This is for convenience to the user. Basically, this is

Re: [R] extracting components of a list

Dimitris, wouldn't this be more precise --- sapply(jj,function(x) which(x$b[1]==4)) [[1]] [1] 1 [[2]] numeric(0) [[3]] [1] 1 John Dimitris wrote --- maybe something like this: jj - list(list(a = 1, b = 4:7), list(a = 5, b = 3:6), list(a = 10, b = 4:5)) ### jj[sapply(jj,

[R] Calling C from Fortran

I would like to call C routines from Fortran as suggested in section 5.6 of the Writing R extensions documentation. I'm familiar with Fortran but not with C. I understand the example provided in Fortran: subroutine testit() double precision normrnd, x call rndstart() x = normrnd() call dblepr(X

[R] Dateticks

Hello. I am having the worst time converting x-axis date ticks to real dates. I have tried several suggestions in online help tips and books to no avail. For example, the x-axis has 0, 50, 100, etc, and I want it to have 6/17/03, 8/6/03 etc. See attached (sample). Can anybody help me with

Re: [R] load ing and saving R objects

On Tue, 14 Jun 2005, Prof Brian Ripley wrote: If your file system does not like 15000 files you can always save in a DBMS. Or, switch to a better/more appropriate file system: http://en.wikipedia.org/wiki/Comparison_of_file_systems ReiserFS would allow you to store up to about 1.2 million

[R] plots

Dear all, Is it possible to change the levels in a mosaic plot, the appearance of the level or the levels size? For instance: A C E B D Thanks for your help Adrián __

[R] KMEANS output...

Using R 2.1.0 on Windows 2 questions: 1. Is there a way to parse the output from kmeans within R? 2. If the answer to 1. is convoluted or impossible, how do you save the output from kmeans in a plain text file for further processing outside R? Example: ktx-kmeans(x,12, nstart = 200) I would

Re: [R] KMEANS output...

help(kmeans) describes the return value: a list with a) the complete classification b) the cluster centers c) the within-cluster sums of squares b) the cluster sizes. Thankfully, you don't need any parsing. Reid Huntsinger -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

[R] problem installing packages with compiled-from-source R.app on Mac OS X - Tiger

Hello all, This may be aimed for r-devel, but I encountered this as an R-user and not an R-developer so I start here (having said that, please direct me to R-devel if you think this is appropriate. I am not cross- posting, as I believe this is bad netiquette). I am a recent, but extremely

Re: [R] Dateticks

try this example: x.1 - strptime(6/17/03,'%m/%d/%y') x.1 [1] 2003-06-17 plot(0:250, xaxt='n') dates - x.1 + c(0,50,100,150,200,250) * 86400 # 'dates' is in seconds, so add the appropriate number of days dates [1] 2003-06-17 00:00:00 EDT 2003-08-06 00:00:00 EDT 2003-09-25 00:00:00 EDT

[R] using forecast() in dse2 with an ARMA model having a trend component

(My apologies if this is a repeated posting. I couldn't find any trace of my previous attempt in the archive.) I'm having trouble with forecast() in the dse2 package. It works fine for me on a model without a trend, but gives me NaN output for the forecast values when using a model with a

[R] any function to calculate the lamda coefficient?

i just write one for myself,but the result is different from the spss's. i anyone can send me a copy of function calc the lamda coefficient?(i know i should google first,but right now i can surf the internet.) lamda function(table){ obs1.y-sum(apply(table,1,max))

[R] Plotting rows (or columns) from a matrix in different graphs, not using par

Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function par and then struggling around with axes and labels etc. Luckily, I would always have rows + columns = even

[R] How to fix false positve rates?

Dear R-users, I have a set of 12000 image samples. I can divide this set into two categories: training and testing. I need to classify the test set into a two qualitative outputs: true or false for some characteristic. To do the classification I'm using the packages SVM in e1071 library and

[R] ANN: BioC2005 Conference scheduled for August in Seattle

== BioC2005 Where Software and Biology Connect == http://bioconductor.org/meeting05/ About BioC2005 == This conference will highlight current developments within and beyond Bioconductor, a world-wide open

[R] Survfit,newdata and continuous time varying covariates

Hi All, I am working with three time varying covariates in a coxph model. I cannot seem to figure out how to use survfit and the newdata argument to provide estimated survival curves for two scenarios of one covariate while holding the other two at the mean value. Is it possible to display how

[R] Plotting quiver vector tensor arrows 2d field data

Hi All, I'd like to plot something like http://www.nawcwpns.navy.mil/~weather/mugu/mesodata/analysis.html Looking through the galleries at http://addictedtor.free.fr/graphiques/allgraph.php http://r-spatial.sourceforge.net/gallery/ http://fawn.unibw-hamburg.de/cgi-bin/Rwiki.pl?GraphGallery

Re: [R] using forecast() in dse2 with an ARMA model having a trend component

Scott This works for me: arma.pred.without.trend$forecast[[1]][,1] [1] 106.0038 105.9789 105.9605 105.9396 105.9224 105.9052 105.8926 105.8849 [9] 105.8812 105.8880 105.9043 105.9240 105.9579 105.9878 105.9901 106.0095 [17] 106.0555 106.0782 106.0644 106.0427 106.0297 106.0072 106.0126

Re: [R] problem installing packages with compiled-from-source R.app on Mac OS X - Tiger

You need to use R-patched, nowadays 2.1.1 beta, not the released R 2.1.0. As far as I know it is fixed there. On Tue, 14 Jun 2005, Constantinos Antoniou wrote: Hello all, This may be aimed for r-devel, but I encountered this as an R-user and not an R-developer so I start here (having said

Re: [R] some suggestion

Shame nobody has included this function in their package! Would be useful to have as a standard function! Sander. Dimitris Rizopoulos wrote: take a look at this function by Kevin Wright RSiteSearch(sort.data.frame) Best, Dimitris Dimitris Rizopoulos Ph.D. Student

[R] bug in rpart?

Dear R-helpers, Can you help me to see why code 1 gives error while code 2 runs fine? The only difference in the data is the distribution of age categories. I am attaching the session after the code. Many thanks. XL library(survival) library(rpart) # code 1 n - 20 age - rep(1:3, c(2, 3,

[R] trivial installation question

Hello, I try to install udunits and RNetCDF packages With the first one after R CMD INSTALL udunits_1.2.tar.gz the script looks for the library, does not find and says to edit udunits_1.0/udunits/src/Makevars.in with links to an include file and location of the library. When I put the correct path

Re: [R] Plotting rows (or columns) from a matrix in different graphs, not using par

On Tue, 2005-06-14 at 20:00 +0200, Camarda, Carlo Giovanni wrote: Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function par and then struggling around with

Re: [R] problem installing packages with compiled-from-source R.app on Mac OS X - Tiger

Thank you. I am compiling it now... Costas On 14 2005, at 10:13 , Prof Brian Ripley wrote: You need to use R-patched, nowadays 2.1.1 beta, not the released R 2.1.0. As far as I know it is fixed there. On Tue, 14 Jun 2005, Constantinos Antoniou wrote: Hello all, This may be aimed

[R] within and between subject calculation

Dear helpers in this forum, I have the following question: Suppose I have the following data set: id x y 023 1 2 023 2 5 023 4 6 023 5 7 412 2 5 412 3 4 412 4 6 412 7 9 220 5 7 220 4 8 220 9 8 .. and i want to calculate sum_{i=1}^k sum_{j=1}^{n_i}x_{ij}*y_{ij} is there a simple way to do

[R] Matrix stability problem

Hello, This is not a problem with R, the calculated results are mathematically correct. This a matrix stability problem. Because of measuring errors, my matrix solution is a bit off. Here is what my equations look like: A11 x11+A12 x12 +A13 x13 = b1 A21 x21+A22 x21 +A23 x23 = b2 A31 x31+A32

Re: [R] Manipulating dates

On 6/14/05, Richard Hillary [EMAIL PROTECTED] wrote: Hello, Given a vector of characters, or factors, denoting the date in the following way: 28/03/2000, is there a method of 1) Computing the earliest of these dates; 2) Using this as a base, then converting all the other dates into

Re: [R] Dateticks

On 6/14/05, Bernard L. Dillard [EMAIL PROTECTED] wrote: Hello. I am having the worst time converting x-axis date ticks to real dates. I have tried several suggestions in online help tips and books to no avail. For example, the x-axis has 0, 50, 100, etc, and I want it to have 6/17/03,

Re: [R] Plotting rows (or columns) from a matrix in different graphs, not using par

On 6/14/05, Camarda, Carlo Giovanni [EMAIL PROTECTED] wrote: Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function par and then struggling around with axes and

Re: [R] Puzzled in utilising summary.lm() to obtain Var(x)

On 6/14/05, Ajay Narottam Shah [EMAIL PROTECTED] wrote: I have a program which is doing a few thousand runs of lm(). Suppose it is a simple model y = a + bx1 + cx2 + e I have the R object d where d - summary(lm(y ~ x1 + x2)) I would like to obtain Var(x2) out of d. How might I do it?

Re: [R] load ing and saving R objects

On 6/14/05, Richard Mott [EMAIL PROTECTED] wrote: Does anyone know a way to do the following: Save a large number of R objects to a file (like load() does) but then read back only a small named subset of them . As far as I can see, load() reads back everything. The context is: I have

Re: [R] How to fix false positve rates?

-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Anderson de Rezende Rocha Sent: Tuesday, June 14, 2005 11:01 AM To: r-help@stat.math.ethz.ch Subject: [R] How to fix false positve rates? Dear R-users, I have a set of 12000 image samples. I can divide

[R] lattice, panel.grid, and scales=list(tick.number=XXX)

I have a Lattice plot in which I want to adjust the number of tick marks used, and I want to have the drawn grid reflect that change. Here is what I'm doing: bwplot(var1 ~ var2, data=df, scales=list(tick.number=10), panel=function(...) { panel.grid(h=0,v=-1,...);

Re: [R] Matrix stability problem

It looks to me like the x measurements are normalized to sum to 1. So whatever measurement error there is gets spread around, so to speak. It would help if you could explain the setting a little more fully? Why is it, for example, that you know the A values from experiment to experiment (they do

Re: [R] Preparing timestamped data for fourier analysis

Milos Zarkovic [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] I believe that FFT is not appropriate. However Lomb-Scargle periodogram could be used. This may interest you: (Preprint of submitted paper) Detecting periodic patterns in unevenly spaced gene expression time series using

Re: [R] lattice, panel.grid, and scales=list(tick.number=XXX)

If you look at the code of panel.grid, you'll see why it doesn't work -- it does not use any of the scale parameters. Moreover the Help page for panel.grid explicitly warns that the h,v=-1 specification may not work, so no promises have been broken. I'm not sure how bwplot handles grids, since

Re: [R] within and between subject calculation

On 6/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Dear helpers in this forum, I have the following question: Suppose I have the following data set: id x y 023 1 2 023 2 5 023 4 6 023 5 7 412 2 5 412 3 4 412 4 6 412 7 9 220 5 7 220 4 8 220 9 8 .. and i want to calculate

[R] anova.lme error

Hi, I am working with R version 2.1.0, and I seem to have run into what looks like a bug. I get the same error message when I run R on Windows as well as when I run it on Linux. When I call anova to do a LR test from inside a function, I get an error. The same call works outside of a

[R] Reducing the FPR (false positive rate)

Hello R-USERS, I think some people didn't understand my question. What I want is to use the training set to minimize the FALSE POSITIVE RATE. I think it is possible. I sacrifice ACCURACY to have less FALSE POSITIVES. I don't want a classifier result with 5% of FPR and, for example, 93% of

Re: [R] Calling C from Fortran

Dear Gilles, You can try this. Save the fortran file as 1.f Save the C file as 2.c Then R CMD SHLIB 1.f 2.c You should obtain a file named 1.so. Then start you R program, then dyn.load('/where/is/your/1.so') .Fortran(testit) You should obtain a random number which is normally distributed.

Re: [R] ordinary polynomial coefficients from orthogonal polynomials?

Prof Brian Ripley wrote: On Tue, 14 Jun 2005, Frank E Harrell Jr wrote: Prof Brian Ripley wrote: On Tue, 14 Jun 2005, James Salsman wrote: How can ordinary polynomial coefficients be calculated from an orthogonal polynomial fit? Why would you want to do that? predict() is perfectly

Re: [R] RGui crashes on wle call

Prof Brian Ripley wrote: On Mon, 13 Jun 2005, Chris Bergstresser wrote: I'm seeing the following commands reliably produce a crash in RGui, version 2.0.1, for both my home and office machine: rm(list = ls(all = TRUE)); load(dataset.R); library(wle); data.wle = wle.lm(abortion ~ year

Re: [R] Puzzled in utilising summary.lm() to obtain Var(x)

I have a program which is doing a few thousand runs of lm(). Suppose it is a simple model y = a + bx1 + cx2 + e I have the R object d where d - summary(lm(y ~ x1 + x2)) I would like to obtain Var(x2) out of d. How might I do it? I can, of course, always do sd(x2). But it

Re: [R] Puzzled in utilising summary.lm() to obtain Var(x)

On 6/14/05, Ajay Narottam Shah [EMAIL PROTECTED] wrote: I have a program which is doing a few thousand runs of lm(). Suppose it is a simple model y = a + bx1 + cx2 + e I have the R object d where d - summary(lm(y ~ x1 + x2)) I would like to obtain Var(x2) out of d. How

Re: [R] Plotting rows (or columns) from a matrix in different graphs, not using par

On 6/14/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 6/14/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 6/14/05, Camarda, Carlo Giovanni [EMAIL PROTECTED] wrote: Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in

Re: [R] Logistic regression with more than two choices

Dear all R-users, I am a new user of R and I am trying to build a discrete choice model (with more than two alternatives A, B, C and D) using logistic regression. I have data that describes the observed choice probabilities and some background information. An example below describes the data: