On 8/4/2005 1:03 PM, Shengzhe Wu wrote:
> Hello,
>
> When I open a new device by "windows()", how to hide the 3 buttons on
> the top-right corner of this window, since I want to make these
> buttons not work, and the window can be closed only by "dev.off()" or
> "graphics.off()".
There are no exp
On 8/4/2005 11:39 AM, Shengzhe Wu wrote:
> and how about R on windows system? since I am currently using R of
> windows version.
Hit "ESC" in the console. If that doesn't work, there's a bug. Report
it to the package maintainer if you're running package code, or here if
it's base code. Be sur
Shengzhe Wu wrote:
> Hello,
>
> When I open a new device by "windows()", how to hide the 3 buttons on
> the top-right corner of this window, since I want to make these
> buttons not work, and the window can be closed only by "dev.off()" or
> "graphics.off()".
You cannot.
Uwe Ligges
>
> Thank
and how about R on windows system? since I am currently using R of
windows version.
Thank you,
Shengzhe
On 8/4/05, Roger D. Peng <[EMAIL PROTECTED]> wrote:
> On Unix you can send one of the 'USR' signals to 'kill'. See
> ?Signals.
>
> -roger
>
> Shengzhe Wu wrote:
> > Hello,
> >
> > How to sto
On Unix you can send one of the 'USR' signals to 'kill'. See
?Signals.
-roger
Shengzhe Wu wrote:
> Hello,
>
> How to stop the process when there is a mishandling making R system
> frozen? If there is a way other than quitting the system when frozen
> occurs?
>
> Thank you,
> Shengzhe
>
> ___
On Wed, 3 Aug 2005, ecoinfo wrote:
> Dear Dr. Bivand and R-users,
> I have a 5 by 5 grid, say, location[1:5,1:5], and I want to know the
> indices of 8 neighbours of each cell. For example, for location[2,2], its
> neighbour coordinates are [1,1:3], [2,1], [2,3] and [3,1:3]. Sometimes I
> also
Rangesh Kunnavakkam wrote:
> I have a large vector of around 12597 elements and I wish to calculate
> p-value for each element using a formula of something like:
> p-value= 1- exp^(kexp^(-labda))
> I was wondering someone could give some ideas how to implement for each
>
if labda is the elements of the vector and you know what kexp is , you can
use apply
apply(your_vector, 1, function(x) 1- exp^(kexp^(-x)))
HTH
Jean
On Wed, 3 Aug 2005, Rangesh Kunnavakkam wrote:
> I have a large vector of around 12597 elements and I wish to calculate
> p-value for each eleme
Wensui Liu wrote:
> Jorge,
>
> My understanding about sas.get is that it actually runs SAS in the
> backend, output the SAS data, and then import it into R again.
He has SAS running locally.
>
> It might be better to convert SAS to csv format and read it into R.
For that, the Hmisc package's sa
Jorge,
My understanding about sas.get is that it actually runs SAS in the
backend, output the SAS data, and then import it into R again.
It might be better to convert SAS to csv format and read it into R.
On 8/2/05, Jorge Sirgo <[EMAIL PROTECTED]> wrote:
> I am running R 2.0.1 on Windows trying
On 7/28/2005 7:26 AM, Shengzhe Wu wrote:
> Hello,
>
> When I use "windows()" to open several windows, e.g. 4 windows, device
> 2, device 3, device 4 and device 5. Normally using "dev.set()" to
> specify the current active window, if there is a way to specify the
> active window by mouse click? Whe
Mark,
Please do not post the same question to both R-help and BioC-help
mailing lists because 1) there are many people who are on both lists and
2) people's replies will be archived in two different places making it
harder to others to search in future.
Please see the responses on BioC-help list.
wu sz writes:
> trainSet = as.data.frame(scale(trainSet, center = T, scale = T))
> trainSet.plsr = mvr(formula, ncomp = 14, data = trainSet, method =
> "kernelpls",
> CV = TRUE, validation = "LOO", model = TRUE, x = TRUE,
> y = TRUE)
[Two side notes here:
On 7/20/05, Ronaldo Reis-Jr. <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I try to put a legend in a xyplot graphic.
>
> xyplot(y~x|g,ylim=c(0,80),xlim=c(10,40),as.table=T,layout=c(2,3), ylab="Número
> de machos capturados",xlab=expression(paste("Temperatura (",degree,"C)")),
> key=list(corner=c(0,0),x=0
On Tue, 19 Jul 2005, wu sz wrote:
> Thanks a lot! But DF[c(4,6)] <- DF[c(6,4)] seems to just exchange
> the data, not the names. If exchanging the columns(both data and
> names) needs two steps?
>
> DF[c(4,6)] <- DF[c(6,4)]
> names(DF)[c(4,6)] <- names(DF)[c(6,4)]
Yes, it does. You did howeve
Thanks a lot! But DF[c(4,6)] <- DF[c(6,4)] seems to just exchange
the data, not the names. If exchanging the columns(both data and
names) needs two steps?
DF[c(4,6)] <- DF[c(6,4)]
names(DF)[c(4,6)] <- names(DF)[c(6,4)]
Shengzhe
__
R-help@stat.math.e
On Tue, 19 Jul 2005, Prof Brian Ripley wrote:
> Just change the names! E.g.
>
> names(DF)[c(4,6)] <- names(DF)[c(6,4)]
>
> Strictly a data frame has names, not column names, hence the use the
> names() and names<-() functions here.
That answered the subject line. If you want to exchange the col
Just change the names! E.g.
names(DF)[c(4,6)] <- names(DF)[c(6,4)]
Strictly a data frame has names, not column names, hence the use the
names() and names<-() functions here.
On Tue, 19 Jul 2005, wu sz wrote:
> I have a data frame with 15 variables, and want to exchange the data
> of 4th colum
Marc Schwartz wrote:
> On Thu, 2005-07-14 at 14:38 +0200, wu sz wrote:
>
>>Hello,
>>
>>I wish to plot some figures in a window in turn, but the size of these
>>figures is different, so how can I change the size of the window by
>>resetting the parameters before each plotting?
>>
>>Thank you,
>>She
On Thu, 2005-07-14 at 14:38 +0200, wu sz wrote:
> Hello,
>
> I wish to plot some figures in a window in turn, but the size of these
> figures is different, so how can I change the size of the window by
> resetting the parameters before each plotting?
>
> Thank you,
> Shengzhe
Other than dragging
Not sure you can do this. You might have to launch separate
graphics windows.
-roger
wu sz wrote:
> Hello,
>
> I wish to plot some figures in a window in turn, but the size of these
> figures is different, so how can I change the size of the window by
> resetting the parameters before each plo
See ?which Hint: arr.ind=TRUE
Simon.
At 09:28 AM 14/07/2005, wu sz wrote:
>Hello,
>
>I have a data set matrix of 1200 * 15. How can I get the position of a
>specific value in the matrix?
>
>I use "seq(along = x)[x > value]" to look for the position of the
>value in the matrix, but "seq" can just
Use which(..., arr.ind=TRUE); e.g.,
> m <- matrix(runif(12), 3, 4)
> which(m > .8, arr.ind=TRUE)
row col
[1,] 1 3
[2,] 2 3
[3,] 3 3
[4,] 3 4
> m
[,1] [,2] [,3] [,4]
[1,] 0.2148183 0.08251853 0.9444718 0.4487148
[2,] 0.5386863 0.49673282 0.8054240 0.51
Hello,
a proposed solution of Bill Venables is archieved on the S-News mailing
list:
http://www.biostat.wustl.edu/archives/html/s-news/2001-07/msg00035.html
and if I remember it correctly (and if the variance matrix is estimated
from the data), another similar way is simply to use the Euclidean
You will find previous discussion about pairs in this list. There are limits on
what has been included within the functionality, but you can write your own
panel functions. Here's a starting point.
x <- runif(100)
dim(x) <- c(20,5)
panel.cor1 <- function(x,y,...){
mycor<- cor(x,y)
if (mycor
wu sz <[EMAIL PROTECTED]> writes:
> Hello,
>
> I use package "tcltk" to do some GUI programming, and want to find a
> function which can do the operation "click a button", just like using
> a mouse to click. If tkevent.generate can do that? I tried it as
> below, but failed. Please give me a hint
Dear Mulholland,
I'm sorry, I was not clearly and objective sufficiently.
Below you can see what I'm was trying to do:
# D2Mah by JCFaria and Gabor Grothiendieck: 10/7/05 20:46:41
D2Mah = function(y, x) {
stopifnot(is.data.frame(y), !missing(x))
stopifnot(dim(y)[1] != dim(x)[1])
y= as.
I think we now have a very good R function here.
Thanks for all Gabor!
JCFaria
Gabor Grothendieck wrote:
> This one adds the labels:
>
>
> D2Mah4 = function(y, x) {
>
> stopifnot(is.data.frame(y), !missing(x))
> stopifnot(dim(y)[1] != dim(x)[1])
> y= as.matrix(y)
> x= as.factor(x)
This one adds the labels:
D2Mah4 = function(y, x) {
stopifnot(is.data.frame(y), !missing(x))
stopifnot(dim(y)[1] != dim(x)[1])
y= as.matrix(y)
x= as.factor(x)
man = manova(y ~ x)
E= summary(man)$SS[2] #Matrix E
S= as.matrix(E$Residuals)/man$df.residual
InvS = solve(S)
m
Indeed, it is very nice Gabor (as always)!
So, a doubt: how to preserve the 'rowname' and 'colname' of D2, like in the
first function? I think it is useful to posterior analyzes (as cluster, for
example).
Regards,
# A small correction (reference to gtools was eliminated)
D2Mah2 = function(y, x
And here is one more simplification using the buildin mahalanobis
function:
D2Mah3 = function(y, x) {
stopifnot(is.data.frame(y), !missing(x))
stopifnot(dim(y)[1] != dim(x)[1])
y= as.matrix(y)
x= as.factor(x)
man = manova(y ~ x)
E= summary(man)$SS[2] #Matrix E
S= as.matrix
I think you could simplify this by replacing everything after the
nObjects = nrow(mds) line with just these two statements.
f <- function(a,b) mapply(function(a,b)
(mds[a,] - mds[b,])%*%InvS%*%(mds[a,] - mds[b,]), a,b)
D2 <- outer(seq(nObjects), seq(nObjects), f)
This also eliminates dep
Well, as I did not get a satisfactory reply to the original question I tried to
make a basic function that, I find, solve the question.
I think it is not the better function, but it is working.
So, perhaps it can be useful to other people.
#
# Calculate the matrix of Mahalanobis Distances betwe
In addition to Gabor's solution, you might be interested in the
combinations function from the gtools package.
library(gtools)
combinations(5,2)
[,1] [,2]
[1,]12
[2,]13
[3,]14
[4,]15
[5,]23
[6,]24
[7,]25
[8,]34
[9,]
Gabor Grothendieck wrote:
> On 7/9/05, Jose Claudio Faria <[EMAIL PROTECTED]> wrote:
>
>>Dear R users,
>>
>>The solution is probably simple but I need someone to point me to it.
>>How can I to generate a matrix from a numeric sequence of 1:10 like 'A' or 'B'
>>below:
>>
>>A)
>>|---
On 7/9/05, Jose Claudio Faria <[EMAIL PROTECTED]> wrote:
> Dear R users,
>
> The solution is probably simple but I need someone to point me to it.
> How can I to generate a matrix from a numeric sequence of 1:10 like 'A' or 'B'
> below:
>
> A)
> ||
> | 1 2 3 4 5 |
> |
Christian Hennig wrote:
> Dear Jose,
>
> normal mixture clustering (mclust) operates on points times variables data
> and not on a distance matrix. Therefore
> it doesn't make sense to compute Mahalanobis distances before using
> mclust.
> Furthermore, cluster analysis based on distance matrices
What class is "newx" when you get the error message? Is it a vector
or a 1xK array? If the former, force it to be an array. (Hint:
"array(1:4, dim=c(2,2))[1,]" is a vector of length 2, while "array(1:4,
dim=c(2,2))[1,,drop=FALSE]" is a 1 x 2 matrix..(
If this does not sol
Dear Jose,
normal mixture clustering (mclust) operates on points times variables data
and not on a distance matrix. Therefore
it doesn't make sense to compute Mahalanobis distances before using
mclust.
Furthermore, cluster analysis based on distance matrices (hclust or pam,
say) operates on a poin
Would this suffice:
n <- 10
NaCon <- matrix(rnorm(n * n, mean=10.77, sd=0.02), 10, 10)
MeanNaCon <- rowMeans(NaCon)
[...]
??
Andy
> From: John Sorkin
>
> I am having trouble with apply. Could someone suggest changes to the
> code below that will eliminate the for loops?
> R 2.1.1 patched
> Win
, June 30, 2005 12:04 PM
To: Ghosh, Sandeep
Cc: Sundar Dorai-Raj; r-help@stat.math.ethz.ch
Subject: Re: [R] Help with stripplot
On Thursday 30 June 2005 11:40, Ghosh, Sandeep wrote:
> Thanks for the prompt reply. I need some more help. I want the datapoints
> for each group to show u
On Thursday 30 June 2005 11:40, Ghosh, Sandeep wrote:
> Thanks for the prompt reply. I need some more help. I want the datapoints
> for each group to show up in a different color, like all wt data points to
> be green, het to be blue and hom to be red. Also I need the median point of
> all the data
>
>
> Ghosh, Sandeep wrote:
>> For the following code is there a way to make the jitter all line up
>> horizontally, instead of them being just randomly spread around a value.
>> So for ex if there are multiple values at 63 for genotype wt then all
>> the values should be plotted on the same y valu
Ghosh, Sandeep wrote:
> For the following code is there a way to make the jitter all line up
> horizontally, instead of them being just randomly spread around a value. So
> for ex if there are multiple values at 63 for genotype wt then all the values
> should be plotted on the same y value of
The general location model is assumed: continuous variables are assumed
normal conditional on the categorical variables. ft.
--
Fernando TUSELLe-mail:
Departamento de Econometría y Estadística [EMAIL PROTECTED]
Facultad de CC.EE. y Empresariales Te
On 6/24/05, Domenico Cozzetto <[EMAIL PROTECTED]> wrote:
> Dear all,
> I have a matrix of dissimilarities and I'd like to get a 2d embedding. I'm
> not an expert of these techniques, so I'm a bit confused... Furthermore I
> was not able to find on the web any satisfactory tutorial. So even though
>
Hi Domenico,
multidimensional scaling methods (functions cmdscale, isoMDS, sammon)
start with a dissimilarity matrix while PCA starts from a data
matrix with variables. Equivalence of PCA and classical MDS refers to
the situation where you use the Euclidean distances between the data
points on whi
On 6/24/05, Domenico Cozzetto <[EMAIL PROTECTED]> wrote:
> Furthermore I
> was not able to find on the web any satisfactory tutorial. So even though
> this may be not the most appropriate place to discuss about this issues, I'd
> be very grateful to those who will reply.
http://ordination.okstate.
Hi,
Try allt.som$visual
Edgar Acuna
UPRM
On Thu, 9 Jun 2005, Ken Termiso wrote:
> Hi all,
>
> I originally posted this to the bioconductor group, but maybe it's better
> suited to the r-help...
>
> I'm using som() to partition samples of gene expression data into clusters.
> The point is to class
This seems to have more to do with NULLs than NAs. For instance:
> a <- data.frame(matrix(1:8, nrow=2))
> a
X1 X2 X3 X4
1 1 3 5 7
2 2 4 6 8
> a[a$X2 == 4,]$X1 <- NULL
> a
X1 X2 X3 X4
1 1 3 5 7
2 4 6 8 4
James
On 8/06/2005 7:15 a.m., Liaw, Andy wrote:
> There's something peculi
There's something peculiar that I do not understand here. However, did you
realize that the thing you are assigning into parts of `a' is NULL? Check
you're my.test.boot.ci.1: It's NULL.
Be that as it may, I get:
> a <- data.frame(matrix(1:4, nrow=2), X3=NA, X4=NA)
> a
X1 X2 X3 X4
1 1 3 NA
"Liaw, Andy" <[EMAIL PROTECTED]> writes:
> > R : Copyright 1999, The R Development Core Team
> > Version 0.64.1 (May 8, 1999)
>
> For starter, try upgrading R to a version from this century.
Or millennium for that matter
> The version
> you have seems very strange: AFAICR by 1999 R ve
> From: luciana
>
> Hello!
> I need help to use the package pls.pcr in R.
>
> I installed R in an IRIX 6.5, using the version of R 0.64.1 from
> sgifreeware(I didn't get to install the newest version using
> make). I
> need to use the package pls.pcr and when I give the command:
>
> # R
Dear all,
i have a question as regards to the read.table /table command:
I read in the following file:
Method1 23
Method2 12
Method3 43
Method4 76
Method1 2
Method3 4
If i use the "table" command I get the following matrix:
V2
V1 2 4 12 23 43 76
Method1 1 0 0 1 0 0
Metho
On 6/3/2005 12:17 PM, Marco Zucchelli wrote:
Hi,
I am trying to open a file with guidlgOpen. on windows. When the file is in a directory whose name contains a space, the name is splitted into 2, i.e.
Afile <- guiDlgOpen(title= " Open file",defaultFile="",defaultDir="", multi=FALSE, filters =
It
is
really
hard
to
read
your
message.
Anyway,
I
guess
you
are
looking
for
?reshape
Uwe Ligges
jose silva wrote:
Dear all:
I have this:
A1 B1 C1 D1 E1
A2 B2 C2 D2 E2
A3 B3 C3 D3 E3
And I want this
A1 E1
This is a VERY special question. Please redirect your question to the
appropriate mailing list. See
http://mailman.csd.univie.ac.at/mailman/listinfo/rcom-l
for more information.
Uwe Ligges
Constant Depièreux wrote:
Dear All,
I am experiencing a problem for which I need some help.
I have cr
On 5/19/05, Ted Harding <[EMAIL PROTECTED]> wrote:
> Hi Gabor,
>
> On 05-May-05 Gabor Grothendieck wrote:
> > On 5/5/05, Ted Harding <[EMAIL PROTECTED]> wrote:
> >> [...]
> >> However, while representing the raw data in such a form is
> >> well supported by R, it seems to me that extracting data
>
Hi Gabor,
On 05-May-05 Gabor Grothendieck wrote:
> On 5/5/05, Ted Harding <[EMAIL PROTECTED]> wrote:
>> [...]
>> However, while representing the raw data in such a form is
>> well supported by R, it seems to me that extracting data
>> in a way adapted to different analyses requires users to
>> cre
The scope of variables within show.a is not affected by the environment
that show.a is called in. The scope of the variables in show.a
is determined by the lexical scope of show.a although you can change
this via:
environment(show.a) <- my.env
show.a() # 200
If you want to create a function
On 5/13/05, Matthieu Cornec <[EMAIL PROTECTED]> wrote:
> Hello,
>
> Does anyone know how to write the files created by the call of
> "texi2dvi" in another directory ?
Here are two possibilities:
1. Copy the input files into whatever directory you want using
file.copy and then setwd to that direc
On Tue, 10 May 2005, Stephen Robert Welsh wrote:
Is there a way to print off a hard copy of your plots when using R on
a Unix system. I can't seem to find any commands in any of the
documentation.
?dev.print.
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statis
Giuseppe Biondi Zoccai wrote:
I would like to perform a bootstrap validation of a backward stepwise
logistic regression analysis, but I am a beginner with R and I am not
sure of how to do it.
Is there anyone that can send me a sample file in tab format (that I
can modify in Excel by pasting my data
On 5/5/05, Ted Harding <[EMAIL PROTECTED]> wrote:
> On 05-May-05 Peter Dalgaard wrote:
> > [...]
> > Both systems are victims of the curse of the rectangular data set to
> > some extent. Prototypically, you record the sex of a rat along with
> > every single measurement on it, as if the rat could c
On Thu, 5 May 2005, Christoph Lehmann wrote:
I heard that 'R' does not do a very good job at handling large
datasets, is this true?
importing huge datasets in a data.frame with e.g. a subsequent step of
conversion of some columns into factors may lead into memory troubles
(probably due to memor
I
heard that 'R' does not do a very good job at handling large datasets, is
this true?
importing huge datasets in a data.frame with e.g. a subsequent step of
conversion of some columns into factors may lead into memory troubles
(probably due to memory overhead when building out factors). But we
"Angus Repper" <[EMAIL PROTECTED]> writes:
> Hello
>
>
>
> I am a long-time SAS user, but am new to R. I specifically am looking for
> information pertaining to generating graphics for web output. I would like
> to create dynamic graphics (in the form of generalized reports) for my web
> site
Ingmar Visser wrote:
Hi Uwe,
Thanks for that answer. However, one thing remains unclear, which latex file
do I use here, ie the one generated from Rdconv, or the one generated as an
intermediate step in Rd2dvi?
I meant the one generated from Rdconv.
Uwe
In the latter case, wouldn't I get a double
Hi Uwe,
Thanks for that answer. However, one thing remains unclear, which latex file
do I use here, ie the one generated from Rdconv, or the one generated as an
intermediate step in Rd2dvi?
In the latter case, wouldn't I get a double \documentclass statement and the
like?
thanks in advance, ingmar
Ingmar Visser wrote:
Hi all,
I'm writing a vignette for my package, and I would like to include some of
the package help files in there as well. Is there an easy way of doing so?
I tried using R CMD Rdconv to generate latex files from .Rd files but I am
not sure how to include these into a .Rnw fil
"Ghosh, Sandeep" <[EMAIL PROTECTED]> writes:
> For the r script below
>
> >datestr <- "01/01/2004"
> >as.POSIXct(as.Date(datestr, "%d/%m/%Y"))
> I get the following output
> "2003-12-31 18:00:00 Central Standard Time"
>
> Why is the date a day before. I guess its something to do with the
> tim
Hi
Deepayan Sarkar wrote:
On Friday 22 April 2005 03:29, Sander Oom wrote:
Dear R users,
If I manage to sort out this graph, it is certainly a candidate for the
new R graph gallery
(http://addictedtor.free.fr/graphiques/displayGallery.php)!
I created the following lattice graph:
library(lattice)
tm
On Fri, 22 Apr 2005, Ghosh, Sandeep wrote:
For the r script below
datestr <- "01/01/2004"
as.POSIXct(as.Date(datestr, "%d/%m/%Y"))
I get the following output
"2003-12-31 18:00:00 Central Standard Time"
Why is the date a day before. I guess its something to do with the time,
but is there a way to g
On Friday 22 April 2005 03:29, Sander Oom wrote:
> Dear R users,
>
> If I manage to sort out this graph, it is certainly a candidate for the
> new R graph gallery
> (http://addictedtor.free.fr/graphiques/displayGallery.php)!
>
> I created the following lattice graph:
>
> library(lattice)
> tmp <- e
Well done !
The result of your work has been uploaded to the gallery (graph 48).
Romain
Le 22.04.2005 22:39, Paul Murrell a écrit :
Hi
Deepayan Sarkar wrote:
On Friday 22 April 2005 03:29, Sander Oom wrote:
Dear R users,
If I manage to sort out this graph, it is certainly a candidate for the
new R
7310
Thanks
Mike
- Original Message -
From: "Ted Harding" <[EMAIL PROTECTED]>
To: "Mike White" <[EMAIL PROTECTED]>
Cc:
Sent: Wednesday, April 20, 2005 8:58 AM
Subject: RE: [R] Help with predict.lm
> Sorry, I was doing this too late last night!
> All st
On 20-Apr-05 Mike White wrote:
> Ted
> Thank you for giving me so much of your time to provide
> an explanation of the likelihood ratio approach to the
> calibration problem. It has certainly provided me with a
> new insight into this problem. I will check out the
> references you mentioned to get
Sorry, I was doing this too late last night!
All stands as before except for the calculation at the end
which is now corrected as follows:
On 19-Apr-05 Ted Harding wrote:
[code repeated for reference]
> The following function implements the above expressions.
> It is a very crude approach to solvi
On 19-Apr-05 Ted Harding wrote:
> On 19-Apr-05 Mike White wrote:
>> Hi
>> I have measured the UV absorbance (abs) of 10 solutions
>> of a substance at known concentrations (conc) and have
>> used a linear model to plot a calibration graph with
>> confidence limits. I now want to predict the concen
root
results<-list(Predicted=xpre, Confidence=conf) ## returns results as list
to avoid warning message
return(results)
}
Thanks
Mike White
- Original Message -
From: "Ted Harding" <[EMAIL PROTECTED]>
To: "Mike White" <[EMAIL PROTECTED]>
Cc:
Sent: Tuesda
On 19-Apr-05 Mike White wrote:
> Hi
> I have measured the UV absorbance (abs) of 10 solutions
> of a substance at known concentrations (conc) and have
> used a linear model to plot a calibration graph with
> confidence limits. I now want to predict the concentration
> of solutions with UV absorban
w.abs) # gives slightly different results than
the calib function
Thanks to all who replied,
Mike White
- Original Message -
From: "Liaw, Andy" <[EMAIL PROTECTED]>
To: "'Mike White'" <[EMAIL PROTECTED]>;
Sent: Tuesday, April 19, 2005 1:53 PM
Subject:
> Hi
> I have measured the UV absorbance (abs) of 10 solutions of a
> substance at known concentrations (conc) and have used a
> linear model to plot a calibration graph with confidence
> limits. I now want to predict the concentration of solutions
> with UV absorbance results given in the ne
Will this help?
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/1202.html
[Found by RSiteSearch("calibration") in R-2.1.0.]
Andy
> From: Mike White
>
> Hi
> I have measured the UV absorbance (abs) of 10 solutions of a
> substance at
> known concentrations (conc) and have used a linear model to
On Sun, 2005-04-17 at 02:38 +0800, wenqing li wrote:
> Hey, all: Do we have a convenient command(s) to extract the variance
> components from a fitted model by lm (actually it's a nexted model)?
>
> e.g.: using the following codes we could get MSA,MSB(A) and MSE. How
> to get the variance componen
Briggs, Meredith M <[EMAIL PROTECTED]> wrote on 2005-Apr-15:
>
> Hello
Hi,
> How do I use function 'MERGE" to combine the FILE A and FILE B below to make
> FILE C?
>
> Thank you
>
>
>
> FILE A
> 140151167
> 30.1 11.4 40
>
> dat1 <- data.frame(var1=c(140, 151, 167), var2=c(30.1, 11.4, 40))
> dat2 <- data.frame(var1=c(140, 167), var3=c(5.7, 30.3))
> merge(dat1, dat2, all=TRUE)
var1 var2 var3
1 140 30.1 5.7
2 151 11.4 NA
3 167 40.0 30.3
Matt Austin
Statistician
Amgen
One Amgen Center Drive
M/S 24-2-C
Thousa
Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of
> michael watson (IAH-C)
> Sent: Wednesday, April 06, 2005 7:18 AM
> To: John Fox
> Cc: r-help
> Subject: RE: [R] Help with three-way anova
>
> Hi John
>
> Thanks for your help, that was a ve
anks
Mick
-Original Message-
From: John Fox [mailto:[EMAIL PROTECTED]
Sent: 06 April 2005 12:52
To: michael watson (IAH-C)
Cc: 'r-help'; [EMAIL PROTECTED]
Subject: RE: [R] Help with three-way anova
Dear Mick,
For a three-way ANOVA, the difference between aov() and lm() is mostly
in th
chael watson (IAH-C)
> Sent: Wednesday, April 06, 2005 4:31 AM
> To: [EMAIL PROTECTED]
> Cc: r-help
> Subject: RE: [R] Help with three-way anova
>
> OK, now I am lost.
>
> I went from using aov(), which I fully understand, to lm()
> which I probably don't. I di
rying to
figure out what contrasts matrix I need to use...
Many thanks for your help
Mick
-Original Message-
From: Federico Calboli [mailto:[EMAIL PROTECTED]
Sent: 06 April 2005 10:15
To: michael watson (IAH-C)
Cc: r-help
Subject: RE: [R] Help with three-way anova
On Wed, 2005-04-06
On Wed, 2005-04-06 at 09:11 +0100, michael watson (IAH-C) wrote:
> OK, so I tried using lm() instead of aov() and they give similar
> results:
>
> My.aov <- aov(IL.4 ~ Infected + Vaccinated + Lesions, data)
> My.lm <- lm(IL.4 ~ Infected + Vaccinated + Lesions, data)
Incidentally, if you want
ferent to
zero? How do I interpret that?
Many thanks
Mick
-Original Message-
From: Federico Calboli [mailto:[EMAIL PROTECTED]
Sent: 05 April 2005 16:33
To: michael watson (IAH-C)
Cc: r-help
Subject: Re: [R] Help with three-way anova
On Tue, 2005-04-05 at 15:51 +0100, michael watson (IAH-C
Dear Michael,
For unbalanced data, you might want to take a look at the Anova()
function in the car package.
As well, it probably makes sense to read something about how linear
models are expressed in R. ?lm and ?formula both have some information
about model formulas; the Introduction to R manua
On Tue, 2005-04-05 at 15:51 +0100, michael watson (IAH-C) wrote:
> So, what I want to know is:
>
> 1) Given my unbalanced experimental design, is it valid to use aov?
I'd say no. Use lm() instead, save your analysis in an object and then
possibly use drop1() to check the analysis
> 2) Have I us
Agnes Gault wrote:
Hello!
I am an 'R beginner'. I am trying to check if my data follow a negative
binomial function.
the command i've typed in is:
> nbdo=rnegbin(58,mu=27.82759,theta=0.7349851)
> ks.test(do$DO,nbdo)
Each time i do that, p given is different
The p-values are different each time
Agnes Gault wrote:
Hello!
I am an 'R beginner'. I am trying to check if my data follow a negative
binomial function.
the command i've typed in is:
> nbdo=rnegbin(58,mu=27.82759,theta=0.7349851)
> ks.test(do$DO,nbdo)
Each time i do that, p given is different and i get this warning message:
'War
you can get the estimated covariance matrix of the random-effects
using:
library(nlme)
m <- lme(Orthodont)
# scaled by the residuals variance
pdMatrix(m$modelStruct$reStruct)
# raw
lapply(pdMatrix(m$modelStruct$reStruct), "*", m$sigma^2)
However, it seems that in your model you use as a fixed-eff
On Tuesday 29 March 2005 19:16, Jeff D. Hamann wrote:
> I'm trying to plot a grid over a levelplot
>
> print( levelplot( var1.pred~x+y,
> data=saw.pred,
> aspect=mapasp(saw.pred),
> col.regions=terrain.colors(80),
> main=main) )
>
I'm not sure I am answereing your question, but here goes
# Create a vector with 2400 items
x <- runif(2400)
# Create a 600 by 4 matrix
y <- matrix(x,ncol = 4)
#If you needed the matrix to be done row by row
#y <- matrix(x,ncol = 4,byrow = T)
# Extract values from a particular part of the matri
601 - 700 of 1084 matches
Mail list logo