Hello,
> require(graphics)
> summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
> THSD <- TukeyHSD(fm1, "tension", ordered = TRUE)
> p.adj <- THSD$tension[,4]
HTH
Pascal
Le 13/12/2012 16:08, raz a écrit :
Hey,
I have this TukeyHSD output from which I would like to extract only t
Hi
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of dada
> Sent: Thursday, December 13, 2012 12:41 AM
> To: r-help@r-project.org
> Subject: [R] neural net
>
> Hi
> I would like to do neural netowrk analysis on my data. It look
Hey,
I have this TukeyHSD output from which I would like to extract only the
P-values (p adj, last number).
The problem is that the test output is a character list.
How can I "break" this sentence to separate the Pv?
Tukey multiple comparisons of means
95% family-wise confidence level
Fit:
I apologize if my "scream" went out to the whole list serve, I just
intended to email the moderator. I have emailed a few times, and am still
on the list. I tried the link as well, but perhaps I am doing something
wrong.
On Wed, Dec 12, 2012 at 11:18 PM, Pascal Oettli wrote:
> Hello,
>
> First,
As the posting guide asked, please ask the package maintainer directly
(Cc:ed here).
On 12/12/2012 19:34, Aspro wrote:
Hi,
I'm currently using EMA package to make clustering and heatmaps.
The online doc concerning the package gives the following example code:
data(marty)
c<-clustering(marty,
Hello,
On the web page, search "To unsubscribe from R-help, get a password
reminder, or change your subscription options enter your subscription
email address:"
Pascal
Le 13/12/2012 13:21, Elizabeth Fuller Bettini a écrit :
I apologize if my "scream" went out to the whole list serve, I jus
Hello,
First, no need to scream.
Second, David Winsemius already responded to you, last Dec.,9 (topic was
"[R] Mean-Centering Question"). Here is a hint:
https://stat.ethz.ch/mailman/listinfo/r-help (bottom of the page)
Pascal
Le 13/12/2012 13:01, Elizabeth Fuller Bettini a écrit :
PLEASE RE
On Dec 12, 2012, at 6:34 PM, Jichun Chan wrote:
Hi,
Does anyone know how to write a command to generate time-to-event
data for Cox's regression?
There are examples in the rms package for the cph function. I also see
example in the help pages for the survival package.
--
David Winsemi
download.file( url = filename , destfile = location.to.save.locally , mode
= 'wb' )
On Wed, Dec 12, 2012 at 2:48 PM, Alemu Tadesse wrote:
> I am trying to download the tar files on the website below
>
> filename<-"
> http://rredc.nrel.gov/solar/old_data/nsrdb/1991-2010/SolarAnywhere/x.tar";
> wh
Dear all
I am now running a MCMC iteration in the R program. But it is always
stucked in some loop. This cause big problems for my research. So I want to
know whether we can skip the current dataset and move to next simulated
data when the iteration is stucked? Alternatively, can the MCMC chain ski
Hi
I would like to do neural netowrk analysis on my data. It look like this:
drugparam1 param2 param3 param4 param5 class
A 111 15 125 40 0.5 1
B 347 13 280 55 3 2
C 335 9 119 89 -40 1
D 477
Dear officer
I have a question concerning running R when I am doing my research. Can you
help me to figure that out?
I am now running a MCMC iteration in the R program. But it is always
stucked in some loop. This cause big problems for my research. So I want to
know whether we can skip the current
Hi Arun,
I have a question about choosing a subset of data.
I have a matrix of 5000 in 3500.
I want to choose specific variables such as proc1 to proc1000 and Lab1 to
Lab1600 and put it into a new matrix to run some correlation analysis.
Since I have 3500 variables, I don't know from what column
Hi there. At first glance it sounded to me as an obvious "no-no" question.
But, for some reason, I ran some trials and results looked pretty
intriguing.
So, I checked 14 genotypes (8 plants from each randomly chosen in the
field) on 4 different dates and measured them under 2 different
temperatur
Thank you so much.
It worked very well:)
Best,Farnoosh Sheikhi
Cc: R help
Sent: Wednesday, December 12, 2012 12:59 PM
Subject: Re: Subset of Data
Hi Farnoosh,
Try this:
set.seed(151)
mat1<-matrix(sample(1:400,100,replace=TRUE),ncol=20)
set.seed(15)
colnam
The mailing list ate the attachments, so here they are again.
R code
https://gist.github.com/4270628
Log
http://pastebin.com/0e49CTsL
Andrew
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Andrew Ziem
Sent: Wednesday, December
Hi Farnoosh,
Try this:
set.seed(151)
mat1<-matrix(sample(1:400,100,replace=TRUE),ncol=20)
set.seed(15)
colnames(mat1)<-paste(sample(c("proc","Lab","other"),20,replace=TRUE),sample(1:45,20,replace=FALSE),sep="")
mat1[,grepl("proc|Lab",colnames(mat1))]
# Lab37 proc35 Lab3 Lab24 proc6 proc36 Lab
Hi,
Does anyone know how to write a command to generate time-to-event data for
Cox's regression?
Scomet
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
Hi,
I'm currently using EMA package to make clustering and heatmaps.
The online doc concerning the package gives the following example code:
data(marty)
c<-clustering(marty, metric="pearson", method="ward")
clustering.plot(c, title="Hierarchical Clustering\nPearson-Ward")
which is working perfec
Hello,
I cannot reproduce the error:
> library(ggplot2)
Attaching package: 'ggplot2'
The following object(s) are masked from 'package:latticeExtra':
layer
> library(cshapes)
Loading required package: maptools
Loading required package: foreign
Checking rgeos availability: TRUE
Loading req
Dear parallel users and developers,
I might have encountered a bug in the function 'mclapply' of package
'parallel'. I construct a matrix using the same input data and code with a
single difference: Once I use mclapply and the other time lapply.
Shockingly the result is NOT the same.
To evaluate
Thank you, Rui! This is incredibly helpful -- anonymous functions
are new to me, and I appreciate being shown how useful they are.
Best regards,
David
On Wed, Dec 12, 2012 at 10:12 AM, Rui Barradas wrote:
> Hello,
>
> As for the first question try
>
> scoreset <- lapply(pcl, function(x) x$scor
soon yi ymail.com> writes:
>
> Hi
>
> I am using ggplot to overlay two regression lines on a scatter plot each
> corresponding to a treatment group.
>
> The default plot gives a different slope for each treatment group. However,
> in some cases i want the lines to be parallel -ie no significan
Something like
tbp <- as.numeric( gsub("(.*-)?([0-9]+)\\+?","\\2", Trade_Price_Band) )
would get you the numbers you want. The '500+' is then treated as 500, though.
One way of getting round that might be
tbp <- as.numeric( gsub("(.*-)?([0-9]+)","\\2", Trade_Price_Band) )
then
tbp[is.na(t
I found that the optim() function does not always reach the real minimum. Is it
because the solution is trapped at a local minimum?
Thanks!
Ken
On Dec 12, 2012, at 2:17 PM, Jeff Newmiller wrote:
> ... origin pro?
>
> Then why are you here?
>
> It is not clear from your message that this ha
One solution that does not require matrix multiplication:
Remember that the steady state vector is in the nullspace of I - T.
Therefore:
require(MASS)
n1 <- Null(t(diag(nrow(T)) - T))
n1 / sum(n1)
On Wed, Dec 12, 2012 at 2:19 AM, annek wrote:
> Hi,
> I have a transition matrix T for which I
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of annek
> Sent: Tuesday, December 11, 2012 11:19 PM
> To: r-help@r-project.org
> Subject: [R] Matrix multiplication
>
> Hi,
> I have a transition matrix T for which I want to find
... origin pro?
Then why are you here?
It is not clear from your message that this has anything to do with R.
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.
I extracting splitting rules from Greg Ridgeway's GBM 1.6-3.2 in R 2.15.2, so I
can run classification in a production system outside of R. I have it working
and verified for a dummy data set with all variable types (numeric, factor,
ordered) and missing values, but in the titanic survivors dat
Hi,
#A lengthy solution for the simple problem.
#dat1 is the dataset
res<-unlist(lapply(strsplit(dat1[,1],split="-"),function(x) if(length(x)==1)
rep(x,2) else x))
res[seq(from=2,to=length(res),by=2)]
#or
gsub("[-]","",unlist(regmatches(dat1[,1],gregexpr("-.*|.*\\+",dat1[,1]
A.K.
- O
Hi,
Try ?merge() or ?join()
dat1<-read.table(text="
ID Null
1 NA 1
2 86.1 2
3 88.9 3
4 100.0 4
5 80.6 5
6 97.2 6
",sep="",header=TRUE)
dat2<-read.table(text="
SCD.TD cluster
1 1 4
2 1 17
3 2 4
4 2
A friend told me about some Big 12.12.12 announcement coming out of
XLSolutions about R-PLUS and they won't respond to my email.
Anyone knows if they will having a free version of R-PLUS ?
Thanks
Eugene
[[alternative HTML version deleted]]
Happy Holidays !
XLSolutions January-February 2013 R courses schedule is now available
online at
9
USA cities for with 13 new courses: *** Suggest a future course
date/city
(1) R-PLUS: A Point-and-Click Approach to R
(2) S-PLUS / R : Programming Essentials.
(3) R/S+ Fundamentals and Programming T
Hello,
As for the first question try
scoreset <- lapply(pcl, function(x) x$scores[, 1])
do.call(cbind, scoreset)
As for the second question, you want to know which columns in 'datasets'
have NA's?
colidx <- apply(datasets, 2, function(x) any(is.na(x)))
datasets[, colidx] # These have NA's
Hi
I am using ggplot to overlay two regression lines on a scatter plot each
corresponding to a treatment group.
The default plot gives a different slope for each treatment group. However,
in some cases i want the lines to be parallel -ie no significant
interaction.
My code:
ggplot(data=df,X,Y,
Hi All:
I have a dataset that when estimated as a dynamic state-space model, appears
that the trend term would be better modeled with long-tails, such as a
t-distribution, rather than Gaussian. I have looked at the manuals for the
packages dlm, KFAS and sspir, and if I can estimate such a mode
Sorry, I just realized I didn't send the message below in plain text.
-David Romano
On Wed, Dec 12, 2012 at 9:14 AM, David Romano wrote:
>
> Hi everyone,
>
> Suppose I have a 3D array of datasets, where say dimension 1 corresponds
> to cases, dimension 2 to datasets, and dimension 3 to observatio
Hi everyone,
Suppose I have a 3D array of datasets, where say dimension 1 corresponds to
cases, dimension 2 to datasets, and dimension 3 to observations within a
dataset. As an example, suppose I do the following:
> x <- sample(1:20, 48, replace=TRUE)
> datasets <- array(x, dim=c(4,3,2))
Here,
Dear All
the problem I have is as follows.
I have attribute data in a number of data.frames identified 1:58 (the
column Null is just there to stop it becoming a list, which caused me
trouble, the data I am interested in is in column 1, although other
data frames have multiple categories/cols). The
Another option (better or worse probably depends on many things) is to use:
b <- head(a, -1)
On Wed, Dec 12, 2012 at 1:14 AM, e-letter wrote:
> On 12/12/2012, Daniel Nordlund wrote:
> >> -Original Message-
> >> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org
>
Chad-
The plot you hope to create is meaningful when the variable temperature is
treated as a continuous variable. Below is some code that treats
temperature as a continuous variable and creates the plot. Note that
temperature enters the model "linearly", and you may want to explore other
possib
Dear Chad,
Did you post your entire dataset? If so:
1) your model is too complex for the amount of data you have. See the quotes
below...
2) There is complete separation, leading to large parameter estimates and fits
very close to 0 and 1 (in terms of probabilities)
3) You fit temperature as a
On 12.12.2012 00:05, Marian Talbert wrote:
I'm trying to use Rprof() to identify bottlenecks and speed up a particullary
slow section of code which reads in a portion of a tif file and compares
each of the values to values of predictors used for model fitting. I've
written up an example that a
On 12.12.2012 16:32, John Kerpel wrote:
Folks:
I keep getting the following error message (I'm on Windows 7, R-2.15.2, and
tried a reboot...). Thx!
Either you do not have write permission on that directory or you have
the package loaded already,
Uwe Ligges
John
install.packages("car
dev.cur()
On 12 December 2012 10:45, Williams wrote:
>
> __
> 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, minima
On Dec 12, 2012, at 7:32 AM, John Kerpel wrote:
Folks:
I keep getting the following error message (I'm on Windows 7,
R-2.15.2, and
tried a reboot...). Thx!
Actually it's only a warning. Those are not the same in R.
John
install.packages("caret")Installing package(s) into C:/Program
Folks:
I keep getting the following error message (I'm on Windows 7, R-2.15.2, and
tried a reboot...). Thx!
John
> install.packages("caret")Installing package(s) into C:/Program
> Files/R/R-2.15.2/library
(as lib is unspecified)trying URL
'http://streaming.stat.iastate.edu/CRAN/bin/window
Dear R Wizards,
After much frustration and days of confusion I have finally broken down and
am asking for help, which I dont like doing, but I just cant figure this
one out on my own. Ive conducted a laboratory experiment testing the
effects of temperature and salinity on whether or not a biol
Hello,
I need plot map using ggplot()
I use such code:
library("ggplot2")
library(cshapes)
cshp.data <- cshp()
map <- cshp.data[cshp.data$COWCODE==369,]
map_mp<- list(
geom_polygon(aes(long, lat, group = group), data = map, fill = "grey70",
colour = "grey60", inherit.aes = FAL
__
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, self-contained, reproducible code.
If your final goal is to write the csv files, and only the csv files, then
then this (not tested) should do it.
for (i in unique(appended$dp) ) {
tmp <- subset(appended, dp==i & sampled==0)
write.table(tmp,
file= file.path('output', paste0('set',i,'.csv')),
sep=',',
row.names=FALSE
Hi Fredrik,
Here is an alternative:
df <- data.frame(cat=LETTERS[1:4],num=rnorm(4))
df$g <- df$cat
levels(df$g) <- c("G1", "G1", "G2", "G2")
ggplot(df, aes(x=cat, y=num)) + geom_point() +
geom_line(aes(group=g))
Best,
Ista
On Wed, Dec 12, 2012 at 3:28 AM, Fredrik Karlsson wrote:
> Dear list
On Tue, Dec 4, 2012 at 2:12 PM, David Romano wrote:
>
>
> On Tue, Dec 4, 2012 at 10:22 AM, Duncan Murdoch
> wrote:
>
>> On 04/12/2012 12:54 PM, David Romano wrote:
>>
>>> Hi everyone,
>>>
>>> I normally include a call to browser() as I'm working out the kinks in my
>>> scripts, and I am always a
assuming you want more than just the first four rows
tpb <- read.table( text = " Trade_Price_Band x
1 0-30 0.6237240
2 101-150 0.6743857
3 151-200 0.6778513
4 201-300 0.6640293
5 301-400 0.6630991
6 31-50 0.6314547
7
Hi Steven,
Not sure if you want to understand regular expressions in general or just the
solution to your particular problem. If it's the former and you'd be willing to
read a book on the subject, I'd recommend "Mastering Regular Expressions" by
Jeffrey Friedl. I'm about halfway through now, an
Le mardi 11 décembre 2012 à 23:39 +0100, Richard Zijdeman a écrit :
> Dear Milan, please see my results inline
>
> On 11 Dec 2012, at 16:58, Milan Bouchet-Valat wrote:
>
> > Le mardi 11 décembre 2012 à 16:41 +0100, Richard Zijdeman a écrit :
> >> Dear Milan,
> >>
> >> thank you for kind suggest
Hi,
I met this problem.
Trade_Price_Band x
1 0-30 0.6237240
2 101-150 0.6743857
3 151-200 0.6778513
4 201-300 0.6640293
5 301-400 0.6630991
6 31-50 0.6314547
7 401-500 0.6776249
8 500+ 0.6557705
On 12/12/2012 09:10, peter dalgaard wrote:
On Dec 12, 2012, at 01:48 , arun wrote:
Hi,
Try this:
pdf("broke.pdf")
palette(rainbow(6))
plot(x=x1,y=y1,col=y1,xlim=c(-10,20))
palette(heat.colors(6))
points(x=x2,y=y2,col=y2)
dev.off()
A.K.
Yep. Notice though that this is not specific
Would this accomplish what you're trying to do?
http://www.r-bloggers.com/define-intermediate-color-steps-for-colorramppalette/
Or some combination of ?colorRampPalette and ?cut or classIntervals() from
package classInt.
On Wed, Dec 12, 2012 at 7:00 AM, jeff6868 wrote:
> Hello everybody,
>
>
Hello
Package 'RColorBrewer' + function 'colorRampPalette'
or
?two.colors
HTH
Pascal
Le 12/12/12 21:00, jeff6868 a écrit :
Hello everybody,
I'm trying to create my own color palette on R, in order to interpolate some
different temperature data on different maps (daily means, seasonal
means
Hello everybody,
I'm trying to create my own color palette on R, in order to interpolate some
different temperature data on different maps (daily means, seasonal
means,...).
I would like to create a color palette which works for each map, so I need a
color palette between -40 and +40°C. Sometimes
On Dec 12, 2012, at 08:19 , annek wrote:
> Hi,
> I have a transition matrix T for which I want to find the steady state matrix
> for. This could be approximated by taking T^n , for large n.
>
> T= [ 0.8797 0.0382 0.0527 0.0008
> 0.02120.8002 0.0041 0.0143
> 0.09810
You need to install the party package and then load the package before you can
use its functions:
install.packages("party")
library(party)
On 11.12.2012, at 16:56, Peterson, Kim - DNR wrote:
> Greetings! I'm trying to use function varimpAUC in the party package
> (party_1.0-3 released Septemb
On Dec 12, 2012, at 01:48 , arun wrote:
> Hi,
> Try this:
> pdf("broke.pdf")
> palette(rainbow(6))
> plot(x=x1,y=y1,col=y1,xlim=c(-10,20))
> palette(heat.colors(6))
> points(x=x2,y=y2,col=y2)
> dev.off()
> A.K.
Yep. Notice though that this is not specific to dev.copy2pdf, the same effect
Try this
install.packages("expm")
library(expm)
T = matrix(c( 0.8797 , 0.0382 , 0.0527 , 0.0008,
0.0212 , 0.8002 , 0.0041 , 0.0143,
0.0981 , 0.0273 , 0.8802 , 0.0527,
0.0010 , 0.1343 , 0.0630 , 0.9322),
nrow = 4, ncol = 4,byr
Dear list,
I've been using plotmeans {gplots} a lot before, and found the
"connect" argument to be quite useful.
I've moved to ggplot for several reasons, but would still like to
connect lines conditionally, somehow.
Is it possible to do?
Small example:
df <- data.frame(cat=LETTERS[1:4],num=rnor
On 12/12/2012, Daniel Nordlund wrote:
>> -Original Message-
>> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
>> On Behalf Of e-letter
>> Sent: Tuesday, December 11, 2012 11:45 PM
>> To: r-help@r-project.org
>> Subject: [R] remove last row of a data frame
>>
>> Re
67 matches
Mail list logo
