Hello Ista,
I would like to thank you for your reply!
Your code is indeed an improvement.
This is my code now and there are two things still missing.
DataToPlot<-matrix(data=runif(9),nrow=3,dimnames=list(seq(1,3),seq(4,6)))
require(reshape)
require(ggplot2)
require(raster)
cols<-colours()
cols
Hi All,
Let's say someone uses an R package (prefer to not disclose which) that comes
with two different distributions,the commercial and non-commercial license. Now
the non commercial license has conditions,bit it never runs out. So if this
person knows they should be using the commercial licen
Q: how to sum columns in R?
A: use a web search engine
Le 21/02/2013 16:33, nandita srivastava a écrit :
how to sum columns in R?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
Le jeudi 21 février 2013 à 13:16 +0400, Lawr Eskin a écrit :
> Hello dear R-help mailing list.
>
>
> Looks like the same issue in Russian:
>
>
>
> library(RCurl)
>
> library(XML)
>
> u = " http://www.cian.ru/cat.php?deal_type=2&obl_id=1&room1=1";
>
> a = getURL(u)
>
> a # Here - the Russia
Hi
?ncol
Regards
Petr
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of nandita srivastava
> Sent: Thursday, February 21, 2013 8:34 AM
> To: r-help@r-project.org
> Subject: [R] (no subject)
>
> how to sum columns in R?
>
>
Hi
I'm trying to measure the effect of a policy intervention (Box and Tiao, 1975).
This query has to do with the coding of the model rather than with the
particulars of my dataset, so I'm not providing the actual dataset (or a
simulated one) in this case, apart from some general description.
Th
On Thu, Feb 21, 2013 at 9:52 AM, Pancho Mulongeni
wrote:
> Hi All,
> Let's say someone uses an R package (prefer to not disclose which) that comes
> with two different distributions,the commercial and non-commercial license.
> Now the non commercial license has conditions,bit it never runs out.
Hello,
You've got a bug, Arun, the parenthesis are not right in
((1/sd1)*sqrt(2*pi)). It should be (1/(sd1*sqrt(2*pi))).
The requsted values can also be computed using ?dnorm.
d1 <- (1/(sd1*sqrt(2*pi)))*exp(-((x1-mean1)^2)/(2*(sd1^2)))
d2 <- dnorm(x1, mean = mean1, sd = sd1)
all.equal(d1, d
?colSums
MW
On Feb 21, 2013, at 11:26 AM, PIKAL Petr wrote:
> Hi
>
> ?ncol
>
> Regards
> Petr
>
>
>> -Original Message-
>> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
>> project.org] On Behalf Of nandita srivastava
>> Sent: Thursday, February 21, 2013 8:34 AM
>> To:
Hi All
Having been pointed the use of events and roots in deSolve, I was able to
implement the Izchikevich model of spiking neurons. However, I'm not too
sure of defining the event. The deSolve documentation says:
An event is triggered when the ball hits the ground (height = 0) Then
velocity (y2)
Hello,
http://lmgtfy.com/?q=sum+columns+in+R
The first hit looks promising.
Hope this helps,
Rui Barradas
Em 21-02-2013 07:33, nandita srivastava escreveu:
how to sum columns in R?
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
Hi,
On Thu, Feb 21, 2013 at 4:12 AM, Alaios wrote:
> Hello Ista,
> I would like to thank you for your reply!
>
> Your code is indeed an improvement.
>
> This is my code now and there are two things still missing.
>
> DataToPlot<-matrix(data=runif(9),nrow=3,dimnames=list(seq(1,3),seq(4,6)))
> requ
Let me try to explain this then.
I have numbers like
DataToPlot
4 5 6
1 0.4454995 0.4462009 0.4286807
2 0.3761550 0.5423205 0.6500785
3 0.3779496 0.4671437 0.1799601
and I want to have a color legend with each color that is mapped to the
following printed values
"1"
http://www.fischeredv.de/zxripc/3bg8jc&mc9en1iqo46=lgyc8dv37ec694oxxqf2enjlgwsr70
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.
Dear R experts,
I have two matrix (seq & mat) & I want to retrieve in a new matrix all the
numbers from mat that =1 (corresponding to the same row/ column position) in
seq, or all the numbers in mat that =-1 in seq. - Replace all the numbers with
NA if it's not 1/-1 in seq. There ar
Le jeudi 21 février 2013 à 18:31 +0400, Lawr Eskin a écrit :
> Hi Milan,
>
> a <- getURL(con, .encoding = "UTF-8")
> Encoding(a)
> > [1] "UTF-8"
> a # Here - the UTF-8 codes looks like fine.
> htmlParse(a, encoding = "UTF-8") ###again same encoding issue
And what if you try this:
a2 <- htmlPars
Hi,
I guess this is what you wanted:
seq1<-matrix(c(1,-1,0,1,1,-1,0,0,-1,1,1,NA),3,4)
seq1
# [,1] [,2] [,3] [,4]
#[1,] 1 1 0 1
#[2,] -1 1 0 1
#[3,] 0 -1 -1 NA
set.seed(15)
mat1<-matrix(rnorm(12),3)
mat1
# [,1] [,2] [,3] [,4]
#[1,]
Hi Zhengyu
Surely not the most elegant piece of coding, but this should do the trick:
> m=rep(NA,nrow(seq)*ncol(seq))
> mat1=matrix(m,3)
> for(i in 1:nrow(seq)){
+ for(j in 1:ncol(seq)){
+ ifelse((seq[i,j]==1 | seq[i,j]==-1),
+ mat1[i,j]<-mat[i,j],mat1[i,j]<-NA)
+ }
+ }
> mat1
e.g. Wi
Hi,
On Thu, Feb 21, 2013 at 9:21 AM, Alaios wrote:
> Let me try to explain this then.
>
> I have numbers like
>
> DataToPlot
> 4 5 6
> 1 0.4454995 0.4462009 0.4286807
> 2 0.3761550 0.5423205 0.6500785
> 3 0.3779496 0.4671437 0.1799601
>
> and I want to have a color legen
Dear all
I prepared some image processing routine which depended on package ReadImages.
Basically I imported image to R, changed rgb to grayscale, normalised an image
made some selections and changed to image by imagematrix.
Recently I learned that this package was removed from CRAN and can not
On 21-02-2013, at 15:39, JiangZhengyu wrote:
> Dear R experts,
>
> I have two matrix (seq & mat) & I want to retrieve in a new matrix all the
> numbers from mat that =1 (corresponding to the same row/ column position) in
> seq, or all the numbers in mat that =-1 in seq. - Replace all the numb
Dear Berend, Mark, Jose, Arun, Great! Thank you so much for all your replies
with different codings. They all work well except one - because NAs in matrix A
need to take care of. Best,Zhengyu
> Subject: Re: [R] Ask for help: find corresponding elements between matrix
> From: b...@xs4all.nl
>
Hi Bert,
Thanks for drawing my attention to "simplify" argument and for the examples. I
understand know.
Thanks.
Dan
-Original Message-
From: Bert Gunter [mailto:gunter.ber...@gene.com]
Sent: Wednesday, February 20, 2013 4:25 PM
To: Lopez, Dan
Cc: R help (r-help@r-project.org)
Subject
Hi: bierens has a paper on modelling beer sales in the netherland ( I'm
pretty sure it's on the net.
if not, I have a copy somewhere I think ) using an ARIMAX. why don't you
take his paper and his data and see if you get the same estimates using R.
That's one way if you'll know if you're doing the
> scs2<-data.frame(lapply(scs2, factor))
Calling data.frame() on the output of lapply() can result in changing column
names
and will drop attributes that the input data.frame may have had. I prefer to
modify
the original data.frame instead of making a new one from scratch to avoid these
proble
to anyone who was mistakenly looking for the armax beer paper using Bierens
as the author. I made a mistake and apologize because the author was
franses and the paper is at the link below:
http://repub.eur.nl/res/pub/2112/eur_franses_MR1.pdf
On Thu, Feb 21, 2013 at 11:23 AM, Paul Bernal wrote:
Thanks Mark.
I had looked up "herman bierens beer netherlands" but nothing came up, except
the paper by Philip Hans Franses which I didn't know so I'll have a look at it.
Thanks for letting others note about this.
I had toyed with the idea of contacting the authors of the TSA package directly
Also something like
## generate zero-inflated Poisson data
n <- 100
x1 <- rnorm(n)
x2 <- rnorm(n)
ind <- rbinom(n,1,1/(1+exp(-1-0.1*x1)))
y <- ifelse(ind,rpois(n,lambda=exp(1+0.2*x2)),0)
There is at least zeroinfl in the CRAN package pscl and fmr in Lindsey's
non-CRAN package gnlm for fitting z
Hi,
Try this:
You can save it as either individual jpeg files or as a single pdf document
containing all the figures.
set.seed(54)
res<- matrix(sample(20:60,25,replace=TRUE),ncol=5)
res1<- res[,-1]
b<- res[,1]
#jpeg
mypath<- file.path("/home/arunksa111/Trial",
paste("myplot_",names1,".jpg",sep=
Hi Milan!
> Encoding(a)[1] "unknown"
2013/2/21 Milan Bouchet-Valat
>
>> Le jeudi 21 février 2013 à 13:16 +0400, Lawr Eskin a écrit :
>> > Hello dear R-help mailing list.
>> >
>> >
>> > Looks like the same issue in Russian:
>> >
>> >
>> >
>> > library(RCurl)
>> >
>> > library(XML)
>> >
>> >
Anybody knows about a library for functional median polish?
Thanks
--
Francisco Montes
Dpt. d'Estadística i I. O.
Universitat de València
E-46100-Burjassot
SPAIN
Tlf: +34 963544306
Fax: +34 963544735
http://www.uv.es/~montes
I have tried removing a parenthesis to the right of poly and I get this
error:
Error in bboxx(bbox(as.points(poly)), nsim = 29, s = seq(5, 80, 5)) :
unused argument(s) (nsim = 29, s = seq(5, 80, 5))
It looks like you have one comma for poly, one for as.points, and one for
bbox. The original e
iconv trued before in various try, same issue and result with encoding =
unknown
now try sub - same issue
2013/2/21 Milan Bouchet-Valat
> Le jeudi 21 février 2013 à 18:31 +0400, Lawr Eskin a écrit :
> > Hi Milan,
> >
> > a <- getURL(con, .encoding = "UTF-8")
> > Encoding(a)
> > > [1] "UTF-8"
>
Hi Milan,
a <- getURL(con, .encoding = "UTF-8")
Encoding(a)
> [1] "UTF-8"
a # Here - the UTF-8 codes looks like fine.
htmlParse(a, encoding = "UTF-8") ###again same encoding issue
>>why didn't getURL() detect and set a's encoding correctly?
I think there are page issue because another sites works
Hi all,
I am using package lavaan and have created a structural equation model with
two exogenous and seven endogenous variables with the following
relationships
#specify the model
m1 = ' # regressions
D ~ ma + hs + b4 + b5 + b15 + b16
ma ~ hs + b4 + b5 + b15 + b16
I only thought one was missing since it kept giving me error that it wasn't
set. But all seems to working OK so I guess I don't need to create anything.
Thank you for your help.
On Wed, Feb 20, 2013 at 7:28 PM, S Ellison wrote:
>
>
> > Should I create an output directory somewhere either in Pro
Thank you. I ended up converting the excel file to a csv file and using
mytable<-read.csv(myfile, header=T)
and it worked so I abandoned using the tab delimited txt file.
On Wed, Feb 20, 2013 at 9:29 PM, Jim Lemon wrote:
> On 02/21/2013 07:10 AM, Joanna Papakonstantinou wrote:
>
>> Some of the n
Dear all,
I have a data frame, which looks like this:
Subject | Block | Trial | Feature1 | Feature2
1 | 1 | 1 | ... | ...
1 | 1 | 2 | ... | ...
1 | 2 | 1 | ... | ...
1 | 2 | 2 | ... | ...
1 | 3 | 1 | ... | ...
...| ...| ...| ... | ...
Can I remove the "Trial" column by averaging all the row
Great point, William. I agree your approach is the one to take to preserve
attributes.
Thanks for following up.
From: William Dunlap
elp (r-help@r-project.org)"
Sent: Thursday, February 21, 2013 11:32 AM
Subject: RE: [R] Having trouble converting a dataf
On Thu, Feb 21, 2013 at 10:46 AM, Joanna Papakonstantinou
wrote:
> I only thought one was missing since it kept giving me error that it wasn't
> set. But all seems to working OK so I guess I don't need to create anything.
>
No. The first two gave a response of 1 which means they failed.
Hello,
Try package sos first:
library(sos)
findFn('functional median') # 'functional median polish' returns nothing
Package fda.usc is first.
findFn('functional anova') # also gets package fda.usc
Hope this helps,
Rui Barradas
Em 21-02-2013 18:15, francisco.mon...@uv.es escreveu:
Anybo
As a follow-up, I am able to achieve what I'd like to do in terms of a set
range of values for amplitude using filled.contour() and specifying the
zlim argument - however, in the end I'd still like to have the
possibilities for layering that ggplot offers available to me; so, please
let me know if
...
The OP might also wish to look at the "Robust" task view on CRAN for
other approaches to dealing with "discrepant" data. I leave it to to
him to define what this means and to determine whether anything there
serves his needs.
Cheers,
Bert
On Thu, Feb 21, 2013 at 12:54 PM, Rui Barradas wrote:
Hi Ben,
How about
lim <- with(input, {
c(min(amplitude),
max(amplitude))})
p <- ggplot(sub1, aes(x=x,y=y,fill=amplitude))
p +
geom_tile() +
scale_fill_gradientn("Custom
Colours",
colours=c("purple",
"blue",
"gre
Hi,
May be this helps:
dat1<- read.table(text="
Subject Block Trial Feature1 Feature2
1 1 1 48 40
1 1 2 62 18
1 2 1 34 43
1 2 2 51 34
1 3 1 64 14
",sep="",header=TRUE)
res1<-do.call(rbind,lapply(split(dat1,dat1$Block),function(x)
data.frame(unique(x[,1:2]),t(colMeans(x[
another option that i think is easy & intuitive :)
library(sqldf)
mydata <- sqldf( "select Subject , Block , avg( Feature1 ) as Feature1 ,
avg( Feature2 ) as Feature2 , [..keep going..] , avg( Feature10 ) as
Feature10 from yourdata group by Subject , Block" )
On Thu, Feb 21, 2013 at 4:45 PM,
Excellent! Now I see that documentation in discrete_scale(). Thanks.
*Ben Caldwell*
Graduate Fellow
University of California, Berkeley
130 Mulford Hall #3114
Berkeley, CA 94720
Office 223 Mulford Hall
(510)859-3358
On Thu, Feb 21, 2013 at 1:26 PM, Ista Zahn wrote:
> Hi Ben,
>
> How about
>
>
Many find the functions in the plyr package more convenient to use than the
do.call(rbind, lapply(split(...),...) business:
> library(plyr)
> ddply(dat1, .(Subject,Block), summarize, MeanFeature1=mean(Feature1),
MeanFeature2=mean(Feature2))
Subject Block MeanFeature1 MeanFeature2
1
Hi Bill,
Great info.
The problem is what was originally given to me looks like DPUT1 below (random
sample of 25).
This is the only format they can give me this in and the data already looks
molten. So I applied reshape2::dcast which resulted in a dataframe made of
character vectors; except fo
#So I tried I tried this which had no effect
keepcols<- grepl("Q1_",names(scs.c2))
levels(scs.c2[,keepcols])<-list(NoResp="",NotImportant="not
important",SomewhatImpt="somewhat
important",Important="important",VeryImpt="very important")
#then this which also failed. It coerced a
On Thu, Feb 21, 2013 at 5:22 PM, Joanna Papakonstantinou
wrote:
> Below is what happened when I ran
>> testInstalledPackages(scope="base",errorsAreFatal=FALSE) Testing examples
>> for package ‘base’ comparing ‘base-Ex.Rout’ to ‘base-Ex.Rout.prev’ ... files
>> differ in number of lines: Testing exa
(apologies for not having edited the to/cc list)
__
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-conta
the world's most minor typo on
https://github.com/hadley/devtools/wiki/Rcpp
optimisiation
helpful page, btw! 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-pr
Below is what happened when I ran
> testInstalledPackages(scope="base",errorsAreFatal=FALSE) Testing examples
for package base comparing base-Ex.Rout to base-Ex.Rout.prev
...files differ in number of lines: Testing
examples for package tools comparing tools-Ex.Rout to
tools-Ex.Rout.save
I want to analyze binary, multinomial, and count outcomes (as well as
the occasional continuous one) for clustered data.
The more I search the less I know, and so I'm hoping the list can
provide me some guidance about which of the many alternatives to choose.
The nlme package seemed the obvious
One more link to look at
http://glmm.wikidot.com/faq
This is the r-sig-mixed-models FAQ.
On Thu, Feb 21, 2013 at 8:53 PM, Ross Boylan wrote:
> I want to analyze binary, multinomial, and count outcomes (as well as the
> occasional continuous one) for clustered data.
> The more I search the less
Dear All,
I've 2 arrays A and B:
dim(A) = 100, 10, 1000
dim(B) = 100, 20, 900
I know there are 5 columns of values common to both arrays. I want to ask
how to merge the 2 arrays. Thanks in advance!
Best Regards,
Ray
[[alternative HTML version deleted]]
I think this specification is insufficient to respond accurately to. Please
make a reproducible subset of your data (or simulated data) and provide it in
dput form, and describe your desired result data set more clearly.
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducib
Thanks, Jeff.
Here is a simplified hypothetical sample (sorry for the clumsy code):
A1 <- matrix(1:5, nrow=5, ncol=1)
A2 <- matrix(6:10, nrow=5, ncol=1)
A3 <- matrix(11:15, nrow=5, ncol=1)
A4 <- matrix(16:20, nrow=5, ncol=1)
A5 <- matrix(21:25, nrow=5, ncol=1)
A6 <- matrix(26:30, nrow=5, ncol=1)
B
Look at the harvestr package for one way to control multiple parallel
simulations and make sure that they have different seeds.
On Wed, Feb 20, 2013 at 4:24 PM, C W wrote:
> Thanks, Greg. I think you are right. I did simulation one right after
> the other, less than 2 seconds.
>
> But still,
Huh? This answer seems to have absolutely no bearing on the question.
cheers,
Rolf Turner
On 02/21/2013 11:26 PM, PIKAL Petr wrote:
Hi
?ncol
Regards
Petr
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
project.org] On Behalf Of nandita
Some possibilities: The EBImage package on Bioconductor; the jpeg and png
packages read jpeg and png images.
On Thu, Feb 21, 2013 at 8:26 AM, PIKAL Petr wrote:
> Dear all
>
> I prepared some image processing routine which depended on package
> ReadImages. Basically I imported image to R, chang
On 22/02/2013 05:44, Greg Snow wrote:
Some possibilities: The EBImage package on Bioconductor; the jpeg and png
packages read jpeg and png images.
Almost everyone using ReadImages was using it to read TIFF images, which
package tiff can do. See the 'Data Import/Export Manual',
http://cran.r
63 matches
Mail list logo
