[R] help

I want to extract the first non-zero digit of any vector, as for example from a-runif(100,0,1). Tried it using grep (...) but didn´t work. Any help is very much appreciated. Thanks, Manica -- visite http://ipimar-iniap.ipimar.pt/neomav/ __ [EMAIL

Re: [R] Correct Computer Modern font in postscript(..) output

family=ComputerModern gives you CMSL10, which is the `correct font'. From src/main/devPS.c: /* Computer Modern as recoded by Brian D'Urso */ { ComputerModern, {CM_regular_10.afm, CM_boldx_10.afm, CM_italic_10.afm, CM_boldx_italic_10.afm, CM_symbol_10.afm} }, and

[R] [R-pkgs] New versions: mvbutils and debug packages

Dear R users New versions of the 'mvbutils' and 'debug' packages are now available on CRAN, both in source form and as precompiled binaries. 'mvbutils' offers the following (as well as many miscellaneous utilities): ? hiearchical, searchable project organization, with workspaces switchable

Re: [R] Sweave and R output: possible to suppress Schunk tags in *.tex file output?

On Fri, 12 Mar 2004 08:38:53 -0800, Waichler, Scott R (WSR) wrote: Hello, I would like to fill the rows of a Latex tabular environment with output from R, as in \begin{table} \caption{Table caption.} \label{tab:events} \begin{tabular}{c r r r r r} \hline

[R] Re: R-help Digest, Vol 13, Issue 14

You could look at @Article{ e1071-papers:meyer+leisch+hornik:2003, author= {David Meyer and Friedrich Leisch and Kurt Hornik}, title = {The Support Vector Machine under Test}, journal = {Neurocomputing}, year = 2003, month = {September},

[R] [R-pkgs] Subselect package - Version 0.7.1

A new version (0.7.1) of package 'subselect' has been uploaded to CRAN. Package 'subselect' provides functions which assess the quality of variable subsets as surrogates for a full data set, in an exploratory data analysis, and search for subsets which are optimal under various criteria. As

Re: [R] nnet classification accuracy vs. other models

My experience is that nnet needs a lot of tuning, not only in terms of numbers of layers, but also in terms of the other parameters. My first results where I kept very much of the default parameter values with nnet have been very bad, as bad as you say. (But as Brian Ripley already wrote, it's not

Re: [R] help

I want to extract the first non-zero digit of any vector, as for example from a-runif(100,0,1). a[min((1:length(a))[a!=0])] _ Simon Wood [EMAIL PROTECTED]www.stats.gla.ac.uk/~simon/ Department of Statistics,

[R] expanding some values in logical vector

Dear all In automatic dropout evaluation function I construct an index (pointer), which will be TRUE at unusual values. Then I need to expand these TRUE values a little bit forward and backward. Example: having span=5, from vector idx-rep(F,10) idx[4]-T idx [1] FALSE FALSE FALSE TRUE

[R] optim() and La.svd() error code 1

Hello everybody. I am using optim() to minimize a function of 19 variables and I repeatably get the following error message: R source(/users/sat/rksh/goldstein/emulator/optimizer.R) sann objective function values initial value 5044.955275 Error in La.svd(x, nu, nv, method) : error code 1

Re: [R] help

Manica wrote: I want to extract the first non-zero digit of any vector, as for example from a-runif(100,0,1). Tried it using grep (...) but didn´t work. Any help is very much appreciated. Thanks, Manica How about: a - runif(10, 0, 1) b - as.character(a) nonzero - regexpr([1-9], b)

RE: [R] help

On 12-Mar-04 Manica wrote: I want to extract the first non-zero digit of any vector, as for example from a-runif(100,0,1). Tried it using grep (...) but didn´t work. Any help is very much appreciated. Thanks, Manica runif will put real numbers in 'a'. It is highly likely that any of them

[R] Test Relation between Time Series? = Cross Correlation Analysis with ccf()

Dear R world, When investigating two time series and applying a cross correlation ccf() function, the results show that correlation is maximal (ccf=0.8) at a lag of 0.1. e.g. ccf(Inv.KBDIn,ts.medNDII, type=c(correlation), na.action=na.omit, main=c(Cross-Correlation of Inverse KBDI against

[R] creating a ps. file

Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the second cycle of the loop. I

RE: [R] help

I guess what you want is: a - abs(a) floor( a / 10^floor( log10( a ) ) ) Bendix -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED]

[R] imputation of sub-threshold values

Is there a good way in R to impute values which exist, but are less than the detection level for an assay? Thanks, Jonathan Williams OPTIMA Radcliffe Infirmary Woodstock Road OXFORD OX2 6HE Tel +1865 (2)24356 __ [EMAIL PROTECTED] mailing list

Re: [R] expanding some values in logical vector

Try: x-rep(FALSE,20); x[c(4,10,15)]-TRUE x [1] FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE TRUE FALSE FALSE [13] FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE x[outer(which(x),-1:1,+)]-T x [1] FALSE FALSE TRUE TRUE TRUE FALSE FALSE FALSE TRUE TRUE TRUE FALSE [13] FALSE

Re: [R] Problems with SJava instalation

Em Dom 14 Mar 2004 07:26, Ronaldo Reis Jr. escreveu: Em Sex 05 Mar 2004 10:29, Jean Coursol escreveu: Hi, I try to execute ldconfig and the SJava dont work. I try to install with gcc-3.3 (was gcc-2.95 in my first instalation). Where is the problem? Thanks Ronaldo Hi, I finally got to

[R] printing graphics whithout save us

hello, it wanted to know as printing a graphics on a screen, without to have save in a file, already use dev.print and give me an error in PS , unable to start device postscript and a warning printing via file = not implemented.Thanks Ruben __ [EMAIL

RE: [R] creating a ps. file

paste is the function you need. Bendix C. -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED] www.biostat.ku.dk/~bxc

Re: [R] No traceback available when using try(...)

The Luke's proposition it works fine. I have extended the code proposed by Luke, to implemtent a myTry() and a mytraceback(). They print the function stack and the variables of each functions, it may helps. -- mytry.R

[R] Simple numeric as.is question

I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is, but the specifics escape me. I have used x-read.table(clipboard, header=F) to import from a txt file. How

Re: [R] creating a ps. file

On Mon, 2004-03-15 at 05:37, [EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the second

Re: [R] creating a ps. file

for (i in 1:n){ dev.print(postcript,file=paste(c:/Rfigures/plot_,i,.ps,sep=)) } Best! A.S. Alessandro Semeria Models and Simulations Laboratory Montecatini Environmental Research Center (Edison Group), Via Ciro Menotti 48, 48023 Marina di Ravenna (RA), Italy

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the

[R] R equiv to proc gremove in maps package

Is there an R equivalent to SAS's proc gremove? You would use this procedure to combine the units on an existing map, for example to build a map of Metropolitan Statistical Areas (MSAs) from the [US] counties dataset where the internal boundries surround the MSAs (which are groups of counties)

RE: [R] imputation of sub-threshold values

-Original Message- From: Jonathan Williams [mailto:[EMAIL PROTECTED] Sent: 15 March 2004 12:48 To: Ethz. Ch Subject: [R] imputation of sub-threshold values Security Warning: If you are not sure an attachment is safe to open contact Andy on x234. There are 0 attachments

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the

Re: [R] expanding some values in logical vector

Thank you very much. I myself thought about something similar but I was not so smart to use outer. Cheers Petr On 15 Mar 2004 at 13:51, Peter Wolf wrote: Try: x-rep(FALSE,20); x[c(4,10,15)]-TRUE x [1] FALSE FALSE FALSE TRUE FALSE FALSE FALSE FALSE FALSE TRUE FALSE FALSE [13]

Re: [R] creating a ps. file

[EMAIL PROTECTED] wrote: Dear all I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps) I overwrite my results with the

RE: [R] Simple numeric as.is question

If read.table() is not recognizing the data as numeric, there's probably something in the content that confused read.table(). You can try using x - scan(clipboard) instead and see if and how it chokes. Andy From: Janet Gannon I am reading a list of numbers from my clipboard, and have

RE: [R] creating a ps. file

paste(c:/Rfigures/plot, i, .ps, sep=) -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of [EMAIL PROTECTED] Sent: 15 March 2004 11:37 To: [EMAIL PROTECTED] Subject: [R] creating a ps. file Dear all I wrote a routine. At the end of each cycle of the

RE: [R] Simple numeric as.is question

Look at the help for: ?as.numeric HTH, Andy -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Janet Gannon Sent: Monday, March 15, 2004 7:04 AM To: [EMAIL PROTECTED] Subject: [R] Simple numeric as.is question I am reading a list of numbers from

[R] norm of complex number

Hi! I want o calc the norm of a complex number ( |c| where c is complex). Does the function abs() this ? Thanks, -- Frank G. Zoellner AG Angewandte Informatik Technische Fakultat Universitat Bielefeld phone: +49(0)521-106-2951 fax: +49(0)521-106-2992 email: [EMAIL PROTECTED]

[R] setting x-y axis at origin

Hello - I'm just getting into 'R' and am having trouble setting up the x-y axis to share (0,0). In the example posted here: http://geocities.com/chris_ciotti/Images/part1.pdf, each axis has a 0 which I do not want. Any help on getting a graph starting at (0,0) would be greatly appreciated.

Re: [R] Simple numeric as.is question

Janet - Try x2 - as.numeric(as.character(x)) hist(x2) I'm not a Windows user, so I can't test this before sending. It might solve the problem, might not. (Flame !! : This is just ONE MORE example of the difficulties caused by the default behavior of read.table() to make things into

Re: [R] creating a ps. file

Before the loop, open a postscript device with onefile=FALSE, and see the advice in the postscript help page for how to specify the file argument. Then each plot will produce a separate file. (Remember to do dev.off() after the loop.) On Mon, 15 Mar 2004 [EMAIL PROTECTED] wrote: I wrote a

Re: [R] creating a ps. file

On Mon, Mar 15, 2004 at 12:37:02PM +0100, [EMAIL PROTECTED] wrote: I wrote a routine. At the end of each cycle of the loop I would like to save the result (plot) in a postcriptfile. Of course if I just use dev.print in the following way: dev.print(device=postcript, 'c:/Rfigures/plot_1.ps)

Re: [R] R equiv to proc gremove in maps package

Fred Rohde wrote: Is there an R equivalent to SAS's proc gremove? You would use this procedure to combine the units on an existing map, for example to build a map of Metropolitan Statistical Areas (MSAs) from the [US] counties dataset where the internal boundries surround the MSAs (which are

Re: [R] setting x-y axis at origin

I don't understand exactly what you want but the following might help: plot(1:11, 1:11, xlim=c(2, 5), ylim=c(3, 7)) Is this what you want? spencer graves christopher ciotti wrote: Hello - I'm just getting into 'R' and am having trouble setting up the x-y axis to share (0,0). In

Re: [R] setting x-y axis at origin

christopher ciotti wrote: Hello - I'm just getting into 'R' and am having trouble setting up the x-y axis to share (0,0). In the example posted here: http://geocities.com/chris_ciotti/Images/part1.pdf, each axis has a 0 which I do not want. Any help on getting a graph starting at (0,0) would

Re: [R] norm of complex number

That's called the modulus, and computed by function Mod(). On Mon, 15 Mar 2004, Frank Gerrit Zoellner wrote: I want o calc the norm of a complex number ( |c| where c is complex). Does the function abs() this ? -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied

Re: [R] norm of complex number

abs also works (in R 1.8.1 under Windows 2000), as we see just by trying it: abs(1+1i) = 1.414214 hope this helps. spencer graves Prof Brian Ripley wrote: That's called the modulus, and computed by function Mod(). On Mon, 15 Mar 2004, Frank Gerrit Zoellner wrote: I want

Re: [R] Simple numeric as.is question

Janet Gannon wrote: I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is, but the specifics escape me. I have used x-read.table(clipboard, header=F) to import

Re: [R] Simple numeric as.is question

On Mon, Mar 15, 2004 at 09:04:05AM -0500, Janet Gannon wrote: I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is, but the specifics escape me. I have used

RE: [R] Correct Computer Modern font in postscript(..) output

Brian, Thanks for the reponse. I apologise if I gave the impression of being accusatory - not my intention. My point was that when I use italic(), in R with ComputerModern, I get output different to the output I get when I use \italic in TeX/LaTex. I'm neither an R nor TeX pro (for want of

Re: [R] imputation of sub-threshold values

On Mon, 15 Mar 2004, Jonathan Williams wrote: Is there a good way in R to impute values which exist, but are less than the detection level for an assay? If there were a good way to do it, it would probably be implementable or implemented in R. If you can persuade the people measuring the

Re: [R] Simple numeric as.is question

Timur Elzhov [EMAIL PROTECTED] writes: On Mon, Mar 15, 2004 at 09:04:05AM -0500, Janet Gannon wrote: I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is,

Re: [R] Simple numeric as.is question

Janet Gannon wrote: I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is, but the specifics escape me. I have used x-read.table(clipboard, header=F) to import from

Re: [R] Simple numeric as.is question

Janet Gannon wrote: I am reading a list of numbers from my clipboard, and have been successful, except I can't make a histogram as R doesn't recognize my variable as numeric. I know I need to use as.is, but the specifics escape me. I have used x-read.table(clipboard, header=F) to import from

Re: [R] imputation of sub-threshold values

Thomas Lumley [EMAIL PROTECTED] writes: On Mon, 15 Mar 2004, Jonathan Williams wrote: Is there a good way in R to impute values which exist, but are less than the detection level for an assay? If there were a good way to do it, it would probably be implementable or implemented in R. If

RE: [R] Simple numeric as.is question

From: Peter Dalgaard Sent: Monday, March 15, 2004 3:51 PM To: Timur Elzhov Cc: [EMAIL PROTECTED]; Janet Gannon Subject: Re: [R] Simple numeric as.is question Timur Elzhov [EMAIL PROTECTED] writes: On Mon, Mar 15, 2004 at 09:04:05AM -0500, Janet Gannon wrote: I am reading a list

Re: [R] printing graphics whithout save us

[EMAIL PROTECTED] wrote: hello, it wanted to know as printing a graphics on a screen, without to have save in a file, already use dev.print and give me an error in PS , unable to start device postscript and a warning printing via file = not implemented.Thanks Ruben Please read ?dev.print

Re: [R] R CMD check errors

On Fri, 12 Mar 2004, Uwe Ligges wrote: Thomas Stabla wrote: Hello again, I tried to isolate the source code, which causes the error messages: * checking S3 generic/method consistency ... WARNING Error in .tryQuietly({ : Error in library(package, lib.loc = lib.loc, character.only =

[R] make check on Solaris 8 fails due to plot

Dear all, I am having trouble trying to install R-1.8.1 on a Sun Solaris 8 (Generic_108528-23 version) machine. The configuration was successful but make check fails. I traced the the problem to the plot() function. 1 + 1 [1] 2 capabilities() jpeg pngtcltk X11GNOME

Re: [R] R CMD check errors

Thomas Stabla wrote: On Fri, 12 Mar 2004, Uwe Ligges wrote: Thomas Stabla wrote: Hello again, I tried to isolate the source code, which causes the error messages: * checking S3 generic/method consistency ... WARNING Error in .tryQuietly({ : Error in library(package, lib.loc = lib.loc,

RE: [R] Simple numeric as.is question

On Mon, 15 Mar 2004, Andy Bunn wrote: Look at the help for: ?as.numeric and make sure you read there that it is the wrong thing (as a column of a data frame will be a factor). Better to look this up in the FAQ! HTH, Andy -Original Message- From: [EMAIL PROTECTED]

Re: [R] norm of complex number

?complex reveals Mod: Mod(1+1i) [1] 1.414214 hope this helps. spencer graves p.s. help.search(complex) reveals complex, etc. Frank Gerrit Zoellner wrote: Hi! I want o calc the norm of a complex number ( |c| where c is complex). Does the function abs() this ? Thanks,

Re: [R] make check on Solaris 8 fails due to plot

Adaikalavan Ramasamy [EMAIL PROTECTED] writes: Dear all, I am having trouble trying to install R-1.8.1 on a Sun Solaris 8 (Generic_108528-23 version) machine. The configuration was successful but make check fails. I traced the the problem to the plot() function. 1 + 1 [1] 2

Re: [R] setting x-y axis at origin

?par, look at xaxs and yaxs, especially style i. On Mon, 15 Mar 2004, christopher ciotti wrote: I'm just getting into 'R' and am having trouble setting up the x-y axis to share (0,0). In the example posted here: http://geocities.com/chris_ciotti/Images/part1.pdf, each axis has a 0 which

Re: [R] setting x-y axis at origin

On Mon, 2004-03-15 at 08:45, christopher ciotti wrote: Hello - I'm just getting into 'R' and am having trouble setting up the x-y axis to share (0,0). In the example posted here: http://geocities.com/chris_ciotti/Images/part1.pdf, each axis has a 0 which I do not want. Any help on

RE: [R] norm of complex number

Not abs(), but Mod(). See ?complex. HTH, Andy From: Frank Gerrit Zoellner Hi! I want o calc the norm of a complex number ( |c| where c is complex). Does the function abs() this ? Thanks, -- Frank G. Zoellner AG Angewandte Informatik Technische Fakultat Universitat Bielefeld

Re: [R] Simple numeric as.is question

On Mon, 15 Mar 2004, Tom Blackwell wrote: Try x2 - as.numeric(as.character(x)) hist(x2) I'm not a Windows user, so I can't test this before sending. It might solve the problem, might not. (Flame !! : This is just ONE MORE example of the difficulties caused by the default behavior

Re: [R] Simple numeric

Regarding the flame, your might be interested in the following: read.table(myfile, as.is=T) which will interpret character columns as character rather than factor and still reads in the numeric columns as numbers. You can also specify specific columns such as as.is=3:4 if you want columns 3

[R] creating vignettes ... ERROR

The package build fails with: * creating vignettes ... ERROR /package/R/R-1.8.1/linux/lib/R/bin/texi2dvi: pdflatex exited with bad status, quitting. /package/R/R-1.8.1/linux/lib/R/bin/texi2dvi: see SamplesSession.log for errors. Error in texi2dvi(file = bft, pdf = TRUE, clean = TRUE, quiet =

Re: [R] norm of complex number

Frank == Frank Gerrit Zoellner [EMAIL PROTECTED] on Mon, 15 Mar 2004 16:31:58 +0100 writes: Frank Hi! I want o calc the norm of a complex number ( |c| Frank where c is complex). Does the function abs() this ? Yes, it does - and help(abs) will tell you so in the next version of R.

[R] effect size

Hi, Having searched google '[R] aov effect size' without any results I wonder if I not completely miss something. Is there any R function that calculates the effect size of an AOV's main effect or interaction effect? It should be related to the F's and the degree of freedom of the error,

[R] spearman rank correlation problem

Hello R gurus, I want to calculate the Spearman rho between two ranked lists. I am getting results with cor.test that differ in comparison to my own spearman function: my.spearman function(l1, l2) { if(length(l1) != length(l2)) stop(lists must have same length) r1 - rank(l1) r2

[R] A beginners question

Hi. I used R years ago at Uni but haven't touched it since, so I'm back to beginner level now. After looking through the manuals for a while I've come to the conclusion I don't know where to start looking for what I'm looking for, so I decided to ask for help. I've got a fairly simple data

[R] rate of change

Hello I am wondering, how do I find if R has a certain funciton to do a given task. do I just type help.search(rate). I am just trying to find a function to calculate the rate of change for a variable. I could come up with one if there isn't any allready builtin. thanks

[R] How to fit a classification model in R?

I have about 20 000 cases with discrete variables (some are counts, some are factors). I'm interested in fitting a series of models outcome ~ 1 outcome ~ sex outcome ~ sex + age outcome ~ age * sex outcome ~ age * sex + location ... I do NOT expect to get any statistical significance out of this

[R] glm questions

Greetings, everybody. Can I ask some glm questions? 1. How do you find out -2*lnL(saturated model)? In the output from glm, I find: Null deviance: which I think is -2[lnL(null) - lnL(saturated)] Residual deviance: -2[lnL(fitted) - lnL(saturated)] The Null model is the one that includes the

Re: [R] How to fit a classification model in R?

Perhaps step(lm.object) helps, what use a backward elimination of attributes , further i remind in library(leaps) some helpfuel things for modellimg selection. christian Am Dienstag, 16. März 2004 05:53 schrieb Richard A. O'Keefe: I have about 20 000 cases with discrete variables (some

Re: [R] make check on Solaris 8 fails due to plot

See the R-admin manual: this happens for certain broken versions of Solaris-sparc gcc. If that is not the cause we would need to know much more about what you used to build R. We do regularly test on Solaris 8, with gcc 3.3.3 and Forte 7. On Mon, 15 Mar 2004, Adaikalavan Ramasamy wrote: