Re: [R] artificial data matrix with 100000 rows

' option Is '' reached getOption(max.print) '' too difficult to read? You *can* increase the 'max.print' option as much as you like, and that's why I said 'not at all unable' above. Regards, Martin PS For instance, m[(nrow(m)-10):nrow(m),] PS [,1] [,2] [,3] PS [1,] NA

Re: [R] FW: variable format

Frank E Harrell Jr wrote: Martin Becker wrote: Dear Cory, I am not familiar with SAS, but is this what you are looking for? divisionTable - matrix(c(1, New England, 2, Middle Atlantic, 3, East North Central, 4

Re: [R] Matlab's lsqnonlin

! Martin Maechler, ETH Zurich KateM KateM Katharine Mullen KateM mail: Department of Physics and Astronomy, Faculty of Sciences KateM Vrije Universiteit Amsterdam, de Boelelaan 1081 KateM 1081 HV Amsterdam, The Netherlands KateM room: T.1.06 KateM tel: +31 205987870

Re: [R] FW: variable format

(factor(a$divisionOld,levels=divisionTable[,1],labels=divisionTable[,2])) a$divisionNew [1] NA New EnglandMiddle Atlantic [4] East North Central West North Central South Atlantic Kind regards, Martin Cory Nissen schrieb: Anybody

[R] negative value for AIC and BIC

Hi all, I obtained negative values for AIC and BIC criteria for a particular model that I have developped... I don't remember to have negative values for these crietria for others applications, so I am a little suprised... Could anyone tell me if something is wrong or his conclusion

Re: [R] integers

CH == Christoph Heibl [EMAIL PROTECTED] on Tue, 4 Sep 2007 11:53:43 +0200 writes: CH Hi list, CH The function is.integer tests if an object is of type integer: CH see e.g.: CH is.integer(12) # FALSE CH is.real(12) # TRUE CH mode(12)# numeric

Re: [R] Importing huge XML-Files

intermediate objects. Martin Alexander Heidrich [EMAIL PROTECTED] writes: Dear all, for my diploma thesis I have to import huge XML-Files into R for statistical processing - huge means a size about 33 MB. I'm using the XML-Package version 1.9 As far as reading the complete file into R via

Re: [R] Incomplete Gamma function

of pbeta(*, log=TRUE) has been dramatically improved in some extreme tail/range cases. -- which benefits pt(), pf(), pbinom() etc in equivalent situations. Martin Maechler, ETH Zurich and R-core AN Anup AN - Got a little couch AN potato

Re: [R] R and Web Applications

Hi Chris -- RWebServices provides a way to expose R functionality as (SOAP-based) web services. http://bioconductor.org/packages/2.0/bioc/html/RWebServices.html This might be more than you are looking for, and has not been tested on MacOS (which seems to be your platform). Martin Chris Parkin

Re: [R] Behaviour of very large numbers

-5.442553e+80i [9] 3.892581e+80-4.389625e+80i 3.136461e+80-3.536955e+80i Regards, Martin # now this: x - seq(-51,-49,length=100) x^B [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN snip is.numeric(x^B) [1] TRUE is.real(x^B

Re: [R] Rmpi and x86

the library and indicates that it is starting lam. If I have already issued lamboot, then library(Rmpi) loads as expected. Where do things go wrong for you? Martin Edna Bell [EMAIL PROTECTED] writes: Dear R Gurus: Is there a problem with Rmpi on x86 with SUSE 10.1, please? I've tried everything

Re: [R] Rmpi and x86

sorting out the shared library issue) Martin Edna Bell [EMAIL PROTECTED] writes: Here is what happens: Note: lam-7.1.4 linux-tw9c:/home/bell/Desktop/R-2.5.1/bin # ./R CMD INSTALL --clean Rmpi_0.5-3.tar.gz * Installing to library '/home/bell/Desktop/R-2.5.1/library' * Installing *source

Re: [R] library(fCalendar) timeDate(12.03.2005, format=%d.%m.%Y)

Martin Maechler wrote: OL == Ola Lindqvist [EMAIL PROTECTED] on Tue, 21 Aug 2007 14:32:19 +0200 writes: OL Thanks! OL Seems to work fine now! Well, for your example. But sorry to say, the patch breaks other cases. I'm investigating further (and will hopefully

[R] R mailinglist ETH website now with trusted certificate

for several R-project.org servers, notably also the subversion (R source) server, and these will hopefully be put in place as well with the next few weeks. Martin Maechler, ETH Zurich R core __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] library(fCalendar) timeDate(12.03.2005, format=%d.%m.%Y)

of fCalendar once Diethelm Wuertz is back from wherever; I've already made more changes) Martin Maechler, ETH Zurich [but different department than D.Wuertz] OL Martin Becker wrote: Dear Ola, I think you spotted a small bug in *package* fCalendar. Explicit specification should

Re: [R] library(fCalendar) timeDate(12.03.2005, format=%d.%m.%Y)

() (which you are asked for in the posting guide). If you are using Windows and don't know how to apply the patch, you can download a patched binary version here: http://www.saar-gate.net/download/fCalendar_251.70.zip Regards, Martin PS: Maybe r-sig-finance is more appropriate

Re: [R] recommended combo of apps for new user?

interacting with R without getting too bogged down in setting up a new database system -- which sounds like a job all on its own. After I understand a bit more about R's capabilities I'll ease into the database part. Cheers, Martin On 8/19/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: Regarding RODBC

Re: [R] recommended combo of apps for new user?

interacting with R without getting too bogged down in setting up a new database system -- which sounds like a job all on its own. After I understand a bit more about R's capabilities I'll ease into the database part. Cheers, Martin __ R-help@stat.math.ethz.ch

Re: [R] Any parser generator / code assistance for R?

schools If you are still interested in R, look out for documents with S4 (or formal methods and classes) and R in the title ;-) Regards, Martin __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] recommended combo of apps for new user?

productive? (for example, a code editor). Thank you for your time, Martin J. Brown Portland, Oregon [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] recommended combo of apps for new user?

recommend? 3) Any other software I need to learn that would make my work in R more productive? (for example, a code editor). Thank you for your time, Martin J. Brown Portland, Oregon __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] function to find coodinates in an array

GaGr == Gabor Grothendieck [EMAIL PROTECTED] on Thu, 16 Aug 2007 23:46:28 -0400 writes: GaGr Get the indices using expand.grid and then reorder GaGr them: set.seed(1); X - array(rnorm(24), 2:4) # input GaGr X # look at X GaGr do.call(expand.grid, sapply(dim(X),

Re: [R] matching elements from two vectors

: '%in%' is trivially defined via match(). Type `%in%` in your R session. Regards, Martin GF On 8/17/07, jim holtman [EMAIL PROTECTED] wrote: Also if you want all the matches x[x %in% y] [1] 2 3 3 3 On 8/17/07, Gon�alo Ferraz [EMAIL PROTECTED] wrote

Re: [R] residual plots for lmer in lme4 package

Hi Margaret, Have a look at qqmath in the lattice package. ?qqmath Hank On Aug 16, 2007, at 2:45 AM, Margaret Gardiner-Garden wrote: Hi, I was wondering if I might be able to ask some advice about doing residual plots for the lmer function in the lme4 package. Our group's aim is to

Re: [R] Possible memory leak with R v.2.5.0

in the right direction, Martin Peter Waltman [EMAIL PROTECTED] writes: I'm working with a very large matrix ( 22k rows x 2k cols) of RNA expression data with R v.2.5.0 on a RedHat Enterprise machine, x86_64 architecture. The relevant code is below, but I call a function

Re: [R] ung�ltige Versionsspezifikation

old to fit to a current version of cairo (rather 'cairoDevice' ??). And BTW: It is a *package*, not a library!!! Martin Maechler, ETH Zurich JK --- Mag. Ferri Leberl [EMAIL PROTECTED] wrote: Dear everybody, excuse me if this question ist trivial, however, I have now looked

Re: [R] getting lapply() to work for a new class

, but that implies that Test 'is' a list). Martin Pijus Virketis [EMAIL PROTECTED] writes: Thank you. When I tried to set as.list() in baseenv(), I learned that its bindings are locked. Does this mean that the thing to do is just to write my own lapply, which does the coercion using my private as.list

Re: [R] Mixture of Normals with Large Data

of the original data; something others have proposed as well. Martin BertG On 8/7/07, Bert Gunter [EMAIL PROTECTED] BertG wrote: Why would anyone want to fit a mixture of normals with 110 million observations?? Any questions about the distribution that you would care to ask can

[R] Secondary axis

, the second one from 0 to 40). The problem is compounded by the fact that the two data sets have very different frequencies: one consists of twice-monthly measurements, the other of hourly measurements. I would be very grateful for advice on how to do this. Thanks in advance, Patrick Martin

Re: [R] deriv for psigamma

is in R*/src/main/deriv.c and patches which implement the above and (there are few more 'FIXME's there ... ;-) against https://svn.r-project.org/R/trunk/src/main/deriv.c are welcome - after useR!2007 Martin Maechler, ETH Zurich __ R-help@stat.math.ethz.ch

[R] italic greek symbols

Hi Folks, I am using R 2.5.1 on a Mac OS X 10.4.9, via ESS. I would like to try to get an italic mu onto a plot axis label. I note that in a previous email, (Thu, 4 May 2006 19:41:41 +0100 (BST)), Brian Ripley wrote, There is no italic symbol font available on most devices. So unless you try to

Re: [R] italic greek symbols

message: OS reports request to set locale to greek cannot be honored in: Sys.setlocale(LC_CTYPE, greek) Cheers, Hank On Jul 31, 2007, at 3:00 PM, Prof Brian Ripley wrote: On Tue, 31 Jul 2007, Martin Henry H. Stevens wrote: Hi Folks, I am using R 2.5.1 on a Mac OS X 10.4.9, via ESS. I would

Re: [R] Slightly OT - use of R

be the number of subscribers (which I'm sure Martin Maechler can find out). Of course, some people will have subscribed under more than one email address, so that would somewhat over-estimate the number of people who subscribe. But it can be traded off (to a somewhat unknown

[R] options(OutDec) etc {was Strange warning in summary.lm}

the computations can use textConnection() etc.. Martin Maechler, ETH Zurich OT == ONKELINX, Thierry [EMAIL PROTECTED] on Wed, 25 Jul 2007 11:04:43 +0200 writes: OT Dear Peter, Uwe and Brian, OT I've found some more problems with options(OutDec = ,). OT 1) as.numeric yields NA

Re: [R] lme or gls prediction intervals

it is subtle and complicated to figure out what the appropriate hypothesis is, but I was hoping to get some feedback. Cheers, Hank On Jul 25, 2007, at 3:21 AM, [EMAIL PROTECTED] wrote: Martin, Have you checked ?intervals.gls This intervals are approximate, but this would be the obvious

Re: [R] Is Rwiki down?

Me too. Hank On Jul 25, 2007, at 7:50 AM, Jin Lo wrote: Dear R-friends, is Rwiki down? I've been trying to login for the past couple of days with no success. Yours sincerely, Jin __ __ wherever

[R] lme or gls prediction intervals

Hi folks, I am trying to generate 95% confidence intervals for a gls model using predict.nlme with R version 2.5.1 (2007-06-27) . nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-83. I have looked in help, and I can do it for lm and glm models, and I can generate

Re: [R] Computing the rank of a matrix.

RV == RAVI VARADHAN [EMAIL PROTECTED] on Sat, 07 Apr 2007 18:39:36 -0400 writes: this is a bit a late reply... better late than never RV Hi Martin, Hi Ravi, thanks for your research. RV I played a bit with rankMat. I will first state the RV

Re: [R] initalizing and checking validity of S4 classes

'invalid' after the base class has finished with initialize). Martin Bojanowski, M.J. (Michal) [EMAIL PROTECTED] writes: Dear useRs and wizaRds, I am currently developing a set of functions using S4 classes. On the way I encountered the problem exemplified with the code below. For some reason

Re: [R] Package design, placement of legacy functions

in *.R). A user of your **installed** package will be able to use system.file(legacy.R, package = pkg) e.g., as source(system.file(legacy.R, package = pkg)) WA Thanks for the guidance. . . you're welcome. Martin Maechler, ETH Zurich

Re: [R] Rmpi installation error

Sunitha -- I think Rmpi finds an MPICH1 installation, rather than MPICH2. Martin Sunithak [EMAIL PROTECTED] writes: Hello everybody, I am trying to install Rmpi on 64-bit HP-Proliant server. I am facing the following error: * Installing *source* package 'Rmpi' ... Try to find mpi.h

Re: [R] nearest correlation to polychoric

RV == Ravi Varadhan [EMAIL PROTECTED] on Fri, 13 Jul 2007 10:52:36 -0400 writes: RV Martin, I sent you the Matlab code for this which I had RV obtained from Nich Higham. RV Cheng, Sheung Hun and Higham, Nick (1998) A Modified RV Cholesky Algorithm Based on a Symmetric

Re: [R] Extending Matrix class

for a future version of R, since the inherited methods don't know what to do with Weighting, but it would already be interesting if all necessary new methods could be defined semi-automatically: 1) call the corresponding Matrix method 2) pass on all my extra slots Regards, Martin

Re: [R] create a matrix from an excel file?

TS == Tavpritesh Sethi [EMAIL PROTECTED] on Sat, 14 Jul 2007 12:59:03 +0530 writes: TS how do you create a matrix from an excel file read into TS R by the command read.delim? data.matrix(.) {is typical a bit more useful than as.matrix(.) }

Re: [R] R file via SSH

You may find screen useful, see e.g. http://tolstoy.newcastle.edu.au/R/e2/help/07/02/10824.html for a short description how to use it with R (reattaching works via screen -r). (Maybe your sysadmin has to install screen first, it's a unix tool.) KR Martin [EMAIL PROTECTED] wrote

Re: [R] nearest correlation to polychoric

Indefinite Factorization; \emph{SIAM J. Matrix Anal.\ Appl.}, \bold{19}, 1097--1110. Jens, could you make your code (mentioned below) available to the community, or even donate to be included as a new method of posdefify() ? Regards, Martin Maechler, ETH Zurich DR Best, Dimitris DR

Re: [R] how to revert to an older limma version?

and provide commented, minimal, self-contained, reproducible code. -- Martin Morgan Bioconductor / Computational Biology http://bioconductor.org __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Problem/bug with smooth.spline and all.knots=T

H == Hubertus [EMAIL PROTECTED] on Wed, 4 Jul 2007 14:58:44 +0200 writes: H Dear list, H if I do H smooth.spline(tmpSec, tmpT, all.knots=T) H with the attached data, Thanks for providing the data via URL (see below) H with the attached data, I get this error-message:

Re: [R] Lookups in R

. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Martin Morgan Bioconductor

[R] speed and looping issues; calculations on big datasets

this function for matrices of size 10^2 x 10^5 or even more. so i would be very thankful for any ideas, suggestions etc to improve this cheers martin __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] align() function missing in R ?

'') than the S-plus one. Note however that you'd typically not be allowed to copy the S-plus implementation. Martin ML == Markus Loecher [EMAIL PROTECTED] on Thu, 28 Jun 2007 11:10:51 -0400 writes: ML Dear list members, I switched from Splus to R a few ML years ago and so far found

Re: [R] Loading problem with XML_1.9

'external pointer', or whether the savings is just in the R representation of the document. Also, depending on how you use the resulting document, you might want to watch out for the memory leak mentioned in http://www.omegahat.org/RSXML/Changes Martin Luo Weijun [EMAIL PROTECTED] writes: Hello

[R] CORRECTION: Re: BioC2007, August 6-7 ** early registration ends July 1 **

BioC2007: Where Software and Biology Connect August 6-7 in Seattle, WA, USA The BioC2007 developer-focused meeting is on August 8, the day *after* the main conference. In addition, I neglected to mention the poster session on Monday evening. Best, Martin Morgan -- Bioconductor / Computational

[R] BioC2007, August 6-7 ** early registration ends July 1 **

is on August 5th. For (preliminary) details, please visit: http://wiki.fhcrc.org/bioc/Seattle_Dev_Meeting_2007 See you in August! Martin Morgan -- Bioconductor / Computational Biology http://bioconductor.org __ R-help@stat.math.ethz.ch mailing list

Re: [R] a string to enviroment or function

be helpful (though AnnotationDbi is only available with R-devel and revmap seems not to be documented). res - revmap(hgu95av2ENTREZID) hgu95av2ENTREZID[[1190_at]] [1] 5800 res[[5800]] [1] 1190_at 32199_at Martin Weiwei Shi [EMAIL PROTECTED] writes: then how to do this f1 - function(mylab

Re: [R] extract index during execution of sapply

, minimal, self-contained, reproducible code. -- Martin Morgan Bioconductor / Computational Biology http://bioconductor.org __ 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] Boxplot issues

(ideally cut-n-pastable) R code examples. Regards, Martin Maechler SE Steve Ellison __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Matrix *package*, CHOLMOD error: problem too large

part of the Math S4 group generic. I'd say that there's only historical reason for them *not* to be part of Math, and I am likely going to propose to change this Martin Maechler DM Duncan Murdoch __ R-help@stat.math.ethz.ch mailing list https

Re: [R] if within a function

for defensive programming where I would not use them in simple one liners, but would use them otherwise ) Martin Maechler, ETH Zurich __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Could not find lmer function in {Matrix} package

) to get to lmer. Regards, Martin SB I have attempted these actions in both MacOSx R Versions 2.4.1 and SB 2.5.0. I have also tried this in Version 2.5.1. beta on a Windows SB machine. SB Thanks in advance for any suggestions. SB Steve

Re: [R] Matrix library error: should never happen; please report

; [and then do report to [EMAIL PROTECTED] JQ Thanks -- Jose Quesada, PhD. Best regards, Martin Maechler, ETH Zurich __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

Re: [R] Odp: Odp: outlying

$rank) (Tcrit - qt(0.995, df = df.resid)) ## Std.error of mean ~= sqrt(1/n Var(X_i)) = 1/sqrt(n) sqrt(Var(X_i)) cc[,1] + c(-1,1) * sqrt(n) * Tcrit * cc[,Std. Error] ## -6.391201 18.555177 --- Martin Maechler, ETH Zurich __ R-help@stat.math.ethz.ch

Re: [R] Augment 'Matrix' matrices

SH == Scott Hyde [EMAIL PROTECTED] on Mon, 18 Jun 2007 16:59:00 -1000 (HST) writes: SH Martin, How does Matrix implement augmented matrices? I SH tried this and got the expected result: {Replying to R-help, since this question has come up several times } V=matrix(1,2,3

Re: [R] How do I avoid a loop?

is about ten times faster than Gabor's on my machine (time to compile the C code excluded!). Regards, Martin On 6/19/07, Feng, Ken [EMAIL PROTECTED] wrote: Hi, I start with an array of booleans: x - c( TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, TRUE, FALSE, TRUE ); I want

Re: [R] Rmpi and rsprng for Windows

Eric - http://www.stats.uwo.ca/faculty/yu/Rmpi/ (the Rmpi package author page) has helpful Windows instructions. Martin Eric Ferreira [EMAIL PROTECTED] writes: Dear f_R_iends, I'm new on parallel programming and trying to use a machine with Windows to access a linux computer cluster. I

Re: [R] How do I avoid a loop?

ten times faster than Gabor's on my machine (time to compile the C code excluded!). Regards, Martin On 6/19/07, Feng, Ken [EMAIL PROTECTED] wrote: Hi, I start with an array of booleans: x - c( TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, TRUE, FALSE, TRUE ); I want to define an y

Re: [R] Responding to a posting in the digest

is TH to start to compose a new message and copy what you need TH from the digest. TH While I've never reveived R-help in digest format TH myself, according to Martin Maechler: TH http://finzi.psych.upenn.edu/R/Rhelp02a/archive/59429.html TH Please open the URL

Re: [R] fSeries - Ox - ver: 240.10068 - Steps to make it work

= TRUE) in .garchOxFit() (fSeries version 240.10068) should read write(x = series, file = OxSeries.csv, ncolumns = 1, append = TRUE) instead. Incorporating the changes for [EMAIL PROTECTED] could be a good occasion to fix this as well :-) Regards, Martin Ian Gregory wrote: -Bugs

Re: [R] fSeries - Ox - ver: 240.10068 - Steps to make it work

' (in your current workspace) instead of 'IanSeries' (which is even worse than aborting...). Regards, Martin Ian Gregory schrieb: I tried the following and it works for me (after the changes to make): library(fSeries) data(dem2gbp) IanSeries = dem2gbp[, 1] garchOxFit(~garch(1,1),IanSeries) Any

Re: [R] Package update announcements

MHHS == Martin Henry H Stevens [EMAIL PROTECTED] on Mon, 11 Jun 2007 17:23:46 -0400 writes: MHHS Hi Folks, I was wondering what everyone thought about MHHS adding a sentence to each package update announcement MHHS that described what the package did. R extensions are MHHS

Re: [R] Package update announcements

Hi Martin, Maybe we should appeal to each others' vanity (Just think how many hundreds---even thousands---of people will use your package once they know about it) Hmmm. Maybe that would have the opposite of the intended effect. Cheers, Hank On Jun 12, 2007, at 3:16 AM, Martin Maechler

Re: [R] {spam?} Nonlinear Regression

Hi tronter, PLEASE do read the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code. Hank On Jun 7, 2007, at 5:50 PM, tronter wrote: Hello I followed the example in page 59, chapter 11 of the 'Introduction to R'

[R] autoregressive spectral density estimate by andrews' plug-in method?

coefficients at frequency zero. However, I have a time series of length 160 and I need to estimate its spectral density via Andrews methodology. Any suggestions will be appreciated. Excuse me if I am asking something obvious. Regards, Martin

Re: [R] Fwd: Using odesolve to produce non-negative solutions

Hi Jeremy, First, setting hmax to a small number could prevent a large step, if you think that is a problem. Second, however, I don't see how you can get a negative population size when using the log trick. I would think that that would prevent completely any negative values of N (i.e.

Re: [R] Fwd: Using odesolve to produce non-negative solutions

Hi Spencer, I have copied Woody Setzer. I have no idea whether lsoda can estimate parameters that could take imaginary values. Hank On Jun 11, 2007, at 12:52 PM, Spencer Graves wrote: in line Martin Henry H. Stevens wrote: Hi Jeremy, First, setting hmax to a small number could prevent

[R] Package update announcements

Hi Folks, I was wondering what everyone thought about adding a sentence to each package update announcement that described what the package did. R extensions are so numerous that it is difficult to keep up with them. Would it be appropriate to ask package developers to add a brief sentence

Re: [R] pnorm how to decide lower-tail true or false

really want q = 1 - pfoo(x) = Prob_{foo}[X x] then you can get this directly via q - pfoo(x, lower.tail = FALSE, ) Simple example with R : pnorm(10) [1] 1 1 - pnorm(10) [1] 0 pnorm(10, lower.tail=FALSE) [1] 7.619853e-24 Regards, Martin Maechler, ETH Zurich

Re: [R] Tools For Preparing Data For Analysis

Is there an example available of this sort of problematic data that requires this kind of data screening and filtering? For many of us, this issue would be nice to learn about, and deal with within R. If a package could be created, that would be optimal for some of us. I would like to

Re: [R] Why is the R mailing list so hard to figure out?

its silly footer (aka advertizement) to all mails sent out, luring people to use it. I (as mailing list manager) have considered more than once to just evaporate such footers as I evaporate most other spam-footers (as e.g. from hotmail, yahoo, gmx, ...). Martin Maechler, ETH Zurich VE

Re: [R] Interaction term in lmer

HI Emine, The parentheses specify a variable as random or fixed. You are using center for random intercepts, but fixed for interactions. You may want mer(cbind( yvect, nvect-yvect) ~ trt + ( 1 | center) + (1 | trt:center), family = binomial, niter = 25, method = Laplace, control = list

Re: [R] opinions please: text editors and reporting/Sweave?

as in free beer not as in free speech. In other words, Winshell is *not* 'Free Software' / 'Software Libre' nor is it Open Source Software. Tinn-R, as R itself, *is* Free Software (and Emacs and ESS and (La)TeX are too). Martin Maechler, ETH Zurich Jared On 5/31/07, Duncan Murdoch [EMAIL PROTECTED

Re: [R] determining a parent function name

Hi sundar -- maybe myerr - function(err) err$call foo - function() stop() tryCatch({ foo() }, error=myerr) foo() suggests a way to catch errors without having to change existing code or re-invent stop? Martin Sundar Dorai-Raj [EMAIL PROTECTED] writes: Hi, Vladimir, Sorry, didn't see

[R] Different fonts on different axes

Hi Folks, How do I get red bold font on my y axis and black standard font on my x axis? plot(runif(10), ylab=Red, Bold?, xlab=Black, standard?) Any pointers or examples would be great. Thanks! Hank Dr. Hank Stevens, Assistant Professor 338 Pearson Hall Botany Department Miami University

Re: [R] Different fonts on different axes

: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Martin Henry H. Stevens Sent: Thursday, May 31, 2007 11:00 AM To: R-Help Subject: [R] Different fonts on different axes Hi Folks, How do I get red bold font on my y axis and black standard font on my x axis? plot(runif

Re: [R] Aggregate to find majority level of a factor

How about tapply? plot - gl(2,3); plot type - letters[c(1,2,2,1,1,1)]; type tapply(type, list(plot), function(x) {tabl - table(x) names(tabl[tabl==max (tabl)])}) Hank On May 31, 2007, at 3:25 PM, Thompson, Jonathan wrote: I want to use the

Re: [R] How to check for existence url from within a function?

returns a logical vector, typically c(FALSE, FALSE) and using that in if(.) gives at least a warning. Martin __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

[R] about the unscaled covariances from a summary.lm object

: solve(t(X)%*%X), but I get a matrix quite different from the matrix given by cov.unscaled. Is it just computational instability, or I am missing something important? Regards, Martin __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

[R] exemples, tutorial on lmer

Hi all, I have some difficulties to work with the function lmer from lme4 Does somebody have a tutorial or different examples to use this function? Thanks, Oliver. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] exemples, tutorial on lmer

Hi Oliver, You could start with R News 2005, no. 1. Also the PDF associated with lme4, Implementation.pdf. Hnak On May 29, 2007, at 10:35 AM, Olivier MARTIN wrote: Hi all, I have some difficulties to work with the function lmer from lme4 Does somebody have a tutorial or different examples

Re: [R] ratio distribution - missing attachment

(Doug Bates uses foo_R.txt) and it should make it to the lists; as a third alternative, just cut paste the corresponding text into your e-mal. I think your R function should make it to R-help (and its archives), so I'd be thankful for a repost. Martin Ravi == Ravi Varadhan [EMAIL PROTECTED

Re: [R] File path expansion

Duncan == Duncan Murdoch [EMAIL PROTECTED] on Fri, 25 May 2007 15:38:54 -0400 writes: Duncan On 5/25/2007 1:09 PM, Prof Brian Ripley wrote: On Fri, 25 May 2007, Martin Maechler wrote: path.expand(~) [1] /home/maechler Yes, but beware that may not do what

Re: [R] normality tests [Broadcast]

: It's a very bad idea ... Martin Maechler, ETH Zurich LuckeJF -Original Message- From: LuckeJF [EMAIL PROTECTED] LuckeJF [mailto:[EMAIL PROTECTED] On Behalf LuckeJF Of Liaw, Andy Sent: Friday, May 25, 2007 12:04 PM LuckeJF To: [EMAIL PROTECTED]; Frank E Harrell Jr Cc

[R] lmer and scale parameter in glmm model

Hi all, I try to fit a glmm model with binomial distribution and I would to verify that the scale parameter is close to 1... the lmer function gives the following result : Estimated scale (compare to 1 ) 0.766783 But I would like to know how this estimation (0.766783) is performed, and I

Re: [R] Make check failure for R-2.4.1

Adam == Adam Witney [EMAIL PROTECTED] on Fri, 25 May 2007 09:38:29 +0100 writes: Adam Thanks for your replies Details inline below: Adam On 24/5/07 17:12, Martin Maechler [EMAIL PROTECTED] wrote: UweL == Uwe Ligges [EMAIL PROTECTED] on Thu, 24 May 2007 17:34:16 +0200

Re: [R] Make check failure for R-2.4.1

) : Adam object digits not found well, that's a typo - I think - you should have been able to fix (I said something like ...). Just do replace the '==' by '=' Martin Adam Thanks again... Adam adam __ R-help@stat.math.ethz.ch mailing list https

Re: [R] File path expansion

path.expand(~) [1] /home/maechler RobMcG == McGehee, Robert [EMAIL PROTECTED] on Fri, 25 May 2007 11:44:27 -0400 writes: RobMcG R-Help, RobMcG I discovered a mis-feature is ghostscript, which is used by the bitmap RobMcG function. It seems that specifying file names in the

Re: [R] windows to unix

*the windows-non-expert ..} Martin Erin thanks in advance, Erin Sincerely, Erin Erin Hodgess Erin Associate Professor Erin Department of Computer and Mathematical Sciences Erin University of Houston - Downtown Erin mailto: [EMAIL PROTECTED

Re: [R] I made some progress on my previous systemfit question but still not quite there

be able to select the equations via [], e.g., eqns=eqn[1:2] for the first and second equation. HTH, Martin I would want demand to be AUD, supply to be CHF, etc. If someone has time to run the code, they would see what I mean. Thanks a lot

Re: [R] Make check failure for R-2.4.1

? Martin UweL Uwe Ligges UweL Adam Witney wrote: I'm trying to install R-2.4.1, everything configure's and make's OK, but the make check fails: running code in 'd-p-q-r-tests.R' ...make[3]: *** [d-p-q-r-tests.Rout] Error 1 make[3]: Leaving directory `/usr

Re: [R] Question about setReplaceMethod

Hi Adschai -- You'll want to return the value whose slot you have modified: setReplaceMethod(setX, foo, function(this,value) { [EMAIL PROTECTED] - value this # add this line }) Martin [EMAIL PROTECTED

  1   2   3   4   5   6   7   8   9   10   >