-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf
Of Samuel Okoye
Sent: Friday, March 14, 2008 8:21 AM
To: [EMAIL PROTECTED]
Subject: [R] SAS data
Hello,
I am trying to read the SAS file MyData.sa7bdat in R! This file is saved
under
D:\data! I
Hello, it seems www.rmetrics.org has been unavailable for some days.
Would anyone know of a mirror? Thanks, Ruud
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
At 15:52 14/03/2008, [EMAIL PROTECTED] wrote:
Hi,
I would like to perform a forward selection procedure on a data set
with 6 observations and 10 predictors. I tried to run it with
regsubsets (I set nvmax=number of observations) but I keep getting
these warning messages:
Warning messages:
1: 5
Hi
I spent hours looking over my formula. Somehow I cant find the reason
why it gives me different answer.
help appreciated.
x =
as.matrix(read.table(http://www.niehs.nih.gov/research/atniehs/core/microarrays/docs/heinloth.txt,1))
x = t(x)#now rows are subjects, cols are genes
x =
On Sat, 2008-03-15 at 12:48 -0500, daniel jupiter wrote:
Hi all,
I apologize for what might be a silly question.
I am interested in doing a one way anova.
This is not too hard in and of itself, either with anova, aov or oneway.test
.
However, I need to
1) get pvalues,
if obj is the
Hi!
I think the results agree:
using simulated data:
set.seed(1)
library(mvtnorm)
x-rmvnorm(44,rep(0,10))
y = x[(26:45)-1,1:10]
x = x[(2:25)-1,1:10]
p = ncol(x); p
nx = nrow(x); nx
ny = nrow(y); ny
n = nx+ny; n
# (t(x)-colMeans(x)) %*% t(t(x)-colMeans(x))
T2 = nx*ny/n *
Hello,
is there already a function in any R package which does
source code formatting of deparsed lists?
Let's create the following list:
x - list(a = round(rnorm(3), 2),
b = round(rnorm(3), 2))
xx -c(aa = round(rnorm(30)), f = function(a) a + b, list(x, x))
Now, I want deparse it
Is this basically what you want?
xx -list(aa = round(rnorm(30),2), f = function(a) a + b, a=x, b=x)
xx
$aa
[1] 1.34 -0.21 -0.18 -0.10 0.71 -0.07 -0.04 -0.68 -0.32 0.06 -0.59
0.53 -1.52 0.31 -1.54 -0.30
[17] -0.53 -0.65 -0.06 -1.91 1.18 -1.66 -0.46 -1.12 -0.75 2.09 0.02
-1.29 -1.64
jim holtman wrote:
Is this basically what you want?
xx -list(aa = round(rnorm(30),2), f = function(a) a + b, a=x, b=x)
xx
$aa
[1] 1.34 -0.21 -0.18 -0.10 0.71 -0.07 -0.04 -0.68 -0.32 0.06 -0.59
0.53 -1.52 0.31 -1.54 -0.30
[17] -0.53 -0.65 -0.06 -1.91 1.18 -1.66 -0.46 -1.12 -0.75
Hi,
Is it possible to supply a new method for the %*% operator? I need to
provide a new method for working on variables of a newly defined class,
ad. I've had no problems overloading +, * etc.., using code such as:
+.ad - function(a,b = NULL)
{
# further code here
}
I've tried to do the
Dear useRs,
I have several problems in using rgl-0.77 (and recent earlier
versions) on Gentoo Linux with a custom-built v. 2.6.22 kernel.
Currently I use R-2.6.1.
When I build rgl,
# R CMD INSTALL /home/liviu/inst/dwn/R/rgl_0.77.tar.gz
or
install.packages(rgl, dependencies=TRUE, method =wget),
Joe Cainey jcainey at gmail.com writes:
Is it possible to supply a new method for the %*% operator?
clipped
I've tried to do the same thing with %*%:
%*%.ad - function(a,b)
{
# further code here
}
However this doesn't work; the new method is never called and the standard
%*%
On 16/03/2008 3:38 PM, Liviu Andronic wrote:
Dear useRs,
I have several problems in using rgl-0.77 (and recent earlier
versions) on Gentoo Linux with a custom-built v. 2.6.22 kernel.
Currently I use R-2.6.1.
When I build rgl,
# R CMD INSTALL /home/liviu/inst/dwn/R/rgl_0.77.tar.gz
or
Hi,
I am trying to use the Fisher scoring method with a geometric distribution,
with p = .07, 100 observations from the geom distrib, and 10 iterations.
I cannot quite get the code to work.
Can anyone see the mistake?
n - 100
p - 0.07
x - rgeom(n, p)
s - sum(x)
f - function(x, p)
On 16 March 2008 at 17:00, Duncan Murdoch wrote:
| On 16/03/2008 3:38 PM, Liviu Andronic wrote:
| Dear useRs,
|
| I have several problems in using rgl-0.77 (and recent earlier
| versions) on Gentoo Linux with a custom-built v. 2.6.22 kernel.
| Currently I use R-2.6.1.
|
| When I build
Thanks, Ken.
1. How can I find S4 methods for a given function given class(es)
of objects? The 'showMethods' function lists available generics for a
given function; showMethods('%*%') just produced for me a list of 52
different signatures for %*%. However, I don't know how to find
I have a problem I cannot get over for a long time. Imagine I have a
data frame with 17 colums. 16 of them are craniometric variables of
Cervus elaphus and one contains age.
The data frame has 83 rows.
I want to write a loop which plots the values of each craniometric
variable against the
Take a look at 'matplot'. If you want to loop, try
for (i in 1:16) plot(df[[paste(V, i, sep=)]] ~ df$AGE)
On 3/16/08, Ing. Michal Kneifl, Ph.D. [EMAIL PROTECTED] wrote:
I have a problem I cannot get over for a long time. Imagine I have a
data frame with 17 colums. 16 of them are craniometric
You can try this also:
sapply(DF[-match(AGE, names(DF))], plot, x=DF$AGE)
On 16/03/2008, Ing. Michal Kneifl, Ph.D. [EMAIL PROTECTED] wrote:
I have a problem I cannot get over for a long time. Imagine I have a
data frame with 17 colums. 16 of them are craniometric variables of
Cervus elaphus
Hi the list
I am fighting with the twins setAs and setIs...
Here are some questions and comments (comments to myself but that migth
be wrong, it is why I am posting them)
1. Very surprising : using setIs define 'is', 'as-' but not 'as' ???
2. Using setAs define 'as', 'as-' but not 'is'...
What
On Mar 16, 2008, at 8:12 PM, Christophe Genolini wrote:
Hi the list
I am fighting with the twins setAs and setIs...
Here are some questions and comments (comments to myself but that
migth
be wrong, it is why I am posting them)
1. Very surprising : using setIs define 'is', 'as-' but not
I am doing a reanalysis of results that have previously been published.
My hope was to demonstrate the value of adoption of more modern
regression methods in preference to the traditional approach of
univariate stratification. I have encountered a puzzle regarding
differences between I
Dear all,
I have a question regarding the use of stepAIC and polynomial (quadratic to be
specific) terms in a binary logistic regression model. I read in McCullagh and
Nelder, (1989, p 89) and as far as I remember from my statistics cources,
higher-degree polynomial effects should not be
Spencer Graves wrote:
Thanks, Ken.
1. How can I find S4 methods for a given function given class(es)
of objects? The 'showMethods' function lists available generics for a
given function; showMethods('%*%') just produced for me a list of 52
different signatures for %*%.
Dear Martin:
This is wonderful. Thank you very much.
It would be a great help if your suggestions could be added to
See Also and Examples for methods.
Thanks again,
Spencer Graves
Martin Morgan wrote:
Spencer Graves wrote:
Thanks, Ken.
1. How can I
At 08:50 PM 3/16/2008, caspar wrote:
Dear all,
I have a question regarding the use of stepAIC and polynomial
(quadratic to be specific) terms in a binary logistic regression
model. I read in McCullagh and Nelder, (1989, p 89) and as far as I
remember from my statistics cources, higher-degree
--- jim holtman [EMAIL PROTECTED] wrote:
Glad I could help. You might want to post it back
to R-Help so that
others can see what was done.
On Sun, Mar 16, 2008 at 6:02 PM, Felipe Carrillo
[EMAIL PROTECTED] wrote:
Jim: I owe you man, this is great,I never thought
that
I could
There is absolutely no reason to remove age altogether. Notice that
typically age and age^2 are highly correlated. To see this, consider 100
people with mean age 35 and 95% tolerance limite between 15 and 55:
age - rnorm(100, 35, 10)
cor(age, age^2)
[1] 0.9898186
So if you use raw age and
Hi r-users,
How do find the mean for each row? Thank you in advance for your help.
1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 Day Totals
10.0 0.0 0.0 0.0 0.6 0.0 8.4 0.0 29.4 0.0 38.4
20.0 0.0 1.8 0.0 22.4 0.0 0.2 0.4 0.8 0.0
How about rowSums(x)/ncol(x), where x is the matrix?
On Mon, Mar 17, 2008 at 1:48 PM, Roslina Zakaria [EMAIL PROTECTED] wrote:
Hi r-users,
How do find the mean for each row? Thank you in advance for your help.
1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 Day Totals
1
30 matches
Mail list logo