Re: [R] Help on GAM

I want to do this model using mgcv caco1.gam-gam(caco1$Pa~s(caco1$T10)+s(caco1$S10)+s(caco1$C10), data=caco1) But I want to indicate the degrees of freedom of each term, how I have to do it? Then, I saw that always the results refered to positive interactions, there are not negative ones?

Re: [R] random uniform sample of points on an ellipsoid (e.g. WG

Following up to Russell Senior's original query: On 21 Feb 2007, Russell Senior wrote: I am interested in making a random sample from a uniform distribution of points over the surface of the earth, using the WGS84 ellipsoid as a model for the earth. I know how to do this for a sphere, but

[R] absent

Je serai absent(e) à partir du 01/03/2007 de retour le 05/03/2007. Bonne année 2007 Si vous n'etes pas destinataires de ce message, merci d'aver...{{dropped}} __ R-help@stat.math.ethz.ch mailing list

[R] HELP GAM

I think it will be better if I explained what I want to do. I have the abundance of Penilia (zooplankter) as y and I have the temperature, salinity and chlorophyll concentration at different depths. I want to know which how the abundance of Penilia is modulate by those environmetal variables and

[R] stl cycle sub-series plot

Dear R-people, I tried to build a seasonal cycle sub-series plot of my time series using: monthplot(nitrat.stl, choice=seasonal) However, I only get the horizontal lines of the mean values and not the vertical lines for each year. Can anybody help me? Sincerely, Anja --

[R] Matrix Library failing to load

Dear R users I am trying to load the package lmer4 but it seems the Matrix library is failing to load and giving me the following error message : library(Matrix) Error in importIntoEnv(impenv, impnames, ns, impvars) : object 'Logic' is not exported by 'namespace:methods' Error:

[R] pheno package

Dear R-people, I tried to use a command seqMK(x) of the pheno package. However, the installation of the package fails with following error message: Error in library(pkg, character.only = TRUE) : 'pheno' is not a valid package -- installed 2.0.0? What is the matter? Can anybody help

[R] How to read in this data format?

Hi, I recieved an ascii file, containing following information: $$ Experiment Number: $$ Associated Data: FUNCTION 1 Scan1 Retention Time 0.017 399.8112184 399.87420 399.9372152 Scan2 Retention Time 0.021 399.8112181 399.8742

Re: [R] how to apply the function cut( ) to many columns in a data.frame?

ahimsa campos-arceiz wrote: Dear useRs, In a data.frame (df) I have several columns (x1, x2, x3xn) containing data as a continuous numerical response: df var x1x2 x3 1143 147 137 2 9393 117 316439 101 4123 11897 5

Re: [R] RDA and trend surface regression

Dear All, PCNM and more general MEM can be found in the package spacemakeR available here (http://biomserv.univ-lyon1.fr/~dray/). As said by Jari, I do not want to put the package on CRAN. I have no no problems to include it in another package (vegan or spdep). It has a vignette which can be

[R] Reshape data

I have two data frames with a number of common variables that I wish to put into long format. The first contains the variables es(Csales) [1] terr Dec.02 Jan.03 Feb.03 Mar.03 Apr.03 May.03 Jun.03 [9] Jul.03 Aug.03 Sep.03 Oct.03 Nov.03 Dec.03 Jan.04 Feb.04 [17] Mar.04 Apr.04 May.04

[R] R (Input from Keyboard): How do I remove error?

Hi, I am trying to read some value from keyboard in one of my R program. I want to do manipulation on data based on this value. But I am facing a problem which is described below with the help of an example: My code is similar to: : : cat(\n,Enter value:: ,\n) y-scan(n=1) a - b*y cat(\n,Enter

[R] function with Multiple Output

Dear R user, I have a simple question for you: I have created a global function that evoke other subsidiary functions. when I run the global function I want to store the outcomes of the of each subsidiary function into a variables. an examples my.fun=function(vector, index){ a=fun.a(vector,

Re: [R] R (Input from Keyboard): How do I remove error?

?readline Note that this gives you a character vector, so you might want to use as.numeric() after that. Petr d. sarthi maheshwari napsal(a): Hi, I am trying to read some value from keyboard in one of my R program. I want to do manipulation on data based on this value. But I am facing a

Re: [R] How to read in this data format?

Well, not extremely elegant, but should work: 1) open your file in some ascii text editor, delete the rubbish at the beginning up to line Scan 1, and replace all spaces in names - e.g. make a mass replace of 'Retention Time' by let say 'RetentionTime'. 2) Use read.table(), matrix() and

[R] object is not subsettable

Dear colleagues, I've just come across a problem with the following command which is a part of the metaOverview.R code file provided as an monography- accompanying file at http://www.bioconductor.org/docs/mogr/metadata: ## R hasChr - eapply(GOTERM, function(x)

Re: [R] interactions and GAM

Don't know if this helps, but... gam in package mgcv will let you set up smooths that interact with factors using the `by' variable mechanism. See ?gam.models, particularly the last example. Prediction is supported. Simon -- Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2

[R] barplot2, gap.barplot

Hello, I try to handle a simple bar-plot, but it turns out to be not as simple as I thought. 1) I have created a .dat-File, e.g. test.dat: DATADATA-SEM 2.2 0.32 6.2 1.30 12.71.61 48.63.08 4.1 0.86 4.5 0.32 1.5 1.13 1.2 1.08 The first row is the data

Re: [R] Help on GAM

On Thursday 01 March 2007 08:38, Dacha Atienza wrote: I want to do this model using mgcv caco1.gam-gam(caco1$Pa~s(caco1$T10)+s(caco1$S10)+s(caco1$C10), data=caco1) But I want to indicate the degrees of freedom of each term, how I have to do it? e.g. s(caco1$T10,k=6,fx=TRUE) will give you a 5

Re: [R] legend question

? par it is the xpd you're looking for. x - seq(-pi, pi, len = 65) par(xpd=TRUE) plot(x, sin(x), type=l, col = 2,xpd=NA) legend(x = 0, y = -1.5, legend text, pch = 1, xjust = 0.5) --- Jenny Barnes [EMAIL PROTECTED] wrote: Hi folks, Do you mind if I ask a related question that I have

[R] need indexing help again

Dear all, i have a problem which deals with indexing again, i guess. I am a beginner and though i read the S Poetry stuff it is not that easy for me to use the examples for my problems, maybe all this needs some practice... Here´s the question. As i want to perform a factor analysis i need

Re: [R] random uniform sample of points on an ellipsoid (e.g. WG

Hello, I do not know if I am completely out of it but . . . if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the equation of an ellipsoid and A = T'T (cholesky) then T.[x,y,z] should be a point in the ellipsoid ? isn't it? Heberto Ghezzo McGill University Montreal Canada

Re: [R] barplot2, gap.barplot

On Thu, 2007-03-01 at 14:56 +0100, Marc A. Rohling wrote: Hello, I try to handle a simple bar-plot, but it turns out to be not as simple as I thought. 1) I have created a .dat-File, e.g. test.dat: DATA DATA-SEM 2.2 0.32 6.2 1.30 12.7 1.61 48.6 3.08 4.1 0.86 4.5 0.32

Re: [R] random uniform sample of points on an ellipsoid (e.g. WG

R Heberto Ghezzo, Dr wrote: I do not know if I am completely out of it but . . . if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the equation of an ellipsoid and A = T'T (cholesky) then T.[x,y,z] should be a point in the ellipsoid ? isn't it? Yes, it's a point _on_ the

[R] count the # of appearances...

Hi there, is there a possibility to count the number of appearances of an element in a vector ? i mean of any given element.. deliver all elements which are exactly xtimes in this vector ? thx in advance !! __ R-help@stat.math.ethz.ch mailing list

[R] question about returning Random Effects' covariance matrix estimate using lme fitting

Dear all I am fitting and analyzing linear mixed-effects models using the R command 'lme'. The following is the results: dental.fit - lme(fixed = distance~age, random = ~age + cluster = ~subject, data = dental) summary(dental.fit) Variance/Covariance Components Estimates:

Re: [R] random uniform sample of points on an ellipsoid (e.g. WG

On 01-Mar-07 R Heberto Ghezzo, Dr wrote: Hello, I do not know if I am completely out of it but . . . if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the equation of an ellipsoid and A = T'T (cholesky) then T.[x,y,z] should be a point in the ellipsoid ? isn't it? Heberto Ghezzo

Re: [R] count the # of appearances...

Hello Bunny/Lautioscrew, sum(your vector == your chosen element) should do what you want... HTH, Ranjan On Thu, 1 Mar 2007 15:20:19 +0100 bunny , lautloscrew.com [EMAIL PROTECTED] wrote: Hi there, is there a possibility to count the number of appearances of an element in a vector ? i

Re: [R] count the # of appearances...

?table Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://med.kuleuven.be/biostat/

Re: [R] count the # of appearances...

Dear Bunny, ? table might be what you wish. Best wishes Wolfgang Huber EBI bunny , lautloscrew.com wrote: Hi there, is there a possibility to count the number of appearances of an element in a vector ? i mean of any given element.. deliver all elements which are exactly xtimes

[R] time series data ploting

hello, I have plots http://uosis.mif.vu.lt/~roka5178/time series.JPG I must to do almost the same. I do not have any data. but i know that this is time series data with parabolic changeable mean. I can not generate tie series data with parabolic mean. how can i do this? and have

Re: [R] count the # of appearances... [Broadcast]

Take a look at table. Hope this helps, Matt -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of bunny , lautloscrew.com Sent: Thursday, March 01, 2007 9:20 AM To: r-help@stat.math.ethz.ch Subject: [R] count the # of appearances... [Broadcast] Hi there, is

[R] question about xtable and Hmisc

I would like to get rid of the row numbers using xtable and latex. The commands d = cbind(1:10,rep(1:2,5)) ans = xtable(d) latex(ans) gives output containing \begin{tabular}{rrr} \hline 1 2 \\ \hline 1 1.00 1.00 \\ 2 2.00 2.00 \\ 3 3.00 1.00 \\ 4 4.00 2.00 \\

Re: [R] count the # of appearances...

On Thu, 2007-03-01 at 15:20 +0100, bunny , lautloscrew.com wrote: Hi there, is there a possibility to count the number of appearances of an element in a vector ? i mean of any given element.. deliver all elements which are exactly xtimes in this vector ? thx in advance !! Vec -

Re: [R] how to apply the function cut( ) to many columns in a data.frame?

From: Chuck Cleland ahimsa campos-arceiz wrote: Dear useRs, In a data.frame (df) I have several columns (x1, x2, x3xn) containing data as a continuous numerical response: df var x1x2 x3 1143 147 137 2 9393 117 316439 101

Re: [R] count the # of appearances...

bunny , lautloscrew.com wrote: Hi there, is there a possibility to count the number of appearances of an element in a vector ? i mean of any given element.. deliver all elements which are exactly xtimes in this vector ? thx in advance !! myvec - c(0, rep(1:6, each = 2), 7, 7, 7)

[R] Another newbie book recommandation question

I hope this question is sufficiently different from the other requests for book recommendations that it's not repetitious. If not, I apologize in advance. I'm curious what standard reference books working statisticians, or biostatisticians, have within easy reach of their desk. I'm a computer

Re: [R] count the # of appearances...

On 01-Mar-07 bunny , lautloscrew.com wrote: Hi there, is there a possibility to count the number of appearances of an element in a vector ? i mean of any given element.. deliver all elements which are exactly xtimes in this vector ? thx in advance !! If it is a specific element

Re: [R] Another newbie book recommandation question

Oh, Boy. This might result in a data dump since each of us has a personal library. Here are the top dozen or so from mine: 1. Agresti, Alan, Categorical Data Analysis, 2nd ed., Wiley, 2002 2. Box, George E. P., William G. Hunter, and J. Stewart Hunter, Statistics for

Re: [R] question about xtable and Hmisc

steve fisk at bowdoin.edu writes: I would like to get rid of the row numbers using xtable and latex. but I don't want the row numbers. Is it possible to get rid of them? Also, if x is a data frame, latex(x) contains the row numbers. Can I get rid of them here as well? Steve Steve,

Re: [R] Another newbie book recommandation question

for the size of your data file, I think R can handle it. of course, it also depends on your hardware. however, it might not be a good idea to do heavy data manipulation work in R. stata has very good routine for survey analysis. i am not sure if R is as good as stata in terms of survey analysis.

Re: [R] How to read in this data format?

You can't expect general-purpose tools like read.table in R to be able to deal with highly specialized file format. Here's what I'd start. It doesn't put data in the format you specified exactly, but I doubt you'll need that. This might be sufficient for your purpose: dat -

[R] Creating a vector of variable bin widths

Dear R users, I am having a little trouble with grouping data. ---Detailed explanation (summary below) A small sample of my data is below (which has already been rounded and grouped a little from the raw data for clarity). I am sampling data from an unknown

Re: [R] PROC TABULATE with R

You can also have a look at structable() in vcd, especially the indexing functions. The problem with OLAP in R is, that you will have to create sth. like a hierarchical factor to handle rollup/drilldown correctly. And you will have to decide whether it's purely memory-based (fast calculations,

[R] Query about data manipulation

Hi Thanks much for the prompt response to my earlier enquiry on packages for regression analyses. Along the same topic(?), I have another question about which I could use some input. I am retreiving data from a MySQL database using RODBC. The table has many BLOB columns and each BLOB column has

[R] Fit Student Copula

Hello everybody, I have a big problem that I do not manage to solve ! I will be very grateful if you can solve this ! I want to fit a t Copula with the copula package : student.cop - ellipCopula(t, param = c(0.5, 0.6, 0.7), dim = 3, dispstr = un,df=5) x-rcopula(student.cop,1000) fit -

[R] HELP on GAM MODELS

I think it will be better if I explained what I want to do. I have the abundance of Penilia (zooplankter) as y and I have the temperature, salinity and chlorophyll concentration at different depths. I want to know which how the abundance of Penilia is modulate by those environmetal variables and

Re: [R] Query about data manipulation

On Thu, 2007-03-01 at 08:34 -0800, lalitha viswanath wrote: Hi Thanks much for the prompt response to my earlier enquiry on packages for regression analyses. Along the same topic(?), I have another question about which I could use some input. I am retreiving data from a MySQL database

Re: [R] Double-banger function names: preferences and suggestions

Definitely not #2. Prefer #1 but #3 is ok as well. Thanks for contributing and inquiring. - Original Message - From: hadley wickham [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Sunday, February 25, 2007 7:44 AM Subject: [R] Double-banger function names: preferences and suggestions

[R] setting font in plots

Dear Reader I am trying to change the font in a plot and after several trials finally came up with the following code. plot(var_a, var_b, pch = 16, font.lab = 10, font = 10) points(var_a, var_c, pch = 3, font = 10) legend(0.1, 0.8, legend = c(“var_b”, “var_c”), pch = c(16,3)) It does change

[R] Where to find ROracle for R 2.2.1

Hello I am looking for ROracle - Module but found only a version that can installed with R Version 2.30. I am running a production system where I can't change the currently installed R. That's why I need a ROracle - Module that is compatible with R version 2.2.1 Thanks Reinhard

[R] matrices dimensions limitation

Hello, I have several questions around matrices size limitation: can i create a matrix with 70 rows? does it depends on the number of columns? if so can Is there a way to ask to R, given the number of columns, the max number of rows I can have? I need in fact to do regression

[R] Row-wise two sample T-test on subsets of a matrix

Hello all, I am trying to run a two sample t-test on a matrix which is a 196002*22 matrix. I want to run the t-test, row-wise, with the first 11 columns being a part of the first group and columns 12-22 being a part of the second group. I tried running something like (temp.matrix being my

[R] R File IO Slow?

Is R file IO slow in general or am I missing something? It takes me 5 minutes to do a load(MYFILE) where MYFILE is a 27 MB Rdata file. Is there any way to speed this up? The one idea I have is having R call a C or Perl routine, reading the file in that language, converting the data in to R

Re: [R] Double-banger function names: preferences and suggestions

I agree with Jason: !2, prefer 1, can accept 3. Clint BowmanINTERNET: [EMAIL PROTECTED] Air Dispersion Modeler INTERNET: [EMAIL PROTECTED] Air Quality Program VOICE: (360) 407-6815 Department of Ecology FAX:(360)

Re: [R] How to read in this data format?

Read in the data using readLines, extract out all desired lines (namely those containing only numbers, dots and spaces or those with the word Time) and remove Retention from all lines so that all remaining lines have two fields. Now that we have desired lines and all lines have two fields read

Re: [R] R File IO Slow?

Just an idea: Two things that can slow down save()/load() is if you save() in ASCII format or a compressed binary format. If this is your case for MYFILE, try to resave in a non-compressed binary format. See ?save for details. /HB On 3/1/07, ramzi abboud [EMAIL PROTECTED] wrote: Is R file IO

[R] GLM problems.

Dear R users I'm new in R management and maybe It's a silly question. I'm working with GLM to obtain predictive models. I have some problesm with the prediction instruction: DatosTotal - read.csv(Var_perdizcsv.csv, sep =;) edvariable - edit(DatosTotal) pre - predict(rlfinal, DatosTotal, type

Re: [R] R File IO Slow?

A 27MB .RData file is relatively big, in may experience. What do you think is slow? Maybe it's your computer that is slow? -roger ramzi abboud wrote: Is R file IO slow in general or am I missing something? It takes me 5 minutes to do a load(MYFILE) where MYFILE is a 27 MB Rdata file. Is

Re: [R] object is not subsettable

Hi Radek, Radek Blatny [EMAIL PROTECTED] writes: Dear colleagues, I've just come across a problem with the following command which is a part of the metaOverview.R code file provided as an monography- accompanying file at http://www.bioconductor.org/docs/mogr/metadata: The appropriate

[R] number of levels for a factor

I have temp list which have 19 data.frame I want to get number of levels for pr in the first dat.frame I do this like this: temp[[1]]$pr just has 1A24 after I do nlevels(temp[[1]]$pr) I expect to get 1, but I get 19 anyone know why? tail(temp[[1]]$pr) [1] 1A24 1A24 1A24 1A24 1A24 1A24 19

Re: [R] R File IO Slow?

It is not slow on my system. The file was 34MB on disk and took about 37 seconds to write out (probably mostly disk I/O on my laptop) and 12 seconds to read in after I flushed the system cache. x - runif(27e6/4) # creates a 34MB file on disk object.size(x) [1] 5424

Re: [R] setting font in plots

On Thu, 2007-03-01 at 23:48 +1000, Mark Allen wrote: Dear Reader I am trying to change the font in a plot and after several trials finally came up with the following code. plot(var_a, var_b, pch = 16, font.lab = 10, font = 10) points(var_a, var_c, pch = 3, font = 10) legend(0.1, 0.8,

Re: [R] setting font in plots

Hi Mark Allen wrote: Dear Reader I am trying to change the font in a plot and after several trials finally came up with the following code. plot(var_a, var_b, pch = 16, font.lab = 10, font = 10) points(var_a, var_c, pch = 3, font = 10) legend(0.1, 0.8, legend = c(“var_b”, “var_c”),

Re: [R] matrices dimensions limitation

On Thu, 2007-03-01 at 18:20 +0100, Bruno C. wrote: Hello, I have several questions around matrices size limitation: can i create a matrix with 70 rows? does it depends on the number of columns? if so can Is there a way to ask to R, given the number of columns, the max

Re: [R] How to read in this data format?

Here is one way of doing it using the reshape package: # test data from email x - $$ Experiment Number: + $$ Associated Data: + + FUNCTION 1 + + Scan 1 + Retention Time 0.017 + + 399.8112 184 + 399.8742 0 + 399.9372 152 + + + Scan 2 + Retention Time 0.021 + + 399.8112 181 + 399.8742

Re: [R] Row-wise two sample T-test on subsets of a matrix

Here's one suggestion: convert the matrix into a three-dimensional array and use apply on it. Ranjan On Thu, 1 Mar 2007 11:51:29 -0600 (CST) Nameeta Lobo [EMAIL PROTECTED] wrote: Hello all, I am trying to run a two sample t-test on a matrix which is a 196002*22 matrix. I want to run the

Re: [R] GLM problems.

Daniel Jimenez wrote: Dear R users I'm new in R management and maybe It's a silly question. I'm working with GLM to obtain predictive models. I have some problesm with the prediction instruction: DatosTotal - read.csv(Var_perdizcsv.csv, sep =;) edvariable - edit(DatosTotal) pre -

Re: [R] R File IO Slow?

On Thu, 2007-03-01 at 09:22 -0800, ramzi abboud wrote: Is R file IO slow in general or am I missing something? It takes me 5 minutes to do a load(MYFILE) where MYFILE is a 27 MB Rdata file. Is there any way to speed this up? The one idea I have is having R call a C or Perl routine,

[R] problem with throwing lines out of matrix

Dear all, again i have a problem with special case of dropping lines out of a matrix. i tried the following, where throwout is a vector of length 18 with simple id values that should be compared to any for (k in 1:length(throwout)) { mymatrix=matrix(mymatrix[-(forfact[k]),],ncol=4)

Re: [R] R File IO Slow?

I decided to run an experiment: just reading in a file which is 78MB in binary format (of ints). It takes less than 30s using a laptop with 512 MB RAM, 2.3 GHz Intel-4 single processor. At that point, I did not notice that Ramzi was talking about a .RData file. For huge files, I usually do not

Re: [R] Row-wise two sample T-test on subsets of a matrix

I'd suggest to use the function mt.teststat in the package multtest or rowttests in the package genefilter. Both can be found athe the bioconductor webpage (www.bioconductor.org) Stef On Thu, Mar 01, 2007 at 12:48:29PM -0600, Ranjan Maitra wrote: RanjanHere's one suggestion: convert the

Re: [R] problem with throwing lines out of matrix [Broadcast]

You didn't' actually tell us what forfact is. You could try something like (untested) My.matrix - my.matrix[!(my.matrix[, 1] %in% remove.values), , drop = FALSE] Also, reading some of the introductory R manuals you can find on CRAN in the documentation section will really give you a leg up on

Re: [R] count the # of appearances...

You can split() the original vector according to its elements into a list, then use lapply() to count the lengths of the list: x - round(runif(100,1,10)) unlist(lapply(split(x,f=x),length)) 1 2 3 4 5 6 7 8 9 10 5 14 9 13 11 8 8 14 12 6 This can be useful if you want to use some

[R] Default par() options

The following question/idea came up on the RKWard development mailing list, but might be of general interest: Is there a nice way to customize the default look of all graphs on all devices? I.e. a way to - for instance - set the following options before each plot: par(bg=light gray, las=2,

Re: [R] question about xtable and Hmisc

Unfortunately, this applies to print.xtable, and not to latex. I want to know how to eliminate them using latex() Greg Johnson wrote: steve fisk at bowdoin.edu writes: I would like to get rid of the row numbers using xtable and latex. but I don't want the row numbers. Is it possible to

Re: [R] Fit Student Copula

set.seed(1) student.cop - ellipCopula(t, param = c(0.5, 0.6, 0.7), dim = 3, dispstr = un,df=5) x-rcopula(student.cop,1000) fit - fitCopula(x, student.cop, c(0.5,0.5,0.5,5)) Warning messages: 1: NaNs produced in: qt(p, df, lower.tail, log.p) 2: NaNs produced in: sqrt((df + d)/(Q + df)) 3:

[R] [R-pkgs] Update of X2R sent to CRAN

A new version of X2R has just been uploaded to CRAN. It should be available at mirrors within a few days. This contains revisions to the For2R component to fix a bug in which data frames were not written correctly if the user did not pass row labels. The new version is supplied as files

Re: [R] Query about data manipulation

Use 'strsplit' to create a vector and then 'rbind' into a matrix or create a dataframe: x - c(id1\tid2\tid3\tid4, 1\t2\t3\t4, 5\t6\t7\t8) x [1] id1\tid2\tid3\tid4 1\t2\t3\t4 5\t6\t7\t8 (y - strsplit(x, \t)) [[1]] [1] id1 id2 id3 id4 [[2]] [1] 1 2 3 4 [[3]] [1] 5 6 7 8

Re: [R] question about xtable and Hmisc

Sorry, meant for this to go to the whole list. On Mar 1, 2007, at 3:29 PM, Charilaos Skiadas wrote: On Mar 1, 2007, at 2:52 PM, steve wrote: Unfortunately, this applies to print.xtable, and not to latex. I want to know how to eliminate them using latex() 1) Why do you need to use

Re: [R] How to read in this data format?

Dear All, thanks for the replies, Jim Holtman has given a solution which fits my needs, but Gabor Grothendieck did the same thing, but it looks like the coding will allow faster processing (should check this out tomorrow on a big datafile). @gabor: I don't understand the use of the grep

Re: [R] Query about data manipulation

Hi Thanks much for that input.It was extremely helpful. I am seeking some input about another stumbling block using RODBC; SQLQuery et al with large BLOB values. Although the following query dataFromDB - sqlQuery(channel, select uncompress(columnName) from tableName where Id=id ); returns just

[R] R code for Statistical Models in S ?

I just acquired a copy of Statistical Models in S, I guess most commonly known as the white book, and realized to my dismay that most of the code is not directly executable in R, and I was wondering if there was a source discussing the things that are different and what the new ways of

Re: [R] Row-wise two sample T-test on subsets of a matrix

Ranjan Maitra napsal(a): Here's one suggestion: convert the matrix into a three-dimensional array and use apply on it. Converting to 3 dims should not be neccessary: m - matrix(rnorm(110),ncol=22) t.list - apply(m,1,function(x){t.test(x[1:11],x[12:22],paired=TRUE)}) However, I have no idea

Re: [R] Default par() options

I am no expert on these topics but currently I am solving a similar issue using the .Rprofile file and the .First function. So maybe it's enough to put .First - function(){ par(whatever you want) further instructions if neccessary } Petr Thomas Friedrichsmeier napsal(a): The following

[R] Preserving order in an intersection

Hi. Here's an odd request that actually arose out of my own bad planning. Suppose I do the following: which(v1 %in% v2) I will get a set of indices for v1, and they will be ordered in the same order that v1 is in. I want the indices of the intersection for v1 ordered according to v2. I do

Re: [R] How to read in this data format?

On 3/1/07, Bart Joosen [EMAIL PROTECTED] wrote: Dear All, thanks for the replies, Jim Holtman has given a solution which fits my needs, but Gabor Grothendieck did the same thing, but it looks like the coding will allow faster processing (should check this out tomorrow on a big datafile).

Re: [R] R code for Statistical Models in S ?

The White Book provides the original S Language Specification. This was what existed at Bell labs way back then. Subsequent implementations, both S-Plus and R, will differ on details. Also, a lot of development effort has flowed over the dam since publication, so both implementations contain lots

[R] Simplest question ever...

Let's say i have a = c(1, 4, 5) b = c(2, 6, 7) and i have matrix m, what's an efficient way of access m[1, 2], m[4, 6], m[5, 7] like of course m[a, b] = is not going to do, but what's an expression that will allow me to have that list? Thanks! -- View this message in context:

Re: [R] R code for Statistical Models in S ?

You might want to start looking at the FAQ's http://cran.r-project.org/faqs.html in particular http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-and-S robert On 3/1/07, Charilaos Skiadas [EMAIL PROTECTED] wrote: I just acquired a copy of Statistical Models in S, I guess most commonly known as

[R] permutation tests on autocorrelations

Dear list, I have this huge array of numbers, of dimensions 67 x 33 x 51 x 6, the 6 being the replications. I wanted to test for evidence of autocorrelation between the 6 replications, marginally. I can calculate the first-order autocorrelation very easily using an appropriately defined

Re: [R] question about xtable and Hmisc

Also, if x is a data frame, latex(x) contains the row numbers. Can I get rid of them here as well? I think you are asking for the rowname=NULL argument. latex(x, rowname=NULL) See ?latex to confirm if that is what you are looking for. Rich __

[R] covariance question which has nothing to do with R

This is a covariance calculation question so nothing to do with R but maybe someone could help me anyway. Suppose, I have two random variables X and Y whose means are both known to be zero and I want to get an estimate of their covariance. I have n sample pairs (X1,Y1) (X2,Y2) . . . . .

[R] passing arguments to multiple functions

Hi, I am writing a function and I want to pass arguments to a function I will call internally without having to specify all the possible arguments to pass. Usually I would use '...' but I want to do this for two functions that I will call, and the two functions do not take the same arguments.

[R] Using R for devices trial

I would like to use R for submissions to FDA/CDRH (the medical device company I work for currently uses only SAS). Previous postings to the list regarding R and 21 CFR 11 compliance have been very helpful. However, reluctance to using open source software for statistical analyses and reporting

Re: [R] Simplest question ever...

I'm not sure this is the most efficient, but how about: diag(m[a,b]) ? On 3/1/07, yoo [EMAIL PROTECTED] wrote: Let's say i have a = c(1, 4, 5) b = c(2, 6, 7) and i have matrix m, what's an efficient way of access m[1, 2], m[4, 6], m[5, 7] like of course m[a, b] = is not going to

[R] FTP download from ftp.sec.gov

All, I have managed to download files from web sites and ftp sites using R, so just for fun I tried to do so from the SEC's ftp site using the following code: ftp - ftp://anonymous:[EMAIL PROTECTED]/edgar/full-index/company.idx download.file(url=ftp, destfile=test.txt) And it does not work. R

Re: [R] Default par() options

Thomas: I am not sure exactly what you are asking for below, but I wonder if your query could be satisfied by the judicious use of the ... argument in a wrapper function to par(), like myPar=function(bg=lightgray, pch=19,...)par(bg=bg,pch=pch,...) or perhaps myX11 - function(width=10,

Re: [R] Query about data manipulation

On Thu, 2007-03-01 at 12:35 -0800, lalitha viswanath wrote: Hi Thanks much for that input.It was extremely helpful. I am seeking some input about another stumbling block using RODBC; SQLQuery et al with large BLOB values. Although the following query dataFromDB - sqlQuery(channel,

Re: [R] Simplest question ever...

try: a = c(1, 4, 5) b = c(2, 6, 7) m - matrix(1:49,7,7) m[cbind(a,b)] # 'array' indexing [1] 8 39 47 m [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,]18 15 22 29 36 43 [2,]29 16 23 30 37 44 [3,]3 10 17 24 31 38 45 [4,]4 11 18 25

  1   2   >