Hi, a solution could be:
# example matrix a:
a <- matrix(1:100, 10, 10)
a[, 1] <- (sample(c("aa","bb" , "ab"), 10, rep=TRUE))
a <- a[order(a[, 1]), ] # order the matrix by row = 1
#subsetting a:
lev <- levels(as.factor(a[, 1]))
subs <- list()
for(i in 1:length(lev)) {
subs[[i]] <- a[a[, 1] %
agroscope.admin.ch> writes:
>
> 2.50% 97.50%
> intercept 29787966
> glands13143611
> damage169 6144
> treatment L1 778 6703
> treatment L4 -3899 5125
> Length-1817 1828
>
plotrix::plotCI(..., err="x")
??
__
The split function does essentially this, but puts the results into a list
rather than using the dangerous and messy assign function. The overall
syntax is simpler as well.
On Thu, Feb 12, 2015 at 3:14 AM, Jim Lemon wrote:
> Hi Samarvir,
> Assuming that you want to generate a separate data fram
On 12.02.2015 16:05, arnabkm2007 wrote:
Hi,
I am trying to run the following model for OpenBUGS and want to use
R2OpenBUGS package. The model specifies weibull distribution for censored
data.
weibull.model <- function()
{
for(i in 1:n)
{
exp.alpha[i] ~ dgamma(a.alpha,
On 12.02.2015 19:07, Li, Yan wrote:
I meant R3.0.3. I need this package working in R3.0.3.
I have 32bit and 64bit dlls both included in the package.
So it should be fine. Maybe you compiled it for a different R version or
on another OS
Best,
Uwe Ligges
-Original Message-
To which I will add, you could have tried
count.fields('test.dat', sep='$' )
which I expect would have given you 3,3,2,3, and hence a pointer to where
the problem is.
count.fields is mentioned in the "See Also" section of ?read.table
--
Don MacQueen
Lawrence Livermore National Laboratory
70
## The next step is to think of this in the analysis of covariance setting.
## this example is based on
library(HH)
demo("ancova", package="HH", ask=FALSE)
item<-rep(letters[1:8], each=18)
day<-rep((0:5)*100, 24)
set<-rep(rep(1:3, each=6), 8)
test<-data.frame(item, day, set)
set.seed(111)
test$
On 2/12/2015 10:22 AM, John Sorkin wrote:
Windows 7, 64-bit.
I am trying to install RStudio. Before installing RStudio, I installed R 3.1.2. During the installation or R, I installled (as per the default) 32- and 64-bit packages. When I tried to install RStudio, I received the message
R does
Hi,
Maybe the following link helps.
http://r.789695.n4.nabble.com/box-cox-td4363944.html
--
View this message in context:
http://r.789695.n4.nabble.com/SIMPLE-question-tp4703176p4703180.html
Sent from the R help mailing list archive at Nabble.com.
It seems that Surv objects do not sort correctly. This seems to be a bug.
Anyone else found this?
>survival.data
[1] 4+ 3 1+ 2 5+
>class(survival.data)
[1] "Surv"
>sort(survival.data)
[1] 2 1+ 4+ 3 5+
An easy work-around is to define a function sort.Surv
>sort.Surv<-function(a){ord<-orde
On 12/02/2015 16:26, Professor Bickis wrote:
It seems that Surv objects do not sort correctly. This seems to be a bug.
Anyone else found this?
This is presumably about Surv() from package survival, not mentioned.
There was a bug, corrected in R-devel (and I will port to R-patched
before 3.
On Thu, Feb 12, 2015 at 7:12 PM, arnaud gaboury
wrote:
> On Thu, Feb 12, 2015 at 3:40 PM, arnaud gaboury
> wrote:
>> I have two df (and dt):
>>
>> df1
>> structure(list(name = c("poisonivy", "poisonivy", "poisonivy",
>> "poisonivy", "poisonivy", "poisonivy", "poisonivy", "poisonivy",
>> "cruzecon
I did not know the SVG file did not come through. I thought SVG should be
able to pass through the filter. Here is a PDF file along with an PNG.
Guess one of them should be able to pass.
祝好,
He who is worthy to receive his days and nights is worthy to receive* all
else* fr
On Thu, Feb 12, 2015 at 3:40 PM, arnaud gaboury
wrote:
> I have two df (and dt):
>
> df1
> structure(list(name = c("poisonivy", "poisonivy", "poisonivy",
> "poisonivy", "poisonivy", "poisonivy", "poisonivy", "poisonivy",
> "cruzecontrol", "agreenmamba", "agreenmamba", "vairis", "vairis",
> "vairis
I meant R3.0.3. I need this package working in R3.0.3.
I have 32bit and 64bit dlls both included in the package.
-Original Message-
From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com]
Sent: Thursday, February 12, 2015 1:01 PM
To: Li, Yan; r-help@r-project.org
Subject: Re: [R] 32bit R
On 12/02/2015 11:08 AM, Li, Yan wrote:
Dear All,
I build a R package which will need to call external C functions. I registered
the C functions in the NAMESPACE file and include 32bit and 64bit dlls in the
packages. If I load the package in 64bit R and calls the external C functions,
it works
Hello,
Try the following.
install.packages("sos") #if not yet
Then,
library(sos)
findFn("boxcox")
There are several hits, maybe you could start with package car
Hope this helps,
Rui Barradas
Em 12-02-2015 16:19, CHIRIBOGA Xavier escreveu:
Hi everybody!
I want to do a boxcox transform
Dear John,
I think you have to specify the directory of R in the R-studio installation
process. Somehow the installation wizard is unable to find the address of R
in your system.
Thanks & Regards,
Arnab
Arnab Kumar Maity
Graduate Teaching Assistant
Division of Statistics
Northern Illinois Uni
It seems that Surv objects do not sort correctly. This seems to be a bug.
Anyone else found this?
> survival.data
[1] 4+ 3 1+ 2 5+
> class(survival.data)
[1] "Surv"
> sort(survival.data)
[1] 2 1+ 4+ 3 5+
An easy work-around is to define a function sort.Surv
sort.Surv<-function(a){ord<-or
2.50% 97.50%
intercept 29787966
glands 13143611
damage 169 6144
treatment L1778 6703
treatment L4-3899 5125
Length -1817 1828
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.
Hi everybody!
I want to do a boxcox transformation, but I got this:
Error: could not find function "boxcox"
What can I do? I am using R studio.
Thanks!!!
Xavier
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://s
Hi,
I am trying to run the following model for OpenBUGS and want to use
R2OpenBUGS package. The model specifies weibull distribution for censored
data.
weibull.model <- function()
{
for(i in 1:n)
{
exp.alpha[i] ~ dgamma(a.alpha, b.alpha)
alpha[i] <- l
You should install Rtools first.
Then you will need to use
install.packages(c(A"","B","etc"), type="source", dependencies=TRUE)
Rich
On Thu, Feb 12, 2015 at 11:40 AM, PIKAL Petr wrote:
> Dear all
>
> I just switched to new version
>
>> version
>_
> platform i386-w64-mingw3
Maybe https://groups.google.com/forum/#!forum/shiny-discuss ?
Hadley
On Thu, Feb 12, 2015 at 9:07 AM, Doran, Harold wrote:
> I found a google user group for shiny, and am curious if there is an SIG as
> well. Didn't see one in my searches, but looking for an active place to ask
> questions and
Time to (re)read the Posting Guide... questions about unreleased versions of R
are off topic here. Go to R-devel.
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#.
Ah, ask in the RStudio support area?
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO
Dear all
I just switched to new version
> version
_
platform i386-w64-mingw32
arch i386
os mingw32
system i386, mingw32
status Under development (unstable)
major 3
minor 2.0
year 2015
month 02
day
Windows 7, 64-bit.
I am trying to install RStudio. Before installing RStudio, I installed R 3.1.2.
During the installation or R, I installled (as per the default) 32- and 64-bit
packages. When I tried to install RStudio, I received the message
R does not appear to be installed. Please install R
Dear All,
I build a R package which will need to call external C functions. I registered
the C functions in the NAMESPACE file and include 32bit and 64bit dlls in the
packages. If I load the package in 64bit R and calls the external C functions,
it works fine. However if I load the package in 3
I'm a bit blind today. I read df as a dput() .
John Kane
Kingston ON Canada
-Original Message-
From: hyil...@gmail.com
Sent: Thu, 12 Feb 2015 23:38:01 +0800
To: jrkrid...@inbox.com
Subject: Re: [R] ggplot2 shifting bars to only overlap in groups
You are most likely simply not running th
On 12/02/2015 13:56, michael.eisenr...@agroscope.admin.ch wrote:
There are a number of options for doing this. You might consider using
forest plots as provided by one of a number of the meta-analysis
packages. See the CRAN MetaAnalysis task view for more details.
--
Michael
http://www.dewe
Thank you very much, Ben. Your answer has been very useful.
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
You are most likely simply not running the whole lines of code: note that
the first line is:
N = 32
Best
,
He who is worthy to receive his days and nights is worthy to receive* all
else* from you (and me).
The Prophet, Gibr
And?
John Kane
Kingston ON Canada
> -Original Message-
> From: michael.eisenr...@agroscope.admin.ch
> Sent: Thu, 12 Feb 2015 13:56:37 +
> To: r-help@r-project.org
> Subject: [R] horizontal bar plots for CI visualization
FR
I am gettting the error"
Error in rep_len(rep.int(seq_len(n), rep.int(k, n)), length) :
object 'N' not found
Also your image did not come through. Try sending it as a pdf file.
when I try to create
df<- data.frame(gender=gl(2,1,N, c("male","female")),
direction=gl(2,2,N, c("up","do
I found a google user group for shiny, and am curious if there is an SIG as
well. Didn't see one in my searches, but looking for an active place to ask
questions and share code.
Thanks.
[[alternative HTML version deleted]]
__
R-help@r-project
Hi Rich
> -Original Message-
> From: Richard M. Heiberger [mailto:r...@temple.edu]
> Sent: Wednesday, February 11, 2015 10:53 PM
> To: PIKAL Petr
> Cc: r-help@r-project.org
> Subject: Re: [R] suggestion for optimal plotting to show significant
> differences
>
> Petr,
>
> My first attempt i
I have two df (and dt):
df1
structure(list(name = c("poisonivy", "poisonivy", "poisonivy",
"poisonivy", "poisonivy", "poisonivy", "poisonivy", "poisonivy",
"cruzecontrol", "agreenmamba", "agreenmamba", "vairis", "vairis",
"vairis", "vairis", "vairis", "vairis", "xaeth"), text = c("ok",
"need items
Hi all,
I have four factors for a continuous time variable along with its
confidence interval. I would like to produce a publication quality error
bar chart that is clear to understand. For now, I used colors, x axis
position, facets and alpha level to distinguish them.
I would like to overlap ea
2.50% 97.50%
intercept 29787966
glands 13143611
damage 169 6144
treatment L1778 6703
treatment L4-3899 5125
Length -1817 1828
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.
> R-ouput
>
> Df Sum SqMean SqF value
> Pr(>F)
> treatment314015 4672 65.021 < 2e-16
> ***
> variety 3 5883 1961 27.295
> 2.95e-12 ***
> treatment:variet
Hello Everyone...
I am fairly new to R. I would like to use the adehabitatHR package to create
kernel density and isopleths from my sea turtle GPS data. I�m running into some
issues�
Basically I am having trouble assigning IDs and XY fields within R to create
KDs and MCPs.
Thank you in advan
I would not say it is fully "solved" since in using Nelder-Mead
you did not get the Hessian.
The issue is almost certainly that there is an implicit bound due to
log() or sqrt() where a parameter gets to be near zero and the finite
difference approximation of derivatives steps over the cliff. Prob
You don't need to do these operations in pieces so the mapply is unnecessary.
Neither dplyr nor data.table can go faster than (assuming your data frame is
called DF):
DF$dtm <- as.POSIXct( as.numeric( DF$ts ), tz="GMT", origin="1970-01-01" )
They can in this case save you from having to retype
> I would like to query that if it is possible for R to terminate a program
> such as
> notepad, word, etc.?
pskill from the tools package says it supports the Windows system call
'TerminateProcess'.
You can get a process's ID using the windows command line 'tasklist'.
I got this to work on my
I am little lost between all the possibilities to apply a function to
a data.frame or data.table.
Here is mine:
structure(list(name = c("poisonivy", "poisonivy", "poisonivy",
"poisonivy", "poisonivy", "poisonivy", "poisonivy", "poisonivy",
"cruzecontrol", "agreenmamba", "agreenmamba", "vairis", "
I done a mistake when I paste the script in the message. An update:
url.csv <-
https://dl.dropboxusercontent.com/u/34009642/cepajabo07_wide_acumulado.csv
data02 <- read.table(url.csv, header=TRUE, sep="\t", dec=",")
head(data02)
timepoints <- 1:5 # 5 days
time <- t(matrix(rep(timepoints, 120),
Hi Samarvir,
Assuming that you want to generate a separate data frame for each
value of "Name",
# name of initial data frame is ssdf
for(nameval in unique(ssdf$Name)) assign(nameval,ssdf[ssdf$Name==nameval,])
This will produce as many data frames as there are unique values of
ssdf$Name, each name
hello,
I am cleaning some large data with 4 million observation and 7 variable.
Of the 7 variables , 1 is name/string
I want to subset data, which have same name
Example-
Name var1 var2 var3 var4 var5 var6
aa- - - - --
ab
bd
ac
ad
af
ba
bd
aa
av
i want
Make use of the plyr and reshape2 package (both on CRAN):
library(plyr)
d<-adply(ArrayDiseaseCor, 1:2)
# adply calls function identity by default
d<-melt(d)
d<-subset(d,value>.5)
head(d)
You will have to rename columns, or adjust arguments in melt/adply.
Note: use set.seed before sampling for rep
see inline
On 12 February 2015 at 09:10, Sven E. Templer wrote:
> Make use of the plyr and reshape2 package (both on CRAN):
>
I forgot:
library(reshape2)
> library(plyr)
> d<-adply(ArrayDiseaseCor, 1:2)
> # adply calls function identity by default
> d<-melt(d)
> d<-subset(d,value>.5)
> head(d)
51 matches
Mail list logo