A new version of the package FRAILTYPACK is now available on CRAN.
-- possibility to fit now a Shared and a Joint Frailty model with a log-normal
distribution for the random effects.
-- possibility to deal with interval-censored data (for a shared frailty model)
-- possibility to fit a joint
Absolutely correct. The documentation wasn't clear on that. Changed for the
next version of the package. And the escalc() function now no longer checks for
sd2i, since it is not needed anyway (for SMCR). For now, you will just have
to set sd2i to something (e.g., 0).
Thanks for bringing this
version.string R version 2.15.2 (2012-10-26)
I am trying to install biOps on MacOS X 10.8.2
First, I have tiff, fftw-3, jpeg
and set paths like
cd /usr/include
sudo ln -s /usr/local/include/fftw3.h
for x in /usr/local/include/j*.h; do sudo ln -s $x; done
for x in /usr/local/include/tiff*.h; do
I have two binary files(rasters) with the same dimensions. The first file is
called `over` and the second is `corr`. I want to replace values in `over`
by `NA` whenever `corr` is greater than 0.5.
to read the two files we can use:
conne - file(C:corr.bin,rb)
over -
See the recent discussion on R-sig-mac (the place to ask questions about
OS X).
On 04/03/2013 11:00, ishi soichi wrote:
version.string R version 2.15.2 (2012-10-26)
I am trying to install biOps on MacOS X 10.8.2
First, I have tiff, fftw-3, jpeg
and set paths like
cd /usr/include
sudo ln -s
HI,
Assume I have the data frame generated from R as the following:
Product Price market_share
A 10010%
B 1109%
C 12020%
D 90 61%
What I want to do is to have this kind of excel user interface by
Hello,
I create a plot from a coxph object called fit.ads4:
plot(survfit(fit.ads4))
plot is located at:
https://www.dropbox.com/s/9jswrzid7mp1u62/survfit%20plot.png
I also create the following survfit statistics:
print(survfit(fit.ads4),print.rmean=T)
Call: survfit(formula = fit.ads4)
Please help!
I made an oa design using the DoE package of R.Then, I created the 18 blocks
containing 4 choice sets and converted these into questionnaire form using
rotation method of the CE support package.However,some of the alternatives in
the choice sets are infeasible/unrealistic and
Hi,
I encountered a segfault, memory not mapped error when using adjbox in
robustbase. In trying to recreate the issue I found that the error
occurs only for large sample size. Here is the code.
require(robustbase)
Loading required package: robustbase
x - rnorm(10)
y - rep(1, 10)
Hi everyone!! I am new in R and I want to create a simple R function for
estimating historical-VaR. In y_IBM returns, there are 2300 observations. For
evaluation I take the next 2000 observations,
then I abandon the latest 300 observations. Firstly, I use the window which has
the fix
length
I'm late to this discussion, but let me try to put it in another context.
Assume that I wanted to know whether kids who live west of their school or east of
their shool are more likely to be early (some hypothesis about walking slower if the sun
is in their eyes). So I create a 0/1 variable
Hi,
I am using R on Ubuntu 12.10 (64bit). This morning, Ubuntu's software
updater automatically installed updates to R's base system (version
2.15.3; via the CRAN PPA). Now R does not work anymore. Here is what I
get when I simply enter R on the shell prompt:
bodenhof FUKUOKA~R
cannot find
Hi,
Please see http://stackoverflow.com/q/15090817, where it was reported
that sudo chmod 644 /etc/R/Renviron fixed this issue. If that doesn't
work for you, consider posting to
https://stat.ethz.ch/mailman/listinfo/r-sig-debian
Best,
Ista
On Mon, Mar 4, 2013 at 9:56 AM, Ulrich Bodenhofer
On 4 March 2013 at 10:01, Ista Zahn wrote:
| Please see http://stackoverflow.com/q/15090817, where it was reported
| that sudo chmod 644 /etc/R/Renviron fixed this issue. If that doesn't
That issue affected only pre-releases of 2.15.3. The 2.15.3 release does not
have that issue. It also should
Yes, it can cause problems. And speaking for myself, I'd say it's not
worth the risk, because it's easy enough to find alternative variable
names that are close enough to the notation of your formulas that
remembering should be no problem. For example, tt, cc, and mmatrix
might do it.
-Don
--
Dear Wolfgang,
Thank you for your clarification - and your excellent package !
Best wishes,
Markus
-Ursprüngliche Nachricht-
Von: Viechtbauer Wolfgang (STAT)
[mailto:wolfgang.viechtba...@maastrichtuniversity.nl]
Gesendet: Montag, 4. März 2013 10:45
An: Markus Kösters;
I was also facing the same problem.
It may be possible that in your computer Java may not be installed, try
installing Java. My problem was solved after installing it.
http://www.java.com/en/download/manual.jsp
Thanks.
Sent from Windows Mail
[[alternative HTML version deleted]]
HI,
Check these links:
http://r.789695.n4.nabble.com/Re-new-question-td4659908.html
http://r.789695.n4.nabble.com/Re-reading-data-td4658705.html#a4659751
It shows reading files (not binary files) from multiple subdirectories and
doing some calculations.
You may need to change it according to
It seems within the last ~3 months Ive been having issues with writing text or
csv files from a R data frame. The problem is multifold and it is hard to
filter out what is going on and where the problem is. So, Im hoping someone
else has come across this and may provide insight.
My
I just upgraded from 2.15.2 to 2.15.3 on our Slackware systems and see the
following when testing that the new version properly loads:
The following object(s) are masked from 'package:stats':
cor
The following object(s) are masked from 'package:base':
as.Date, as.Date.numeric
What,
Hello everyone!
Does anyone of you know how I could solve the following problem.
I guess, it is not a very difficult question, but I simply lack of the
right idea:
I have a dataset containing data of convictions. This dataset contains 4
columns:
- personId: individual number that identifies the
What Jim said separately is correct, and I would suggest following his
advice.
But there are some points worth looking at in your method.
See this example:
item1 - item2 - item3 - item4 - 1:4
matrix1-cbind(item1, item2, item3, item4)
z - c(TRUE,TRUE,FALSE,TRUE)
matrix2 -
On 04/03/2013 10:09 AM, Kerry wrote:
It seems within the last ~3 months Ive been having issues with writing text or
csv files from a R data frame. The problem is multifold and it is hard to
filter out what is going on and where the problem is. So, Im hoping someone
else has come across
Hi Rich,
Immediately after you see one of those messages, do, e.g.,
find('cor')
It should tell you that you have more than one object named 'cor' in your
search path, and where they all are. Then you can decide if it's what you
want (probably not, but can't say from here).
-Don
--
Don
Here's the first 5 lines of my dataset:
structure(list(id = structure(c(1L, 1L, 1L, 1L, 1L), .Label = c(F07001,
F07002, F07003, F07004, F07005, F07006, F07008, F07009,
F07010, F07011, F07014, F07015, F07017, F07018, F07019,
F07020, F07021, F07022, F07023, F07024, F10001, F10002,
F10004,
On 04/03/2013 10:52 AM, Kerry wrote:
Here's the first 5 lines of my dataset:
structure(list(id = structure(c(1L, 1L, 1L, 1L, 1L), .Label = c(F07001,
F07002, F07003, F07004, F07005, F07006, F07008, F07009,
F07010, F07011, F07014, F07015, F07017, F07018, F07019,
F07020, F07021, F07022, F07023,
Here is one way. There will be many ways to do it; I offer this one
because it is very general.
-Don
tmp - split(testdata, testdata$personId)
myfun - function(df) {
dfo - df
if (any(df$law=='SVG')) dfo$svg - 1 else dfo$svg - 0
dfo
}
tmpo - lapply(tmp,myfun)
testout - do.call('rbind',
I realize my command code is not writing a CSV file, I already pointed out that
CSV file seems to work OK but not in the TXT format. Regardless of that, there
should be no problems in ArcGIS in reading the date column - in ArcGIS it will
simply recognize it as a text field by default. As I
Have a look at cast() from the reshape package.
library(reshape)
cast(personId ~ law, data = testdata, value = article, fun = length)
cast(personId ~ law, data = testdata, value = article, fun = function(x){1 *
(length(x) 0)})
Van:
There’s more than one way to skin a cat, here is another
mm-model.matrix(~personId+law+0,testdata)
merge(testdata,aggregate(mm[,-1],list(personId=mm[,personId]),max))
cheers
Am 04.03.2013 16:44, schrieb David Studer:
Hello everyone!
Does anyone of you know how I could solve the following
B == Baan baanba...@gmail.com
on Mon, 4 Mar 2013 15:02:02 +0530 writes:
B Hi, I encountered a segfault, memory not mapped error
B when using adjbox in robustbase. In trying to recreate
B the issue I found that the error occurs only for large
B sample size. Here is the code.
On Mon, 4 Mar 2013, MacQueen, Don wrote:
Immediately after you see one of those messages, do, e.g.,
find('cor')
It should tell you that you have more than one object named 'cor' in your
search path, and where they all are. Then you can decide if it's what you
want (probably not, but can't
Thanks a lot!!!
On Mon, Mar 4, 2013 at 12:46 PM, MacQueen, Don macque...@llnl.gov wrote:
What Jim said separately is correct, and I would suggest following his
advice.
But there are some points worth looking at in your method.
See this example:
item1 - item2 - item3 - item4 - 1:4
On 04/03/2013 11:15 AM, Kerry wrote:
I realize my command code is not writing a CSV file, I already pointed
out that CSV file seems to work OK but not in the TXT format.
Sorry. In that case, I think you really do have to go to R-sig-Geo to
find someone who knows about ArcGIS.
Duncan
Your description of diagnosis uses non-R software (off topic here). Please
either describe the difference in the files (you may need a hex editor or the
hexbin package to detect the differences) or supply the files that behave
differently (this may require some alternate route than this mailing
Hi List,
Â
Iâ m analysing the selectivity of
resting site use by forest carnivores through mixed modelling techniques and I
wonder which will be the best r package to deal with several aspects
simultaneously:
-Â Â Â Â Â Â Â Â Â binomial
variable response;
-Â Â Â Â Â Â Â Â Â possible
spatial
Thank you Martin. Look forward to the fix.
Regards
Baan
On Monday 04 March 2013 10:19 PM, Martin Maechler wrote:
B == Baan baanba...@gmail.com
on Mon, 4 Mar 2013 15:02:02 +0530 writes:
B Hi, I encountered a segfault, memory not mapped error
B when using adjbox in robustbase.
Does it just not work or does it not do the right thing? The reason it doesn't
work is that you are writing 'T = length(returns) x_foc = vector(length=n) N =
T-(n+1)' on one line instead of using three lines. However, your description of
what you want to do also doesn't seem to correspond to
1. You are right with your statements.
2. You should address all further discussion to the r-sig-mixed-models
list, not here.
-- Bert
On Mon, Mar 4, 2013 at 9:06 AM, Filipe Carvalho
filipescpcarva...@yahoo.com wrote:
Hi List,
I’ m analysing the selectivity of
resting site use by forest
Hello again,
Let say I have following vector:
set.seed(1)
Vec - sample(LETTERS[1:5], 10, replace = TRUE)
Vec
Now with each repeated letter, I like to add suffix programatically.
Therefore I want to get following vector:
c(B, B1, C, E, B2, E1, E2, D, D1, A)
Can somebody tell me how to achieve
On Mar 4, 2013, at 8:28 AM, Rich Shepard wrote:
On Mon, 4 Mar 2013, MacQueen, Don wrote:
Immediately after you see one of those messages, do, e.g.,
find('cor')
It should tell you that you have more than one object named 'cor' in your
search path, and where they all are. Then you can
On Mar 4, 2013, at 5:39 AM, Tammy Ma wrote:
HI,
Assume I have the data frame generated from R as the following:
Product Price market_share
A 10010%
B 1109%
C 12020%
D 90 61%
What I want
Use ave(). The following tries to save a little time and simplify things
by only processing the duplicated entries.
Vec
[1] B B C E B E E D D A
f - function(v) {
+ d - duplicated(v)
+ v[d] - ave(v[d], v[d], FUN=function(vdi)paste0(vdi, seq_along(vdi)))
+ v
+ }
Christofer:
This reminds me of those IQ puzzles I took so many years ago as a kid:
Given the numbers 7, 22, 43, 5, 26,...what are the next 3 numbers
in this series? I don't recall having a clue, and when I got older and
more mathematical, generally came to the conclusion that it could
Hi
It sounds like you just want to control the size of the external raster
image that you are creating. For example ...
library(grid)
im - matrix(0:1, ncol=3, nrow=3)
# grid.newpage()
grid.raster(im, interpolate=FALSE)
png(ident.png, 3, 3)
grid.raster(im, interpolate=FALSE)
dev.off()
I am trying to create contingency tables (to evaluate prior to performing
Pearson's Chi-Squared test for independence). I would like to see column
and row totals as well as expected and observed values and cell counts.
I tried to use the package contingency. tables but get the following
warning:
From time to time since the release of R 2.15.3 was announced
I have tried to upgrade from 2.15.2 by doing
sudo yum install R
I keep getting
Package R-2.15.2-1.fc17.x86_64 already installed and latest version
Nothing to do
Has the upgrade really not propagated to the Fedora
Hi Elisa,
I am not sure about the output you wanted.
dat1- read.csv(rightest.csv,sep=,,header=TRUE,check.names=FALSE)
dat2- as.dist(dat1[,-1],upper=F,diag=F)
vec1- as.vector(dat2)
label1=c(0-25,25-50,50-75)
Count1- as.data.frame(table(cut(vec1,breaks=seq(0,75,25),labels=label1)))
#Overall count
I never saw the original note nor its resubmission. Nor could I find it
in the R-help archives.
Is it just me?
Not that it really matters a damn --- I don't use SAS either --- I'm just
curious, and find the situation mysterious.
cheers,
Rolf
On 03/04/2013 04:32 PM, Frank
On Mar 4, 2013, at 1:48 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
From time to time since the release of R 2.15.3 was announced
I have tried to upgrade from 2.15.2 by doing
sudo yum install R
I keep getting
Package R-2.15.2-1.fc17.x86_64 already installed and latest version
On Mon, Mar 4, 2013 at 2:57 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
I never saw the original note nor its resubmission. Nor could I find it
in the R-help archives.
Is it just me?
No, same here.
Not that it really matters a damn --- I don't use SAS either --- I'm just
curious, and
Hi,
You could use:
res-unsplit(lapply(split(Vec,Vec),function(x) if(length(x)1)
c(head(x,1),paste0(head(x,-1),seq_along(head(x,-1 else x),Vec)
res
# [1] B B1 C E B2 E1 E2 D D1 A
A.K.
- Original Message -
From: Christofer Bogaso bogaso.christo...@gmail.com
To: r-help
On Mar 4, 2013, at 11:31 AM, William Dunlap wrote:
Use ave(). The following tries to save a little time and simplify things
by only processing the duplicated entries.
Vec
[1] B B C E B E E D D A
f - function(v) {
+ d - duplicated(v)
+ v[d] - ave(v[d], v[d],
Le lundi 04 mars 2013 à 12:37 -0600, Joanna Papakonstantinou a écrit :
I am trying to create contingency tables (to evaluate prior to performing
Pearson's Chi-Squared test for independence). I would like to see column
and row totals as well as expected and observed values and cell counts.
I
Thanks Marc. All very clear, cogent and useful. Sorry for not
sending my original post to r-sig-fedora; I *knew* about that,
but had (in my dotage?) forgotten.
cheers,
Rolf
On 03/05/2013 09:01 AM, Marc Schwartz wrote:
On Mar 4, 2013, at 1:48 PM, Rolf Turner
Le mardi 05 mars 2013 à 08:57 +1300, Rolf Turner a écrit :
I never saw the original note nor its resubmission. Nor could I find it
in the R-help archives.
Is it just me?
Not that it really matters a damn --- I don't use SAS either --- I'm just
curious, and find the situation mysterious.
On 2013-03-04 12:04, Ista Zahn wrote:
On Mon, Mar 4, 2013 at 2:57 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
I never saw the original note nor its resubmission. Nor could I find it
in the R-help archives.
Is it just me?
No, same here.
Think Nabble!!!
(where there are often waits for
Hi
An alternative is the 'gridSVG' package
http://www.stat.auckland.ac.nz/~paul/gridSVG/gridsvg.pdf
Paul
On 02/03/13 06:12, Joseph Clark wrote:
Hi all:I found some great demonstrations of interactive presentation graphics generated in R with the SVGAnnotation package,
Hi,
We have comment questions from a survey that we need to categorize. What
package and functions can I use in R to help do this?
Daniel Lopez
Lawrence Livermore Labs
SHRM
[[alternative HTML version deleted]]
__
R-help@r-project.org
On Mar 4, 2013, at 12:18 PM, Peter Ehlers wrote:
On 2013-03-04 12:04, Ista Zahn wrote:
On Mon, Mar 4, 2013 at 2:57 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
I never saw the original note nor its resubmission. Nor could I find it
in the R-help archives.
Is it just me?
No, same
Hi,
Sorry, I forgot names(which()
Names2-apply(res1,1,function(x) paste(names(which(x!=0)),collapse=,))
res2- data.frame(Frequency=apply(res1,1,function(x) sum(1*(x!=0))),
stations=Names2,stringsAsFactors=FALSE)
res2
# Frequency
#0-25 121
#25-50 122
#50-75 76
Daniel,
Looking at CRAN, the following might be useful:
memiscTools for Management of Survey Data, Graphics, Programming,
Statistics, and Simulation
questionr Functions to make surveys processing easier
surveyanalysis of complex survey samples
surveydataTools to
Hi,
dat1- read.csv(rightest.csv,sep=,,header=TRUE,check.names=FALSE)
dat2- as.dist(dat1[,-1],upper=F,diag=F)
vec1- as.vector(dat2)
label1=c(0-25,25-50,50-75)
Name1-unlist(lapply(0:123,function(i) rep(i+1,i)))
dat3-data.frame(Name1,vec1)
res-t(aggregate(.~Name1,data=dat3,function(x)
Thank you.
I actually ended up using:
CrossTable(mdt)
Cell Contents
|-|
| N |
| Chi-square contribution |
| N / Row Total |
| N / Col Total |
| N / Table Total |
|-|
Total Observations in
On 2013-03-04 12:30, David Winsemius wrote:
On Mar 4, 2013, at 12:18 PM, Peter Ehlers wrote:
On 2013-03-04 12:04, Ista Zahn wrote:
On Mon, Mar 4, 2013 at 2:57 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
I never saw the original note nor its resubmission. Nor could I find it
in the
On 03/05/2013 09:28 AM, Lopez, Dan wrote:
Hi,
We have comment questions from a survey that we need to categorize. What
package and functions can I use in R to help do this?
There was a fairly recent posting about a package called qdap that might
be relevant. See:
R 2.15
OS X and Windows
Colleagues,
I am running a lengthy R script ( 20K lines of code) within a GUI build in
RealBasic. There are rare instances in which an error is triggered in R. I am
trying to trap these based on their syntax. I have identified three general
formats for error
On Mar 4, 2013, at 3:09 PM, Fisher Dennis wrote:
R 2.15
OS X and Windows
Colleagues,
I am running a lengthy R script ( 20K lines of code) within a GUI build in
RealBasic. There are rare instances in which an error is triggered in R. I
am trying to trap these based on their syntax.
There is something that I wish I had that I think would help me a lot to be a
better R programmer, that I think would probably help many others as well.
I put the wish out there in the hopes that someone might think it was worth
doing at some point.
I wish I had the code of some substantial,
Hi,
I am attaching a plot where you can see there are a few jumps (plots 1, 4,
5 and 6), due to incidents with the measuring sensors (basically someone
touching the sensor). I need to revert those changes to have a plot without
unreal measurements, so make those fragments go back to its original
On Mar 4, 2013, at 3:42 PM, andrewH wrote:
There is something that I wish I had that I think would help me a lot to be a
better R programmer, that I think would probably help many others as well.
I put the wish out there in the hopes that someone might think it was worth
doing at some
Hi,
If you meant to get the array index from `res1`, then this will do it. But, if
it is from 'dat3`, it will be huge as each index is unique.
indx-which(apply(res1,1,function(x) x!=0) ,arr.ind=TRUE)
Names2-do.call(rbind,lapply(1:3,function(i) {x-indx[indx[,2]==i,];
Thanks for the illustrative example. In my project actually my supervisor
wanted to estimate the probabilities using a conditional MLE approach,
which happens to be the case that *uses clogit() while trying to achieve
aim b in your words*.
I learned that clogit() is based on the sufficient
Hi all, I have a couple of questions.
1) Is there a more effective way to remove the following pattern
(patternrm) from the station_id_Timeseries (see below)
patternrm - c(_2000_4_OND, _2001_1_JFM, _2001_2_AMJ, _2001_3_JAS,
_2001_4_OND, _2002_1_JFM, _2002_2_AMJ, _2002_3_JAS, _2002_4_OND,
Hi Cesar
Not sure what you actually want to accomplish
?rle may give you some ideas eg (I have added some to return to the
good section)
x = c(246,251,250,255,5987,5991,5994,599,255,259,262,267)
xdiff = diff(x)
xdiff
[1] 5-1 5 5732 4 3 -5395 -344 4 3 5
Hi,
May be you can try:
res1- gsub(\\_.*\\_.*\\_.*,,x)
res2-sub((\\w)(\\w)(\\w)(\\w)(\\w)(\\w)(\\w)(\\w)(\\w)(\\w)(\\w)$, ,
x,perl=TRUE)
identical(res1,res2)
#[1] TRUE
A.K.
- Original Message -
From: Irucka Embry iruc...@mail2world.com
To: r-help@r-project.org
Cc:
Sent: Monday,
The book that helped me break into R and more advanced texts was
Crawley's Statistics: An Introduction with R. Very light read that
assumes no prior knowledge with stats or R. I am using it to teach my
fellow grad students R and all agree it was worth scrimping pennies to
get. He also has a
hello sir
myself Jyoti Sharma, and i am working as preoject fellow in IGIB Delhi.
I need your help to know how to create a package as well as how to
post that package to the CRAN mirror, for public use. i have used
package.skeleton() command but its not working properly.
So, please suggest me
Hi Andrew: Not that I've gone through it all yet but the draft of hadley's
book at https://github.com/hadley/devtools/wiki/Introduction has a lot if
not all of the commands you refer to and all of their gory details along
with many examples. No matter what you're budget, given that the book will
On Mar 4, 2013, at 7:42 PM, Jyoti Sharma wrote:
hello sir
myself Jyoti Sharma, and i am working as preoject fellow in IGIB Delhi.
I need your help to know how to create a package as well as how to
post that package to the CRAN mirror, for public use. i have used
package.skeleton()
Hi Arun, thank you for your assistance. The code for part 1 worked
perfectly.
I am including here a list of all of the site names for part 2. Due to
the nature of the station ids it has been difficult for me to find a
good pattern to use with regular expressions.
Thank you all in advance for
curiosity question: I was wondering whether the R binaries and BLAS
libraries for ubuntu linux are compiled using SSE4 and AVX support.
this probably can go a long way towards a unified memory bus GPGPU
substitute.
Ivo Welch (ivo.we...@gmail.com)
http://www.ivo-welch.info/
Hi,
I'm trying to use the data which I generate within R to make images in .bmp
format to be lateron printed by a printer.
My first thought was the RImageJ package, but this seems to be discontinued.
What I am currently doing is generating a matrix of grey values, which needs to
be parsed
83 matches
Mail list logo