summary.lm tells you F, df, and p value.
try:
summary(regression9)
On 2/6/07, Jason R. Finley [EMAIL PROTECTED] wrote:
Hello,
I have spent a good deal of time searching for an answer to this but
have come up empty-handed; I apologize if I missed something that is
common knowledge.
I am
Hi,
# Dimensions(in characters) of the internal pager.
pgrows = 48
pgcolumns = 128
Thanks, that worked.
Regards,
John.
-Original Message-
From: talepanda [mailto:[EMAIL PROTECTED]
Sent: 06 February 2007 06:37
To: Gavin, John
Cc: R-help@stat.math.ethz.ch
Subject: Re: [R]
Jason R. Finley wrote:
Hello,
I have spent a good deal of time searching for an answer to this but
have come up empty-handed; I apologize if I missed something that is
common knowledge.
I am trying to figure out how to get an ANOVA table that shows the
sum of squares. degrees of
Please discuss this on R-devel: it is way off-topic for R-help.
'11' is the number of arguments, it is not 'magic'.
On Tue, 6 Feb 2007, Geoff Russell wrote:
Hi All,
I've started work on a MetaPost device driver (please don't hold your
breath).
I've copied the XFig driver and renamed
On Tue, 6 Feb 2007, talepanda wrote:
For size, maybe:
# Dimensions(in characters) of the internal pager.
#pgrows = 25
#pgcolumns = 80
pgrows = 48
pgcolumns = 128
in Rconsole, but location cannot be handled.
For the very good reason that you can have multiple pagers and I at least
do
On Mon, 5 Feb 2007, Leeds, Mark (IED) wrote:
I had a couple of questions about the ar function that i was hoping
someone could answer.
[...]
I call the ar function with aic=TRUE as below so that it picks the order
of the ar model based on BIC.
It actually does as it says it does (BIC !=
For the very good reason that you can have multiple pagers and I at least
do not want them exactly on top of one another.
Yes. and I knew for the first time now that location of new pager is
randomly determined.
I think it is cool way.
On 2/6/07, Prof Brian Ripley [EMAIL PROTECTED] wrote:
On
Hi,
I have a set of covariance matrices but not the original data. I want to carry
out some exploratory factor analysis. So, I am trying to construct a covariance
matrix list as the input for factanal. I can construct a list which includes
the cov, the centers, and the n.obs. But it doesn't
I want to do moran's I test in R language. I try to use gearymoran in
Package ade4,moran in Package spdep, and Moran.I in Package ape. But
I do not know how to do it because data format is different.
My data:
xy dbh
111.0310.7 7
118.110.28 1.2
Hi Ansel,
thank you for the response
generally speaking, wavelets are known to be good at
extracting signal from noisy data and are adaptive but
I am not familiar with any R implementation of wavelets.
I used wavelets before.
1. wavelets is just a particular case of orthogonal function
Hello
I am currently translating some Matlab code into R and would like to know
what people would recommend as the equivalent in R for spline(x,y,xx) in
matlab. The latter performs a cubic spline interpolation where x,y are the
given data vectors and xx is the new abscissa vector.
It allows for
--- Stefan Grosse [EMAIL PROTECTED] wrote:
You are posting in the wrong mailing list, mail your
question to the lyx
users list: http://www.lyx.org/internet/mailing.php
OOPs I wondered why I wasn't getting a reply. I use
the same account for both lists. My appologies to the
list
John
Peter Dalgaard wrote:
Jason R. Finley wrote:
Hello,
I have spent a good deal of time searching for an answer to this but
have come up empty-handed; I apologize if I missed something that is
common knowledge.
I am trying to figure out how to get an ANOVA table that shows the
sum of
First, I use lattice data rather than point pattern data, so this may not be
applicable.
The following is an example of the Moran's I test using the south dataset,
analyzing the variable HR60 for dependence.
I use GeoDa to create my spatial weight (https://www.geoda.uiuc.edu/), and
I call the
Hi,
I've got some data in a data frame arranged like this:
x y z othervariables
0.1 0.2 1.7 0.01
0.2 0.2 1.3 0.23
0.2 0.3 1.1 0.43
etc
I'd like to plot a heatmap of this data, with x and y as the (x,y) co-
ords
Hi all,
for the data below I would like to
1. generate a dummy variable for each group gr of the same composition by
people, then save each portion in a separate file,
2. compute the frequency of 1's in x for each person by group
gr. So, mike will have freq=2/3, as he has two 1 and one 0 in
yyan liu zhliur at yahoo.com writes:
I used the function lme in R and Splus-7. With the same dataset and same
argument for the function, I got
quite different estimation results from these two software. Anyone has this
experience before?
From the FAQ:
R by default uses treatment contrasts
You can use
plot(y~x,col=color.index.in.palette.defined.from(z),pch=20,type=p)
where
color.index.in.palette.defined.from(z)
is a function or an expression, returning either a color index in a
predefined palette or any other color representation, suitable for R. This
is described in ?par.
I
The help page says
covmat: A covariance matrix, or a covariance list as returned by
'cov.wt'. Of course, correlation matrices are covariance
matrices.
and there is an example of a covariance list (ability.cov).
factanal(factors = 2, covmat = ability.cov)
Hi
I have table1 which has the foll. columns
id age rate
and table2 which has the foll. columns
id count
I wish to get data from table1 for all the ids which
are persent in table2 and where the rate is not equal
to 999.
The ids in table2 are a subset of those in table1 and
every id in table2 has
I installed R (R-2.4.1-win32.exe) on a new computer with Windows Vista
and a 64 bit operating system (hp dv9000 with intel core t7200). The
base R runs fine, but I can not get any of the packages to load. From
within R I choose install packages choose a site then a package. I
tried installing
Hi
I have table1 which has the foll. columns
id age rate
and table2 which has the foll. columns
id count
I wish to get data from table1 for all the ids which
are persent in table2 and where the rate is not equal
to 999.
The ids in table2 are a subset of those in table1 and
every id in table2 has
subset(table1, rate != 999id == table2$id)
On 2/6/07, lalitha viswanath [EMAIL PROTECTED] wrote:
Hi
I have table1 which has the foll. columns
id age rate
and table2 which has the foll. columns
id count
I wish to get data from table1 for all the ids which
are persent in table2 and where
Wensui Liu wrote:
subset(table1, rate != 999id == table2$id)
On 2/6/07, lalitha viswanath [EMAIL PROTECTED] wrote:
Hi
I have table1 which has the foll. columns
id age rate
and table2 which has the foll. columns
id count
I wish to get data from table1 for all the ids which
are
Hi,
I would like to make a conditional sum for certain columns in a dataframe.
This is how my dataframe looks like:
Clinic Rep ColM1 ColM2 ColM3 ... ColM40
A 1 10 01
B 1 0-1 00
C 1 0-1 1-1
A 2 11 -1
On 2/6/2007 10:33 AM, Daniel O'Shea wrote:
I installed R (R-2.4.1-win32.exe) on a new computer with Windows Vista
and a 64 bit operating system (hp dv9000 with intel core t7200). The
base R runs fine, but I can not get any of the packages to load. From
within R I choose install packages
Daniel O'Shea wrote:
I installed R (R-2.4.1-win32.exe) on a new computer with Windows Vista
and a 64 bit operating system (hp dv9000 with intel core t7200). The
base R runs fine, but I can not get any of the packages to load. From
within R I choose install packages choose a site then a
Yes, consider:
df=data.frame(Clinic=rep(c(A, B, C), 2),
Rep=c(1, 1, 1, 2, 2, 2),
colM1=c(1, 0, 0, 1, -1, 1),
ColM2=c(0, -1, -1, 1, 0, 0),
ColM3=c(0, 0, 1, -1, 0, 1),
ColM40=c(1, 0, -1, 0, 1, -1)
)
Clinic=1; Rep=2
Hi Seth
I'm running R on WinXP. On the first time I've installed Rdbi and RdbiPgSQL
from zip files not downloaded from Bioconductor. No problems during
installation.
Now I tried to install from biocLite and a message was showed for RdbiPgSQL.
biocLite(c(Rdbi, RdbiPgSQL))
Running getBioC version
The other day, CNN had a story on working at Google. Out of curiosity, I
went to the Google employment web site (I'm not looking, but just
curious). In perusing their job posts for statisticians, preference is
given to those who use R and python. Other languages, S-Plus and
something called SAS
On Tue, 6 Feb 2007, Jian Zhang wrote:
I want to do moran's I test in R language. I try to use gearymoran in
Package ade4,moran in Package spdep, and Moran.I in Package ape. But
I do not know how to do it because data format is different.
Please consider posting on the R-sig-geo mailing list.
Eduardo Dutra de Armas [EMAIL PROTECTED] writes:
Hi Seth
I'm running R on WinXP. On the first time I've installed Rdbi and RdbiPgSQL
from zip files not downloaded from Bioconductor. No problems during
installation.
Now I tried to install from biocLite and a message was showed for RdbiPgSQL.
Dear all,
Here is my questions:
1- Under a WINDOWS installation of R-2.4.1, can
we change the naming of a new ploting device open
by the command windows()?. Instead of the
default name e.g.: Device 2 I would like to use
something like Density plot or whatever!
2- Under a MAC OS X
I would like the option to specify alternative scale parameters when
using the gamma family, log link glm. In particular I would like the
option to specify any of the following:
1. maximum likelihood estimate
2. moment estimator/Pearson's
3. total deviance estimator
Is this easy? Possible?
Confirmation Of Results For Electronic Email Winners
MICROWORLD INTERNATIONAL PROMOTIONS
EUROPEAN CENTRAL OFFICE
162 ZUIDOOST 1109
AMSTERDARM THE NETHERLANDS
We happily announce to you the draw of the Microworld International
promotions organized for email users all over the World as part of
I think you mean 'shape parameter'. If so, see the MASS package and
?gamma.shape.
glm() _is_ providing you with the MLE of the scale parameter, but really
no estimate of the shape (although summary.glm makes use of one).
On Tue, 6 Feb 2007, WILLIE, JILL wrote:
I would like the option to
Hello,
I'm using lattice and opened an X11 device with the following call
X11(width=5,height=5,pointsize=1)
I then ran the following code
library(lattice)
x-rnorm(30,sd=2)
y-runif(30)
xyplot(y~x,pch=.,col=black,cex=1)
If i
On Tue, 6 Feb 2007, Prof Brian Ripley wrote:
I think you mean 'shape parameter'. If so, see the MASS package and
?gamma.shape.
Also http://www.stats.ox.ac.uk/pub/MASS4/#Complements
leads to several pages of discussion.
glm() _is_ providing you with the MLE of the scale parameter, but
How do you know the points are not being plotted?
All you can tell is that your graphics card / monitor / X11 driver is not
displaying them, not at all the same thing.
My guess is that the rectangle is being displayed only if the differences
in its discretized x and y coordinates are both
Thank you. You are correct, the shape parameter is what I need to
change I think I see how to use the MASS package to do it...or if not,
at least I have enough now to figure it out.
A question to reconcile terminology which will speed me up, if you have
time to help me a bit more: phi = 'scale
On Feb 6, 2007, at 4:05 PM, Prof Brian Ripley wrote:
How do you know the points are not being plotted?
All you can tell is that your graphics card / monitor / X11 driver
is not
displaying them, not at all the same thing.
Yes, that is what I meant but I suppose I conveyed it wrongly.
On Tue, 6 Feb 2007, WILLIE, JILL wrote:
Thank you. You are correct, the shape parameter is what I need to
change I think I see how to use the MASS package to do it...or if not,
at least I have enough now to figure it out.
A question to reconcile terminology which will speed me up, if you
As someone who has (reluctantly) sent job postings to R Help, I think
that a SIG would be a good idea.
Max
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Doran, Harold
Sent: Tuesday, February 06, 2007 2:08 PM
To: R-help@stat.math.ethz.ch
Subject: [R] R
Our (financial) firm has struggled immensely with finding appropriate
candidates and I believe this would be *tremendously* useful. I really hope
that your suggestion is implemented. Importantly, I believe it would help
create a commercial presence for R, and secondly reduce use of products like
... two main drawbacks of R at our firm (as viewed by our IT dept) are lack
of
guaranteed support as well as the difficulty in finding candidates.
-- Just an aside: lack of guaranteed support -- absolutely true in theory,
absolutely false in practice. I doubt that the voluntary support found on
I have an object, ARestopt, that is the result of a call to the arima
function. I want to calculate test statistics ( null is coeffs are zero
) for each of the estimated coefficients.
My coefficents are in a vector called ARestopt$coef and my covariance
matrix is ARestopt$var.coef.
I thought
On 2/6/07, Serguei Kaniovski [EMAIL PROTECTED] wrote:
Hi all,
for the data below I would like to
1. generate a dummy variable for each group gr of the same composition by
people, then save each portion in a separate file,
This gives back a list that you can go through and write out each
I wanted to print the first and last rows of some dataframes in Sweave
using dots in columns to separate the two parts. Head and tail almost
work, but I have problems with factors and row names.
z-data.frame(id=letters[1:26], x=sample(1:26,26))
rbind(head(z,3), ., tail(z,1))
id x
1
That sounds like a good idea. The name R makes it especially hard to
find job postings, resumes or do any other type of search. Googling
resume+sas or job opening+sas is quick and fairly effective (less a
few airline jobs). Doing that with R is of course futile. At the risk of
getting flamed, it's
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce
our*** R/S Advanced Programming ***course in San Francisco on March
15-16 and
*** Data Mining: Practical Tools and Techniques in R/Splus *** course in
Salt Lake City on March 26-27
Ask for group discount and reserve your
Thanks for that Brian,
I have worked through the examples. They work because the covmat were produced
by the cov.wt which provides output as a list object. I am trying to construct
my own list object to use as the covmat. There are no obvious instructions on
how to do this.
So, here is what I
Chris,
Not precisely what you want (does both rows and columns), but can be
specialized to your case easily.
Use with xtable in a chunk with 'results=tex' as in
xtable( dot.matrix( diag(10) ) )
Chuck
dot.matrix - function(x , rws=NULL, cls=NULL, in.math=FALSE )
{
## Purpose:
Given two columns of type character in a dataframe of the form:
col1col2
31* 66
0 0*
102*66
71* 80
31 2*
66 31*
47 38*
How do I generate the following dataframe? Ie. col1 contains row item
with * and col2 contains row member without *
col1col2
I've been looking for job that allows me to use R/S+ since I got out
of graduate school 2 years ago but with no success. I am wondering if
there is something that can be done to promote the use of R in
industry.
It's been very frustrating to see people doing statistics using
excel/spss and even
Alistair,
I have worked through the examples. They work because the covmat
were produced by the cov.wt which provides output as a list object.
I am trying to construct my own list object to use as the covmat.
There are no obvious instructions on how to do this.
So, here is what I have done
I think several ways can do that.
my code is:
data.frame(t(apply(dat,1,function(x)as.numeric(sub(\\*,,x[sort(1:2,grep(\\*,x)==2)])
HTH
On 2/7/07, Dale Steele [EMAIL PROTECTED] wrote:
Given two columns of type character in a dataframe of the form:
col1col2
31* 66
0 0*
I'm sorry, I forgot my definition of dat.
dat in code of previous post is:
dat-data.frame(col1=c(31*,0,102*,71*,31,66,47),col2=c(66,0*,66,80,2*,31*,38*))
dat
col1 col2
1 31* 66
20 0*
3 102* 66
4 71* 80
5 31 2*
6 66 31*
7 47 38*
On 2/7/07, talepanda [EMAIL PROTECTED]
Hi again Jasjeet, Walter,
I have a further question about an error message I get when running
multinomRob. I am simulating a dataset where I look at the effect of
making a previous categorical choice on the probability of making the
same choice later on. Given the following code:
n - 20
Try this:
DF - data.frame(
col1 = factor(c(31*, 0, 102*, 71*, 31, 66, 47)),
col2 = factor(c(66, 0*, 66, 80, 2*, 31*, 38))
)
replace(DF, TRUE, as.numeric(sub(*, , as.matrix(DF), fixed = TRUE)))
On 2/6/07, Dale Steele [EMAIL PROTECTED] wrote:
Given two columns of type character
Bob,
Far from flaming you, I think you made a good point - one that I
imagine most people who use R have come across. The name R is a big
impediment to effective online searches. As a check, I entered R
software, SAS software, SPSS software, and S+ software into
google. The R 'hit rate' was only
Roger,
Error in if (logliklambda loglik) bvec - blambda :
missing value where TRUE/FALSE needed
In addition: Warning message:
NaNs produced in: sqrt(sigma2GN)
That message comes from the Newton algorithm (defined in source file
multinomMLE.R). It would be better if we
Hi,
Viragra 1,80
Ciralis 3,00
Lervitra 3,35
http://officecop.printeryml-com
Replace - with . in the above link
--
would bring him here, beyond the reach of Dumbledores help and
protection, and into my waiting arms. And here he is ...the boy you all
believed had been my downfall...
Dear expeRts,
I'd like to picture histogram of ones.
For example,
hist(rep(1, 100), col=lightblue, border=black)
A bin is pictured ranging on [0,1] and range of x axis is also [0,1].
How can I adjust the bin range from [0,1] to [0.9, 1] ?
I am looking forward to hearing from expeRt.
If I correctly understand,
hist(rep(1, 100), col=lightblue, border=black,breaks=0:10*0.1)
see:
breaks arg in ?hist
On 2/7/07, Dong H. Oh [EMAIL PROTECTED] wrote:
Dear expeRts,
I'd like to picture histogram of ones.
For example,
hist(rep(1, 100), col=lightblue, border=black)
A bin is
I need to make weighted boxplots. I found that ggplot makes them. I
would however like to label them with the boxplot statistics (the
median, q1 and q3). In the boxplot function in r-base, I could output
boxplot statistics and then write a text on the plot to place the
labels. How would one do it
I know that two of the search engine for R is available:
http://www.rseek.org/
http://www.dangoldstein.com/search_r.html
--
Ýbrahim Mutlay
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
66 matches
Mail list logo