[R] box priece Q statistic

2007-04-19 Thread zahid khan
when we applied box priece Q statistic in stat library how does it work? suppose we have series Yt then either it apply on residuals obtained by regressing yt on Yt-1 .or apply direct on the series Yt - [[alternative HTML version deleted

Re: [R] Accessing R applications from Web

2007-04-19 Thread hadley wickham
> I learned to do this by studying Alfredo Pontillo and Angelo Mineo's > R-php (http://dssm.unipa.it/R-php). Specifically, look at the file > pages/result/index.php They have done lots of other good stuff in > their project, like filtering dangerous commands - you would not want > to let a user ju

Re: [R] Problem with ?curve

2007-04-19 Thread Ron Michael
Dear all R users, Thank you Greg for having time in my problem. Your suggestion is working pretty fine for my problem. But here I am involved for some teaching assignments in a local school. If at least I could use sliders for changing values, it would be very fine. Have anyone any suggestion

[R] glm problem

2007-04-19 Thread rach.s
Hie, I've used GLM function, n try to get a convergence by using a while loop and break it when my absolute value between the 2 values is less than 1e-09. But i got this warning message, what does this mean? how should i correct my algorithm? my start value is 0.5. Error in if (abs((phi_0) - (phi

Re: [R] Windows Vista issue (was Error en R)

2007-04-19 Thread Prof Brian Ripley
This is not an error in R. It is a problem of your permissions under Vista, and was discussed here: https://stat.ethz.ch/pipermail/r-help/2007-April/129661.html a few days ago. On Thu, 19 Apr 2007, Roberto Mario Domínguez Olivero wrote: Mrs and Mr My operative system is "Windows Vista" and

[R] Error en R

2007-04-19 Thread Roberto Mario Domínguez Olivero
Mrs and Mr My operative system is "Windows Vista" and can't update or install pakages of CRAN, the message for all pakages is for example: Erro en zip.unpack(pkg, tmpDir) : no fue posible abrir el archivo 'C:/Program Files/R/R-2.4.1pat/library/file8f653c7/MASS/chtml/MASS.chm' Help me thank yo

[R] sequential for loop

2007-04-19 Thread Michael Toews
Hi all, I'm usually comfortable using the *apply functions for vectorizing loops in R. However, my particular problem now is using it in a sequential operation, which uses values evaluated in an offset of the loop vector. Here is my example using a for loop approach: dat <- data.frame(year=rep

Re: [R] Character coerced to factor and I cannot get it back

2007-04-19 Thread Stephen Tucker
You can also set this option globally with options(stringsAsFactors = TRUE) I believe this was added in R 2.4.0. --- Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > Try this: > > DF <- data.frame(let = letters[1:3], num = 1:3, stringsAsFactors = FALSE) > str(DF) > > > On 4/19/07, John Kane

Re: [R] how to convert the lower triangle of a matrix to a symmetric matrix

2007-04-19 Thread Charles C. Berry
On Thu, 19 Apr 2007, Ranjan Maitra wrote: > Hi, > > I have a vector of p*(p+1)/2 elements, essentially the lower triangle of > a symmetric matrix. I was wondering if there is an easy way to make it > fill a symmetric matrix. I have to do it several times, hence some > efficient approach would b

Re: [R] Problem installing packages

2007-04-19 Thread Vincent Goulet
Le 07-04-19 à 19:02, Peter Dalgaard a écrit : > Dirk Eddelbuettel wrote: >> On Thu, Apr 19, 2007 at 10:51:58PM +0200, Peter Dalgaard wrote: >> >>> Deepankar Basu wrote: >>> I tried to install the package "tseries" but could not do so. Any suggestions are welcome. I am using Ubuntu Linux

Re: [R] Computing an ordering on subsets of a data frame

2007-04-19 Thread Steven McKinney
Hi Lukas, Using by() or its cousins tapply() etc. is tricky, as you need to properly merge results back into X. You can do that by adding a key ID variable to X, and carrying along that key ID variable in calls to by() etc., though I haven't tested out a method. You can also create a new column

Re: [R] Problem installing packages

2007-04-19 Thread Dirk Eddelbuettel
Deepankar, On 19 April 2007 at 21:32, DEEPANKAR BASU wrote: | Thanks Dirk, that worked. | | $ sudo apt-get install r-base-dev r-cran-tseries Glad to be of help, and thanks for reporting back. Now that you have r-base-dev and what it pulled it, your next attempts at installing from CRAN (under s

Re: [R] A particular shuffling on a vector

2007-04-19 Thread Peter Alspach
Emmanuel One option which appears to work: emmanuel <- c(1,1,1,2,2,3,3,3) #emmanuel <- c(1,1,2,3,4,5,6,6,6,6,6,6,6,6,7,8) runs <- rle(emmanuel)[[1]] shuffle <- sample(1:length(runs)) newEmm <- rep(emmanuel[cumsum(runs)[shuffle]], runs[shuffle]) startPos <- sample(1:length(emmanuel), 1) if (start

Re: [R] Problem installing packages

2007-04-19 Thread DEEPANKAR BASU
Thanks Dirk, that worked. $ sudo apt-get install r-base-dev r-cran-tseries Deepankar - Original Message - From: Dirk Eddelbuettel <[EMAIL PROTECTED]> Date: Thursday, April 19, 2007 5:32 pm Subject: Re: [R] Problem installing packages > On Thu, Apr 19, 2007 at 05:19:37PM -0400, Deepanka

[R] how to convert the lower triangle of a matrix to a symmetric matrix

2007-04-19 Thread Ranjan Maitra
Hi, I have a vector of p*(p+1)/2 elements, essentially the lower triangle of a symmetric matrix. I was wondering if there is an easy way to make it fill a symmetric matrix. I have to do it several times, hence some efficient approach would be very useful. Many thanks and best wishes, Ranjan _

Re: [R] R in cron job: X problems

2007-04-19 Thread Liaw, Andy
This is in the FAQ, if I remember correctly... However, alternatively: As Jeff Horner recently pointed out on the list, the Cairo package is a good way of generating png without needing an X display. You may want to look into that. I've just installed cairo on our CentOS boxes and the Cairo pac

Re: [R] R in cron job: X problems

2007-04-19 Thread Henrik Bengtsson
Try using png2() in R.utils, which immitates png() but uses bitmap() and ghostscript to create the PNG file. You need to set 'R_GSCMD' to tell R where ghostscript is located - you can use System$findGhostscript() at startup to let R try to locate ghostscript for you. /H On 4/19/07, Mark Liberman

[R] A particular shuffling on a vector

2007-04-19 Thread Emmanuel Levy
Hello, I was wondering if anyone can think of a straightforward way (without loops) to do the following shuffling: Let's imagine a vector: c(1,1,1,2,2,3,3,3) I would like to derive shuffled vectors __where the same digits are never separated__, although they can be at both ends (periodicity). S

[R] erratic behavior of match()?

2007-04-19 Thread Alan Zaslavsky
>Is this a consequence of machine error or something else? >Could this be overcome? (It works correctly when integers are used in >the sequences as well as in many other circumstances) The usual solution for testing a==b with floating-point round-off error is abs(a-b) X1=seq(0,1,len=11) > X2=seq(

[R] R in cron job: X problems

2007-04-19 Thread Mark Liberman
I'd like to use an R CMD BATCH script as part of a chron job that is set up to run every hour. The trouble is that the script creates a graphical output in a file via png(), and apparently this in turn works through X. When cron invokes the job, no X server is available -- I suppose that the DISP

[R] Loading Dynamic Libraries on OS X

2007-04-19 Thread Saptarshi Guha
Hello, I am playing around with the R_X11.c source code and wanted to load the compiled DLL via dyn.load. It loads properly (and masks the old X11 device) when i load it in the console R (version: version.string R version 2.4.1 (2006-12-18) ) However the same command doe

Re: [R] Filtering

2007-04-19 Thread Gabor Grothendieck
There are some examples here: library(zoo) ?aggregate.zoo and also in vignette("zoo") On 4/19/07, Soare Marcian-Alin <[EMAIL PROTECTED]> wrote: > Hello Everybody, > > How can I filter a dataset? > > Im trying to filter the dataset EuStockMarkets monthly and quarter. > > data(EuStockMarkets) >

Re: [R] Matlab import

2007-04-19 Thread Henrik Bengtsson
On 4/19/07, Spencer Graves <[EMAIL PROTECTED]> wrote: > Have you tried 'getVariable' from Matlab? I also recently failed > to get 'readMat' to work for me -- probably because I hadn't saved the > files using the options Henrik suggested below. > > Fortunately, I was able to get somethi

Re: [R] Character coerced to factor and I cannot get it back

2007-04-19 Thread Gabor Grothendieck
Try this: DF <- data.frame(let = letters[1:3], num = 1:3, stringsAsFactors = FALSE) str(DF) On 4/19/07, John Kane <[EMAIL PROTECTED]> wrote: > > --- Tyler Smith <[EMAIL PROTECTED]> wrote: > > > I really need to sit down with the manual and sort > > factors and classes > > properly. In your case,

[R] Filtering

2007-04-19 Thread Soare Marcian-Alin
Hello Everybody, How can I filter a dataset? Im trying to filter the dataset EuStockMarkets monthly and quarter. data(EuStockMarkets) ftse = EuStockMarkets[,4] Alin [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing lis

Re: [R] Problem installing packages

2007-04-19 Thread Peter Dalgaard
Dirk Eddelbuettel wrote: > On Thu, Apr 19, 2007 at 10:51:58PM +0200, Peter Dalgaard wrote: > >> Deepankar Basu wrote: >> >>> I tried to install the package "tseries" but could not do so. Any >>> suggestions are welcome. I am using Ubuntu Linux and here is what I did: >>> >>> $ sudo R >>>

[R] Package loading shows missing libraries in Mac OS X

2007-04-19 Thread David Daniel
I've been using R actively for only the past few months under Mac OS X, but have found a few situations where I cannot get a package to load. This seems like the same situation reported recently on the thread "R-2.4.1 for MacOS X - languageR, acepack, Hmisc": > > library(Hmisc) > Error in d

Re: [R] Matlab import

2007-04-19 Thread Spencer Graves
Have you tried 'getVariable' from Matlab? I also recently failed to get 'readMat' to work for me -- probably because I hadn't saved the files using the options Henrik suggested below. Fortunately, I was able to get something like the following to work: # Start Matlab server Matla

Re: [R] rbind() of factors in data.frame

2007-04-19 Thread Prof Brian Ripley
On Thu, 19 Apr 2007, Albrecht, Dr. Stefan (AZ Private Equity Partner) wrote: > I would like to inquire, if it is a desired feature that the combination > with rbind() of two data frames with factors columns does not sort the > factors levels of the combined data frame. Yes, and a documented one

Re: [R] Problem installing packages

2007-04-19 Thread Dirk Eddelbuettel
On Thu, Apr 19, 2007 at 05:19:37PM -0400, Deepankar Basu wrote: > While installing R (version 2.4.0) on my Dell desktop (Optiplex GX 270) > on a Ubuntu Linux platform, I had followed the advice contained in the > following link: > http://cran.r-project.org/bin/linux/debian/README > > Briefly, what

Re: [R] erratic behavior of match()?

2007-04-19 Thread Henrik Bengtsson
...and so say google [http://www.google.com/search?q=1%250.1]: "1 modulo 0.1 = 0.1", so end of discussion ;) In bit of a food coma now, but the following is interesting: r = a %% b <=> r = (b*a/b) %% (b*b/b) <=> r = b*((a/b) %% 1) > modulo <- function(a, b) { b * ((a/b) %% 1) }

Re: [R] Problem installing packages

2007-04-19 Thread Dirk Eddelbuettel
On Thu, Apr 19, 2007 at 09:53:04PM +0100, Prof Brian Ripley wrote: > 'grid' is an integral part of R, so this suggests your R installation is > broken. (You could have found that out from the R FAQ Q5.1.1.) Just for the record, it is close to impossible to detach grid from R. Both are in the r-b

Re: [R] inconsistent output using 'round'

2007-04-19 Thread Bert Gunter
It has nothing to do with round() -- it's the "digits" argument of the print method that controls the number of digits in the output, print.default in this case. And the documentation from print.default says for the digits argument: "digits: a non-null value for digits specifies the minimum number

Re: [R] Problem installing packages

2007-04-19 Thread Dirk Eddelbuettel
On Thu, Apr 19, 2007 at 10:51:58PM +0200, Peter Dalgaard wrote: > Deepankar Basu wrote: > > I tried to install the package "tseries" but could not do so. Any > > suggestions are welcome. I am using Ubuntu Linux and here is what I did: > > > > $ sudo R > > > > Then once in R, I did > > > > > >>

Re: [R] general question about plotting multiple regression results

2007-04-19 Thread John Fox
Dear Thomas and Simon, On Thu, 19 Apr 2007 07:37:28 -0700 (PDT) Thomas Lumley <[EMAIL PROTECTED]> wrote: > On Thu, 19 Apr 2007, Simon Pickett wrote: > > > Hi all, > > > > I have been bumbling around with r for years now and still havent > come up > > with a solution for plotting reliable graphs

Re: [R] Problem installing packages

2007-04-19 Thread Deepankar Basu
While installing R (version 2.4.0) on my Dell desktop (Optiplex GX 270) on a Ubuntu Linux platform, I had followed the advice contained in the following link: http://cran.r-project.org/bin/linux/debian/README Briefly, what I had done is the following: I had added "deb http://cran.R-project.org/bi

Re: [R] inconsistent output using 'round'

2007-04-19 Thread Bob Green
Peter, Many thanks. I have never seen a confidence interval from 0.000 to 6265941604681544800.000 - this is a worry. I am also still puzzled why use of digits = 3, produced output which includes 2, 3 and 4 decimal points as per below. The two decimal point values for the coef should h

Re: [R] erratic behavior of match()?

2007-04-19 Thread Duncan Murdoch
On 4/19/2007 4:29 PM, Bernhard Klingenberg wrote: > Thank you! Is floating point arithmetic also the reason why > > 1 %% 0.1 > > gives the "surprising" answer 0.1 (because 0.1 cannot be written as a > fraction with denominator a power of 2, e.g. 1%%0.5 correctly gives 0). > > This seems to go

Re: [R] rbind() of factors in data.frame [Broadcast]

2007-04-19 Thread Raubertas, Richard
Please, no. It is already annoying enough that levels are sorted when creating a factor. Don't compound it by extending this to other functions. In concept the order of the levels of a factor is irrelevant (although in practice it makes a big difference, e.g. when plotting). If so, then why i

Re: [R] Matlab import

2007-04-19 Thread Henrik Bengtsson
Hi, as already mentioned, do not save MAT files in ASCII format but save to binary formats, i.e. do *not* use "-ascii". Moreover, from ?readMat, you find that: "From Matlab v7, _compressed_ MAT version 5 files are used by default [3]. These are not supported. Use 'save -V6' in Matlab t

Re: [R] Error loading libraries in MAC

2007-04-19 Thread Matthew Keller
HI Mayte & Gregory, This is probably a question needing to be posted to r-sig-mac. A search for this problem on that forum turns up lots of hits. I think everyone is having these problems (which makes me pause about whether to switch to Mac, given how much I use R). Below is a message from the r-s

Re: [R] Problem installing packages

2007-04-19 Thread Prof Brian Ripley
'grid' is an integral part of R, so this suggests your R installation is broken. (You could have found that out from the R FAQ Q5.1.1.) You have told us very little (please study the posting guide), not including which version of R and how you installed it. Please try reinstalling R, from th

Re: [R] Problem installing packages

2007-04-19 Thread Peter Dalgaard
Deepankar Basu wrote: > I tried to install the package "tseries" but could not do so. Any > suggestions are welcome. I am using Ubuntu Linux and here is what I did: > > $ sudo R > > Then once in R, I did > > >> install.packages("tseries", dep=TRUE) >> > > which gave a long list of errors.

Re: [R] erratic behavior of match()?

2007-04-19 Thread Alberto Monteiro
Bernhard Klingenberg wrote: > > Thank you! Is floating point arithmetic also the reason why > > 1 %% 0.1 > > gives the "surprising" answer 0.1 > Floating point arithmetic is one of the most Evil things computer science imposed against the mathematicians; any reasonable person, back when comput

Re: [R] erratic behavior of match()?

2007-04-19 Thread Bernhard Klingenberg
Thank you! Is floating point arithmetic also the reason why 1 %% 0.1 gives the "surprising" answer 0.1 (because 0.1 cannot be written as a fraction with denominator a power of 2, e.g. 1%%0.5 correctly gives 0). This seems to go a bit against the statement in the help for '%%', which states "F

[R] k-means and PAM

2007-04-19 Thread nathaniel Grey
Hi, I would like to compare results from a k-means analysis to that done with PAM. However, when I ran the k-means on four clusters it gave me the means for each result with the number of each cluster e.g. variable1 1 33.33 2 12.00 3 109.00 4 67.98 So I know for instance cluster 1 has a

[R] Problem installing packages

2007-04-19 Thread Deepankar Basu
I tried to install the package "tseries" but could not do so. Any suggestions are welcome. I am using Ubuntu Linux and here is what I did: $ sudo R Then once in R, I did > install.packages("tseries", dep=TRUE) which gave a long list of errors. Going through them, I saw that the basic problem w

[R] Fwd: RE: Character coerced to factor and I cannot get it back

2007-04-19 Thread John Kane
t;,"cc","st") > > dd1 > > typeof(dd1$st); class(dd1$st) > > > > dd2 <- cbind(dd, as.character(st)) > > names(dd2) <- c("aa","bb","cc","st&quo

Re: [R] Character coerced to factor and I cannot get it back

2007-04-19 Thread John Kane
--- Tyler Smith <[EMAIL PROTECTED]> wrote: > I really need to sit down with the manual and sort > factors and classes > properly. In your case, I think the problem has > something to do with > the way a list behaves? I'm not sure, but if you > convert your list to > a dataframe it seems to work

Re: [R] Character coerced to factor and I cannot get it back

2007-04-19 Thread Tyler Smith
I really need to sit down with the manual and sort factors and classes properly. In your case, I think the problem has something to do with the way a list behaves? I'm not sure, but if you convert your list to a dataframe it seems to work ok: > dd3 <- as.data.frame(dd1) > typeof(dd3$st) [1] "inte

Re: [R] rJava help

2007-04-19 Thread Prof Brian Ripley
On Thu, 19 Apr 2007, tkobayas wrote: > I use Redhat Linux for R. Now I am trying to use the rJava library. I > downloaded JRE beta 2 from Sun. I typed I am not sure why, as JRE 1.6.0 was released a while back. > [EMAIL PROTECTED] R CMD javareconf In what R version? > Java interpreter : /usr/bi

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Deepankar Basu
Ravi, Thanks a lot for that clarification. Now I see that the code works. Deepankar On Thu, 2007-04-19 at 14:01 -0400, Ravi Varadhan wrote: > Hi Deepankar, > > Dimitris' code works just fine. Your problem is that the output of optim > does not have a corresponding "summary" method. Instead yo

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Ravi Varadhan
Hi Deepankar, Dimitris' code works just fine. Your problem is that the output of optim does not have a corresponding "summary" method. Instead you should simply type the name of the object returned by "optim" to look at the results. > out <- optim(mu.start, mlogl, method = "CG", y = women$J, X

Re: [R] rJava help

2007-04-19 Thread tkobayas
Hi, I use Redhat Linux for R. Now I am trying to use the rJava library. I downloaded JRE beta 2 from Sun. I typed [EMAIL PROTECTED] R CMD javareconf Java interpreter : /usr/bin/java Java version : 1.6.0-beta2 Java home path : /usr/java/jre1.6.0 Java library path: $(JAVA_HOME)/lib/amd64/se

Re: [R] Histogram with uneven bins

2007-04-19 Thread jim holtman
will this work for you? x <- runif(1000, 0, 40) x.c <- cut(x, breaks=c(0, 1, 3, 6, 10, Inf)) barplot(table(x.c)) On 4/19/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > Hi R-helpers > > I would like to produce a histogram with uneven bins (e.g., 0, 1-2, 3-5, > 6-10, 10-20, >20) but I would li

Re: [R] Problem with ?curve

2007-04-19 Thread Greg Snow
Something along these lines is on my todo list for the TeachingDemos package, so any thoughts on what you want the user interface to be would be helpful (do you want sliders to set the x-value? Do you want to click on the plot to see an update? Or just moving the mouse over the plot without clicki

Re: [R] division of decimal number

2007-04-19 Thread Prof Brian Ripley
On Thu, 19 Apr 2007, Gavin Simpson wrote: > On Wed, 2007-04-18 at 16:06 +0100, Barry Rowlingson wrote: >> Schmitt, Corinna wrote: >>> Dear R-Experts, >>> >>> how can I divide the number 0.285 with 2. I need a function. >>> Result: 0.285 / 2 = 0.1425 >> >> Just get the / operator: >> >> > divide

Re: [R] Problems in programming a simple likelihood

2007-04-19 Thread Deepankar Basu
Dimitris, Thanks a lot for your suggestion and also for suggestions that others have provided. I am learning fast and with the help of the R community will be able to get this going pretty soon. Of course, right now I am just trying to learn the language; so I am trying to program a standard probi

Re: [R] Accessing R applications from Web

2007-04-19 Thread Bos, Roger
I guess you would have to do all those step for a Internet based portal, but if you just want to have a web interface on a _Intranet_ portal where you trust your users, you can use Rpad which is much easier in that you don't really have to write any code, just do the install procedure. There is a

Re: [R] Error with strptime

2007-04-19 Thread Stephen Tucker
you have to use "POSIXct" classes to include date-time objects into data frames. strptime() returns an object of class of "POSIXlt". when you do the cbind(), it automatically converts test2 into "POSIXct" you probably want bsamp$spltime<-as.POSIXct(strptime(test,format="%d-%B-%y %H:%M")) (but ple

Re: [R] Conditional power, predictive power

2007-04-19 Thread francogrex
Ok I am replying to my own message! I wrote a "function", it works well but it's a bit twisted because you will have to edit the last file in excel or other. This is to analyze the bayesian predictive power in an analysis where treatment x is compared to treatment y. Example: Total final subject n

[R] Character coerced to factor and I cannot get it back

2007-04-19 Thread John Kane
Something probably obivous but I don't see it. I needed to find the first 1 or two digits of some 5 and 6 digit numbers since they identified research sites while the rest of the number was the plot id. I converted the numbers to characters, got the first 1 or 2 characters as appropriate and w

Re: [R] macro sous R

2007-04-19 Thread Julien Barnier
Hi, > hello I'd like know if it's possible to write macro variables under > R as we know do it with SAS thanks For what I know of both softwares, R and SAS are two very different languages. You don't have "Macros" in R, because you don't need them : you use functions, which are much more powerful

Re: [R] division of decimal number

2007-04-19 Thread Gavin Simpson
On Wed, 2007-04-18 at 16:06 +0100, Barry Rowlingson wrote: > Schmitt, Corinna wrote: > > Dear R-Experts, > > > > how can I divide the number 0.285 with 2. I need a function. > > Result: 0.285 / 2 = 0.1425 > > Just get the / operator: > > > divide = get("/") > > > > divide(0.285,2) > [1

Re: [R] Accessing R applications from Web

2007-04-19 Thread John McNeil
Sarthi, I have done what you describe using PHP and R on a linux server. Basically, you follow these steps: - Create a form to gather input - Parse input from the form with PHP into a file that R can read - Run R commands from PHP with exec(). - Convert the postscript output to an image a web br

[R] rbind() of factors in data.frame

2007-04-19 Thread Albrecht, Dr. Stefan (AZ Private Equity Partner)
Dear all, I would like to inquire, if it is a desired feature that the combination with rbind() of two data frames with factors columns does not sort the factors levels of the combined data frame. > str(rbind(data.frame(a = factor(c(4, 3))), data.frame(a = factor(c(2, 1) 'data.frame': 4

[R] How does max.lag work?

2007-04-19 Thread yoooooo
Let's say I have some data DIFF, then I do > a <- acf(DIFF, lag.max=2) > a$acf , , 1 [,1] [1,] 1.000 [2,] -0.2085724 [3,] -0.1597789 > a <- acf(DIFF) > a$acf , , 1 [,1] [1,] 1. [2,] -0.20857236 [3,] -0.15977892 How come they are different? I thought the

Re: [R] A structured output from for(i...)?

2007-04-19 Thread francogrex
OK it works thanks. francogrex wrote: > Hello I am using the "for (i...)" and a sink() into a file. But the output > I > am having is not arranged in either a vector or any other good structure. > I > would like to have the output in a file directly as a vector so that I do > not have to edit th

[R] MNP,multinomial probit model

2007-04-19 Thread Ozlem Ipekci
I am an undergraduate student in Isik University, in Istanbul/TURKEY. For my term project, I am studying on the relationship between the prices and market shares of 3 automobile models. I need to find estimates for the parameters so that I can continue my study. I found about your MNP package a

Re: [R] read.table with more than one sep

2007-04-19 Thread Gabor Grothendieck
Try this: # test data Input <- "4547;1970.01.01 00:00-1970.01.01 01:00; noData 4547;1970.01.01 00:00-1970.01.01 01:00; noData" # replace next line with Lines <- readLines("myfile.dat") Lines <- readLines(textConnection(Input)) Lines <- gsub("[;-]", " ", Lines) read.table(textConnection(Lines

[R] Histogram with uneven bins

2007-04-19 Thread mtb954
Hi R-helpers I would like to produce a histogram with uneven bins (e.g., 0, 1-2, 3-5, 6-10, 10-20, >20) but I would like the resulting bars to be the same width ( i.e., a bar's width would not be proportional to its corresponding bin size). Also, the x and y axes of my histograms frequently (almo

Re: [R] Using "mean" if two values are identical

2007-04-19 Thread Stephen Tucker
## making data up # make matrix with some equal values > mat <- cbind(x=rnorm(10),y=rnorm(10),z=rnorm(10)) > mat[c(8,9),"y"] <- mat[c(1,7),"x"] > mat x y z [1,] 0.26116849 0.5823529 -0.96924020 [2,] -0.21415406 0.1085396 2.00542549 [3,] 0.56890081 -1.25263

Re: [R] macro sous R

2007-04-19 Thread John Kane
Introduction to R Ch. 10. Writing your own functions. --- elyakhlifi mustapha <[EMAIL PROTECTED]> wrote: > hello I'd like know if it's possible to write macro > variables under R as we know do it with SAS > thanks __ R-help@stat.math.ethz.ch mailing

[R] Error with strptime

2007-04-19 Thread Jean-Louis Abitbol
Dear All, I am trying to convert to POSIXct after pasting a date and a time in character format with strptime. It is probably obvious but I don't understand why I get an error message after bsamp$spltime<-strptime(test,format="%d-%B-%y %H:%M") whereas I can get what I want if I do it in 2 steps

[R] macro sous R

2007-04-19 Thread elyakhlifi mustapha
hello I'd like know if it's possible to write macro variables under R as we know do it with SAS thanks - [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat

Re: [R] Using "mean" if two values are identical

2007-04-19 Thread Julien Barnier
Hi, > I've got a matrix (mail end) with the colnames x, y, z. In this matrix > are different measurements. x and y are risign coordinates. > > My question. Always, if the "x" AND "y" coordinates are the same, I want to > get the mean of their z values. Maybe you could try with the aggregate func

Re: [R] importing excel-file

2007-04-19 Thread Alberto Monteiro
Hans-Peter wrote: > >> This works in most cases. > > Huu? (the package is supposed to work in *all* cases...!). Do you use > the newest version (v1.3.1 or v1.3.2)? If there are any bugs/issues, > please report them to me and they - most likely - will get fixed. > Here a

Re: [R] general question about plotting multiple regression results

2007-04-19 Thread Thomas Lumley
On Thu, 19 Apr 2007, Simon Pickett wrote: > Hi all, > > I have been bumbling around with r for years now and still havent come up > with a solution for plotting reliable graphs of relationships from a > linear regression. termplot() does this for a range of regression models (without interaction

Re: [R] Accessing R applications from Web

2007-04-19 Thread Jeffrey Horner
d. sarthi maheshwari wrote: > Hi > > I am trying to provide web interface to my R application. My requirements > are simple and can be written down as follows : > > 1) User must provide all the input parameters on web interface. > 2) Based on the input values, R will do some computations. > 3) Pr

[R] read.table with more than one sep

2007-04-19 Thread Jan . Schwanbeck
Dear List, somebody knows, if the following operation can be done in an easier way? The data rows which should be read into R look like that: 4547;1970.01.01 00:00-1970.01.01 01:00; noData Unitil now we are doing this procedure: 1. Dividing columns which are separated by ";" => temp <- re

Re: [R] convergence

2007-04-19 Thread Ravi Varadhan
No. --- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Health Division of Geriatric Medicine and Gerontology Johns Hopkins University Ph: (410) 502-2619 Fax: (410) 614-9625 Email: [EMAIL PRO

Re: [R] hello

2007-04-19 Thread Dimitris Rizopoulos
look at R FAQ 7.10 http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-do-I-convert-factors-to-numeric_003f I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgiu

Re: [R] hello

2007-04-19 Thread rolf
[EMAIL PROTECTED] wrote: > ?as.numeric in response to > ok how can I do to cancel le function factor for the moment ican't do > this and I'm stuck with my problem for exemple ... The use of as.numeric() to convert a factor to a numeric vector is fraught wi

Re: [R] plotting command trouble

2007-04-19 Thread Schmitt, Corinna
Thanks. I will try it. Perhaps it works. Corinna -Ursprüngliche Nachricht- Von: Jim Lemon [mailto:[EMAIL PROTECTED] Gesendet: Donnerstag, 19. April 2007 15:16 An: Schmitt, Corinna Betreff: Re: [R] plotting command trouble Schmitt, Corinna wrote: > I know, but it is needed. Any other i

Re: [R] convergence

2007-04-19 Thread Alberto Monteiro
Ted Harding wrote: > > There are various ways round this, but a 'for' loop with > a fixed number of iterations is not usully one of them! > > The simplest is to use while(). A possibly strategy is > > Y.old <- initial.Y > while(TRUE){ > Y <- compute.Y(Y.old, ...) > if(abs(Y - Y.old)

[R] Cohen's d on Adjusted Means (UNCLASSIFIED)

2007-04-19 Thread Bliese, Paul D LTC USAMH
Classification: UNCLASSIFIED Caveats: NONE I have a large group randomized trial (pre-post design) where the randomization was marginally successful. Given the pre-existing differences among groups, it makes sense to report adjusted means (aka least squares means though I estimated them via pre

[R] Using "mean" if two values are identical

2007-04-19 Thread Felix Wave
Hello, I have got a question. I've got a matrix (mail end) with the colnames x, y, z. In this matrix are different measurements. x and y are risign coordinates. My question. Always, if the "x" AND "y" coordinates are the same, I want to get the mean of their z values. e.q. " x" AND "y" in line

Re: [R] hello

2007-04-19 Thread MORNEAU François
?as.numeric -Message d'origine- De : [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] De la part de elyakhlifi mustapha Envoyé : jeudi 19 avril 2007 15:30 À : R-help@stat.math.ethz.ch Objet : [R] hello ok how can I do to cancel le function factor for the moment ican't do this and I'm stuck w

Re: [R] plotting command trouble

2007-04-19 Thread Philipp Pagel
On Thu, Apr 19, 2007 at 02:10:37PM +0200, Schmitt, Corinna wrote: >> Hi, the length of the coordinates are different. >> >> Try: >> plot(seq(0,60, l=100), seq(0,0.896, l=100), type="n", xlab="Zeit [min]", >> ylab="Absorptionsmessung >>bei 600nm",main="Zellwandstabilit?t" ) > I know, but it

Re: [R] matrix building with two variables

2007-04-19 Thread Hans-Peter
2007/4/19, Schmitt, Corinna <[EMAIL PROTECTED]>: > Thanks. > > I did read everything I could but could not understand everything. Hopefully > with more programming practice it will become more less. Here is a great book: Uwe Ligges: Programmieren mit R. worth every "Rappen" :-) _

[R] hello

2007-04-19 Thread elyakhlifi mustapha
ok how can I do to cancel le function factor for the moment ican't do this and I'm stuck with my problem for exemple Id_Rep Id_Geno Val_O Id_TrT1 1 1 64238 145 2 2 1 64238 1 1 3 1 64238 2 1 4 1 64238 1 1 5 1 64

Re: [R] "tree-ID" in any segmentation package available?

2007-04-19 Thread Florian Koller-Meinfelder
Thanks Torsten, That seems to fit the bill nicely! Best, Florian Florian Koller-Meinfelder Research Consulting & Development __ GfK Fernsehforschung GmbH Nordwestring 101 90319 Nürnberg Tel +49 (0)911 395-3554 Fax +49 (0)911 395-4130 www.gfk.com/gfkfernseh

Re: [R] inconsistent output using 'round'

2007-04-19 Thread Peter Dalgaard
Bob Green wrote: > I am hoping for some advice regarding limiting decimal points to 3. > 'Round' produces the desired results except for the 97.5% confidence interval. > Any advice as to how I modify the code to obtain output to 3 decimal > points for all ouput is appreciated, > The rounding

Re: [R] "tree-ID" in any segmentation package available?

2007-04-19 Thread Torsten Hothorn
On Thu, 19 Apr 2007, Florian Koller-Meinfelder wrote: Dear R-helpers, I am looking for a segmentation package that gives some "tree identifier" as output for every observation in the data set (my response variable is binary). I have skimmed through "rpart", "ada" and "adabag": The output "tree

Re: [R] convergence

2007-04-19 Thread Jim Lemon
rach.s wrote: > hie.. > how can i write a loop that makes algorithm keeps repeating until a solution > is converged?do i use a for loop? i know that we can use for loop to ask for > a number of repetitions, but how to use it to ask the algorithm to keep > repeating until a solution is converged? g

Re: [R] convergence

2007-04-19 Thread Doran, Harold
No, use a while loop. Something like change <- 1 while(abs(change) > .001 ){ do stuff change <- updated change } > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of rach.s > Sent: Thursday, April 19, 2007 8:00 AM > To: r-help@stat.math.

Re: [R] convergence

2007-04-19 Thread Ted Harding
On 19-Apr-07 12:00:17, rach.s wrote: > > hie.. > how can i write a loop that makes algorithm keeps repeating > until a solution is converged?do i use a for loop? i know > that we can use for loop to ask for a number of repetitions, > but how to use it to ask the algorithm to keep repeating > until

Re: [R] general question about plotting multiple regression results [Broadcast]

2007-04-19 Thread Liaw, Andy
I suspect you'll greatly benefit a read of Prof. Fox's book(s) on regression models, as well as making use of his "car" package. You may want to read up on partial residual plots and partial regression plots. Andy From: Simon Pickett > > Hi all, > > I have been bumbling around with r for years

Re: [R] is there a function to give significance to correlation?

2007-04-19 Thread Ranjan Maitra
Hi Jenny, On Thu, 19 Apr 2007 12:48:50 +0100 (BST) Jenny Barnes <[EMAIL PROTECTED]> wrote: > Dear R-Help, > > I am trying to find a function that will give me the significance of the > correlation of 2 variables (in the same dimension arrays) correcting for > serial > autocorrelation. I am n

Re: [R] matrix building with two variables

2007-04-19 Thread Julien Barnier
Hi, > I did read everything I could but could not understand everything. Hopefully with more programming > practice it will become more less. Then maybe you read everything you could a bit too fast. Because the answer to your question is in the first document to read, "An introduction to R", sect

  1   2   >