Hi, just a general question: when we do hierarchical clustering, should we
compute the dissimilarity matrix based on scaled dataset or non-scaled dataset?
daisy() in cluster package allow standardizing the variables before calculating
dissimilarity matrix; but dist() doesn't have that option at
Hi, I am using read.ssd() from foreign package to read some SAS datasets. I
have
2 types of SAS datasets, one with sas7bdat extension, the other with ssd01
extension. I have no problem with the first dataset type, but got the following
error message with the 2nd dataset type (with ssd01
Looks like the log file is not appropriately attached. Here it is again. Thanks
for any suggestions.
John
- Original Message
From: array chip arrayprof...@yahoo.com
To: r-help@r-project.org
Sent: Mon, August 2, 2010 2:18:32 PM
Subject: [R] read SAS dataset using read.ssd()
Hi, I am
HI, I was trying to install xlsx package for reading in Excel 2007 files. The
installation went smoothly. But when I tried to load the library, I got the
following error message:
library(xlsx)
Loading required package: xlsxjars
Loading required package: rJava
Error : .onLoad failed in
Problem solved. I need to set Environment variable in windows, not inside R.
1. Click through: start control panel system
2. Click the Advanced tab.
3. Click Environment variables.
4. Click New (under 'system') to add a new variable name and value.
- Original Message
From: array
Hi, I would like to a draw a scatterplot of x1 and x2 (plot (x1, x2)), and also
want to draw a sort of regression line across the data points. But x1 and x2
are just 2 independent variables, so in this case a regression of x1 over x2,
or
vice versa, is not appropriate per se. What would be an
Hi, I have a simple plot by plot(x,y, log='xy'). However, due to the large
range
of values of x, the x-axis annotation is printed as 2e+02 1e+03 5e+03 2e+04
1e+05 instead of 200 1000 5000 2 10. How can I make it printed
as in the later one?
Thanks
John
Thank you David!
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Sent: Mon, August 9, 2010 12:42:06 PM
Subject: Re: [R] x-axis annotation
On Aug 9, 2010, at 3:14 PM, array chip wrote:
Hi, I have a simple
Hi how can print x-axis labels in 45 degree in boxplot() (or plot in general)?
I
can use las=2 to print in 90 degree, but it looks ugly. Is there a simple
option
to do 45 degree easily?
Thanks
John
__
R-help@r-project.org mailing list
Than you Marc.
John
- Original Message
From: Marc Schwartz marc_schwa...@me.com
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Sent: Thu, August 12, 2010 12:17:12 PM
Subject: Re: [R] x-axis label print in 45 degree
On Aug 12, 2010, at 2:14 PM, array chip wrote:
Hi
I searched with print x-axis label in 45 degree which didn't return useful
links. Apparently I used poor search keywords.
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: Marc Schwartz marc_schwa...@me.com
Cc: array chip arrayprof...@yahoo.com; r-help@r
Hi,
When I plot, the axis ticks are printed as 50.00 25.00 10.00 1.00 0.05
0.01, is there any way to print them as 50 25 10 1 0.05 0.01 instead?
Thanks
John
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
suggestions?
Many thanks
John
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc: R-Help r-h...@stat.math.ethz.ch
Sent: Mon, August 16, 2010 1:19:57 PM
Subject: Re: [R] print numbers
On Aug 16, 2010, at 3:53 PM, array chip wrote:
Hi
Hi, I am really puzzled by this. hope someone can help me
I have a 2 small data frames a and b derived from a larger data frames.
They
look exactly the same to me, but identical() always returns FALSE.
a
a b
2 10011048 L
4 10011048 R
6 10011049 L
8 10011049 R
b
a b
1
Thank you all for the suggestions. They all worked. Now, those numbers on the
ticks are produced by default plot function bxp(), i.e. they are different in
each plot, is there any functions I can retrieve them so I can use formatC() or
prettyNum() etc?
Thanks again,
John
- Original
Oops, I overlooked the row names. Sorry for my carelessness.
Thanks
- Original Message
From: Erik Iverson er...@ccbr.umn.edu
To: array chip arrayprof...@yahoo.com
Cc: r-h...@stat.math.ethz.ch
Sent: Mon, August 16, 2010 2:53:37 PM
Subject: Re: [R] identical()
Hello,
array chip
Marc, this works perfectly!
Thanks
John
- Original Message
From: Marc Schwartz marc_schwa...@me.com
To: William Dunlap wdun...@tibco.com
Cc: array chip arrayprof...@yahoo.com; r-help@r-project.org
Sent: Mon, August 16, 2010 3:23:33 PM
Subject: Re: [R] print numbers
Bill et al,
See
Hi all, let me give a simple example:
b-20
I would like to print ylab as P20 where P is printed in Italic font. When I
do the following:
plot(1, ylab=expression(paste(italic(P),b,sep=)))
I got y axis label printed as Pb instead of P20. What is the best solution
to print platmath symbols with
Thanks, yes it worked!
What about if I want to print as P2, A where A is just letter A and 2 is from
variable b.
John
- Original Message
From: baptiste auguie baptiste.aug...@googlemail.com
To: array chip arrayprof...@yahoo.com
Cc: r-h...@stat.math.ethz.ch
Sent: Thu, August 19, 2010
Thanks David!
John
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc: baptiste auguie baptiste.aug...@googlemail.com; r-h...@stat.math.ethz.ch
Sent: Thu, August 19, 2010 11:34:07 AM
Subject: Re: [R] plotmath question
On Aug 19
David, yes, I now see how it worked.
Thanks again,
John
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc: baptiste auguie baptiste.aug...@googlemail.com; r-h...@stat.math.ethz.ch
Sent: Thu, August 19, 2010 12:12:46 PM
Subject
Hi, I am a beginner of xyplot() (or lattice package). On one hand, I
immediately
realized it's a very powerful utility. On the other hand, there are too many
things for me to learn. Still haven't figure out a generalization of the syntax
and usage under many different circumstances.
Let me
What a simple way to do what I want to do! Thanks.
But if there is missing data in variable y, then the averaged lines is broken
where the missing data is present. For example:
dat$y[c(10,185,200,400,450)]-NA
then using type='a' option will result in broken lines.
xyplot(y~day|sex, groups=trt,
Hi, is there anyway I can retrieve the user coordinates for the region of the
heatmap (only the heatmap, not include dendrogram, x- y- axis annotation). I
found that par(usr) didn't give the user coordinates that I want. I want
those
user coordinates to add some additional information to the
Hi all,
I asked this before the holiday, didn't get any response. So would like to
resend the message, hope to get any fresh attention. Since this is not purely
lme technical question, so I also cc-ed R general mailing list, hope to get
some
suggestions from there as well.
I asked some
Hi, suppose I have a data frame as below:
dat-cbind(expand.grid(id=c(1,2,3),time=c(0,3,6),mode=c('R','L'),rep=(1:3)),y=rnorm(54))
I kind of want to squeeze the data frame into a new one with averaged y
over
rep for the same id, time and mode. taking average is easy with tapply:
tapply(dat$y,
R 3 1.175192908 -0.21924078
But I only need one row:
id time modeavg.y
110 R-0.21924078
Thanks again,
John
- Original Message
From: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc
Thank you Peter, yes this is what I need!
John
- Original Message
From: Peter Alspach peter.alsp...@plantandfood.co.nz
To: array chip arrayprof...@yahoo.com; David Winsemius
dwinsem...@comcast.net
Cc: r-help@r-project.org r-help@r-project.org
Sent: Thu, September 9, 2010 4:26:53 PM
from the output, I think it's both.
- Original Message
From: John Sorkin jsor...@grecc.umaryland.edu
To: r-help@r-project.org
Sent: Fri, September 10, 2010 5:25:44 AM
Subject: [R] lme, groupedData, random intercept and slope
Windows Vista
R 2.10.1
Does the following use of
But as far as I know, profile() seems to be de-activated in the lme4 package.
- Original Message
From: Gavin Simpson gavin.simp...@ucl.ac.uk
To: John Sorkin jsor...@grecc.umaryland.edu
Cc: r-help@r-project.org; Bert Gunter gunter.ber...@gene.com
Sent: Fri, September 10, 2010 2:05:37 AM
splines_2.11.1 stats4_2.11.1 tools_2.11.1
Any suggestions would be appreciated.
John
- Original Message
From: Gavin Simpson gavin.simp...@ucl.ac.uk
To: array chip arrayprof...@yahoo.com
Cc: John Sorkin jsor...@grecc.umaryland.edu; r-help@r-project.org; Bert
Gunter
gunter.ber
there is
pretty slow,..
Thanks,
John
- Original Message
From: Gavin Simpson gavin.simp...@ucl.ac.uk
To: array chip arrayprof...@yahoo.com
Cc: John Sorkin jsor...@grecc.umaryland.edu; r-help@r-project.org; Bert
Gunter
gunter.ber...@gene.com
Sent: Fri, September 10, 2010 10:46:16 AM
Hi all, I
When I plot both lines and points using type=c('l', 'p') in xyplot(), if I want
to include in legend both of them using keys=list(lines=list(col=1:3),
points=list(pch=1:3)), the lines and points are plotted side by side in legend.
Is there anyway to plot the points in the middle of
Hi, another question: is there any argument that controls the line width of
axis
box of xyplot()? I tried lwd=2 or lwd.axis=2 in xyplot() or within
scales=list()
argument, without success.
Thanks
John
__
R-help@r-project.org mailing list
: David Winsemius dwinsem...@comcast.net
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Sent: Mon, September 13, 2010 4:05:04 PM
Subject: Re: [R] xyplot legends
On Sep 13, 2010, at 6:25 PM, array chip wrote:
Hi all, I
When I plot both lines and points using type=c('l', 'p
: Daisy Englert Duursma daisy.duur...@gmail.com
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Sent: Mon, September 13, 2010 8:05:53 PM
Subject: Re: [R] xyplot axis line width
check out ?par for all the details on plotting
‘mgp’ The margin line (in ‘mex’ units) for the axis title
Thank you Deepayan. This is exactly what I needed.
John
- Original Message
From: Deepayan Sarkar deepayan.sar...@gmail.com
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Sent: Tue, September 14, 2010 4:52:29 AM
Subject: Re: [R] xyplot axis line width
On Tue, Sep 14
Hi, I asked this on mixed model mailing list, but that list is not very active,
so I'd like to try the general R mailing list. Sorry if anyone receives the
double post.
Hi, I have a dataset of animals receiving some eye treatments. There are 8
treatments, each animal's right and left eye was
) is the same
between aov, lmer and lme. But I am not sure how to relate results between aov
and lmer/lme for the other 2 variance components (id and eye%in%id).
Thanks
John
- Original Message
From: Peter Dalgaard pda...@gmail.com
To: array chip arrayprof...@yahoo.com
Cc: r-help@r
from glht() reasonable/meaningful? If not, will the suggested
1-way ANOVA used together with glht() give us correct post-hoc multiple
comparsion results?
Thank you very much!
John
- Original Message
From: Peter Dalgaard pda...@gmail.com
To: array chip arrayprof...@yahoo.com
Cc: r
Thank you Peter and Ben for your comments.
John
- Original Message
From: Peter Dalgaard pda...@gmail.com
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org; r-sig-mixed-mod...@r-project.org
Sent: Mon, September 20, 2010 12:28:43 PM
Subject: Re: [R] lmer() vs. lme() gave
Is there anyway to make plotting point character being thicker in xyplot? I
mean
not larger which can achieved by cex=2, but thicker. I tried lwd=2, but it
didn't work. I know lwd works in regular plot() not only for lines, but also
for points. For example
plot(1:10, lwd=2)
Thanks
John
Thank you Greg. I also got it work by using panel.points (lwd=2) instead of
using panel.xyplot()
- Original Message
From: Greg Snow greg.s...@imail.org
To: array chip arrayprof...@yahoo.com; r-help@r-project.org
r-help@r-project.org
Sent: Thu, September 23, 2010 2:48:06 PM
Subject
,pch=0:1,cex=2,lwd=2),
text=list(lab=c('A','B'),cex=1.5,font=2)))
Any suggestions?
Thanks
John
- Original Message
From: array chip arrayprof...@yahoo.com
To: Greg Snow greg.s...@imail.org
Cc: r-help@r-project.org
Sent: Thu, September 23, 2010 4:03:00 PM
Subject: Re: [R] how to make
Yes, it does what I want. Thank you Peter! Just wondering what else grid.pars
controls? not just the symbol in legend, right?
John
- Original Message
From: Peter Ehlers ehl...@ucalgary.ca
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org r-help@r-project.org
Sent: Thu
Hi I am using xyplot() to plot on the log scale by using scale=list(log=T)
argument. For example:
xyplot(1:10~1:10, scales=list(log=T))
But the axis labels are printed as scientific notation (10^0.0, etc), instead
of
fixed notation. How can I change that to fixed notation?
options(scipen=4)
Thanks for the suggestion. But my example is just an example, I would prefer to
have some generalized solution, like what options(scipen=4) does in general
graphics, which usually gave pretty axis labels as well.
Any suggestions?
Thanks
John
From: Henrique
may not be lists when relation = same
Syntax problem in this last command?
Thanks
On 27-Sep-10, at 12:16 PM, Henrique Dallazuanna wrote:
Try this:
xyplot(1:10~1:10, scales=list(log = T, labels = round(log(1:10), 4)))
On Mon, Sep 27, 2010 at 4:10 PM, array chip arrayprof...@yahoo.com
, 10)
ans$left$labels$labels - as.character(tick.at)
ans }
- Original Message
From: array chip arrayprof...@yahoo.com
To: Don McKenzie d...@u.washington.edu; Henrique Dallazuanna
www...@gmail.com
Cc: R-help Forum r-help@r-project.org
Sent: Mon, September 27, 2010 3:09:20 PM
Thank you Dr. Sarkar. yscale.components.log10.3 is pretty good choice.
John
- Original Message
From: Deepayan Sarkar deepayan.sar...@gmail.com
To: array chip arrayprof...@yahoo.com
Cc: Don McKenzie d...@u.washington.edu; Henrique Dallazuanna
www...@gmail.com; R-help Forum r-help@r
Hi, I am wondering if anyone can propose a simple/best way to do the following:
Let's say I have a data frame
dat -
cbind(expand.grid(mode=c('right','left'),time=0:3,id=c('p1','p2','p3')),y=c(3,5,rep(4,6),6,2,rep(3,6),4,4,rep(2,6)))
dat
mode time id y
1 right0 p1 3
2 left0 p1
))
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
spec...@stat.berkeley.edu
On Fri, 1 Oct 2010, array chip wrote:
Hi, I am wondering if anyone can propose a simple/best way to do
Hi, is there a way to retrieve the extremes of the user coordinates of the
plotting region, like what par(usr) does in general graphics? I'd like to use
them to print additional texts at certain place inside each panel. Thanks
John
[[alternative HTML version deleted]]
Hi, how can I make the point characters thicker (NOT larger) in xyplot when
groups= argument is used?
dat-data.frame(x=1:100,y=1:100,group=rep(LETTERS[1:5],each=20))
### lwd=2 doesn't work here
xyplot(y~x,groups=group,data=dat,col=1:4,pch=1:4,lwd=2)
### lwd=2 works with panel.points(),
Hi, is there a package for getting type II or type III tests on mixed models
(lme or lmer), just like what Anova() in car package does for aov, lm, etc.?
Thanks
John
[[alternative HTML version deleted]]
__
R-help@r-project.org
Dear R users, I am not asking questions specifically on R, but I know there are
many statistical experts here in the R community, so here it goes my questions:
Freedman (1982) propose an approximation of sample size/power calculation based
on log-rank test using the formula below (This is what
%
significance level with 80% power
any comments are appreciated.
Â
John
--- On Thu, 5/6/10, Joris Meys jorism...@gmail.com wrote:
From: Joris Meys jorism...@gmail.com
Subject: Re: [R] sample size for survival curves
To: array chip arrayprof...@yahoo.com
Date: Thursday, May 6, 2010, 8:12 PM
Thanks Kevin. I thought the time t is at the end of follow-up (length of
follow-up)?
John
--- On Thu, 5/6/10, Kevin E. Thorpe kevin.tho...@utoronto.ca wrote:
From: Kevin E. Thorpe kevin.tho...@utoronto.ca
Subject: Re: [R] sample size for survival curves
To: array chip arrayprof
Hi, I would like to use the curve() function to draw the predicted curve from
an nls() object. for example:
dd-read.table(dd.txt,sep='\t',header=T,row.names=1)
obj-nls(y~c+(d-c)/(1+(x/e)^b),data=dd,start=list(b=-1, c=0, d=100, e=150))
coef(obj)
b c d e
Subject: Re: [R] automate curve drawing on nls() object
To: array chip arrayprof...@yahoo.com, r-help@r-project.org
Date: Tuesday, May 18, 2010, 7:42 PM
I can't directly answer your question
regarding 'expression', but can you just replace b, c,d, and
e with coef(obj)[1], coef(obj)[2], ...
etc
18, 2010, 7:33 PM
array chip arrayprofile
at yahoo.com writes:
Hi, I would like to use the curve() function to draw
the
predicted curve from an nls() object. for example:
Is there a reason you don't want to use plot() and
predict() ???
dd-read.table(dd.txt,sep='\t',header
I do know it, but wanted to spend a little more effort to make generalized
function for this type of plot.
Thanks
--- On Tue, 5/18/10, Shi, Tao shida...@yahoo.com wrote:
From: Shi, Tao shida...@yahoo.com
Subject: Re: [R] automate curve drawing on nls() object
To: array chip arrayprof
Hi, I am currently trying to do dose-response curves
using weighted 4-parameter model (4PL). The weighting was based on
1/(expected variance) derived from historical data. I tried both drm() from drc
package, and nls(), found very different
results derived from drm() vs. nls() using weights=
Hi, I am currently trying to do dose-response curves
using weighted 4-parameter model (4PL). The weighting was based on
1/(expected variance) derived from historical data. I tried both drm() from drc
package, and nls(), found very different
results derived from drm() vs. nls() using weights=
Hi, is there a way to create one image file (like using win.metafile(), bmp(),
etc) that contained multiple pages of plots, just like what postscript() does
in creating PDF file?
Thanks
John
__
R-help@r-project.org mailing list
Hi, I am wondering if there is a way to extract studentized residues from a
nls() object? I searched archive, someone posted a similar question before
(http://tolstoy.newcastle.edu.au/R/e6/help/09/04/10845.html), but didn't get an
answer.
Thanks for any suggestions,
John
Hi, I am wondering how I can specify no intercept in a mixed model using
lmer().
Here is an example dataset attached (test.txt). There are 3 workers, in 5
days, measured a response variable y on independent variable x. I want to
use a quadratic term (x2 in the dataset) to model the
Hi, is it possible to specify a constant intercept (based on prior knowledge)
in linear regression using lm()?
Thanks
John
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Hi, I asked this before, but haven't got any response. So would like to have
another try. thanks for help. Also tried twice to join the model mailing list
so that I can ask question there, but still haven't got permission to join that
list yet.
===
Hi, I am wondering how I can
Hi, this is not R technical question per se. I know there are many excellent
statisticians in this list, so here my questions: I have dataset with ~1800
observations and 50 independent variables, so there are about 35 samples per
variable. Is it wise to build a stable multiple logistic model
on interactions. No point in trying I'd say. If you really need all
that information, you might want to take a look at some dimension
reduction methods first.
Cheers
Joris
On Mon, Jun 14, 2010 at 2:55 PM, Marc Schwartz marc_schwa...@me.com wrote:
On Jun 13, 2010, at 10:20 PM, array chip
Hi, I am learning xyplot. I have an example dataset attached.
plotdata-read.table(plotdata.txt,sep='\t',header=T,row.names=1)
head(plotdata,n=4)
y x type
1 -4.309601 -0.7448405A
2 -4.715421 0.7875994A
3 -2.310638 0.5455310A
4 -2.685803 10.4116868A
+0.1117*x-0.0006357*x*x)
})
type='b' doesn't give me a smooth line, so it didn't produce what I want.
Thanks
John
- Original Message
From: Bert Gunter gunter.ber...@gene.com
To: array chip arrayprof...@yahoo.com; r-help@r-project.org
Sent: Thu, July 8, 2010 4:25:51 PM
Subject: RE: [R
Hi, I have been playing with graphic parameters to
meet my specific plot requirement. I would like to
draw a plot in the top right region of the plot
region, I can do that with a number of parameters such
as plt, mar or oma. For example setting the
parameters:
par(plt=c(0.28,0.956,0.25,0.9))
or
Hi is there a package for diagonal linear discriminant
analysis (diagonal LDA)?
Thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
Hi I am using forestplot() in rmeta package on a
dataset of 45 point estimates with corresponding
confidence intervals. The resulting plot was just too
large and plotted out of the graphic window that I can
not see whole picture. Is there anyway to fix this
problem?
Thanks
Hi, it seems that the fastbw() in the Design package
only works with variable of class factor according
to the help page if I understand correctly. Is there
any R function/package that do stepwise variable
selection for a Cox model with continuous independent
variables?
Thank you
John
Rob Tibshirani propose to use lasso with Cox
regression for variable selection in his 1997 paper
The lasso method for variable selection in the Cox
model published in Statistics In Medicine 16:385. I
understand the lars() function in lars package
implemented lasso, but it does not do lasso with
Hi,
When I install new packages from CRAN, I frequently
find that some packages were missing from the download
queue. For example, on one of my computer with R2.6.2,
I can not find package glmpath from the download
queue. On my other computer with R2.5.1, I could still
find that particular
Hi,
I am new to model selection by coefficient shrinkage
method such as lasso. And I became particularly
interested in variable selection in Cox regression by
lasso. I became aware of the coxpath() in R package
glmpath does lasso on Cox model. I have tried the
sample script on the help page of
Hi, I am wondering how can I add a common title to the
top of a page that contains multiple plots using,
e.g., par(mfrow=c(2,3))?
thanks
John Zhang
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of
array chip
Sent: Wednesday, 5 March 2008 10:48 a.m.
To: [EMAIL PROTECTED]
Subject: [R] multiple plots with a title
Hi, I am wondering how can I add a common title to
the top of
a page
Hi I got an error message using datadist() from Design package:
library(Design,T)
dd - datadist(beta.final)
options(datadist=dd)
lrm(Disease ~ gsct+apcct+rarct, x=TRUE, y=TRUE)
Error in eval(expr, envir, enclos) : object Disease not found
All variables inclduing response variable
Hi, another question about validate() in Design library. The arugment B of
this function is number of repetition for method=bootstrap, which is easy to
understand; but for method=crossvalidation, B is the number of groups of
omitted observations. This is confusing, I don't understand what it
:
From: Jorge Ivan Velez jorgeivanve...@gmail.com
Subject: Re: [R] datadist() in Design library
To: array chip arrayprof...@yahoo.com
Cc: R mailing list r-help@r-project.org
Date: Thursday, July 9, 2009, 6:55 PM
Dear John,
Have you tried it specifying the 'data'
argument as suggested in lrm
Hi, I am wondering if there is a way to specify the prevalence of events in
logistic regression using lrm() from Design package? Linear Discriminant
Analysis using lda() from MASS library has an argument prior= that we can use
to specify the prevalent of events when the actual dataset being
Hi, I would like to rbind 2 data frames. They both some common column names,
but also some unique column names each, is there any simple function that rbind
these 2 data frames with filling NAs for those columns of unique names?
thanks
__
Hi, I am just wondering if there is a test available for testing if a linear
fit of an independent variable in a Cox regression is enough? Thanks for any
suggestions.
John Zhang
[[elided Yahoo spam]]
Hi, is there a function to calculate fractional ranks?
Thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal,
Hi, all, I found that the smooth.spline() function produces different results
between R and S-Plus. I was trying to play different parameters of the function
without any success. The script of the function contains Fortran code, so it
seems impossible to port the code from S-Plus to R (or I
PROTECTED] wrote:
supclust
On 11/9/07, array chip [EMAIL PROTECTED]
wrote:
Hi is there a package for diagonal linear
discriminant
analysis (diagonal LDA)?
Thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
Hi,
The code round(0.0002234,4) will return 2e-04, how can
I return 0.0002 instead? I want to print 0.0002
instead of 2e-04 for formating purpose.
thanks
Looking for last minute shopping deals?
Hi,
I am trying to test some contrasts, using glht() in
multcomp package on fixed effects in a linear mixed
model fitted with lme() in nlme package. The command I
used is:
## a simple randomized block design,
## type is fixed effect
## batch is random effect
## model with
Hi, I have been trying glht() from multcomp package
and contrast() from contrast package to test a
contrast that I am interested in.
With the following simulated dataset (fixed effect
type with 3 levels (b, m, t), and random effect
batch of 4 levels, a randomized block design with
interaction),
Hi, I am running a simple one-way ANOVA with an
independent factot variable treat (3 levels: a, b
and c) and a response variable y. I want to test a
linear relationship of the response among the 3 levels
of the variable treat (ordered a-b-c). I used
glht() from multcomp package. Later I found out
to be ordered then the
default contrasts
will be poly()nomials.
Ross Darnell
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
On Behalf Of array chip
Sent: Friday, 18 January 2008 10:22 AM
To: [EMAIL PROTECTED]
Subject: [R] how to specify a particular contrast
Hi
Hi, I am trying to use F test or Chi-square test to test if 2 5-parameter (A,
B, xmid, scal and H) logistic curves are parallel based on residual sum of
squares.
What's usually done is to first fit the 2 curves using a constraint (or global)
model where all parameters are kept the same except
()
To: array chip arrayprof...@yahoo.com
Cc: r-help@r-project.org
Date: Saturday, March 13, 2010, 9:12 AM
Hi John,
you may have a look at the drc-package (or look at the
maintainers web
site www.bioassay.dk)
I think it does exactly what you want.
HTH
Eik
array chip schrieb:
Hi, I am
Hi, I have very simple balanced randomized block design where I total have 48
observations of a measure of weights of a product, the product was manufactured
at 4 sites, so each site has 12 observations. I want to use lme() from nlme
package to estimate the standard error of the product weight.
Thanks Dennis for the thorough explanation and correction on the design.
John
--- On Thu, 4/1/10, Dennis Murphy djmu...@gmail.com wrote:
From: Dennis Murphy djmu...@gmail.com
Subject: Re: [R] trying to understand lme() results
To: array chip arrayprof...@yahoo.com
Cc: R-help@r-project.org
Date
1 - 100 of 277 matches
Mail list logo
