Re: [R] probability-probability plot

There are PP plots and QQ plots. I've tried PP plots and didn't get much from them. The normal QQ plot [qqnorm, in R], however, I use for all kinds of things. In a data mining situation, I'll get, e.g., 500 p.values, all uniformily distributed under the null hypothesis. I'll trans

Re: [R] jagged array

Access list items using "[[". "res[1]" is a list with only one attribute. "res[1][1]" is the same list with only one attributre. "res[[1]]" is the first attribute of that list, so res[[1]][1] in your example should be 3.125. spencer graves Tolga Uzuner wrote: > Berton Gu

Re: [R] plot question

Trying characters and expressions variously it seems that font.lab applies to character strings but not to expressions so if you want to use an expression just use bold (or whatever) explicitly on the expression. One gotcha is that bold will not work as one might have expected on numbers so they m

Re: [R] jagged array

Berton Gunter wrote: >Please read "An Introduction to R" and the help file for "[", where this is >explained. If you want to use R, you need to expend effort to learn it. > >-- Bert Gunter > >-Original Message- >From: [EMAIL PROTECTED] >[mailto:[EMAIL PROTECTED] On Behalf Of Tolga Uzuner

Re: [R] jagged array

Please read "An Introduction to R" and the help file for "[", where this is explained. If you want to use R, you need to expend effort to learn it. -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Tolga Uzuner Sent: Saturday, July 02, 2005

[R] jagged array

Hi, I have a "jagged array" which looks like this: > res ... [[996]] [1] 1.375 3.375 4.125 4.625 4.875 4.875 6.625 7.125 8.875 [[997]] [1] 1.875 5.125 6.875 7.875 9.875 [[998]] [1] 1.875 5.375 6.625 6.875 8.125 9.375 9.625 [[999]] [1] 1.875 6.875 9.875 [[1000]] [1] 1.875 6.125 6.875 9.375 9.

[R] how to set the position and size of the select.list window

Hello, I use "select.list" to obtain a window of select items, but how can I set the position and size of this window? Thank you! shengzhe __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] Specifying the minimum and maximum of x and y-axis in plot

see ?plot.default for example: plot(1:10,1:10,ylim=c(0,10),xlim=c(0,20)) - Original Message - From: "Tolga Uzuner" <[EMAIL PROTECTED]> To: Sent: Saturday, July 02, 2005 9:39 PM Subject: [R] Specifying the minimum and maximum of x and y-axis in plot > Hi, > > How do I specify the mini

[R] Specifying the minimum and maximum of x and y-axis in plot

Hi, How do I specify the minimum and maximum of the x and y-axis when using plot ? Thanks, Tolga __ 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-g

Re: [R] probability-probability plot

On 07/02/05 15:41, Yulei He wrote: Hi, there. Is there any function in R to plot the probability-probability plot (PP plot)? Suppose I am testing some data against normal. qqnorm might be what you want, or lead to it. -- Jonathan Baron, Professor of Psychology, University of Pennsylvania Ho

[R] plot question

dear list: in the following plot: plot(rnorm(10),rnorm(10),xlab="year",ylab=expression (paste('M x'*10^{3},)),font.lab=2) font.lab=2, but xlab and ylab are different. I want both labels in the same way. help? a.d. - Email Enviado utilizando o se

[R] probability-probability plot

Hi, there. Is there any function in R to plot the probability-probability plot (PP plot)? Suppose I am testing some data against normal. Thanks. Yulei $$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 [EMAIL PROTECTED] 734-647-0305(H) 734-763-04

Re: [R] Is it possible to use glm() with 30 observations?

I agree with Ted: "in model-fitting terms, it is a resounding success!" With any data set having at least one point with a binomial yield of 0 or 100%, you can get this phenomenon by adding series of random numbers sequentially to a model. Eventually, you will add enough variables t

Re: [R] interrupted Y axis

On Sat, 2005-07-02 at 06:09 -0700, Nick Drew wrote: > I did not find an answer to my question after a quick > search using the R > search engine so thought I'd ask away: > > Does any know if there's a function exists to create > an interrupted Y axis? > What I mean by interrupted Y axis is that pa

Re: [R] interrupted Y axis

On 7/2/05, Nick Drew <[EMAIL PROTECTED]> wrote: > I did not find an answer to my question after a quick > search using the R > search engine so thought I'd ask away: > > Does any know if there's a function exists to create > an interrupted Y axis? > What I mean by interrupted Y axis is that part o

[R] Cluster of iris data set from Mahalanobis distances

Dear R list, My question: I'm trying to calculate Mahalanobis distances for 'Species' from the iris data set as below: cat('\n') cat('Cluster analyse of iris data\n') cat('from Mahalanobis distance obtained of SAS\n') cat('\n') n = 3 dat = c(0, 89.86419,

[R] interrupted Y axis

I did not find an answer to my question after a quick search using the R search engine so thought I'd ask away: Does any know if there's a function exists to create an interrupted Y axis? What I mean by interrupted Y axis is that part of the Y axis has been removed or excised to permit one to see

Re: [R] Generating correlated data from uniform distribution

On 02-Jul-05 Peter Dalgaard wrote: > "Jim Brennan" <[EMAIL PROTECTED]> writes: > >> OK now I am skeptical especially when you say in a weird way:-) >> This may be OK but look at plot(x,y) and I am suspicious. Is it still >> alright with this kind of relationship? > ... >> N <- 1 >> rho <- .6 >

Re: [R] Is it possible to use glm() with 30 observations?

On 02-Jul-05 Kerry Bush wrote: > I have a very simple problem. When using glm to fit > binary logistic regression model, sometimes I receive > the following warning: > > Warning messages: > 1: fitted probabilities numerically 0 or 1 occurred > in: glm.fit(x = X, y = Y, weights = weights, start = >

Re: [R] Is it possible to use glm() with 30 observations?

On 2 Jul 2005, at 06:01, Spencer Graves wrote: > The issue is not 30 observations but whether it is possible to > perfectly separate the two possible outcomes. Consider the following: > > tst.glm <- data.frame(x=1:3, y=c(0, 1, 0)) > glm(y~x, family=binomial, data=tst.glm) > > tst2.glm <-

Re: [R] Generating correlated data from uniform distribution

"Jim Brennan" <[EMAIL PROTECTED]> writes: > OK now I am skeptical especially when you say in a weird way:-) > This may be OK but look at plot(x,y) and I am suspicious. Is it still > alright with this kind of relationship? ... > N <- 1 > rho <- .6 > x <- runif(N, -.5,.5) > y <- x * sample(c(1,-

[R] how to call sas in R

Hello all, I would like to know how to call sas code in R. Since I simulate data in R and I need to use sas code (garch-t,egarch and gjr) to estimate it. I need to simulate 500 times with 2000 obs. How I can call that code in R.Also, how I can keep the parameters from the estimate. j=1:500 i=1:200

[R] Generating correlated data from uniform distribution

While you are looking at weird distributions, here is one that we have used in experiments on noise masking to explore the bandwidth of visual mechanisms D'Zmura, M., & Knoblauch, K. (1998). Spectral bandwidths for the detection of color. Vision Research, 20, 3117-28 and G. Monaci, G. Meneg