It looks like no one tested the beta of 2.1.0 for Windows in an East Asian
or cyrillic locale. In most cases example() will fail there.
If this is likely to affect you, please use R-patched instead.
We really do need help from the users in the beta-test period, especially
in languages not used
What does str(tsx) give?
So you are feeding in something the function has no idea about.
Try rsiTA([EMAIL PROTECTED],14)
Tom
> -Original Message-
> From: Neuro LeSuperHéros [mailto:[EMAIL PROTECTED]
> Sent: Wednesday, 20 April 2005 9:42 AM
> To: Mulholland, Tom; r-help@stat.math.ethz.ch
On Wed, 2005-04-20 at 02:08 +0200, Werner Wernersen wrote:
> Hi!
>
> I am trying to get the tick / label under a stacked
> boxplot with only a single
> data row. With >=2 rows it works, but with a single
> one the tick resp. my class
> name is not printed below the boxplot. Can anybody
> point m
Dimitris Rizopoulos wrote:
Hi Brendan,
according to
@Book{pinheiro.bates:00,
author= {J. Pinheiro and D. Bates},
title = {Mixed-Effects Models in S and S-PLUS},
year = {2000},
address = {New York},
publisher = {Springer-Verlag}
}
Section 2.2.8, the optimization procedure fo
Andrew Criswell asks:
Hello All:
Should I conclude from this discussion that there is no
practical
means by which nested generalized mixed models can be compared
from
output produced through glmmPQL or GLMM?
[WNV] The picture is, in my view, not as bleak as this,
Nope, its not that.
I only have one column in there. So
rsiTA(tsx[,2],14)
Error in "[.timeSeries"(tsx, , 2) : subscript out of bounds
I get the orginal error with:
rsiTA(tsx[,1],14)
Error in "[.timeSeries"(close, 1:(length(close) - 1)) :
only 0's may be mixed with negative subscripts
Thanks
Should you be using rsiTA(tsx[,2],14). If you look at the function you will see
it is expecting just the values you want in the calculation.
If you give it a matrix it treats the whole matrix as being price data.
Tom
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTE
Hello All:
Should I conclude from this discussion that there is no practical
means by which nested generalized mixed models can be compared from
output produced through glmmPQL or GLMM? What is one then to do???
Andrew
On Sun, 17 Apr 2005, Deepayan Sarkar wrote:
> On Sunday 17 April 2005 08:39
While you may want to jump in at the deep end I would suggest that you need to
take the time to fully digest some of the basics. This is one of them. There
are a number of ways that can be used. You tell me you have a table, but I
asume that you have a data.frame. I normally associate table with
fSeries Technical Analysis rsiTA problem
Hello,
Im trying to use the rsiTA() function but keep getting this error:
rsiTA(tsx,14)
Error in "[.timeSeries"(close, 1:(length(close) - 1)) :
only 0's may be mixed with negative subscripts
Heres is the first three lines of my data:
tsx[1:3,]
On Tuesday 19 April 2005 18:30, Glen Luckjiff wrote:
> On linux/X11 how do I stop paging of multiple page graphics.
> I've tried par(ask = T). This works in S+6 but not R.
> Thanks very much.
grid::grid.prompt(TRUE)
Deepayan
__
R-help@stat.math.ethz.c
Hi!
I am trying to get the tick / label under a stacked
boxplot with only a single
data row. With >=2 rows it works, but with a single
one the tick resp. my class
name is not printed below the boxplot. Can anybody
point me to what am I doing
wrong?
For example:
boxplot(data.frame(c(3,4,5)),nam
On linux/X11 how do I stop paging of multiple page graphics.
I've tried par(ask = T). This works in S+6 but not R.
Thanks very much.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Already we have seen several questioners confused by the new features of
2.1.0. Resources which may help:
- Read the R-admin manual. It has changed a lot, and covers setting up
locales and languages as well as changes to the configure options.
Everyone should read chapter 6 (Internationalizat
Thank you Prof. Ripley.
The installation is complete and R works great.
Best regards,
Cristian
Quoting Prof Brian Ripley <[EMAIL PROTECTED]>:
> On Tue, 19 Apr 2005 [EMAIL PROTECTED] wrote:
>
> > I am trying to install R-2.1.0 under a Sun platform.
> >
> > I used tar and then when I type
> >
In all likelihood, you'll get an overly optimistic estimate of performance
that way.
Andy
> From: array chip
>
> Hi all,
>
> In Tibshirani's PNAS paper about nearest shrunken
> centroid analysis of microarrays (PNAS vol 99:6567),
> they used cross validation to choose the amount of
> shrinkage
Hi all,
In Tibshirani's PNAS paper about nearest shrunken
centroid analysis of microarrays (PNAS vol 99:6567),
they used cross validation to choose the amount of
shrinkage used in the model, and then test the
performance of the model with the cross-validated
shrinkage in separate independent testi
Thanks you very much.
Yes, this is exactly what I want.
Best wishes,
Jan Sabee
On 4/19/05, Liaw, Andy <[EMAIL PROTECTED]> wrote:
> Is something like this what you're looking for?
>
> fixSome <- function(..., fixed) {
> fList <- list(...)
> for (i in fixed) fList[[i]] <- fList[[i]][1]
>
On Tue, 19 Apr 2005, Ivan Alves wrote:
Dear all,
I obtained and successfully compiled R 2.1.0 on Mac OSX 10.3.9. The
flags used for configuration are:
./configure --enable-utf8 --enable-R-shlib --with-blas='-framework
vecLib' --with-lapack --with-aqua
--with-tcl-config=/Library/Frameworks/Tcl.fram
On Tue, 19 Apr 2005, someone with the silly handle "Neuro LeSuperHéros"
wrote:
Hello, I must be doing something wrong that's very obvious. But I just don't
see it.
Nor do we. There are no such functions in R itself. Looks like this
might be package fBasics, although you did not mention it.
Disregard that last message. It seems to be a temperamental problem.
Sometimes it works sometimes it doesn't.
Jason Fisher
On 4/19/05, Jason Fisher <[EMAIL PROTECTED]> wrote:
> OK, I decided to try something new. Rather than installing RODBC via:
>
> install.packages("RODBC", .Library, repos=
Hi Chris,
You can get the desired effect by using the "at" argument to boxplot,
and by setting up the plot dimensions manually.
men <- data.frame(grp=rep(letters[c(1, 2, 4, 6)], each=10),
response=rnorm(40))
women <- data.frame(grp=rep(letters[c(2:5, 7)], each=10),
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
OK, I decided to try something new. Rather than installing RODBC via:
install.packages("RODBC", .Library, repos="http://cran.r-project.org";,
method="internal", destdir=getwd())
I instead downloaded the RODBC_1.1-3.zip file from
http://www.bioconductor.org/CRAN/ and installed th
Werner Bier <[EMAIL PROTECTED]> writes:
> Hi there,
>
> The optim(par, fn, ...parscale...) function in R requires 'parscale' which is
> defined as:
>
> "A vector of scaling values for the parameters. Optimisation is performed on
> 'par/parscale' and these should be comparable in the sense th
Works for me. Did you update RODBC with update.packages(checkBuilt=TRUE)?
The version I downloaded says in library(help=RODBC)
Built:R 2.1.0; i386-pc-mingw32; 2005-04-09 19:24:26;
windows
and for definiteness
pedump -e /R/library/RODBC/lib/RODBC.dll | grep Clo
v 2.0.1 (sh old!) on Win2k
I think I know the answer to this but I can hope ...
I have data for continuous variables (measures of residents) by a
categorical variable in range (1,22), the units in which they live.
I want to plot these data with a pair of boxplots one above another
with s
Jason Fisher wrote:
Hello...
After installing the precompiled version of R 2.1.0 (congratulations
to the R Development Core Team) for Windows XP (Service Pack 2), I'm
having problems with the "odbcCloseAll" and "odbcClose" functions
within the "RODBC" package. I get pretty much the same error mess
John,
thanks for your response. The fixed effect is estimated at the innermost level.
Andrew
--
Andrew Robinson Ph: 208 885 7115
Department of Forest Resources Fa: 208 885 6226
University of Idaho E : [EMAIL PROTECTED]
PO Box 441133
Dear all,
I obtained and successfully compiled R 2.1.0 on Mac OSX 10.3.9. The
flags used for configuration are:
./configure --enable-utf8 --enable-R-shlib --with-blas='-framework
vecLib' --with-lapack --with-aqua
--with-tcl-config=/Library/Frameworks/Tcl.framework/tclConfig.sh
--with-tk-conf
Hello...
After installing the precompiled version of R 2.1.0 (congratulations
to the R Development Core Team) for Windows XP (Service Pack 2), I'm
having problems with the "odbcCloseAll" and "odbcClose" functions
within the "RODBC" package. I get pretty much the same error message
for both functi
Hello, I must be doing something wrong that's very obvious. But I just
don't see it.
I changed my Windows Time Zone to GMT and my Financial Center to Montreal.
But I still get several warnings.
Sys.timezone()
[1] "GMT Daylight Time"
myFinCenter
[1] "Montreal"
Sys.timeDate()
[1] "Montreal"
[1]
On Wed, 20 Apr 2005 01:36:02 +0800 ronggui wrote:
> it seems that all the existing prop test assume 2 independent or
> matched sample.but in the real world, many situations are not as we
> assume.for example,i do a research on the voter's prefernce through a
> random sampling.and the sample shows
Sorry about the earlier confusion. I'll start fresh here.
I'm getting these error messages when I run "make".
begin installing recommended package foreign
make[2]: *** [foreign.ts] Error 1
make[2]: Leaving directory
`/opt/net/source/R-2.0.1/src/library/Recommended'
make[1]: ***
On Tue, 19 Apr 2005 [EMAIL PROTECTED] wrote:
I am trying to install R-2.1.0 under a Sun platform.
I used tar and then when I type
./configure MAKE=gmake I get the following error:
configure: error: --with-readline=yes (default) and headers/libs are
not available
I would appreciate any help and sug
it seems that all the existing prop test assume 2 independent or matched
sample.but in the real world, many situations are not as we assume.for
example,i do a research on the voter's prefernce through a random sampling.and
the sample shows that 23% of the sample choose A,28% choose B,the others
On Tuesday 19 April 2005 12:23, Barry Rowlingson wrote:
> Deepayan Sarkar wrote:
> > Try (with a sufficiently recent version of R)
> >
> > xyplot(y~x, data=xy, aspect="iso")
> >
> > Is that what you want?
>
> Yes, that's the badger, as we say round here.
>
> Of course my two-point scatterplot i
Has the algorithm used by batchSOM for hexagonal topology changed in some
way between R 1.9.1 and 2.0.1? In comparing between the two versions, I
get identical results for rectangular topology, but very different results
for hexagonal topology. (I tried to find release notes for the VR bundle
On Tue, 2005-04-19 at 12:11 -0500, Deepayan Sarkar wrote:
> On Tuesday 19 April 2005 11:40, Barry Rowlingson wrote:
> > Marc Schwartz wrote:
> > > # Set height value as you require
> > > height <- 3
> > > X11(width = aspect * height, height = height)
> > >
> > > Does that get you what you need
Deepayan Sarkar wrote:
Try (with a sufficiently recent version of R)
xyplot(y~x, data=xy, aspect="iso")
Is that what you want?
Yes, that's the badger, as we say round here.
Of course my two-point scatterplot is a simplification of the actual
data, which consists of nine lines something like a ti
On Tue, 2005-04-19 at 17:40 +0100, Barry Rowlingson wrote:
> Marc Schwartz wrote:
>
> > # Set height value as you require
> > height <- 3
> > X11(width = aspect * height, height = height)
> >
> > Does that get you what you need?
>
> I'm guessing that will work, but it just seems unsatisfa
On Tuesday 19 April 2005 11:40, Barry Rowlingson wrote:
> Marc Schwartz wrote:
> > # Set height value as you require
> > height <- 3
> > X11(width = aspect * height, height = height)
> >
> > Does that get you what you need?
>
> I'm guessing that will work, but it just seems unsatisfactory,
>
Hi.
I am trying to install R-2.1.0 under a Sun platform.
I used tar and then when I type
./configure MAKE=gmake I get the following error:
configure: error: --with-readline=yes (default) and headers/libs are
not available
I would appreciate any help and suggestions.
Best regards,
Cristian Pas
Marc Schwartz wrote:
# Set height value as you require
height <- 3
X11(width = aspect * height, height = height)
Does that get you what you need?
I'm guessing that will work, but it just seems unsatisfactory, since
if I resize my graphics window I mess it up.
If I set par(pty='s') I can h
See FAQ Q7.22
On Tue, 19 Apr 2005, Michael Jerosch-Herold wrote:
I have a problem with making pdf files. I use the code below:
pdf(file="resindex.pdf",width= 6, height= 7, family = "Helvetica", title = "MBF
Data",bg="white")
bwplot(...
dev.off()
When I run a script with the above code in it, by us
> From: Rolf Turner
>
> In response to my question Andy Liaw wrote:
>
> > Here's one way:
>
>
>
> > > dat$rank <- with(dat, ave(x, f, FUN=rank))
>
>
>
> > Don't know how this rates on the slick-quick-sexy scale...
>
> Thanks Andy. That's definitely a 10 on the slick-quick-sexy
In response to my question Andy Liaw wrote:
> Here's one way:
> > dat$rank <- with(dat, ave(x, f, FUN=rank))
> Don't know how this rates on the slick-quick-sexy scale...
Thanks Andy. That's definitely a 10 on the slick-quick-sexy scale.
I (blush!) didn't know about ``ave()'
Michael Jerosch-Herold wrote:
I have a problem with making pdf files. I use the code below:
pdf(file="resindex.pdf",width= 6, height= 7, family = "Helvetica", title = "MBF
Data",bg="white")
bwplot(...
You need to print() lattice graphics!
This is a FAQ.
Uwe Ligges
dev.off()
When I run a script w
On Tue, 2005-04-19 at 16:28 +0100, Barry Rowlingson wrote:
> Suppose I have the following data I want to scatterplot:
>
> > xy
>[,1] [,2]
> [1,]00
> [2,] 214
>
> I start up a graphics window and fire away:
>
> > plot(xy)
>
> - but because the graphics window is s
Hi Thomas,
thanks for the help, it does work.
Regards,
Firas.
On Tue, 19 Apr 2005, Thomas Lumley wrote:
> On Tue, 19 Apr 2005, Firas Swidan wrote:
>
> > Hi,
> > I have a small R question: how to print a single "\" character? I have the
> > following results:
> >> print("\") does not work
> >
> >>
Firas Swidan wrote:
Hi,
I have a small R question: how to print a single "\" character? I have the
following results:
print("\") does not work
print("\\")
[1] "\\"
You can use the lower level function cat() to print exactly what you want:
> cat("\\")
\>
Notice that not even a newline was printed a
On Tue, 19 Apr 2005, Firas Swidan wrote:
Hi,
I have a small R question: how to print a single "\" character? I have the
following results:
print("\") does not work
print("\\")
[1] "\\"
Use cat("\\").
I need to make the following substitution as well, but it does not work
either:
sub("_","\_","g_g"
Here's one way:
> dat <- data.frame(x=runif(10), f=sample(1:2, 10, replace=TRUE))
> dat
x f
1 0.06928792 1
2 0.55972448 1
3 0.39285280 1
4 0.72547126 1
5 0.69027628 1
6 0.51918727 1
7 0.62048040 2
8 0.42240036 1
9 0.90938702 1
10 0.16013661 2
> dat$rank <- with(dat, ave(x, f,
I have a problem with making pdf files. I use the code below:
pdf(file="resindex.pdf",width= 6, height= 7, family = "Helvetica", title = "MBF
Data",bg="white")
bwplot(...
dev.off()
When I run a script with the above code in it, by using "source", it produces a
PDF file which I can not open w
Ted
I agree that with the example data the "inverse regression" approach would
be adequate, but it would be useful to have a function to predict the
results and confidence intervals correctly. I look forward to your solution
to the calibration problem.
In the mean time, I have looked at the calib
Hi,
I have a small R question: how to print a single "\" character? I have the
following results:
> print("\") does not work
> print("\\")
[1] "\\"
I need to make the following substitution as well, but it does not work
either:
> sub("_","\_","g_g")
[1] "g_g"
Thanks in advance,
Firas.
_
On Tue, 19 Apr 2005, Shields, Rusty (IMS) wrote:
Sorry. I posted the contents of foreign.ts.out from an attempted "make"
that included the -m64 options for the gcc, etc.
We haven't see the one that actually failed, as yet. Please start again
from scratch, keep careful records and send us the for
Suppose I have a data frame with two columns ``district'' and
``score'' --- score is numeric; district may be considered
categorical.
I wish to append to this data frame a third column whose entries are
the ranks of ``score'' ***within*** district.
I've tried fiddling about with tapply() and by(
Barry Rowlingson wrote:
Suppose I have the following data I want to scatterplot:
Is there a right way to do this? The only way I can think is to do the
plot without box and axes, and then add them afterwards, but that could
get very messy. Have I missed something obvious?
Package MASS has eqsc
Hi Deepayan.
Thanks for the clarification. For a next release of Lattice, it would
be helpful if you could add a blurb in the help page for densityplot
explicitly stating that the vectors in a "~ a + b" situation must have
the same length. While that is obvious for xyplot (a and b must have
the
Suppose I have the following data I want to scatterplot:
> xy
[,1] [,2]
[1,]00
[2,] 214
I start up a graphics window and fire away:
> plot(xy)
- but because the graphics window is square, the aspect ratio is
wrong. So I add:
> plot(xy, asp=1)
- now the aspect ratio is c
Sorry. I posted the contents of foreign.ts.out from an attempted "make"
that included the -m64 options for the gcc, etc.
I'm running
./configure --prefix=/opt/net/r
My config.site contains:
CC="gcc"
CPPFLAGS="-I/usr/include"
F77="g77"
LDFLAGS="-L/opt/net/utils/gcc3.4.
Sec.FACEA wrote:
Dear all,
I need some help.
I have some problems trying to use the packages for linear programming.
What are "the packages for linear programming"?
I've loaded the packages to R but when I try to use the function R sends
Have you installed it prior to loading it? Any warnings or
Folks-
On R-help, I have seen several message asking for a simple
guide to making packages under Windows. I recently made a
simple package and I decided to take my notes and make them
into a short guide.
http://gsbwww.uchicago.edu/fac/peter.rossi/research/bayes%20book/bayesm/Making%20R%20Packa
On Tuesday 19 April 2005 02:19, Adam Lyon wrote:
> Hi,
>
> I am trying to do a simple superimposing of two density plots and am
> seeing strange behavior (I haven't seen this reported elsewhere).
> This little snipit will make it apparent...
>
> a = rnorm(10)
> b = rnorm(100, mean=-2, sd=0.5)
> den
Dear all,
I need some help.
I have some problems trying to use the packages for linear programming.
I've loaded the packages to R but when I try to use the function R sends
an error warning.
I don't know if something is missing
This is what it shows
# Set up problem: maximize
# x1 + 9 x2 + x3 s
Hi there,
The optim(par, fn, ...parscale...) function in R requires 'parscale' which is
defined as:
"A vector of scaling values for the parameters. Optimisation is performed on
'par/parscale' and these should be comparable in the sense that a unit change
in any element (??) produces a unit c
On Tue, 19 Apr 2005 16:18:14 +0200 Faouzi LYAZRHI wrote:
> Hi,
> I would like to select a few cases (for example cases corresponding to
> sex=male) to make summary for another variable.
Look at ?tapply. E.g.
tapply(x, group, summary)
Z
> How can I do this.
> Thanks for your help
> Fawtzy
>
>
try this litte example:
save the code below in a file test.Rnw and call then
Sweave("test.Rnw", driver = RweaveHTML())
--
testing r2html
look at this: here you can write some text
.
<>=
summary(data.frame(c(1,2,3), c(3,4,5)))
@
insert some graphics
<>=
print(plot(c(1:30), c(1:30)))
@
--
hope it
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
The Tuesday update script came back to get me! I knew it would.
update.packages has changed (for the better) with this release. Look at the
NEWS file:
The 'CRAN' argument to update.packages(), old.packages(),
new.packages(), download.packages() and install.packages() is
de
subset(your.data.frame, your.data.frame$sex == 'male')
cheers
c
Faouzi LYAZRHI wrote:
Hi,
I would like to select a few cases (for example cases corresponding to
sex=male) to make summary for another variable.
How can I do this.
Thanks for your help
Fawtzy
_
Andy
Thanks, the link was very helpful. I havn't checked the code for the calib
function but it appears to work with my data.
[NB the return function generates a warning message with later verions of R
and needs to be amended to return a list]
In the last line of my code, calibration.lm should ha
presuming that you have a data.frame look at "?subset()", e.g.,
# dat is your data.frame
subset(dat, sex=="male")
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Be
Hi,
I would like to select a few cases (for example cases corresponding to
sex=male) to make summary for another variable.
How can I do this.
Thanks for your help
Fawtzy
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-he
Hello,
Running R 2.1.0 in a Win XP
I put the snipped to automatically update my libraries on Tuesdays that was
presented to the list some time back, It worked with no problem for R 2.0.1 but
now that I installed R 2.1.0 and copy the old Rprofile to the new r/etc I get
an error.
This is my Rprofile:
Thanks,
EJ
Prof Brian Ripley wrote:
Please do this for the *sources* and not the installed package.
If I do, it works for me.
On Tue, 19 Apr 2005, Ernesto Jardim wrote:
Hi,
I want to compile the manual for "methods" so that I can have the
usual indexed man pages in a single document with a ToC. I'
I recently learnt how to use Sweave which is a wonderful tool
After which i also tried to use R2HTML as it would allow many of my
colleagues who dont use latex to be able to use and edit my work. I was
unable to make it work and couldnt find a way to implement it. I got some
errors.
I wonder if y
William Valdar wrote:
From: Brian Ripley
As a first shot, use lm with a matrix response. That fits them all at
once with one QR-decomposition. No analogue for glm or lmer, though,
since for those the iterative fits run do depend on the response.
Thanks Brian, that's very helpful. Also thanks
Is something like this what you're looking for?
fixSome <- function(..., fixed) {
fList <- list(...)
for (i in fixed) fList[[i]] <- fList[[i]][1]
do.call("expand.grid", fList)
}
> age <- c("young","mid","old")
> married <- c("no","yes")
> income<- c("low","high","medium")
Please do this for the *sources* and not the installed package.
If I do, it works for me.
On Tue, 19 Apr 2005, Ernesto Jardim wrote:
Hi,
I want to compile the manual for "methods" so that I can have the usual
indexed man pages in a single document with a ToC. I've tried with
R CMD Rd2dvi --pdf /u
Hi,
I want to compile the manual for "methods" so that I can have the usual
indexed man pages in a single document with a ToC. I've tried with
R CMD Rd2dvi --pdf /usr/local/lib/R/library/methods
but the result is odd, the ToC shows only "as" and "setOldClass" and it
does not compile all man page
OK, the hint by Dimitris applied I just do very simple:
test <- array(c(1:16), dim = c(3,4,3))
test
## I call some enries using an index array
test.ind <- array(rbind(c(1,2,1), c(3,3,2)), dim = c(2,3))
test[test.ind]
## suppose I want all values in the 2nd row and 4th col over
## all three 3rd dime
Dear R-user,
I have a data like this below,
age <- c("young","mid","old")
married <- c("no","yes")
income<- c("low","high","medium")
gender<- c("female","male")
I want to make some of combination data like these,
age.income.dat <- expand.grid(age, married[-c(2)], income
> 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
in the specific example you could consider something like this:
test <- array(c(1:16), dim = c(3,4,3))
##
test.ind <- rbind(c(2,4,1), c(2,4,2), c(2,4,3))
test[test.ind]
test[2,4,]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Ca
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 new.abs
data.frame
On Tue, 19 Apr 2005, Dr. Peter Schlattmann wrote:
when performing a meta analysis I have two results obtained with logLik
and lme, which I do not quite understand.
The first result concerns logLik
m0<-glm(or~1,family=gaussian(),data=temp,weights=1./var)
logLik(m0)
`log Lik.' -7.10697 (df=2)
For co
Hi
I have the following array:
test <- array(c(1:16), dim = c(3,4,3))
test
## I call some enries using an index array
test.ind <- array(rbind(c(1,2,1), c(3,3,2)), dim = c(2,3))
test[test.ind]
## suppose I want all values in the 2nd row and 4th col over
## all three 3rd dimensions
test[2,4,]
how to
I'm sorry if this is too off-topic--feel free to ignore. I am
interested in using pl/R, an amazing "plugin" for the postgresql
database. As is typical of these types of applications, pl/R needs to
link against a shared library. However, it appears that the MacOS R
binary does not build a sta
Shields, Rusty (IMS) wrote:
Did I post this to the right mail list? Is there another list that is
more appropriate for this type of question?
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Shields, Rusty
(IMS)
Sent: Monday, April 18, 2005 12:09 PM
To: r-h
Is the fixed effect estimated at the innermost level? If not,
plots of residuals at that level are surely of limited interest.
qqplots, to be relevant, surely need to assess normality of
effects (rather than residuals) at the level that matters for
the intended inferences.
If the fixed effect is e
Did I post this to the right mail list? Is there another list that is
more appropriate for this type of question?
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Shields, Rusty
(IMS)
Sent: Monday, April 18, 2005 12:09 PM
To: r-help@stat.math.ethz.ch
Subje
Dear all,
when performing a meta analysis I have two results obtained with logLik
and lme, which I do not quite understand.
The results are based on these data:
study or var
1 0.10436 0.299111
2 -0.03046 0.121392
3 0.76547 0.319547
4 -0.19845 0.025400
5 -0.10536 0.025041
6 -0.1165
On 19-Apr-05 Ralf Strobl wrote:
> Dear List,
> can anyone explain me this result (Windows XP, R 2.0.1):
>
> > (0.2-0.1)==0.1
> [1] TRUE
> > (0.3-0.2)==0.1
> [1] FALSE
>
> Regards,
> Ralf Strobl
It is a consequence of the finite length of the binary
expression of decimal fractions, which is not
From: Brian Ripley
As a first shot, use lm with a matrix response. That fits them all at once
with one QR-decomposition. No analogue for glm or lmer, though, since for
those the iterative fits run do depend on the response.
Thanks Brian, that's very helpful. Also thanks to Kevin Wright who
su
> >> example(lm)
> > Error in switch(x[2], "1250" = return("ISO 8859-2"), "1251" =
> > return("KOI8-U"), :
> >argument is missing, with no default
--- src/library/utils/R/iconv.R.orig2005-04-18 19:18:58.0
+0900
+++ src/library/utils/R/iconv.R 2005-04-19 17:37:32.0 +0
On Tue, 19 Apr 2005, ronggui wrote:
Prof Brian Ripley,
I found the error by chance. the error is found in R2.1.0 too.
As far as I can see, if anyone had run any example() in R on a Windows
machine in your locale it would have failed. But no one reported a
problem. I don't even have that local
Look at "?Comparison", especially in the "Note" section:
...
For numerical values, remember == and != do not allow for the finite
representation of fractions, nor for rounding error. Using all.equal
with identical is almost always preferable. See the examples.
...
Best,
Dimitris
Dimitris Ri
1 - 100 of 104 matches
Mail list logo