I have recently acquired a new laptop and a new OS (Ubuntu 16.04) and
have encountered a problem when trying to install rgdal.
The initial hiccups were overcome by installing "libgdal-dev" and
"libproj-dev" on my system, and then re-installing the package "sp" (as
advised by postings that I
While running the vglm() with specification of cumulative(link=cloglog)
with a more complicated model, the following errors appeared. How to go
about solving it?
Error in if ((temp <- sum(wz[, 1:M, drop = FALSE] < wzepsilon)))
warning(paste(temp, :
argument is not interpretable as logical
In
Thank you very much Peter - very enlightening to do it from first principles!
Regards,
Paul
On 4 October 2016 at 10:23:52 am, peter dalgaard (pda...@gmail.com) wrote:
> On 04 Oct 2016, at 00:30 , Paul Sanfilippo wrote:
>
> Hi,
>
> I am trying to replicate a test in the Hosmer - Appl
> On 04 Oct 2016, at 00:30 , Paul Sanfilippo wrote:
>
> Hi,
>
> I am trying to replicate a test in the Hosmer - Applied Logistic regression
> text (pp 289, 3rd ed) that uses a Multivariable Wald test to test the
> equality of coefficients across the 2 logits of a 3 category response
> multin
Thanks Bert,
I realise that and am not distilling it down to the p value. I am primarily
considering the issue of collapsing down in the larger context of the added
value/information that the extra categories give. However, I am curious to see
(as I said from a statistical standpoint) whether t
See inline.
-- Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along
and sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Mon, Oct 3, 2016 at 3:30 PM, Paul Sanfilippo wrote:
> Hi,
>
> I am trying to replicate a t
Hi,
I am trying to replicate a test in the Hosmer - Applied Logistic regression
text (pp 289, 3rd ed) that uses a Multivariable Wald test to test the equality
of coefficients across the 2 logits of a 3 category response multinomial model.
I’d like to see whether (from a statistical standpoint)
Hi Syela,
Are the values in ASFR monotonically increasing with year?
Jim
On Tue, Oct 4, 2016 at 4:23 AM, Syela Mohd Noor wrote:
> Hi all, I had a problem with the parameter estimation of the Brass Gompertz
> model for my dissertation. I run the code for several times based on
> different years
Hi Maryam,
Your labels have been "greeked" as the font is too small to be
displayed properly. If you must use PNG format, specify your image
file at least twice as high.
png("pheatmap.png",width=1254,height=5000)
PDF would be a better choice as you can just zoom in and scroll down.
Jim
On Mon,
> On Oct 3, 2016, at 4:43 AM, maryam moazam wrote:
>
> Hi all,
>
> I'm trying to show about 190 items using "pheatmap" package. I'm working
> with R 3.3.1 on windows 7. I used the following code:
>
> pheatmap(mat_data, cellwidth=20, cluster_rows = FALSE, cluster_cols =
> FALSE, fontsize=5, fon
On Mon, 3 Oct 2016, Frank S. wrote:
Dear R users,
[deleted]
I want to get a list of "k" data tables (or data frames) so that each
contains those individuals who for the first time are at least 65,
looping on each of the dates of vector "v". Let's consider the following
example with 5 ind
Get a SLURM batch file from someone there and make the appropriate changes
to it. Best way to learn.
On Mon, Oct 3, 2016 at 2:10 PM, Dominik Schneider
wrote:
> I typically call Rscript inside an sbatch file.
>
> *batch_r.sh*
> #! /bin/bash
>
> cd /home//aso_regression_project #make sure you c
Hi Frank,
How about
library(lubridate)
dtf <- merge(dt, expand.grid(id = dt$id, refdate = v), by = "id")
dtf[, gt65 := as.period(interval(fborn, refdate), unit = "years") > years(65)]
dtf <- dtf[gt65 == TRUE,][, .SD[refdate == min(refdate)], by = id]
Best,
Ista
On Mon, Oct 3, 2016 at 1:17 PM, F
You could try
dbListTables()
and do the subsetting afterwards.
Also try
dbGetQuery()
instead of dbSendQuery()
(There's no obvious reason that I can see to use dbSendQuery() followed by
fetch. It would be interesting if the two ways return different lists of
tables.)
--
Don MacQueen
Lawrenc
I typically call Rscript inside an sbatch file.
*batch_r.sh*
#! /bin/bash
cd /home//aso_regression_project #make sure you change to your
correct working directory
Rscript scripts/run_splitsample-modeling.R
and on the commandline of the login node:
sbatch batch_r.sh
There are lots of slurm con
Dear Authorized Sir / Madam,
I have an R script file in which it includes this lines:
How can i to submit this R jobs via SLURM? Thanks in advance.
*testscript.R*
data=read.table("seqDist.50", header=FALSE)[-1]
attach(data)
d=as.matrix(data)
library(cluster)
cluster.pam = pam(d,6)
table(cluster.
Hi all, I had a problem with the parameter estimation of the Brass Gompertz
model for my dissertation. I run the code for several times based on
different years and it was fine until the seventh year onward where I got
negative values for the parameter (a) which did not make sense since it
represen
Hi all,
I'm trying to show about 190 items using "pheatmap" package. I'm working
with R 3.3.1 on windows 7. I used the following code:
pheatmap(mat_data, cellwidth=20, cluster_rows = FALSE, cluster_cols =
FALSE, fontsize=5, fontsize_row=5, margins=c(5,10))
The resulted image is attached. Please
John M. Chambers Statistical Software Award 2017
The Statistical Computing Section of the American Statistical
Association announces the competition for the John M. Chambers
Statistical Software Award. In 1998 the Association for Computing
Machinery presented its Software System Award to John
Hi R help,
I am trying to do some AIC model averaging on a CLM model in R and keep getting
the error:
Error in aictab.default(cand.set = Cand.model0, modnames = Modnames0, :
Function not yet defined for this object class
The MuMIn package says that the functions should work for clm and clmm m
Hello,
I've just ran your code and it all went well.
So my doubt is: if you have 1269 rows why choose only 100 and
bootstrap? It doesn't seem to make much sense to me.
Try to run the entire df through DataSummary and compare the results
with the bootstrap results.
Rui Barradas
Citando Bry
Hi all,
Here is the first six rows of my data. In total I have 1269 rows.
My goal is to get conduct nonparametric bootstrap and case resampling.
I would like to randomly select 100 out of the 1269 After that, I wish to
bootstrap that randomly selected 100 out of 1269.
I assume I need to set t
If this is a one-time change, it will be simpler to use the Windows File
explorer (I assume Windows since you are using backslashes). If it needs to be
done using R functions, this seems to work:
> dir.create("A/B/C/D/E", recursive=TRUE)
> list.dirs("A")
[1] "A" "A/B" "A/B/C" "
On Mon, Oct 3, 2016 at 11:13 AM, Vivek Singh wrote:
> Hi All:
>
> I am using R with MySQL, following is my code:
>
> library(RMySQL)
> con <- dbConnect(MySQL(), user="user", password="password", dbname="db",
> host="localhost")
> rs1 <- dbSendQuery(con, "show tables like \"%padded\"")
> all_table
Hi All:
I am using R with MySQL, following is my code:
library(RMySQL)
con <- dbConnect(MySQL(), user="user", password="password", dbname="db",
host="localhost")
rs1 <- dbSendQuery(con, "show tables like \"%padded\"")
all_tables3<- fetch(rs1, n=-1)
There are 25 tables in the database when I exec
- ndjdon : Wicked Fast ndjson Reader
Reads in ndjson significantly faster than jsonlite::stream_in(), flattens
each
JSON record and returns a data.table.
https://cran.r-project.org/web/packages/ndjson/index.html
- htmltidy : Clean Up or Pretty Print Gnarly HTML and XHTML
C-backed packa
Why not use the effects package -- designed for this task.
library(effects)
plot(allEffects(mfit))
?plot.eff # for details
On 10/3/2016 3:15 AM, Faradj Koliev wrote:
Dear all,
I need a little help with plotting predicted probabilities (values). Consider
the following example
**
data(”mtca
Dear all,
I need a little help with plotting predicted probabilities (values). Consider
the following example
**
data(”mtcars”)
mfit = lm(mpg ~ vs + disp + cyl, data=mtcars)
newcar=data.frame(vs=c(0,1), disp=230, cyl=6.188)
Pmodel<–predict(mfit, newcar)
**
I want to plot the effect of ”vs
28 matches
Mail list logo
