What I wondered is if I pass the 'upper=~.' ,
it seems step() thinks the full model is current one. Not adding
anymore. If this is the right answer, is there a better way than
creating fmla argument in the above?
Yes, that is exactly what the help page for step says it means. (So why
are
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.
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
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,-1), N,
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 -
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 =
start,
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
x - runif(N,
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
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,
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 of the
Y
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 part of
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
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)
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
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
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: r-help@stat.math.ethz.ch
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
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
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
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
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
Sent:
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
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 Gunter
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
23 matches
Mail list logo