No. If B is singular, it's impossible to find a matrix A such that A%*%B
be the identity matrix (unless you can find a number x such that x*0=1).
Cheers,
Jerome
On August 14, 2003 10:02 am, Feng Zhang wrote:
Thank, Jerome
The question is if this generalized inverse can make
their product
see ?par and ?plot
plot(data[1,], xlab=Xlab, ylab=Ylab)
thomas
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Jim Lemon replied directly to me with a helpful function that needs to
be adapted. Yet it almost does exactly what I want. I'm posting it so
that it appears in the thread as it could be useful to others.
Thank you Jim
Tito
Le jeu 14/08/2003 à 15:28, Jim Lemon a écrit :
If I understand this,
Undoubtedly a simple question:
I've looked at order() and sort() in the help pages for
R1.7.1. It doesn't appear that these functions are immediately
suited to doing the same thing as
PROC SORT DATA = BLAH;
BY X Y Z;
RUN;
in SAS. I have also checked Frank Harrell's Hmisc library.
Does anyone know of some methods already programmed up ( ie freeware : ) ) in
R for the post hoc methods in anova. Particularly Scheffe's method or
Tukey's?
Anna
__
[EMAIL PROTECTED] mailing list
As a newbie to R, I need to learn my way around (no previous experience
of S).What books, doc for R are recommended? I'm interested primarly in
non linear regression and process modelling (and have downloaded the R
documentation from the site).
Thanks for all tips
Anne
Is there a way to print complex, multi-line text and numeric output
with
good control of the field width and precision of numeric values? In
the
first example in the R Documentation for the sprintf function, there
appears
to be a linefeed code \n, but it is printed literally.
sprintf(%s is %f
Uwe, you suggest
So you almost got it, just read ?shell carefully enough:
filenames - shell(dir D:\\tmp\\*.sasb7dat /b, intern = TRUE)
Unfortunately, this does not work and is why I wrote to Rhelp (OK I used
system instead of shell, but I tried shell first); here filenames is
assigned
I believe it is in the Rcmdr package, which requires the car library to be loaded.
You can also perform an ANOVA using the absolute value of the deviations from each
respective group mean, which is what Levene's Test does.
--
Harold C. Doran
Director of Research and Evaluation
New
Thanks. I made it. For those who have the same problem, here is my data:
corr2
12_at 13_at 14_at 15_at
12_at 1. -0.08452797 -0.2977039 -0.10522948
13_at -0.08452797 1. -0.2955658 0.02353673
14_at -0.29770388 -0.29556577 1.000
On Wed, 2003-08-06 at 16:30, Dirk Eddelbuettel wrote:
I am displaying several series in one plot, and would like to make
them distinct without having to employ a legend.
I managed to color tickmarks, but have been unsuccessful with either one of
a) making tickmarks thicker (without
1e100 is just one example of much bigger number that is exactly
represented
(in floating point). But of course
1e100+1 - 1e100
[1] 0
You mean the biggest number such that adding one changes the result?
I should be extremely careful with
print( 9007199254740994, digits=20)
[1] 9007199254740994
I suggested
sapply(1:length(lis), function (i) {v - lis[[i]]; v[which(v==next)+1]})
Of course that was really dumb. It can be simplified, because the index i
is only used to select a list element, which sapply() wants to do for me
anyway. It should be
sapply(lis,
Dear,
I have XYZ data available in a MySQL database. I get it out, can plot
the data with the plot() function, load it into a geoR datastructure.
But what I actually would like to do is a simple contouring of the data
based on a no Kriging interpolation such as TIN based.
I know the first
Dear Jens,
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
Sent: 13 August 2003 10:01
To: [EMAIL PROTECTED]
Subject: [R] placing labels in polygon center ?
Dear all,
is there any function to calculate the center of a polygon mass in R?
Actually I need to
Hi All,
I have a 8000*8000 matrix and I want to print out a file with the row name,
column name and the value for those point with values satisfying a condition.
I tried using a for loop, however, it took me forever to get the result. Is
there a fast way to do this? Thanks!
Bing
do.call(rbind, LIST)
does exactly
DF - rbind(LIST$X1, LIST$X2, LIST$X3, ..., LIST$XN)
See the help on do.call.
Not that rbind.data.frame is `very fast'. If your data frames are
all of exactly the same type (no coercions needed, factors with the same
set of levels, ...) it may be faster to
Is there a way to print complex, multi-line text and numeric output with
good control of the field width and precision of numeric values? In the
first example in the R Documentation for the sprintf function, there appears
to be a linefeed code \n, but it is printed literally.
sprintf(%s is %f
Hi,
I did answer the question, but didn't send it to the R mailing
list.
I'm working on it, any questions may be sent to me, I'm still working
on the functionallity right now... Examples will be included later.
Kurt.
--
Mail from Uwe Ligges
sent on Tuesday August 12 2003 at 16:52 (GMT+0200):
I've used Anova() from the car package to get marginal (aka type II)
sum-of-squares and tests for linear models with categorical
variables. Is it possible to get marginal SSs also for continuous
variables, when the model includes powers of the continuous variables?
For instance, if A and B are
On Wed, 2003-08-13 at 09:39, Haynes, Maurice (NIH/NICHD) wrote:
Has the Levene test of homogeneity of variance been implemented in any
library in R?
Thanks,
Maurice Haynes
National Institute of Child Health and Human Development
Child and Family Research Section
6705 Rockledge Drive
Hedderik van Rijn [EMAIL PROTECTED] writes:
Some more information about my read.spss issues: on a cleanly
installed Debian Woody machine:
~/TMP % R
R : Copyright 2003, The R Development Core Team
Version 1.7.1 (2003-06-16)
[...]
library(foreign); x - read.spss(dataDef.sav)
On Tue, 12 Aug 2003, [iso-8859-1] Bjørn-Helge Mevik wrote:
Prof Brian Ripley [EMAIL PROTECTED] writes:
drop1 is the part of R that does type II sum of squares, and it works in
your example. So does Anova in the current car:
I'm sorry, I should have included an example to clarify what I
On 12 Aug 2003, Peter Dalgaard BSA wrote:
Simon Fear [EMAIL PROTECTED] writes:
system(COMMAND.COM /c dir /b, intern=T, show.output.on.console=T)
character(0)
print(system(COMMAND.COM /c dir /b, intern=F,
show.output.on.console=T))
[1] 0
In both cases here a DOS window opens
(Ted Harding) [EMAIL PROTECTED] writes:
With a bit of experimentation I have determined (I think)
that on my R implementation the largest positive integer
that is exactly represented is (2^53 - 1), based on
(((2^53)-1)+1) - ((2^53)-1)
[1] 1
((2^53)+1) - (2^53)
[1] 0
Those integer
Hi!
I used the following commands
mydata-read.table(dataALLAMLtrain.txt, header=TRUE, sep
=\t,row.names=NULL)
It reads data without any error
Now if I use
edit(mydata)
It shows only 3916 entries, whereas the actual file contains 7129 entries)
[...]
So it seems R is truncating the
On Thu, 14 Aug 2003 09:47:49 +0200, you wrote:
Yes, I have. I am sorry if I am missing some very basic stuff, but both
order and sort will not only put the NAs at the end (with na.last=TRUE) BUT
also sort in ascending or descending order the rest of the elements and that
is not what I want.
Hi all,
I'm working with a dataset from 10 treatments, each
treatment with 30 subjects, each subject measured 5
times. The plot of the dataset suggests that a
3-parameter logistic could be a reasonable function to
describe the data. When I try to fit the model using
gnls I got the message 'Step
You can use the well-known Taylor series approximation to the
variance of an arbitrary function:
Var( f(X) ) ~= Sum( s[i]^2*D2[i] ) + 2*Sum( Sum( s[i,j]*D[i]*D[j] ) )
where D2[i] is the second partial derivative of f(x) with respect
to the ith parameter and D[j] is the first partial derivative
Hi.
How can i select interval of values from a variables?
I can ordinate my variable and use es.group1-varx[1:12]
if the limit values are linked to first and 12th
observation.
There is an easier way?
Thank you.
__
[EMAIL PROTECTED] mailing list
The answer is certainly not unique. Your email doesn't say
whether you are asking about principal components or simply
Gram-Schmidt orthogonalization.
- tom blackwell - u michigan medical school - ann arbor -
On Wed, 13 Aug 2003, Feng Zhang wrote:
I have a question about determining
Uwe Ligges wrote:
Angel wrote:
Why do
x-b%*%ginv(A)
Why should it be the same???
Here, you are calculating (A+ is G-Inverse)
x = b * A+
Let me add: More exactly, you are calculating
x = b' * A+
because b * A+ doesn't fit and R is somehow intelligent here and
transposes b for you instead
This seems to identify a possible bug in R 1.7.1 under Windows 2000:
tstDf - data.frame(y = 1:11, x=1:11)
fit - nls(y~a/x, data=tstDf, start=list(a=1))
predict(fit, se.fit=TRUE)
[1] 7.0601879 3.5300939 2.3533960 1.7650470 1.4120376 1.1766980 1.0085983
[8] 0.8825235 0.7844653 0.7060188
Hi
On 13 Aug 2003 at 12:25, Ronnen Levinson wrote:
Hello.
Is there an easy-to-use contour plot function analogous to
scatterplot3d that can draw handle a dataset of arbitrary (x,y,z)
triplets? That is, say x, y, and z are each measured quanties, and
exhibit neither order nor regularity.
Frank E Harrell Jr wrote:
sort(c(3,1,NA))
[1] 1 3
Shouldn't NAs be retained by default?
help(sort)
sort(c(3,1,NA),na.last=TRUE)
Cheers
Jason
--
Indigo Industrial Controls Ltd.
64-21-343-545
[EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
It is explained, with an example, on the help page. That and an example
in MASS is all the documentation there is, apart from the source code.
On Tue, 5 Aug 2003, Cecília Shiraiwa wrote:
Could someone help me to understand how to use the argument split on
summary.aov?
I would like to split
On Mon, 11 Aug 2003, Jorge A Ahumada wrote:
I usually startup R from within several directories (usually where the
programs or data of interest are). I like to use a small web browser
called dillo to browse help files. However, by default R looks for
mozilla, so I have to type every
I didn't study your code, but regarding a function to rotate a
vector: Multiplication by an orthogonal matrix does that. You may know
that an orthogonal matrix is a matrix whose transpose is its inverse.
Thus, A is orthogonal if and only if (A %*% t(A)) = identity. One of
the simplest
without seeing the file its hard to tell but one possibility that comes to
mind is that there is a # character in your
file. read.table considers this a comment character.
use the argurment comment.char= and see what happens...
On Wed, 13 Aug 2003, Dibakar Ray wrote:
I am very new to R. I
[EMAIL PROTECTED] writes:
Hello everybody,
could anybody give me a hint, who I can use rbind on a list of data.frames, please?
I have a list with a large number of data.frames of the same structure, like:
LIST - list(X1=data.frame(a=1,b=2), X2=data.frame(a=3,b=4),
Hello,
I try to use CGIwithR on Debian, apache web server. For one reason or
another, I get error messages such as (trying trivial.html/trivial.R),
that's when I look in the server log-files:
Mon Aug 11 16:07:17 2003] [error] [client xxx.xxx.xxx.xxx] Premature end of script
headers:
Why do
x-b%*%ginv(A)
and
x-solve(A,b)
give different results?. It seems that I am missing some basic feature of
matrix indexing.
e.g.:
A-matrix(c(0,-4,4,0),nrow=2,ncol=2)
b-c(-16,0)
x-b%*%ginv(A);x
x-solve(A,b);x
Thanks in advance,
Angel
__
[EMAIL
Regarding the accuracy of the Taylor series approximation, my favorite
reference is Bates Watts (1988) Nonlinear Regression Analysis and Its
Applications (Wiley, esp. pp. 255-260). Recently, Brian Ripley also
Chambers Hastie (1992) Statistical Models in S (Wadsworth, ch. 10) and
Venables
Hi
given an ARMA process and the AR and MA coefficients I need the residuals.
arima() calculates the residuals together with the best AR and MA
coefficients, but I need the coefficients to take known values.
In S-PLUS there is a function arima.filt(). Is there something similar in
R?
Thanks for
Perhaps you could specify the factor order you want using
the factor() function:
my.factor - factor(something, levels=c(A,C,B),
labels=c(A,C,B), ordered=TRUE)
boxplot(my.vector ~ my.factor)
Regards,
Andrew C. Ward
CAPE Centre
Department of Chemical Engineering
The
I was looking at the contrasts function in the search for planned
contrast functions. aov take a contrast variable but I don't see how
to use it. Is there a planned contrast capability in aov?
__
[EMAIL PROTECTED] mailing list
On 14-Aug-03 Feng Zhang wrote:
Thank, Jerome
The question is if this generalized inverse can make
their product to be identity matrix?
On August 14, 2003 09:24 am, Feng Zhang wrote:
Dear R-listers,
I have a dxr matrix Z, where d r.
And the product Z*Z' is a singular square matrix.
On 08/08/03 16:42, Al Piszcz wrote:
Is there a list of all contributed R libraries available through CRAN?
Ideally it would include a one or two line description.
I am looking for a packages() command similar to library()
but that would access the CRAN repository and provide a
listing of the
Barry Rowlingson [EMAIL PROTECTED] provided
functions PolygonArea and PolygonCenterOfMass.
As an exercise in R programming, I thought why don't I vectorise these
and then see if it makes a practical difference.
Here are my versions of his functions. Somehow I ended up with a sign
error when I
Simon Fear [EMAIL PROTECTED] writes:
Peter, thanks for quick reply, but I can't drop it:
system(dir /b, intern=T, show.output.on.console=T)
Error in system(dir /b, intern = T, show.output.on.console = T) :
dir not found
Oh, dir is a built-in, right?... Better wait for the real
Singular matrices are not invertible. However you can calculate the
generalized inverse with the function ginv() from package MASS.
HTH,
Jerome
On August 14, 2003 09:24 am, Feng Zhang wrote:
Dear R-listers,
I have a dxr matrix Z, where d r.
And the product Z*Z' is a singular square
On 13-Aug-03 Barry Rowlingson wrote:
Thompson, Trevor wrote:
I didn't see anything in the help file about . being some kind
of special character. Any idea why R is treating a decimal this
way in these functions? Any suggestions how to get around this?
'.' is the regexpr character for
Martin Maechler [EMAIL PROTECTED] writes:
jane == jane murray [EMAIL PROTECTED]
on Wed, 06 Aug 2003 09:16:48 +0100 writes:
jane Hi Can anyone help with the technique of obtaining
jane leverages from a conditional logistic regression
jane model? The code lm.influence does
Dear R list,
I have a sequence of weekly observations of number of adults and larvae
in various size classes from a butterfly population living in a
subtropical area with pronounced wet and dry seasons. Wet and dry
seasons are each defined 26 weeks long with fixed start and end dates.
The data
Have you read the rw-FAQ and set the amount of virtual memory available
to the R process via --memory-size?
[I doubt it! Please do read the FAQs.]
On Sun, 10 Aug 2003, gowuban wrote:
I have trouble with my cluster analysis using package cluster. diana
and agnes both seem to try to allocate
[EMAIL PROTECTED] wrote:
One possibility is to use a DBMS like MySQL or Postgresql, and RODBC to
connect to these. Search the archives for previous postings about these, have a
look at the first R-Newsletter and at Data Import-Export manual.
If you use PostgreSQL, you might want to try PL/R;
Simon Fear wrote:
How can I best achieve the following (works in Splus):
filenames - dos(dir *.sasb7dat /b)
What I am asking, more generically, is: how can I capture the output of
a
DOS command in R?
I have tried using
system(COMMAND.COM /c dir /b, intern=T, show.output.on.console=T)
where
On Tue, 2003-08-12 at 15:03, Gerhard Prade wrote:
Hello all,
i think i am to silly. I have installed R 1.7.1 (2003-06-16). Installed
some packages like xtables ore xml. I tried out this to installing
packages. Then i tried to make a crosstable like i know it from spss.
They say in this
Colin,
On Tue, Aug 12, 2003 at 10:05:46AM +0100, Colin Gillespie wrote:
I am trying to install R using debian's apt-get system. When I do 'apt-get install
r-base', it appears to install everything correctly. However, when I start up R I
get the error message:
That's how it should work.
Huan -
The difference between the empirical (bootstrap') result and the
theoretical results shows evidence for autocorrelation in the time
series data.
- tom blackwell - u michigan medical school - ann arbor -
On Wed, 6 Aug 2003 [EMAIL PROTECTED] wrote:
This is more a statistical
Dear Maurice,
Brian Ripley once posted a function for Levene's test to the R-help list. I
incorporated a slightly modified version in the Rcmdr package. Since the
function is very simple, I'll just list it here (it doesn't make sense to
use the Rcmdr package just for this function):
Chris Knight wrote:
I am tying myself in knots over subscripts when applied to lists
I have a list along the lines of:
lis-list(c(a,b,next,want1,c),c(d, next, want2, a))
From which I want to extract the values following next in each
member of the list, i.e. something along the lines of
Dear Paul,
Use order() to get the indices of the ordered rows.
See ?order
blah - data.frame(X = rep(3:1, each=4), Y = rep(c(2,1,2,1), 3), Z = rep(2:1, 6))
blah
X Y Z
1 3 2 2
2 3 1 1
3 3 2 2
4 3 1 1
5 2 2 2
6 2 1 1
7 2 2 2
8 2 1 1
9 1 2 2
10 1 1 1
11 1 2 2
12 1 1 1
blah -
I should clarify that `withdrawing support' does not mean that it will no
longer work with Win95/98/ME, just that we will no longer attempt to
ensure that it does.
On Tue, 12 Aug 2003, Simon Fear wrote:
Prof Brian Ripley writes:
If you insist on using a long-obselete OS, please help us
Greetings-
I tried to install 'XML' packages in R and get following error message:
install.packages(XML)
trying URL `http://cran.r-project.org/src/contrib/PACKAGES'
Content type `text/plain; charset=iso-8859-1'
Believe me, using Win98 is NOT my choice. I would only ever use
Unix/Linux.
Are there other statisticians similarly limited?
Like Simon I am very often limited with choice of software and OS when
working for clients. I've been waiting for about 5 months to get an editor
installed on my machine.
I am very new to R. I was trying to load some publicly available Expression
data in to R.
I used the following commands
mydata-read.table(dataALLAMLtrain.txt, header=TRUE, sep
=\t,row.names=NULL)
It reads data without any error
Now if I use
edit(mydata)
It shows only 3916 entries, whereas the
Hi,
Sorry if this is already known...
contour() sometimes draws contour lines that intersect.
Is there a temporary fix?
A dataset which causes problems is at
http://www.maths.uwa.edu.au/~adrian/dumpdata.R
If you try just
source(dumpdata.R)
image(huh)
contour(huh)
the 100 x 100
After installing RSPerl from Omegahat I tried to run the tests. The first
test is test.pl in /usr/lib/R/library/RSPerl/examples on Redhat Linux 8.
For some reason perl can't load R. Has anyone had this behavior? Any
suggestions?
Thanks,
--
Ari
http://damon.ib3.gmu.edu/~kahn
Dreams come true;
Hello
Is there a way to set an option in Rprofile to open a new graphic device
number (and save the old one)? Right now I am embedding a statement
such as if(dev.cur() == 3) get(getOption(device))(4) in my script each
time that I want a new graphic's device opened?
I notice there is an option
On Tue, 12 Aug 2003 20:07:35 -0700, you wrote:
Maybe we'll make = electric -- one press becomes -, two presses
become ==. Emacs is wonderful, sometimes.
Just watch out for named parameters:
plot(1:10,1:10, type - 'l')
probably isn't what you want!
Duncan
I have a matrix for which each row has 12 elements that represent the xyz coordinates
of 4 points.
So each row of M is (x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4). Some points have NA as z
values.
I want another matrix to be the same as M but with the coordinates of those points
with z=NA placed
help.search(cross validation) returns 'cv.glm(boot)'
and the boot package provides many other utilities for this.
- tom blackwell - u michigan medical school - ann arbor -
On Thu, 14 Aug 2003, array chip wrote:
Hi, is there a package for performing leave-one-out
cross validation in
Chris Knight [EMAIL PROTECTED] has
lis-list(c(a,b,next,want1,c),c(d, next, want2, a))
and wants c(want1,want2)
Step 1:
inx - sapply(lis, function(x) which(x == next)) + 1
== 4 3
Step 2:
sapply(1:length(lis), function(i) lis[[i]][inx[i]])
== want1 want2
Think
Simon Fear [EMAIL PROTECTED] writes:
Prof Brian Ripley writes:
If you insist on using a long-obselete OS, please help us continue ot
support it by supplying patches and workarounds. I think we (like
Microsoft) are going to have to think seriously about withdrawing
support
from
Try regexpr(\\., Female.Alabama)
-Original Message-
From: Thompson, Trevor [mailto:[EMAIL PROTECTED]
Sent: 13 August 2003 15:47
To: [EMAIL PROTECTED]
Subject: [R] Regexpr with .
I'm trying to use the regexpr function to locate the decimal in a
character string. Regardless of the
Yes, I have. I am sorry if I am missing some very basic stuff, but both
order and sort will not only put the NAs at the end (with na.last=TRUE) BUT
also sort in ascending or descending order the rest of the elements and that
is not what I want. And with order I would only get the z NAs at the end
I am trying to integrate OpenOffice Calc (OO) with R. One of the steps I am
taking is to create a very simple command line C program with the R shared
lib embedded in it. It would simply pass an R/S expression to R, print the
result, and catch errors. I know this is overkill, but it is a step
Simon Fear [EMAIL PROTECTED] writes:
The following code works, to gsub single quotes to double quotes:
line - gsub(', '', line)
(that's a single quote within doubles then a double within singles if
your
viewer's font is not good).
But The R Language Manual tells me that
Quotes and
Pinaud David [EMAIL PROTECTED] writes:
I'm importing some data (with read.table, file .txt) but I have some
problem: R truncates my numeric data at 2 digits during the
importation... It's the first time that I see that.
Are you sure that R is truncating the data itself or just truncating
the
I'm afraid that's not what I wanted. Let's see if taking your example data I
can explain better.
Starting from:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12]
[1,] 11.1 11.2 11.3 12.1 12.2 12.3 13.1 13.2 13.3 14.1 14.2 14.3
[2,] 21.1 21.2 NA 22.1 22.2 22.3 23.1 23.2
On 12 Aug 2003, Douglas Bates wrote:
nels.tomlinson.1 [EMAIL PROTECTED] writes:
I'd like to propose that now that the underscore-as-assignment-operator
is to be removed from R (good thing, too), that the Windows GUI should
replace the underscore ``_'' with the proper assignment operator
The rgl package offers 3-D rendering (and Gouraud shading) capabilities
through an OpenGL.
I don't know why it is not in cran, I found it at:
http://wsopuppenkiste.wiso.uni-goettingen.de/~dadler/rgl/
In any case that's not what I wanted (although I still have to explore in
more depth this
Hi
On 14 Aug 2003 at 9:47, Angel wrote:
Yes, I have. I am sorry if I am missing some very basic stuff, but
both order and sort will not only put the NAs at the end (with
na.last=TRUE) BUT also sort in ascending or descending order the rest
of the elements and that is not what I want. And
On Thu, 14 Aug 2003 12:43:25 -0400, Duncan Murdoch [EMAIL PROTECTED]
wrote :
Perhaps the description below of what se.fit is supposed to do should
be modified.
I've done that now in the development version (to become 1.8.0).
Err, I mean in the patch version (but it should still end up in
Kevin S. Van Horn [EMAIL PROTECTED] writes:
I'm just starting to learn to use R, and although I'm seeing lots of
functions aimed at doing orthodox statistical analyses, I don't see
the same for Bayesian analyses. What support does R have for Bayesian
statistics?
Bayesian = Orthodox,
On Tue, 12 Aug 2003, Bjørn-Helge Mevik wrote:
Prof Brian D Ripley [EMAIL PROTECTED] writes:
On Tue, 12 Aug 2003, [iso-8859-1] Bjørn-Helge Mevik wrote:
Why should I(x^2) be regarded as subservient to x?
In polynomial regression, it is usual to first consider a linear
model, then a
On Wednesday 13 August 2003 07:43, Brunschwig, Hadassa {PDMM~Basel} wrote:
Hello everyone
I have data grouped by subjects and i used the trellis plot to show the
curve for each subject. I have another vector a which gives me the
highest points the curves for each subjects can achieve ( the
Hello everybody,
could anybody give me a hint, who I can use rbind on a list of data.frames, please?
I have a list with a large number of data.frames of the same structure, like:
LIST- list(X1=data.frame(a=1,b=2), X2=data.frame(a=3,b=4),
X3=data.frame(a=5,b=6), , XN=data.frame(a=i,b=k))
On Thu, 7 Aug 2003, orkun wrote:
[quoting me without attribution]
Those are not predicted values, they are fitted values. Try predicting on
the same set of variables as you printed.
Precisely! From ?predict.glm
newdata: optionally, a new data frame from which to make the
Hello, all,
I'd like to propose that now that the underscore-as-assignment-operator
is to be removed from R (good thing, too), that the Windows GUI should
replace the underscore ``_'' with the proper assignment operator ``-''
when you type in the underscore character.
This is the current
I have uploaded to CRAN a new package named 'its' (Irregular Time-Series).
It
implements irregular time-series as an S4 class, extending the matrix class,
and records the time-stamp of each row in the matrix using POSIX. Print,
plot,
extraction, append, and related functionality are available.
Hi, Ana
try these source:
# Scheffe test, according Zar, J.H. (1984) Biostatistical Analysis.
# Englewood Cliffs: Prentice Hall.
# -
Scheffe1- function(mxa,mxb,ms,na,nb,k,F) {
# multiple comparision
# DMSa=| meanXa- meanXb |
# DMSc=
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Alle 03:25, venerdì 8 agosto 2003, [EMAIL PROTECTED] ha scritto:
If I run the same code as a normal user the following line fails with a
segmentation fault con -
dbConnect(m,host=acomputer,dbname=atable,user=au_ser,password=a_pass
word)
I receive
On Tue, 12 Aug 2003, Bing Zhang wrote:
Hi All,
I have a 8000*8000 matrix and I want to print out a file with the row name,
column name and the value for those point with values satisfying a condition.
I tried using a for loop, however, it took me forever to get the result. Is
there a
In addition to the previous replies, see the r-help thread on pretty
onscreen plots:
http://finzi.psych.upenn.edu/R/Rhelp02/archive/12254.html
HTH,
Jerome
On August 6, 2003 08:47 am, Yao, Minghua wrote:
All,
Anybody can tell how to export a R plot onto Word (or Power Point)?
Many thanks
Haynes, Maurice (NIH/NICHD) wrote:
Has the Levene test of homogeneity of variance been implemented in any
library in R?
levene.test() is in Rcmdr.
Uwe Ligges
Thanks,
Maurice Haynes
National Institute of Child Health and Human Development
Child and Family Research Section
6705 Rockledge Drive
On Thu, 14 Aug 2003 12:43:25 -0400, Duncan Murdoch [EMAIL PROTECTED]
wrote :
Perhaps the description below of what se.fit is supposed to do should
be modified.
I've done that now in the development version (to become 1.8.0).
Duncan Murdoch
__
[EMAIL
Hi Angel,
On 14-Aug-03 Angel wrote:
Yes, I have. I am sorry if I am missing some very basic stuff, but both
order and sort will not only put the NAs at the end (with na.last=TRUE)
BUT also sort in ascending or descending order the rest of the
elements and that is not what I want. And with
I am writing to announce the availability of a new package on CRAN:
boolean, which permits users to run Boolean logit and probit analyses.
BACKGROUND
The impact of independent probabilistic causal processes is often
thought to cumulate in a manner consistent with Boolean logic --
101 - 200 of 280 matches
Mail list logo