On Wed, 10 Jun 2020 15:36:11 -0500
Ana Marija wrote:
> Hello,
>
> I have a data frame like this:
>
> > head(tmp1)
> CHR BP PoldPnew
> 1 1 785989 0.9521 0.09278
> 2 1 1130727 0.4750 0.19010
> 3 1 1156131 0.5289 0.48520
> 4 1 1158631 0.2554 0.18140
> 5 1 1211292 0.2954
Hi Saudi,
I can only make a guess, but that is that a variable having a unique
value for each participant has been read in as a factor. I assume that
"better" is some combination of "hum" and "cul" and exactly what is
WF?
Jim
On Thu, Jun 11, 2020 at 5:27 AM Saudi Sadiq wrote:
>
> Dear
On 2020-06-10 18:01 -0700, Jeff Newmiller wrote:
> On June 10, 2020 5:29:10 PM PDT, Jinsong Zhao wrote:
> >
> > (3,1), (5,1), (5,2), (4,2), (4,3), (1,3), (1,2), (3,2)
>
> M <- matrix(c(2,2,rep(1,12), 2), nrow = 5,byrow = FALSE)
> ix <- expand.grid( r = seq.int( nrow( M ) )
> , c
M <- matrix(c(2,2,rep(1,12), 2), nrow = 5,byrow = FALSE)
ix <- expand.grid( r = seq.int( nrow( M ) )
, c = seq.int( ncol( M ) )
)
ix[ 1 == c(M), ]
On June 10, 2020 5:29:10 PM PDT, Jinsong Zhao wrote:
>Hi there,
>
>I have a matrix similar as:
>
>M <-
Hi there,
I have a matrix similar as:
M <- matrix(c(2,2,rep(1,12), 2), nrow = 5,byrow = FALSE)
I hope to get the border subscript of the block with value 1. In the
above example, I hope to get:
(3,1), (5,1), (5,2), (4,2), (4,3), (1,3), (1,2), (3,2)
Is there any function can do that? or any
There are more than one way to do it, and it would help if you
provided some sample data.
But here is an example for one way to do it:
examp.dat <- as.data.frame(matrix(sample(1:5, 100*6, replace=TRUE), ncol=6)
tmp.count <- apply(examp.dat, 1, function(x) sum(x>=3))
examp2.dat <-
What did you expect?
I'm assuming two plots (based on the subject) and side by side based on
the code (nrow =1)
But you are getting several graphs (facets) on the row and only expected
2?
What is in CHR? i.e. summary(tmp1$CHR)
I'm assuming its not a factor with 2 elements...?
> Hello,
Your best chance to get some interest is to adapt an existing package
such as linprog or lpSolve to use your algorithm. Then there will be
sufficient structure to allow R users and developers to see your
ideas working, even if they are not efficiently programmed. It's
always easier to start with
Estimado Manuel Mendoza
Habría que verlo, en forma rápida y con un posible error de mi parte, puede
radicar en que la primer ventana tiene los parámetros y las segundas no, o
dicho de otra forma, algún índice se crea en la primera y en las sucesivas
no se crea nuevamente este índice, se me ocurre
forgot to mention, the training and testing dataframes are composed of
4 IVs (one double numeric IV and three factor IVs) and one DV
(dichotomous factor, i.e. true or false).
The training dataframe consists of 48819 rows and test dataframe
consists of 24408 rows.
Thanks again.
Hi everyone. I'm using the kernlab ksvm function with the rbfdot
kernel for a binary classification problem and getting a strange
result back. The predictions seem to be very accurate judging by the
training results provided by the algorithm, but I'm unable to generate
a confusion
Hi Mozhgan,
This is pretty obscure. If you have entered all of the variables in
your attached list, You can start by looking at the summary of the
output. This may help:
https://stat.ethz.ch/R-manual/R-devel/library/mgcv/html/summary.gam.html
Jim
On Thu, Jun 11, 2020 at 4:47 AM Mozhgan
HI Jim
I run it like this:
Rscript --no-save Manhattan_plot.R
and just in case I added: stringsAsFactors=FALSE
so the script looks like this:
library(qqman)
results_log <- read.table("logistic_results_M3.assoc.logistic.C",
head=TRUE,stringsAsFactors=FALSE)
jpeg("M3.assoc.logistic.jpeg")
Hi Ana,
The problem may be that the JPEG device doesn't handle transparency.
Perhaps PNG?
Jim
On Thu, Jun 11, 2020 at 6:48 AM Ana Marija wrote:
>
> Hello,
>
> I have a data frame like this:
>
> > head(tmp1)
> CHR BP PoldPnew
> 1 1 785989 0.9521 0.09278
> 2 1 1130727 0.4750
Hi Ana,
I don't have the qqman package, but is your "P" column in
"M3.assoc.logistic.C" numeric or has it been read in as a factor?
Jim
On Thu, Jun 11, 2020 at 4:13 AM Ana Marija wrote:
>
> Hello,
>
> I do have a file like this:
> head M3.assoc.logistic.C
> CHR SNP BP P
> 1 1:785989:T:C 785989
?dput
We cannot tell how the columns are being stored in memory from your head()
output.
On June 10, 2020 1:36:11 PM PDT, Ana Marija wrote:
>Hello,
>
>I have a data frame like this:
>
>> head(tmp1)
> CHR BP PoldPnew
>1 1 785989 0.9521 0.09278
>2 1 1130727 0.4750 0.19010
>3 1
?do.call -- takes a list of arguments to a function
... as in
do.call(merge, yourlist) ## or similar perhaps
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
If that is an exact copy of the issued command than it should be
throwing an error related to the use of "smart quotes".
Furthermore it was not clear if that command was the proximate cause of
the error or perhaps it was encountered when you tried to load (rather
than install) the Rserve
Dear Sir/Madam,
Hope everyone is safe and sound. I appreciate your help a lot.
I am evaluating two Arabic subtitles of a humorous English scene and asked
263 participants (part) to evaluate the two subtitles (named Standard
Arabic, SA, and Egyptian Arabic, EA) via a questionnaire that asked them
Hello there,
I hope this email finds you well.
I’m just having a difficulty running RServe, which I’m trying to get
communicate with Tableau. I have this error message every time I try to do this:
Error: long vectors not supported yet: qap_encode.c:36
Fatal error: unable to initialize the JIT
Muy buenas, ¿a ver si hay alguien que sepa por qué en este loop, si hago,
p.e., i = 1 y corro las 2 filas de dentro, me abre una ventana y me hace el
PDP de frg, es decir, lo hace bien, pero si corro todo el loop me abre las
ventanas pero las deja vacías?
predictores <-
hi, I am trying to fuse (cbind, merge... NOT rbind) several dataframes with
different numbers of rows, all df included in a list, and using the code
extract shown below. The function merge() works well with two df but not
more than two...I have 85 dataframes to join in this way (85 df in the
Hello
thank you for reading email
please help me,
I have used mgcv package for GAM model.
To get contributions of each variable, Which code should I write in the R
software?
please see attached file
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE
Hello,
I'm analyzing a 6-question questionnaire. I need to exclude participants
who answered more than 4 questions with a grade of 3 or higher.
How do I write such a code?
Thank you!
[[alternative HTML version deleted]]
__
R-help@r-project.org
Friends:
i am a retired Professor -
not having any access to the resources (human/financial/business/whatever)
that may be required -
therefore i am seeking implementation of my algorithm 'spdspds' -
a novel algorithm for solving Linear Programming Problems with O(L^1.5)
computational complexity -
Hello,
I do have a file like this:
head M3.assoc.logistic.C
CHR SNP BP P
1 1:785989:T:C 785989 0.4544
1 1:785989:T:C 785989 0.689
1 1:1130727:A:C 1130727 0.05068
1 1:1130727:A:C 1130727 0.07381
1 1:1156131:T:C 1156131 0.6008
1 1:1156131:T:C 1156131 0.8685
...
And I don't have any "NA" or "inf"
I'm not sure this is really a statistical problem, in the sense of looking for
a convenient but arbitrary statistical function; to do it well is more of a
physicochemical modelling problem.
I can't give you an answer but maybe a direction I'd consider if I wanted to
take it seriously ...
You
On Wed, 10 Jun 2020, Jeff Newmiller wrote:
Fix your format specification?
?strptime
I have been trying to convert European short dates formatted as dd/mm/yy
into the ISO 8601 but the function as.Dates interprets them as American
ones (mm/dd/yy), thus I get:
Look at Hadley Wickham's
Fix your format specification?
?strptime
On June 10, 2020 1:20:01 AM PDT, Luigi Marongiu
wrote:
>Hello,
>I have been trying to convert European short dates formatted as
>dd/mm/yy into the ISO 8601 but the function as.Dates interprets them
>as American ones (mm/dd/yy), thus I get:
>
>```
Note: My warning was for "stepwise" regression, which is what *you wrote*,
not "segmented".
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 Wed, Jun 10, 2020
The package incidence provides the function fit_optim_split that
returns a predicted doubling time in days. Would this be serial
interval described in the epidemiology manuals?
Thank you
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,
Dear all
To make my problem more on topic I would like to ask about weird results
from segmented fit, despite of Bert's warning.
Here is my data
temp <- structure(list(V1 = c(0L, 15L, 30L, 45L, 60L, 75L, 90L, 105L, 120L,
135L, 150L, 165L, 180L, 195L, 210L, 225L, 240L, 255L, 270L, 285L, 300L,
maybe
isoDates <- as.Date(oriDates, format = "%d/%m/%y")
Heinz
Luigi Marongiu wrote/hat geschrieben on/am 10.06.2020 10:20:
Hello,
I have been trying to convert European short dates formatted as
dd/mm/yy into the ISO 8601 but the function as.Dates interprets them
as American ones (mm/dd/yy),
Luigi,
Try format = "%d/%m/%y"
Berend Hasselman
> On 10 Jun 2020, at 10:20, Luigi Marongiu wrote:
>
> ISO 8601
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Thank you!
On Wed, Jun 10, 2020 at 10:29 AM Daniel Nordlund wrote:
>
> On 6/10/2020 1:20 AM, Luigi Marongiu wrote:
> > isoDates = as.Date(oriDates, format = "%m/%d/%y")
>
> You need to use the format for European short dates.
>
> isoDates = as.Date(oriDates, format = "%d/%m/%y")
>
>
> Hope this
Hi!
Should it be:
as.Date(oriDates, format="%d/%m/%y") # See the order of %d and %m!!
This command seems to work for me, here the output:
[1] "2020-01-23" "2020-01-24" "2020-01-25" "2020-01-26" "2020-01-27"
"2020-01-28" "2020-01-29" "2020-01-30"
[9] "2020-01-31" "2020-02-01" "2020-02-02"
On 6/10/2020 1:20 AM, Luigi Marongiu wrote:
isoDates = as.Date(oriDates, format = "%m/%d/%y")
You need to use the format for European short dates.
isoDates = as.Date(oriDates, format = "%d/%m/%y")
Hope this is helpful,
Dan
--
Daniel Nordlund
Port Townsend, WA USA
On Wed, 10 Jun 2020 10:20:01 +0200
Luigi Marongiu wrote:
> the function as.Dates interprets them as American ones (mm/dd/yy)
> isoDates = as.Date(oriDates, format = "%m/%d/%y")
> How can I convert properly?
Pass the correct format? (Swap m and d in the format string.)
--
Best regards,
Ivan
Hello,
I have been trying to convert European short dates formatted as
dd/mm/yy into the ISO 8601 but the function as.Dates interprets them
as American ones (mm/dd/yy), thus I get:
```
oriDates = c("23/01/20", "24/01/20", "25/01/20", "26/01/20",
"27/01/20", "28/01/20", "29/01/20", "30/01/20",
Hi
External heating. Normally I would use TA instrumentation but for technical
reasons it is impossible. And other complicating factor is that temperature
rise is from beginning almost parabolic (it's derivation is straight line).
Therefore I started with double exponential fit, which is
40 matches
Mail list logo
