Re: [R] Need help

2008-10-27 Thread Tal Galili
Hello John.After a simple search in google, I came up with this: http://www.ats.ucla.edu/stat/R/seminars/Repeated_Measures/repeated_measures.htm Enjoy, Tal. On Mon, Oct 27, 2008 at 6:15 PM, JohnLi <[EMAIL PROTECTED]> wrote: > Hi All, > > I am new comer here and for Statistics R. I want to kn

Re: [R] ttkcombobox

2008-10-27 Thread Sundar Dorai-Raj
Hi, Greg, Yes, that helps immensely! Thanks, --sundar Greg Snow said the following on 10/27/2008 7:14 PM: Here is one example of using ttkcombobox, hopefully it helps: library(tcltk) have.ttk <- function () { # from Prof. Ripley as.character(tcl("info", "tclversion")) >= "8.5" } testf

[R] repeated measure one way anova

2008-10-27 Thread JohnLi
Hi all, I am a new comer for Statistics R. I am using R for one way repeated measure anova, for example, on the following data consisting of three groups. I appreciate your help if anybody help me show the procedure in R. John c2 c3 c4 85.83 75.86 84.19 85.91 73.18 85.9 86.12

[R] Kernel regression with random design

2008-10-27 Thread tzeshian
Hi All, I would like to find out which package in R that can handle nadaraya-watson and local linear estimators with random design or unequally spaced x's variables. Appreciate the help. Thank you. -- View this message in context: http://www.nabble.com/Kernel-regression-with-random-design-tp

Re: [R] read multi TXT file in R

2008-10-27 Thread jim holtman
for (i in 2:4){ input <- read.delim(paste('vegetation_', i, '.txt', sep='')) process the file } On Mon, Oct 27, 2008 at 6:53 PM, Alessandro <[EMAIL PROTECTED]> wrote: > Hi All, > > > > sorry for the question. I have a problem to understand how I could read > several txt file in R

[R] help on package or code for simutaneous equation probit(logit) model

2008-10-27 Thread Yong Wang
Dear List I am trying to fit a simutaneous equation logit model. i.e., the response variables of the structured equations are binomial, I am not sure if systemfit can do this job. A google search doesn't yield too much helpful information. Your knowledge on any other packages or codes are appreciat

Re: [R] Simple plotting question

2008-10-27 Thread Greg Snow
Creating your own plot from scratch will probably not be that hard. You can plot using type='n' to set up the plot. Then the symbols function will plot squares or rectangles with the length of the sides specified by you. The segments function can then add the lines. If you want a shape that

Re: [R] ttkcombobox

2008-10-27 Thread Greg Snow
Here is one example of using ttkcombobox, hopefully it helps: library(tcltk) have.ttk <- function () { # from Prof. Ripley as.character(tcl("info", "tclversion")) >= "8.5" } testfunc <- function() { tt <- tktoplevel() if( !have.ttk() ) stop( "Need version 8.5 or greater of t

Re: [R] License Question

2008-10-27 Thread Duncan Murdoch
[EMAIL PROTECTED] wrote: Katrin, R is completely free - no License. We've built a supported user friendly version: R-PLUS This statement is not correct. R is licensed under the GPL. This is a "free software" license, but it is not the same as having no license. There are oblig

Re: [R] Exponential regression (Y = exp(a*X)) and standard error of Ŷi

2008-10-27 Thread David Winsemius
I am reasonably sure that the log link with a zero intercept duplicates your efforts with nls(). Take a look at this: > txt <- "X Y + 2 2.927064 + 5 14.6582 + 4 8.567706 + 3 5.007817 + 1 1.710867 + 6 25.07823 + 4 8.567706 + 7

[R] Package "evd"

2008-10-27 Thread Manjunath B G
As I tried to draw random sample from Hüsler-Reiss density with the following code, library(evd) margins <- cbind(0, 1, seq(-0.5,0.5,0.1)) x <- rbvevd(101, dep = 1.7, model = "hr", mar1 = margins) [,1] [,2] [1,] -0.56662298 0.8448505 [2,] 0.55824918 1.0217279 [3,]

Re: [R] Installing gnomeGUI in ubuntu

2008-10-27 Thread Kjetil Halvorsen
Sorry. I forgot to say I am on ubuntu intrepid 8.10 RC, in case that matters. Kjetil On Mon, Oct 27, 2008 at 8:28 PM, Kjetil Halvorsen < [EMAIL PROTECTED]> wrote: > Hola! > > I am trying to compile gnomeGUI with R 2.8.0, using > inside R (started with sudo R) > >install.packages("gnomeGUI", dep=

[R] Installing gnomeGUI in ubuntu

2008-10-27 Thread Kjetil Halvorsen
Hola! I am trying to compile gnomeGUI with R 2.8.0, using inside R (started with sudo R) >install.packages("gnomeGUI", dep=TRUE) this fails with: checking for string.h... (cached) yes checking for strings.h... (cached) yes checking for gnome-config... no checking for gnomeConf.sh file in /usr/lo

[R] Simple plotting question

2008-10-27 Thread mcolaco
I'm trying to make a plot very similar to a forest plot, with points of different sizes with means, standard deviations, and a reference line. However, instead of using the standard deviations to determine the size of the data points, I would like to choose my own weightings. It seems like it sh

Re: [R] License Question

2008-10-27 Thread rplus
Katrin, R is completely free - no License. We've built a supported user friendly version: R-PLUS You can freely evaluate R-PLUS - the current free trial has no expiration date. At a click, import data from different formats, use Excel-like spreadsheet for manipulating y

[R] SciViews / Komodo installation problem

2008-10-27 Thread Basti-Fantasti
Dear Forum, I have a problem with the installation of SciViews / Komodo under OSX. After installing, SciViews and Komodo I get the following error message in R: "29:37: execution error: File Komodo wasn’t found. (-43)" I also set up the symbolic link, that is required. The connection between K

Re: [R] License Question

2008-10-27 Thread Tom Backer Johnsen
Freiberger, Katrin wrote: Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issues that need to be taken in consideration, any fees to pay by the compan

[R] read multi TXT file in R

2008-10-27 Thread Alessandro
Hi All, sorry for the question. I have a problem to understand how I could read several txt file in R with read.delim file: vegetation_2.txt, vegetation_3.txt, vegetation_4.txt, etc etc I tried with: for (i in 2:4) {file=paste("vegetation_",i,sep="") <- read.delim (file=paste("vege

Re: [R] Transferring results from R to MS Word2

2008-10-27 Thread Tom Backer Johnsen
Felipe Carrillo wrote: I've responded to this one in the original thread. Tom __ R-help@r-project.org 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 co

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Tom Backer Johnsen
Felipe Carrillo wrote: Are you still can't get the data into R? I sent you this -mail last week, did you read it? Yes I did read it, and I have installed the package and looked at the documentation. What I (after an admittently superficial inspection due to limited available time this weeke

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Tom Backer Johnsen
Greg Snow wrote: Another alternative (since you asked) is the RExcel project (http://rcom.univie.ac.at). I don't know if it will work better for your projects or not, but it may be worth a look. The basic idea is that it uses Excel as a front end and R as the background computational engine

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Tom Backer Johnsen
Greg Snow wrote: Any timing is of course your decision, but my observation (being one of my company's test cases for updating office) is that transitioning from MSWord 2003 to OpenOffice Writer will be less work/stress than transitioning from MSWord 2003 to MSWord 2007. That agrees with my o

Re: [R] Problem in installation of "Rgraphviz" package

2008-10-27 Thread Martin Morgan
Rgraphviz is a Bioconductor package so help can be found on the Bioconductor mailing list. Actually, this post, from Saturday, might be what you are looking for... https://stat.ethz.ch/pipermail/bioconductor/2008-October/024747.html Best, Martin ram basnet wrote: Dear R users, I am not so

Re: [R] Problem in installation of "Rgraphviz" package

2008-10-27 Thread stephen sefick
http://www.bioconductor.org/packages/release/bioc/html/Rgraphviz.html try going here downloading the binary and then installing it from your desktop look at ?install.packages hope this helps On Mon, Oct 27, 2008 at 5:21 PM, ram basnet <[EMAIL PROTECTED]> wrote: > Dear R users, > > I am not so use

[R] Problem in installation of "Rgraphviz" package

2008-10-27 Thread ram basnet
Dear R users,   I am not so used to this R software. I have to use the package " Rgraphviz" but found some problem in the installation process. I download this package and store in R library but i am not getting this package in R installation list. I made review in google search net and use the

[R] About Pareto distribution

2008-10-27 Thread Carlos López
Hi again :-) I finally was able to fix the program, thank you all very much for your help :-) Now I have a problem and I don´t know if it is possible to solve it with R, I have a data set, and because it is data from salaries I am suspecting it comes from a Pareto distribution, my questions a

[R] ttkcombobox

2008-10-27 Thread Sundar Dorai-Raj
Hi, all, (sessionInfo at the end) I've been struggling with the tcltk package and can't seem to get the ttkcombobox to work. Here's an example: library(tcltk) p <- tktoplevel() l <- tclVar() ## I don't know if I'm even calling it correctly cb <- ttkcombobox(p, values = letters[1:4], textvaria

Re: [R] License Question

2008-10-27 Thread dsmith
Katrin, R is free and you get great support from this mailling list. Check out R-PLUS by XLSolutions Corporation. The free R-PLUS trial has no expiration date. At a click, import data from different formats, use Excel-like spreadsheet for manipulating your data, create pub

Re: [R] Stuck with FFT

2008-10-27 Thread stephen sefick
have you read in the data? Have you made it a time series? Did you call it z? ?fft ?spectrum On Mon, Oct 27, 2008 at 2:01 PM, Rthoughts <[EMAIL PROTECTED]> wrote: > > Dear all, > > Before I can get into serious Fourier analysis of Radon time-series I am > practising with 24hour and 24.8hour sin

[R] Stuck with FFT

2008-10-27 Thread Rthoughts
Dear all, Before I can get into serious Fourier analysis of Radon time-series I am practising with 24hour and 24.8hour sinusoids to assist with my interpretation of signals ittributed to tidal input to Radon time series. I am stuck. I have tried researching this to no avail. I am awating a book

[R] dotted lines for branches in ape's plot.phylo?

2008-10-27 Thread Janet Young
Hi, I'm very much enjoying using the ape package to produce phylogenetic trees with colored branches (using edge.color). Is it also possible to specify that some branches should be drawn as dotted lines? That would be really useful. I've tried messing around with edge.width, but that does

Re: [R] Arrays of Trellis plots

2008-10-27 Thread Richard . Cotton
> the example below does not work. (i know it's not supposed, but it makes it > clear what i'm trying to achieve) > > par(mfrow=c(2,1)) > xyplot(y~x2|x1,data=dataframe1,pch=20) > xyplot(y~x2|x1,data=dataframe2,pch=20) > > i know i could probably merge the two datasets and do something like > xyp

Re: [R] Arrays of Trellis plots

2008-10-27 Thread erwann rogard
Thank you all, I've got everything I need now. On Mon, Oct 27, 2008 at 2:29 PM, <[EMAIL PROTECTED]> wrote: > > the example below does not work. (i know it's not supposed, but it makes > it > > clear what i'm trying to achieve) > > > > par(mfrow=c(2,1)) > > xyplot(y~x2|x1,data=dataframe1,pch=20) >

Re: [R] Arrays of Trellis plots

2008-10-27 Thread baptiste auguie
Hi, You could use the grid package to place treillis objects in any custom layout you want, for example (inspired by Paul Murrell's R graphics book < http://www.stat.auckland.ac.nz/~paul/RGraphics/rgraphics.html > fig 5.22), library(grid) library(lattice) df <- data.frame(x=rnorm(100),

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Liviu Andronic
On Mon, Oct 27, 2008 at 8:16 PM, Liviu Andronic <[EMAIL PROTECTED]> wrote: > On Mon, Oct 27, 2008 at 7:32 PM, Greg Snow <[EMAIL PROTECTED]> wrote: >> Another alternative (since you asked) is the RExcel project >> (http://rcom.univie.ac.at). >> > I didn't perform much search, but is there anything

Re: [R] Arrays of Trellis plots

2008-10-27 Thread Deepayan Sarkar
On 10/27/08, erwann rogard <[EMAIL PROTECTED]> wrote: > hello, > > the example below does not work. (i know it's not supposed, but it makes it > clear what i'm trying to achieve) > > par(mfrow=c(2,1)) > xyplot(y~x2|x1,data=dataframe1,pch=20) > xyplot(y~x2|x1,data=dataframe2,pch=20) > > i know

Re: [R] License Question

2008-10-27 Thread David Smith
As others have pointed out, you can download and use R for commercial purposes for free. This is allowed under the terms of the GPL under which R is licensed. If your company requires software to be professionally supported, you might want to consider purchasing a version of R that includes a sup

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Liviu Andronic
Hello, On Mon, Oct 27, 2008 at 7:32 PM, Greg Snow <[EMAIL PROTECTED]> wrote: > Another alternative (since you asked) is the RExcel project > (http://rcom.univie.ac.at). > I didn't perform much search, but is there anything similar for OpenOffice or Gnumeric? Basically, is there a cross-platform R

[R] Arrays of Trellis plots

2008-10-27 Thread erwann rogard
hello, the example below does not work. (i know it's not supposed, but it makes it clear what i'm trying to achieve) par(mfrow=c(2,1)) xyplot(y~x2|x1,data=dataframe1,pch=20) xyplot(y~x2|x1,data=dataframe2,pch=20) i know i could probably merge the two datasets and do something like xyplot(y~x2|x1

Re: [R] License Question

2008-10-27 Thread Marc Schwartz
on 10/27/2008 11:26 AM Duncan Murdoch wrote: > On 10/27/2008 9:04 AM, Freiberger, Katrin wrote: >> Dear All, >> >> I learned about R during my studies at Cologne University of Applied >> Science. Now I work at Allianz Dresdner Bauspar AG and I would like to >> install R here too. Is there any licen

[R] Console size/location on start-up in SDI mode

2008-10-27 Thread Folkes, Michael
Does anybody have a recommended code/method for setting start-up location and window size of the console when working in SDI mode? I use tinn-r, which is already splitting my desktop in half for tinn-r and the console. I'd like to see the console in the bottom, right corner to allow for a graphi

Re: [R] Computational problems in R

2008-10-27 Thread A.Noufaily
Many thanks for your suggestions... I am still checking which one is the most useful for my simulations. Concerning using logs, this might be very helpful, but I am not sure if I can use the following: A+B = A*(1 + B/A) = exp(log(A) + log(1 + B/A)) because unfortunately B can be negative. H

Re: [R] License Question

2008-10-27 Thread Duncan Murdoch
On 10/27/2008 1:19 PM, Barry Rowlingson wrote: 2008/10/27 Freiberger, Katrin <[EMAIL PROTECTED]>: Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issu

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Greg Snow
Another alternative (since you asked) is the RExcel project (http://rcom.univie.ac.at). I don't know if it will work better for your projects or not, but it may be worth a look. The basic idea is that it uses Excel as a front end and R as the background computational engine as an Excel plugin

Re: [R] License Question

2008-10-27 Thread Peter Dalgaard
Duncan Murdoch wrote: On 10/27/2008 9:04 AM, Freiberger, Katrin wrote: Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issues that need to be taken i

Re: [R] License Question

2008-10-27 Thread Barry Rowlingson
2008/10/27 Freiberger, Katrin <[EMAIL PROTECTED]>: > Dear All, > > I learned about R during my studies at Cologne University of Applied Science. > Now I work at Allianz Dresdner Bauspar AG and I would like to install R here > too. Is there any license issues that need to be taken in consideration

[R] Error: C stack usage is too close to the limit

2008-10-27 Thread John Poulsen
Hello, Does anyone know in general what causes the following error: Error: C stack usage is too close to the limit I rewrote some R code that worked fine to be more streamlined, and now I get this error. I don't really know where to start to fix the problem. Thanks, John __

Re: [R] Transferring results from R to MS Word

2008-10-27 Thread Greg Snow
Any timing is of course your decision, but my observation (being one of my company's test cases for updating office) is that transitioning from MSWord 2003 to OpenOffice Writer will be less work/stress than transitioning from MSWord 2003 to MSWord 2007. -- Gregory (Greg) L. Snow Ph.D. Statistic

Re: [R] Docbook equivalent of Sweave

2008-10-27 Thread Greg Snow
The odfWeave package does sweave using OpenOffice files instead of LaTeX files. The .odt files are XML based I believe. If you still want to expand to docbook, looking at the code of odfWeave should at least get you started. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

Re: [R] Methods for showing statistics over space

2008-10-27 Thread Greg Snow
Here are a few ideas. The symbols function has a few different symbols that it can add to a plot (circles as you mentioned), one of my favorites is the thermometers, like a barplot for each county, but with a surrounding rectangle for reference so they are easier to compare while not being line

[R] maximum dimension of SVM in e1701

2008-10-27 Thread Aimin Yan
Does anyone know the maximum dimension that svm can deal with in R? I am working on support vector regression. My data set is 1721*41030 1721 samples 41030 predictors Some information about my machine: cat /proc/cpuinfo processor : 0 vendor_id : GenuineIntel cpu family : 15

Re: [R] License Question

2008-10-27 Thread Duncan Murdoch
On 10/27/2008 9:04 AM, Freiberger, Katrin wrote: Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issues that need to be taken in consideration, any fe

[R] Displaying number of Y/N affected by tree in rule form RE: R question/request on rules from rpart

2008-10-27 Thread Sharma, Dhruv
Hi Prof. Williams, thanks for your suggestion. The updated code is below. It turns out it was a matter of displaying the second column in yval to get the number of N and subtracting it from the n column in the frame to get the number of Y remaining in a binary example. onc

Re: [R] create list of data frames

2008-10-27 Thread Benoit Boulinguiez
Hi Ken, Your command works fine and your comment was right as well, here is the final command (inspired from yours) directly on DATA1.x and DATA1.y foo3<-list() for (i in 1:1000) foo3[[i]]<-nls( Qe~LgmFla(a,b,Ce), data=list(Ce=DATA1.x[,i],Qe=DATA1.y[,i]), start=c(a=0.5,b=0.1) ) Thanks

Re: [R] how to plot chi-square distribution in the graph

2008-10-27 Thread Greg Snow
Here is an example of adding the distribution curve to a histogram: x <- rchisq(100, 5) hist(x, prob=TRUE) curve( dchisq(x, df=5), col='green', add=TRUE) curve( dchisq(x, df=10), col='red', add=TRUE ) It may be easier to compare the therotical curve to a density estimate rather than the histogra

[R] Need help

2008-10-27 Thread JohnLi
Hi All, I am new comer here and for Statistics R. I want to know how to use R to calculate one way repeated measure ANOVA, for example, on three groups of data. I appreciate if someone can help me solve this problem. John __ R-help@r-project.or

Re: [R] create list of data frames

2008-10-27 Thread Kenn Konstabel
Hi, > To realize the data frame I've tried this > > for (i in 1:1000) > { > foo<-list(c(foo[],data.frame( Ce=DATA1.x[,i],Qe=DATA1.y[,i]))) > } I think the following would do it: foo <- list() for(i in 1:1000) foo[[i]] <- data.frame(Ce = DATA1.x[,i], Qe=DATA1.y[,i]) But then again, do you re

Re: [R] counting run lengths

2008-10-27 Thread Mario Lavezzi
Dear Richard, Martin, Dimitris and Domenico thank you very much for your help. I must say that the fastest procedure appears to be the one suggested by Richard This runs pretty quickly: unSpells <- nrow(Atr) - apply(Atr,2,function(x) max(which(x==1))) #c(4,0,7,0) If I may abuse of your kin

[R] create list of data frames

2008-10-27 Thread Benoit Boulinguiez
Hi all, I need to realize nonlinear regression on a thousand data sets. I guess the lapply function would help me on that thus I'd like to create a list of data frames, each data frame containing the data as follows: CeQe 1 1.849147 0.1958672 2 10.054250 0.5771036 3 18.077246

[R] Exponential regression (Y = exp(a*X)) and standard error of Ŷi

2008-10-27 Thread alemu gonsamo
[EMAIL PROTECTED]   Hello   First I want to implement exponential regression in R, with out constant for the following formula. Y = exp(a*X) ‘a’ is coefficient I wanted to determine. That I could do also in SPSS but my question is rather to estimate the ‘standard error of  Ŷi  at each Xi. This

Re: [R] How to see source code for na.omit?

2008-10-27 Thread jim holtman
Type 'na.omit' and see where it leads. In this case you see it is an object class and there are several methods, all of which are not visible, so you use 'getAnywhere' to see them: > na.omit function (object, ...) UseMethod("na.omit") > methods(na.omit) [1] na.omit.data.frame* na.omit.default*

Re: [R] replace a few strings in a text file

2008-10-27 Thread baptiste auguie
First of all, let me apologize for the 200 files that would result from running my "minimal example". I've just realised this on my Desktop; it seems i overlooked the arguments passed to sapply(). Thank you for your reply Duncan, I realize now my approach was quite silly. It is indeed much

[R] License Question

2008-10-27 Thread Freiberger, Katrin
Dear All, I learned about R during my studies at Cologne University of Applied Science. Now I work at Allianz Dresdner Bauspar AG and I would like to install R here too. Is there any license issues that need to be taken in consideration, any fees to pay by the company? I know there are answers

[R] How to see source code for na.omit?

2008-10-27 Thread useR
Hi R helpers, I'd like to see source code for some of built-in R functions... for example, I would like too see how "na.omit" was implemented? Thanks? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] counting run lengths

2008-10-27 Thread Martin Morgan
Hi Mario -- This function f <- function(m) { ## next 2 lines due to Bill Dunlap ## http://tolstoy.newcastle.edu.au/R/e4/devel/08/04/1206.html csum <- cumsum(!m) crun <- csum - cummax(m * csum) matrix(ifelse(crun > 0, (crun-1) %% nrow(m) + 1, 0), nrow=nrow(m)) } ret

Re: [R] counting run lengths

2008-10-27 Thread Richard . Cotton
> >> It works, but the for (i in ...) loop slows down the simulation a lot. > >> > >> Any suggestion on how to avoid this loop? (or in general, to speed up > >> this part of the simulation) > Actually, I have not specified the following: i want to consider only > the "most recent" sequence of

Re: [R] Process dataframes from list

2008-10-27 Thread Gábor Csárdi
On Mon, Oct 27, 2008 at 2:24 PM, Wade Wall <[EMAIL PROTECTED]> wrote: > Hi Gabor, > > Thanks for the code, it's what I was looking for. The route I was > thinking about was to use, based on your code, the names from > my.data.frames as a pointer to the actually dataframes. > > I was envisioning so

Re: [R] Process dataframes from list

2008-10-27 Thread Wade Wall
Hi Gabor, Thanks for the code, it's what I was looking for. The route I was thinking about was to use, based on your code, the names from my.data.frames as a pointer to the actually dataframes. I was envisioning something like a <- data.frame(1:10) b <- data.frame(10:1) c <- data.frame(letters[

Re: [R] Question of "Quantile Regression for Longitudinal Data".

2008-10-27 Thread roger koenker
This is really NOT an R-help question. You just need to read the code more carefully and everything will be revealed there. This isn't like the Pisan Cantos. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 21

Re: [R] How do I sign off this server

2008-10-27 Thread Bos, Roger
John, If you really do like the R-help list, you may want to see if your email program has folders and auto-rules. Since all R-help emails have [R] in the subject line, you can put all those email in a separate folder so everything doesn't get mixed up. That is what I do. HTH, Roger -Origi

Re: [R] Pkg rgl: installation fails because of x11

2008-10-27 Thread Ben Bolker
Duncan Murdoch stats.uwo.ca> writes: > > Matthieu Stigler wrote: > > Hello > > > > Im trying to install package rgl in a freshly installed Ubuntu 8.04 > > system. I have a problem (actually is has been reported three times on > > the R-list but the answers were too complicated for me) when > > i

[R] foreign

2008-10-27 Thread Paloma Main
I have problems to read (file read error) Minitab Portable Files (Minitab13) with read.mtp(foreign) for R versions posterior to R.6.0. Thanks P. Main __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] replace a few strings in a text file

2008-10-27 Thread Duncan Murdoch
baptiste auguie wrote: Dear all, I wrote a wrapper to a FORTRAN program using R. The main program uses a text file (~200 lines) as an input describing the simulation to be run. I typically generate the file once with the right parameters using a combination of file(), paste(), cat(). This

Re: [R] Mallows' distance or Earth Mover's distance in R?

2008-10-27 Thread Rainer M Krug
On Sun, Oct 26, 2008 at 4:01 PM, Matthias Kohl <[EMAIL PROTECTED]> wrote: > Hi Rainer, > > we have not implemented the Mallow's distance so far. Out package distrEx > includes implementations of the Hellinger, the Kolmogorov, the total > variation and the Cramer von Mises distance. Combined with ou

Re: [R] Pkg rgl: installation fails because of x11

2008-10-27 Thread Duncan Murdoch
Matthieu Stigler wrote: Hello Im trying to install package rgl in a freshly installed Ubuntu 8.04 system. I have a problem (actually is has been reported three times on the R-list but the answers were too complicated for me) when installing: configure: using libpng dynamic linkage checking for

[R] replace a few strings in a text file

2008-10-27 Thread baptiste auguie
Dear all, I wrote a wrapper to a FORTRAN program using R. The main program uses a text file (~200 lines) as an input describing the simulation to be run. I typically generate the file once with the right parameters using a combination of file(), paste(), cat(). This is fine, and it works

Re: [R] counting run lengths

2008-10-27 Thread Dimitris Rizopoulos
then try the following: Atr <- cbind(rep(1:0, each = 4), 1, c(1, rep(0, 7)), 1) Atr <- rbind(c(0, 1, 0, 1), Atr) apply(Atr, 2, function (x) { rr <- rle(x) if (tail(rr$values, 1) == 0) tail(rr$length, 1) else 0 }) I hope this what you're looking for. Best, Dimitris Mario Lavezzi wrote

Re: [R] counting run lengths

2008-10-27 Thread Mario Lavezzi
Hi Dimitris, thank you very much. Actually, I have not specified the following: i want to consider only the "most recent" sequence of zeros, that is the last part of the time series. That is, If I have: [,1] [,2] [,3] [,4] [1,]0101 [2,]1111 [3,]11

Re: [R] counting run lengths

2008-10-27 Thread Domenico Vistocco
Try this: unSpells[tail(Atr,1)==0] <- apply(Atr,2,function(x)sum(x==0))[tail(Atr,1)==0] Or (if you don't have to preserve the value in the unSpells vector): unSpells <- apply(Atr,2,function(x)sum(x==0)) But in this case you have 0 instead of 1 in the second and fourth position. Ciao, domenico

Re: [R] counting run lengths

2008-10-27 Thread Dimitris Rizopoulos
it's not totally clear to me what exactly do you need in this case, but have a look at the following: Atr <- cbind(rep(1:0, each = 4), 1, c(1, rep(0, 7)), 1) unSpells <- colSums(Atr == 0) unSpells[unSpells == 0] <- 1 unSpells I hope it helps. Best, Dimitris Mario Lavezzi wrote: Hello, I ha

[R] counting run lengths

2008-10-27 Thread Mario Lavezzi
Hello, I have the following problem. I am running simulations on possible states of a set of agents (1=employed, 0=unemployed). I store these simulated time series in a matrix like the following, where rows indicates time periods, columns the number of agents (4 agents and 8 periods in this

Re: [R] Odp: Request: Most repeated sequence considering combinations at each row

2008-10-27 Thread Muhammad Azam
Dear Petr Thanks for the response. Hope it will now help me to proceed. best regards M.Azam From: Petr PIKAL <[EMAIL PROTECTED]> Cc: R Help ; [EMAIL PROTECTED] Sent: Monday, October 27, 2008 8:43:27 AM Subject: Odp: [R] Request: Most repeated sequence considerin

Re: [R] odd behaviour of identical

2008-10-27 Thread Patrick Burns
Wacek Kusnierczyk wrote: smells bad design. Nonsense. One of the key design features of R is that it hides implementation details from users. They are free to think about the substantive issues with their data rather than worrying about computational trivia. There may have been some, bu

Re: [R] weird behavior with the 3rd root....

2008-10-27 Thread Robin Hankin
This comes up from time to time. The problem is that one needs complex numbers to address taking the third root: there are three cube roots for any nonzero number (real or complex). To wit: > (-0.084121928394+0i)^(1/3) [1] 0.2190818+0.3794609i > (-0.084121928394-0i)^(1/3) [1] 0.21908

Re: [R] pattern matching

2008-10-27 Thread John Lande
On Sun, Oct 26, 2008 at 8:06 PM, Duncan Murdoch <[EMAIL PROTECTED]>wrote: > On 26/10/2008 11:54 AM, John Lande wrote: > >> dear all, >> >> I have a little problem I am doing a loop, witha grep function. sometimes >> it >> happens that have the following results >> >> tmp <- grep("x", y) >>> tmp

[R] Odp: Request: Most repeated sequence considering combinations at each row

2008-10-27 Thread Petr PIKAL
Hi not sure if this is what you want. It does not do fuzzy matching but make a exact evaluation equal row sums of arrays. rle(do.call("c",lapply(lapply(l, rowSums), function(x) paste(x, collapse="" Maybe something similar can be done without conversion to character. Regards Petr [EMAIL

[R] Pkg rgl: installation fails because of x11

2008-10-27 Thread Matthieu Stigler
Hello Im trying to install package rgl in a freshly installed Ubuntu 8.04 system. I have a problem (actually is has been reported three times on the R-list but the answers were too complicated for me) when installing: configure: using libpng dynamic linkage checking for X... no configure: error:

Re: [R] Algo. for matrix inverse

2008-10-27 Thread Peter Dalgaard
megh wrote: I am looking fpr a algo to find matrix inverse. Till time I am aware of Gauss-Jordan Elimination procedure to find the same. Are there any other algo. as well? What does R use to find the inverse? There are many algorithms, depending on matrix structure and properties. R interfaces

[R] Question of "Quantile Regression for Longitudinal Data".

2008-10-27 Thread Helen Chen
Hi, I am trying to estimate a quantile regression using panel data. I am trying to use the model that is described in Dr. Koenker's article. So I use the code the that is posted in the following link: http://www.econ.uiuc.edu/~roger/research/panel/rq.fit.panel.R This code run perfectly. Then

Re: [R] help: selection of data

2008-10-27 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 25.10.2008 19:19:45: > #is this what you want? > x <- c(-10, -5, 0, 5, 10, 15, 20) > y <- c(10, 10, 10, -5, -6, -7, 10) > data <- data.frame(x, y) > subset(data, x>0 & y==min(y)) Not exactly. But maybe y.val=which(data$y==min(data$y)) x.val=which(data$x>0) data[x

[R] Bootstrap Panel Data

2008-10-27 Thread Andreas Klein
Hello. I have serveral problems obtaining percentile intervals via bootstrap in panel data: Assuming a TxN panel with T observations across N groups and T>>N. I want to apply the two step Fama/MacBeth procedure (FM), so I have to look, if the mean of the estimated slopes over time is significa