Dear useRs,Is there a command or package in R to do heirarchical clustering by
"modified ward method"??regardseliza
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat
On Sep 29, 2012, at 7:10 AM, Josh Browning wrote:
> Hi useRs,
>
> I'm experiencing something quite weird with glm() and weights, and
> maybe someone can explain what I'm doing wrong. I have a dataset
> where each row represents a single case, and I run
> glm(...,family="binomial") and get my co
On Sep 29, 2012, at 12:46 PM, David Winsemius wrote:
>
> On Sep 29, 2012, at 12:09 PM, Meredith Ballard LaBeau wrote:
>
>> David:
>> I set this up and it still seems to plot the same:
>>
>> maumnf_plot<-boxplot(log_loads~ind,data=nfmaum, horizontal=TRUE, notch=T,
>> outline=FALSE, whisker=0,
Hello,
Note that some graphics parameters of boxplot and bxp are not the same.
For instance 'col' becomes 'boxfill'.
boxp <- function(x, ...){
bp <- boxplot(x, ..., plot = FALSE)
bp$stats[1, ] <- bp$stats[2, ]
bxp(bp, ...)
}
x <- rnorm(1000)
boxplot(x, col ="blue", notch = TRUE)
On Sep 29, 2012, at 12:09 PM, Meredith Ballard LaBeau wrote:
> David:
> I set this up and it still seems to plot the same:
>
> maumnf_plot<-boxplot(log_loads~ind,data=nfmaum, horizontal=TRUE, notch=T,
> outline=FALSE, whisker=0, main="Maumee River Near Future Climate Scenarios",
> ylab="Log L
On Sep 29, 2012, at 10:23 AM, Meredith Ballard LaBeau wrote:
> Good Afternoon-
> I was wanting to alter the boxplot to remove the lower whisker, both the
> whisker line and staple just on the lower end. Is there a way to do this?
> As my code is currently:
> boxplot(log_loads~ind,data=nfmaum, ho
Good Afternoon-
I was wanting to alter the boxplot to remove the lower whisker, both the
whisker line and staple just on the lower end. Is there a way to do this?
As my code is currently:
boxplot(log_loads~ind,data=nfmaum, horizontal=TRUE, notch=T, outline=FALSE,
whisker=0, main="Maumee River Nea
On Sep 29, 2012, at 6:33 AM, Sabine Dietmann wrote:
> Hi All,
>
> I wonder whether it is possible to use the plot function and reduce
> the scales on both +/- quadrants, i.e.
Those are not quadrants, but rather half-planes.
> plot the data from -15. -8 and 8..15 only.
There is a broken axis f
On 29-09-2012, at 15:33, Sabine Dietmann wrote:
> Hi All,
>
> I wonder whether it is possible to use the plot function and reduce
> the scales on both +/- quadrants, i.e.
> plot the data from -15. -8 and 8..15 only.
>
Have a look at
?plot.default
and try
plot(x,y,xlim=c(-15,-8),ylim=c(8,15
HI,
Try this:
source("arrayconvert.txt")
#dat2 is the data frame
dat4<-split(dat2,dat2$p)
res1<-do.call(rbind,lapply(lapply(lapply(dat4,`[`,5),function(x)
matrix(x,dim=list(20,1))),unlist))
dim(res1)
#[1] 50 20
head(res1,3)
# [,1] [,2] [,3] [,4] [,5] [,6]
Hi All,
I wonder whether it is possible to use the plot function and reduce
the scales on both +/- quadrants, i.e.
plot the data from -15. -8 and 8..15 only.
Many thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
On Thu, Sep 27, 2012 at 4:03 AM, CALEF ALEJANDRO RODRIGUEZ CUEVAS
wrote:
> Hi everybody.
>
> I use to work with monthly data, however this is my first time working with
> weekly data.
>
> I have 4 data bases with monthly data and 1 with weekly data, so I want to
> convert my weekly to monthly data
On 12-09-28 11:51 PM, Qiang Wang wrote:
In the following test, the variable x is supposed to be modified within the
function test.save.x. It works as expected util using the save command.
# initialize x.file
x <- 1
save(x, file="x.file")
rm(x) # Remove x from R_GlobalEnv
exists("x") # FALSE
# d
On Sat, Sep 29, 2012 at 7:32 AM, ilai wrote:
> On Fri, Sep 28, 2012 at 11:03 PM, Elaine Kuo wrote:
>
>> Hello
>>
>> Thanks again.
>>
>> I got the attached graph
>> Unsure why the color is still inconsistent.
>> Please kindly share with your R version.
>>
>
> It is not my version you need to worry
On Fri, Sep 28, 2012 at 11:03 PM, Elaine Kuo wrote:
> Hello
>
> Thanks again.
>
> I got the attached graph
> Unsure why the color is still inconsistent.
> Please kindly share with your R version.
>
It is not my version you need to worry about, but your own...
Post the following outputs (without t
On Sat, Sep 29, 2012 at 8:17 AM, wasss wrote:
> Hello,
> I have some problems regarding the function optim. In my simulations I
> generate negative binomial data and get estimates of the parameter using
> the likelihood and the package optim. At some point there appears the
> warning:
>
> - non-
Hi Katja,
try fitting the original model using ML (not REML) with the parameter
method = "ML":
PModell1 <-lme(sqrt(Earthwormsm.2)~
Treatment+Pflanzenfrischmasse+aBodenfeuchte+bBodenfeuchte+Gfrischmasse+Ltrockenmasseanteil+KCN+I+Eindringtiefe,
random=~1|Block/Treatment/Cluster/Patch, data=Test
Hi Z,
I guess it should work.
The apply(geno1,1,...) applies on rows 1:nrow(geno1). In the example dataset,
there were 5 rows, the result vector (tmp) is also of length 5.
A.K.
From: JiangZhengyu
To: smartpink...@yahoo.com
Sent: Friday, September 28,
Dear help community,
I'm a R-beginner and use it for my master thesis.
I've got a mixed model and want to analyse it with lme. There are a lot
Cofactors that coult be relevant. To extract the important ones I want to do
the stepAIC, but always get an error warning.
Structure of my data:
data.fram
HI,
You can try melt(). Not sure how it will perform in large datasets.
library(reshape)
n<-1e3
set.seed(1)
mat1<-matrix(rnorm(n,15),ncol=20,nrow=50)
dat1<-data.frame(mat1)
dat2<-data.frame(value=melt(dat1)[,2])
dim(dat2)
#[1] 1000 1
head(dat2,6)
# value
#1 14.37355
#2 15.18364
#3 14.164
In the following test, the variable x is supposed to be modified within the
function test.save.x. It works as expected util using the save command.
# initialize x.file
x <- 1
save(x, file="x.file")
rm(x) # Remove x from R_GlobalEnv
exists("x") # FALSE
# define test
test.save.x <- function(){
save
Hi,
Try this:
geno<-read.table(text="
P1 P2 P3 P4 P5
1 2 2 3 2
2 2 2 1 1
1 2 1 2 NA
NA 2 3 4 5
1 1 3 1 3
",sep="",header=TRUE,stringsAsFactors=FALSE)
geno1<-as.matrix(geno)
geno1[is.na(geno1)]<-0
tmp<-apply(geno1,1,function(x) ifelse((sum(x!=2)>3) & (sum(x==1)>=1) &
(sum
Hello,
I have some problems regarding the function optim. In my simulations I
generate negative binomial data and get estimates of the parameter using
the likelihood and the package optim. At some point there appears the
warning:
- non-finite finite-difference value []
I think it occurs, when t
The author of the package, Rob J. Hyndman, along with Yeasmin Khandakar,
published the paper "Automatic Time Series Forecasting: The forecast Package
for R" in the Journal of Statistical Software, July 2008, Volume 27, Issue 3.,
which is freely available online. You'll find the answer there.
Bes
On Fri, 28 Sep 2012, David Winsemius wrote:
On Sep 28, 2012, at 4:52 PM, David Winsemius wrote:
On Sep 28, 2012, at 3:16 PM, Nick Fankhauser wrote:
Hello R-Users!
I'm using a heatmap to visualize a matrix of values between -1 and 3.
How can I set the colors so that white is zero, below ze
25 matches
Mail list logo
