[R] Group a dinamic number of vectors in a data.frame

2006-02-20 Thread Daniele Medri
Hi all, I need to create a data.frame from a variable number of vectors. The number of these vectors could change so I need a dinamic way to group all in a data.frame. The number is length(abc). e.g. vectors in my workspace N1 -c(1,2,3,4) N2 -c(1,2,3,4) N3 -c(1,2,3,4) abc -c(1,2,3) the

Re: [R] Group a dinamic number of vectors in a data.frame

2006-02-20 Thread Dimitris Rizopoulos
one way is: N1 - rnorm(4) N2 - rnorm(4) N3 - rnorm(4) N4 - rnorm(4) X1 - LETTERS[1:4] ### nams - c(paste(N, 1:4, sep = ), X1) dat - data.frame(lapply(nams, get)) names(dat) - nams dat # check also sapply(dat, data.class) I hope it helps. Best, Dimitris Dimitris

Re: [R] Group a dinamic number of vectors in a data.frame

2006-02-20 Thread Gabor Grothendieck
data.frame(OneVector, TwoVector, sapply(apropos(^N[0-9]*$), get, simplify = FALSE)) On 2/20/06, Daniele Medri [EMAIL PROTECTED] wrote: Hi all, I need to create a data.frame from a variable number of vectors. The number of these vectors could change so I need a dinamic way to group all in

[R] Plots: displaying mathematical symbols in specific fonts

2006-02-20 Thread Ed Merkle
Dear SavioRs, I am doing some research where characters in different microsoft fonts serve as experimental stimuli. Hence, in plot labels, I would like to display the characters in specific microsoft fonts. I have figured out how to display letters and numbers, but I am having trouble with

[R] Time on x-axis

2006-02-20 Thread Sangeetha Swaminathan
Hello, I just started using the GNU R. I am having trouble plotting my data. I am tryin to plot the following data: TIMESTAMPLOGIN-TIME (hh:mm:ss) (s) 23:55:03 0.990972 23:55:03 0.990972 23:50:04 0.878968 23:45:04 0.969271 23:40:03 0.868848

Re: [R] Time on x-axis

2006-02-20 Thread Gabor Grothendieck
Try this: # test data tt - c(23:05:02, 23:10:02, 23:15:03, 23:20:03, 23:25:03, 23:30:03, 23:35:03, 23:40:03, 23:45:04, 23:50:04, 23:55:03, 23:55:03) x - c(0.575764, 0.738379, 0.567414, 0.663436, 0.599834, 0.679571, 0.88141, 0.868848, 0.969271, 0.878968, 0.990972, 0.990972) library(zoo)

[R] doubt on periodogram

2006-02-20 Thread Pradeep Gunda
hi all, I am using the periodogram function provided by the GeneTS library in R. The periodogram function takes a time-series as input and returns a spec(a vector of values) and their frequencies. My requirement is, I have to loop through the values in this spec and then remove some of these,

[R] Nested AIC

2006-02-20 Thread Aaron MacNeil
Greetings, I have recently come into some confusion over weather or not AIC results for comparing among models requires that they be nested. Reading Burnham Anderson (2002) they are explicit that nested models are not required, but other respected statisticians have suggested that

[R] mva.pairs

2006-02-20 Thread Vijay A Raghavan
Hello, I am using the following code to plot an MVA plot. library(affy) library(Biobase) library(limma) library(gcrma) pd-read.phenoData(Clk.targets.2.txt,header=TRUE, row.names=1,as.is=TRUE,sep=\t) Data - ReadAffy(filenames=pData(pd)$FileName,phenoData=pd)

Re: [R] Linux Distribution Choice

2006-02-20 Thread Ulises M. Alvarez
Ubuntu is a good choice : ) First, I will recommend you to take a look at: http://ubuntuguide.org/ Specially... http://ubuntuguide.org/#extrarepositories It is slightly out of date, but still is useful. Once you are done with that, installing R is quit simple. From a terminal -available from

[R] glob2rx function not working

2006-02-20 Thread mtb954
Dear R users, Inspired by previous list discussion of the glob2rxc function, I am attempting to create a new vector called TOTAL by summing all vectors whose names begin with ABC: TOTAL = sum(list = ls(pattern = glob2rx(ABC*))) I'm running R 2.2.1 on Windows XP. Can anyone say what I'm missing?

Re: [R] Plots: displaying mathematical symbols in specific fonts

2006-02-20 Thread Prof Brian Ripley
On Mon, 20 Feb 2006, Ed Merkle wrote: Dear SavioRs, I am doing some research where characters in different microsoft fonts serve as experimental stimuli. Hence, in plot labels, I would like to display the characters in specific microsoft fonts. I have figured out how to display letters

Re: [R] glob2rx function not working

2006-02-20 Thread Gabor Grothendieck
Your expression is trying to sum a character vector containing the names of the variables that begin with ABC. Try this and try executing each portion of it to understand it better: ABC1 - ABC2 - 1 do.call(sum, lapply(apropos(glob2rx(ABC*)), get)) # 2 On 2/20/06, [EMAIL PROTECTED] [EMAIL

Re: [R] help on dyn.load()

2006-02-20 Thread Feng Tai
Thanks a lot! It works now. Feng From: Roger Bivand [EMAIL PROTECTED] Reply-To: [EMAIL PROTECTED] To: Feng Tai [EMAIL PROTECTED] CC: r-help@stat.math.ethz.ch Subject: Re: [R] help on dyn.load() Date: Mon, 20 Feb 2006 08:13:49 +0100 (CET) On Mon, 20 Feb 2006, Feng Tai wrote: Hi, I used

Re: [R] glob2rx function not working

2006-02-20 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: Dear R users, Inspired by previous list discussion of the glob2rxc function, I am attempting to create a new vector called TOTAL by summing all vectors whose names begin with ABC: TOTAL = sum(list = ls(pattern = glob2rx(ABC*))) I'm running R 2.2.1 on Windows

Re: [R] Linux Distribution Choice

2006-02-20 Thread Dirk Eddelbuettel
Ulises, Thanks for the helpful post but allow me to add one or two corrections: On 20 February 2006 at 11:40, Ulises M. Alvarez wrote: | Ubuntu is a good choice : ) | | First, I will recommend you to take a look at: | http://ubuntuguide.org/ | | Specially... |

Re: [R] doubt on periodogram

2006-02-20 Thread Ben Bolker
Pradeep Gunda pradeep.gunda at gmail.com writes: hi all, I am using the periodogram function provided by the GeneTS library in R. myts - ts(c(1,2,3,40,3,2,5,40,2,3,1,40), frequency=1); myper - periodogram(myts); myper something like myper$spec - myper$spec[myper$spec4]

[R] get predicted values of models fit by lmer()?

2006-02-20 Thread Mingyu Feng
Greetings. I'm having trouble using lmer beyond a first step. After fitting my data using lmer(), I can use summary(model) to see a list of statistics and summary information. But I also need the predicted values, values of parameters for fixed and random effects and plots, etc. I don't how to

[R] Boxplot Help for Neophyte

2006-02-20 Thread Alex Park
R helpers I am getting to grips with R but came across a small problem today that I could not fix by myself. I have 3 text files, each with a single column of data. I read them in using: myData1-scan(C:/Program Files/R/myData1.txt) myData2-scan(C:/Program Files/R/myData2.txt)

Re: [R] Boxplot Help for Neophyte

2006-02-20 Thread Nolwenn LeMeur
Hi Alex, how about, myData-data.frame(myData1=myData1,myData2=myData2,myData3=myData3) boxplot(myData) Nolwenn ** Nolwenn Le Meur, PhD Fred Hutchinson Cancer Research Center Computational Biology 1100 Fairview Ave. N., M2-B876 P.O. Box 19024 Seattle, WA

Re: [R] Boxplot Help for Neophyte

2006-02-20 Thread Marc Schwartz (via MN)
On Mon, 2006-02-20 at 20:27 +, Alex Park wrote: R helpers I am getting to grips with R but came across a small problem today that I could not fix by myself. I have 3 text files, each with a single column of data. I read them in using: myData1-scan(C:/Program Files/R/myData1.txt)

Re: [R] Need frequency distribution for x,y coordinates

2006-02-20 Thread Jim Lemon
mark shanks wrote: Hi, I have a set of data in x,y coordinates across the range of -5 to 5 in each dimension. I would like to obtain the frequency distribution of the different points, and then graph them so you can see which of the points are the most frequently occurring.

[R] panel-VAR

2006-02-20 Thread tang5
Hi, Anybody knows if there is a package for panel VAR ? thanks __ 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] Tobit Regression (residual Assumption)

2006-02-20 Thread Jim Lemon
Hi folks, I exchanged a couple of emails privately about tobit regression and advised that he (bambang pramono) ask the list for pointers to a basic, widely available, introductory reference. I think it would be a great help if someone would kindly suggest such a reference. Thanks, Jim

Re: [R] lattice: calling functions

2006-02-20 Thread Deepayan Sarkar
On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 2/14/2006 9:38 AM, Gabor Grothendieck wrote: On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 2/14/2006 8:56 AM, Wolfram Fischer wrote: I defined three functions: fun0 - function( x=1:5, y=1:5, ... ) xyplot( y ~ x, ... )

[R] need a R-code formatter?

2006-02-20 Thread Michael
Hi all, I am using Tin-R as my editor; I use it because it allows me to send several selected lines to R-console and execute them... In some sense, this is my line-by-line debugger. But it doesn't have a code formatter, I have to layout the indention myself -- when there are many { } blocks

[R] power law

2006-02-20 Thread Glazko, Galina
Dear list, Does anyone know how to fit the power law distribution? I have the empirical distribution and would like to check whether it fits power law (with the power estimated from the data). Any hints are appreciated Best regards Galina [[alternative HTML version

[R] how to close all windows?

2006-02-20 Thread Michael
Hi all, I have 10+ graphic windows and 10+ R help window opening now... How do I close them all at once? Thanks a lot! [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Reanding a windows file

2006-02-20 Thread Alexandra R. M. de Almeida
Dear R users Some one knows how to read a windows file with the adress like D:\CSV\Work\test.csv without change the \ by / by myself?! There is some function or parameter that read like the adress windows form. I tried to search somethig to substitute the \ by /, but I don't find...

[R] SVM weights

2006-02-20 Thread Houssam Nassif
Hello I am using e1071, svm. After tuning the parameters, I would like to get the weights assigned to the different features of my feature vector. This is to perform a feature reduction step by eliminating features with low discrimination values. How can i get the weight vector? best

[R] standardizing data

2006-02-20 Thread Christian Jones
Hello R team, IŽm looking for a way to standardize (z transformation= standard deviation 1 and mean 0) a row of x y coordinates in order to conduct a trend analysis. Does anyone know the command in R? many thanks for help in advance Christian __

[R] Group a dinamic number of vectors in a data.frame

2006-02-20 Thread Daniele Medri
Hi all, I need to create a data.frame from a variable number of vectors. The number of these vectors could change so I need a dinamic way to group all in a data.frame. The number is length(abc). e.g. vectors in my workspace N1 -c(1,2,3,4) N2 -c(1,2,3,4) N3 -c(1,2,3,4) abc -c(1,2,3) the

[R] How to read more than 1 table?

2006-02-20 Thread Atte Tenkanen
Question 1) I want to read many csv-tables and run the same commands for all of them. Is there some simpler solution than this below to solve this problem? for (g in 1:6) { if (g==1){k=kt1_0057} if (g==2){k=kt1_0101} if (g==3){k=kt1_0613} if (g==4){k=staten} if (g==5){k=tenpenny} if

[R] Using LAPACK in C-code to be loaded in R - getting-started-help

2006-02-20 Thread Søren Højsgaard
I want to speed up computations (involving matrices) by writing some C-code to be loaded. In the C-code, I need to invert matrices etc. As I've understood the writing R extensions doc, I can use use #include R_ext/Linpack.h in my .c-file and get access to linpack's facilities within my C-code.

[R] [R-pkgs] Zigurrat updated

2006-02-20 Thread Bob Wheeler
I've updated the Ziggurat normal and exponential generator in SuppDists in accordance with Leong, et.al (2005). A comment on the implementation of the Ziggurat Method. Jour. Stat. Sci. 12-7, 1-4. The fix takes care of problems for very large sets of random values. I tested it a bit on small

Re: [R] Reanding a windows file

2006-02-20 Thread Peter Dalgaard
Alexandra R. M. de Almeida [EMAIL PROTECTED] writes: Dear R users Some one knows how to read a windows file with the adress like D:\CSV\Work\test.csv without change the \ by / by myself?! There is some function or parameter that read like the adress windows form. I tried to search

Re: [R] standardizing data

2006-02-20 Thread Peter Dalgaard
Christian Jones [EMAIL PROTECTED] writes: Hello R team, I´m looking for a way to standardize (z transformation= standard deviation 1 and mean 0) a row of x y coordinates in order to conduct a trend analysis. Does anyone know the command in R? many thanks for help in advance Christian

Re: [R] standardizing data

2006-02-20 Thread Will Terry
Hey there, Did you try scale(x)? This centers and scales the data in a matrix x. Will On Feb 19, 2006, at 1:05 PM, Christian Jones wrote: Hello R team, IŽm looking for a way to standardize (z transformation= standard deviation 1 and mean 0) a row of x y coordinates in order to conduct

[R] var-covar matrices comparison:

2006-02-20 Thread Aldi Kraja
Hi, Using package gclus in R, I have created some graphs that show the trends within subgroups of data and correlations among 9 variables (v1-v9). Being interested for more details on these data I have produced also the var-covar matrices. Question: From a pair of two subsets of data (with 9

Re: [R] Reanding a windows file

2006-02-20 Thread Duncan Murdoch
On 2/16/2006 12:39 PM, Alexandra R. M. de Almeida wrote: Dear R users Some one knows how to read a windows file with the adress like D:\CSV\Work\test.csv without change the \ by / by myself?! There is some function or parameter that read like the adress windows form. I tried to search

Re: [R] How to read more than 1 table?

2006-02-20 Thread Kjetil Brinchmann Halvorsen
Atte Tenkanen wrote: Question 1) I want to read many csv-tables and run the same commands for all of them. Is there some simpler solution than this below to solve this problem? for (g in 1:6) { if (g==1){k=kt1_0057} if (g==2){k=kt1_0101} if (g==3){k=kt1_0613} if (g==4){k=staten} if

Re: [R] need a R-code formatter?

2006-02-20 Thread Søren Højsgaard
You can use emacs with ESS - and/or try to encourage (kindly) the creator of Tinn-R to create such a code-formatter... Best Søren Fra: [EMAIL PROTECTED] på vegne af Michael Sendt: on 15-02-2006 07:05 Til: R-help@stat.math.ethz.ch Emne: [R] need a R-code

Re: [R] How to read more than 1 table?

2006-02-20 Thread Søren Højsgaard
for (k in list.files()){ ... } best Søren Fra: [EMAIL PROTECTED] på vegne af Atte Tenkanen Sendt: ma 20-02-2006 22:40 Til: r-help@stat.math.ethz.ch Emne: [R] How to read more than 1 table? Question 1) I want to read many csv-tables and run the same commands

Re: [R] Using LAPACK in C-code to be loaded in R - getting-started-help

2006-02-20 Thread Prof Brian Ripley
Are you talking about LAPACK (title) or LINPACK (body)? On Mon, 20 Feb 2006, Søren Højsgaard wrote: I want to speed up computations (involving matrices) by writing some C-code to be loaded. In the C-code, I need to invert matrices etc. As I've understood the writing R extensions doc, I can

[R] Unable to configure R 2.2.1 on Solaris 5.10 x86_64

2006-02-20 Thread Aric Gregson
Hello, Apologies for the post here. I have read the R-Admin (learned a lot!) and searched the web for days, but still fail at compiling R on my Ultra 20 running solaris 10 x86 1/06. This is the tail of './configure' output: . checking for dlopen in -ldl... yes checking readline/history.h

[R] linear discriminant analysis in MASS

2006-02-20 Thread Alain Paquette
Hello R people I now know how to run my discriminant analysis with the lda function in MASS: lda.alain=lda(Groupes ~ Ht.D0 + Lc.Dc + Ram + IDF, gr, CV = FALSE) and it works fine. But I am missing a test and cannot find any help on how to get it, if it exist. The S equivalent:

Re: [R] bVar slot of lmer objects and standard errors

2006-02-20 Thread Ulrich Keller
Hello, I'm sorry to resurrect this thread that I started almost two months ago. I've been pretty busy since I posted my question and the issue is not that high on my priority list. Thanks to all those who replied, and I hope I can tickle your interest again. As a reminder, my question was how

Re: [R] Linux Distribution Choice

2006-02-20 Thread Neil Shephard
I thought I'd add my vote for Gentoo, which has been my distro of choice for some time (although Slackware is very good as well, and teaches you a lot about *NIX in general). The main advantage of gentoo is that updating is very easy using the portage system (as pointed out by jon butchar), but

[R] How to sum values across multiple variables using a wildcard?

2006-02-20 Thread mtb954
I have a dataframe called data with 5 records (in rows) each of which has been scored on each of many variables (in columns). Five of the variables are named var1, var2, var3, var4, var5 using headers. The other variables are named using other conventions. I can create a new variable called var6

[R] do.call, browser and traceback

2006-02-20 Thread hadley wickham
A problem that I've encountered when using do.call a lot is very large stack traces, eg: f - function(x) stop() do.call(error, mtcars) traceback() f - function(x) browser() do.call(f, mtcars) I have hacked together my own version of traceback to fix this by limiting the length of each line to 80

Re: [R] Tobit Regression (residual Assumption)

2006-02-20 Thread Gray Calhoun
The second edition of Jeff Wooldridge's Introductory Econometrics, a modern approach has a subsection on tobit regression and should be widely available. It was written for undergraduate economics majors. --Gray On 2/13/06, Jim Lemon [EMAIL PROTECTED] wrote: Hi folks, I exchanged a couple of

Re: [R] do.call, browser and traceback

2006-02-20 Thread Gabor Grothendieck
If f is long then you can get some savings like this: do.call(f, mtcars) # note: used f rather than f This does not solve the whole problem but its a step. On 2/20/06, hadley wickham [EMAIL PROTECTED] wrote: A problem that I've encountered when using do.call a lot is very large stack

Re: [R] How to sum values across multiple variables using a wildcard?

2006-02-20 Thread Gabor Grothendieck
See: ?rowSums On 2/20/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: I have a dataframe called data with 5 records (in rows) each of which has been scored on each of many variables (in columns). Five of the variables are named var1, var2, var3, var4, var5 using headers. The other

[R] R and packages

2006-02-20 Thread Erin Hodgess
Dear R People: Here is yet another strange problem. I'm using R in one of my classes. However, the computer lab has something called Deep Freeze and the students cannot save anything to the hard drive. I had R installed and things were working well. They would save their .Rdata files to

Re: [R] How to sum values across multiple variables using a wildcard?

2006-02-20 Thread Simon Blomberg
data - data.frame(var1=c(1,2,3), var2=c(3,4,5), var3=c(4,5,6), foo = c(100,200,300)) # sum rows with var in their name rowSums(data[, grep(var, names(data))]) 1 2 3 8 11 14 [EMAIL PROTECTED] wrote: I have a dataframe called data with 5 records (in rows) each of which has been scored on

Re: [R] How to sum values across multiple variables using a wildcard?

2006-02-20 Thread Marc Schwartz
On Mon, 2006-02-20 at 18:41 -0600, [EMAIL PROTECTED] wrote: I have a dataframe called data with 5 records (in rows) each of which has been scored on each of many variables (in columns). Five of the variables are named var1, var2, var3, var4, var5 using headers. The other variables are named

Re: [R] How to sum values across multiple variables using a wildcard?

2006-02-20 Thread mtb954
Thanks Gabor, Simon and Marc...I got this to work with the grep() and rowSums examples you provided. Mark On 2/20/06, Marc Schwartz [EMAIL PROTECTED] wrote: On Mon, 2006-02-20 at 18:41 -0600, [EMAIL PROTECTED] wrote: I have a dataframe called data with 5 records (in rows) each of which

Re: [R] R and packages

2006-02-20 Thread Duncan Murdoch
On 2/20/2006 8:07 PM, Erin Hodgess wrote: Dear R People: Here is yet another strange problem. I'm using R in one of my classes. However, the computer lab has something called Deep Freeze and the students cannot save anything to the hard drive. I had R installed and things were

Re: [R] R and packages

2006-02-20 Thread Duncan Murdoch
On 2/20/2006 8:07 PM, Erin Hodgess wrote: Dear R People: Here is yet another strange problem. I'm using R in one of my classes. However, the computer lab has something called Deep Freeze and the students cannot save anything to the hard drive. I had R installed and things were

Re: [R] need a R-code formatter?

2006-02-20 Thread Uwe Ligges
Søren Højsgaard wrote: You can use emacs with ESS - and/or try to encourage (kindly) the creator of Tinn-R to create such a code-formatter... ... or contribute the code formatter to the Tinn-R project which is one of the great ideas of Open Source software. The Tinn-R developers will be

Re: [R] how to close all windows?

2006-02-20 Thread Uwe Ligges
Michael wrote: Hi all, I have 10+ graphic windows and 10+ R help window opening now... How do I close them all at once? Thanks a lot! [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

[R] indexing within panels in xyplot

2006-02-20 Thread Sebastian Luque
Dear R-helpers, I need to show a linear fit through a subset of the data within each combination of levels of two factors. So I prepared an xyplot with different panels for each level of one of the factors, and different symbols within each panel for the levels of the second factor. My problem

Re: [R] Group a dinamic number of vectors in a data.frame

2006-02-20 Thread Uwe Ligges
Daniele Medri wrote: Hi all, I need to create a data.frame from a variable number of vectors. The number of these vectors could change so I need a dinamic way to group all in a data.frame. The number is length(abc). e.g. vectors in my workspace N1 -c(1,2,3,4) N2 -c(1,2,3,4) N3

Re: [R] R and packages

2006-02-20 Thread J Dougherty
You want to take the sysad by the hairy of his chinny-chin-chin and explain the issue in short sentences. He or she ought to be able produce a solution without difficulty, and should not have any problems about doing so. It IS their job. JD On Monday 20 February 2006 17:07, Erin Hodgess

Re: [R] Heckman regression / adjustment for standard errors?

2006-02-20 Thread Arne Henningsen
Hi Brian, On Friday 17 February 2006 22:30, Brian Perron wrote: Hello folks, I am trying to estimate the two-step Heckman regression model. I would like to make an adjustment for intragroup correlations. Stata can implement this with the cluster option, but I am really hoping to stick

Re: [R] indexing within panels in xyplot

2006-02-20 Thread Frede Aakmann Tøgersen
Based on your two first sentences I think the solution is to use xyplot(y ~ x | facA, groups = facB, data = toydf,type=c(p,r)) Try it and see if this is what you want. Best regards Frede Aakmann Tøgersen Scientist Danish Institute of Agricultural Sciences Research Centre Foulum Dept. of