Re: [R] Strange variable names in factor regression

2024-05-09 Thread Michael Dewey
Comment in in-line below On 09/05/2024 13:09, Naresh Gurbuxani wrote: On converting character variables to ordered factors, regression result has strange names. Is it possible to obtain same variable names with and without intercept? Thanks, Naresh mydf <- data.frame(date =

Re: [R] Function with large nested list

2023-12-18 Thread Michael Dewey
Dear Emily Comment in-line On 18/12/2023 09:56, Emily Bakker wrote: Hello list, I want to make a large rulebased algorithm, to provide decision support for drug prescriptions. I have defined the algorithm in a function, with a for loop and many if statements. The structure should be as

Re: [R] [External] Re: Error in setwd(dir) when initializing R

2023-11-21 Thread Michael Dewey
Dear Ana When you installed R it gave you an option to put a shortcut on the desktop. If you did that then you can start the R GUI by clicking on it. You may need to customise it for your preferences. Specifically it starts R with an option cd-user-docs which I delete and I then alter the

Re: [R] Question about R software and output

2023-10-03 Thread Michael Dewey
Dear Charity Since your organisation is a member of King's Health Partners you might like to ask colleagues in KCL for local support. Michael On 02/10/2023 08:48, Ferguson Charity (CEMINFERGUSON) wrote: To whom it may concern, My understanding is that the R software is downloaded from a

Re: [R] How to fix this problem

2023-09-25 Thread Michael Dewey
It looks here as though the E coli column has commas in it so will be treated as character. Michael On 25/09/2023 15:45, avi.e.gr...@gmail.com wrote: David, This may just be the same as your earlier problem. When the type of a column is guessed by looking at the early entries, any

Re: [R] Odd result

2023-09-24 Thread Michael Dewey
Dear David To get the first 46 rows just do KurtzData[1:43,] However really you want to find out why it happened. It looks as though the .csv file you read has lots of blank lines at the end. I would open it in an editor to check that. Michael On 23/09/2023 23:55, Parkhurst, David wrote:

Re: [R] prop.trend.test

2023-09-07 Thread Michael Dewey
Dear Thomas Are you looking for more than smokers / patients Michael On 07/09/2023 14:23, Thomas Subia via R-help wrote: Colleagues  Consider smokers  <- c( 83, 90, 129, 70 ) patients <- c( 86, 93, 136, 82 )  prop.trend.test(smokers, patients)  Output: Chi-squared Test for

Re: [R] To UNSUBSCRIBE

2023-06-25 Thread Michael Dewey
You can find instructions on how to do this on the bottom of e-mailas from the list. On 25/06/2023 15:34, Siddharth Arun wrote: Please unsubscribe me from your mailing list. -- Michael http://www.dewey.myzen.co.uk/home.html __

Re: [R] Error on using Confint function to calculate 95% CI

2023-04-26 Thread Michael Dewey
I am afraid your post is more or less unreadable since you posted in HTML and this is a plain text list. It might also help if you gave more context like the full results of your model. There is a list dedicated to mixed models https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models which

Re: [R] nlme varFixed

2023-03-04 Thread Michael Dewey
Dear Franz Only attachments of specific types are allowed on the R lists and yours did not come through. It might be necessary to show your code as well. Have you thought of trying the mailing list dedicated to mixed models? https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models Michael

Re: [R] My forest plot is not fit to windows in R software

2023-01-11 Thread Michael Dewey
sity On Tuesday, January 10, 2023 at 09:46:40 PM GMT+3:30, Michael Dewey wrote: Dear Fahimeh It would help if you also told us what package you are using as there are several which perform forest plots. I assume it is meta. Your plot is cropped on three side as far as I can see. Wha

Re: [R] My forest plot is not fit to windows in R software

2023-01-10 Thread Michael Dewey
Dear Fahimeh It would help if you also told us what package you are using as there are several which perform forest plots. I assume it is meta. Your plot is cropped on three side as far as I can see. What parameters did you give to the call of the graphics device? What happens if you use a

Re: [R] Select dataframe row containing a digit

2022-11-30 Thread Michael Dewey
I suspect [[:digit:]] might have done what you want. Michael On 30/11/2022 13:04, Luigi Marongiu wrote: Thank you, I have been trying with [:digit:] but did not work. It worked with `df$val[grepl('[0-9]', df$val)] = "NUM"` On Wed, Nov 30, 2022 at 2:02 PM Ivan Krylov wrote: В Wed, 30 Nov

Re: [R] Sensitivity and Specificity

2022-10-24 Thread Michael Dewey
edict_testing[1, 1] : incorrect number of dimensions > specificity function (data, ...) {     UseMethod("specificity") } On Mon, Oct 24, 2022 at 10:45 AM Michael Dewey <mailto:li...@dewey.myzen.co.uk>> wrote: Rather hard to know without seeing what output you expected an

Re: [R] Sensitivity and Specificity

2022-10-24 Thread Michael Dewey
Rather hard to know without seeing what output you expected and what error message you got if any but did you mean to summarise your variable predict before doing anything with it? Michael On 24/10/2022 16:17, greg holly wrote: Hi all R-Help , After partitioning my data to testing and

Re: [R] Forestplot, grid graphics Viewplot grid.arange

2022-10-14 Thread Michael Dewey
There are several other CRAN packages which provide forest plots (see CRAN Task View for details) and they do not all use grip graphics which I think the forestplot package does. It might be worth swapping to one of them. Michael On 14/10/2022 04:34, Jim Lemon wrote: Hi Mary, I didn't see

Re: [R] High concurvity/ collinearity between time and temperature in GAM predicting deaths but low ACF. Does this matter?

2022-06-10 Thread Michael Dewey
days and 1 for heaping days, but was told that meant that I was assuming the effect was the same for all heaping days. If I coded with 12 or 14 levels, wouldn't that leave a lot of noise in the data? Jade On Thu, 9 Jun 2022 at 14:52, Michael Dewey wrote: Dear Jade Do you really need to fit a separate

Re: [R] High concurvity/ collinearity between time and temperature in GAM predicting deaths but low ACF. Does this matter?

2022-06-09 Thread Michael Dewey
Dear Jade Do you really need to fit a separate parameter for each heaping day? Can you not just make it a binary predictor or a categorical one with fewer levels, perhaps 14 (for heaping in each year) or 12 (for each calendar month). I have no idea whether that would help but it seems worth a

Re: [R] Shiny app Deployment not Working

2022-05-27 Thread Michael Dewey
I think people are going to need more details here. 1 where does GeneBook come from if not CRAN? 2 what is the app you are hoping to deploy, is it part of GeneBook? 3 what exactly do you mean by not showing? On 27/05/2022 07:29, ABHIRUP MOITRA wrote: I want to deploy a shiny app where I have

Re: [R] Restoration of "rite" package of R as R-script editor or the like

2022-05-23 Thread Michael Dewey
Notepad++ is fine as a light weight editor with syntax highlighting but the plug-in for R does not seem to have been updated for years. Michael On 23/05/2022 17:37, Rui Barradas wrote: Hello, I was also thinking about Notepad++. And about an editor no one mentinoed, vim. I've learned vi on

Re: [R] Issue about eRM package

2021-09-28 Thread Michael Dewey
I do not know that package so cannot help with that but the dataset did not come through. This mailing list is quite restrictive in what sorts of attachment it allows so I suggest trying something like a .csv file or a .txt file. Michael On 28/09/2021 10:37, 宋启发 wrote: I am using your

Re: [R] What are the pros and cons of the log.p parameter in (p|q)norm and similar?

2021-08-06 Thread Michael Dewey
11:53 AM Duncan Murdoch <mailto:murdoch.dun...@gmail.com>> wrote: On 03/08/2021 12:20 p.m., Michael Dewey wrote: > Short version > > Apart from the ability to work with values of p too small to be of much > practical use what are th

[R] What are the pros and cons of the log.p parameter in (p|q)norm and similar?

2021-08-03 Thread Michael Dewey
Short version Apart from the ability to work with values of p too small to be of much practical use what are the advantages and disadvantages of setting this to TRUE? Longer version I am contemplating upgrading various functions in one of my packages to use this and as far as I can see it

Re: [R] if statement and for loop question

2021-06-01 Thread Michael Dewey
Dear Kai When you ask again it is best to tell us what your input is and what output you were hoping for and what you actually got. If you can make a small data-set which shows all that then your post will be much more likely to get a helpful response. If you want to transfer the data-set to

Re: [R] calculating area of ellipse

2021-05-11 Thread Michael Dewey
Dear Stephen In that application the axes would be sensitivity and specificity (or their inverses) or some transformation of them like logits so the units would be the same. Whether the area has any scientific meaning I am not sure. Michael On 11/05/2021 15:20, Stephen Ellison wrote: In

Re: [R] readline in function call with space in prompt.

2021-02-09 Thread Michael Dewey
The function test as defined below by Jeremie works as I would have expected for me on Windows so I am unable to replicate the problem there. R version 4.0.3 (2020-10-10) Platform: x86_64-w64-mingw32/x64 (64-bit) Running under: Windows 10 x64 (build 19041) Matrix products: default locale: [1]

Re: [R] make new collumns with conditions

2021-01-25 Thread Michael Dewey
Dear Krissie I think you misunderstood Rui's response. He was generating some fake data to test the code not suggesting you rebuild your data frame. Michael On 25/01/2021 16:01, krissievdh wrote: Hi, Thanks for your response. I do get what you're doing. However, the table I sent is just a

Re: [R] make new collumns with conditions

2021-01-25 Thread Michael Dewey
Dear Krissie I think you may be looking for the aggregate command. Note that this is a plain text list so if you post in HTML we do not see what you see. In this case we did not see any red numbers. Michael On 25/01/2021 13:57, krissievdh wrote: Hi, I have a dataset (d_vigi)with this kind

Re: [R] [SPAM] Re: Different results on running Wilcoxon Rank Sum test in R and SPSS

2021-01-19 Thread Michael Dewey
January, 2021, 03:53:27 pm IST, Michael Dewey wrote: Unfortunately your data did not come through. Try using dput() and then pasting that into the body of your e-mail message. On 18/01/2021 17:26, bharat rawlley via R-help wrote: > Hello, > On running the Wilcoxon Rank Sum test in R a

Re: [R] Different results on running Wilcoxon Rank Sum test in R and SPSS

2021-01-19 Thread Michael Dewey
Unfortunately your data did not come through. Try using dput() and then pasting that into the body of your e-mail message. On 18/01/2021 17:26, bharat rawlley via R-help wrote: Hello, On running the Wilcoxon Rank Sum test in R and SPSS, I am getting the following discrepancies which I am

Re: [R] Help with connection issue for R (just joined, leading R for our agency)

2020-12-14 Thread Michael Dewey
Just to add to Petr's comment there are other basic editors with syntax highlighting like Notepad++ which are also OK if you want a fairly minimalist approach. Michael On 14/12/2020 08:16, PIKAL Petr wrote: Hallo Alejandra Although RStudio and ESS could help with some automation (each with

Re: [R] Find the ideal cluster

2020-12-12 Thread Michael Dewey
Dear Jovani If you cross-post on CrossValidated as well as here it is polite to give a link so people do not answer here when someone has already answered there, or vice versa. Michael On 12/12/2020 15:27, Jovani T. de Souza wrote: So, I and some other colleagues developed a hierarchical

Re: [R] scott-knot ESD effect size test

2020-12-02 Thread Michael Dewey
There seems to be a package on CRAN dedicated exclusively to this test. It is called ScottKnotESD rather unoriginally. Michael. On 02/12/2020 10:32, Neha gupta wrote: I have the following data from resample svm= svm$resample$RMSE nn= nn$resample$RMSE we perform the statistical tests like

Re: [R] WG: Packegs for mathematics

2020-11-17 Thread Michael Dewey
Dear Mihaela As well as the search tips you have been given you may want to look at the CRAN Task Views. There is one for Multivariate https://cran.r-project.org/view=Multivariate and one for NumericalMathematics https://cran.r-project.org/view=NumericalMathematics and several others which may

Re: [R] Error message when using "revtools" package

2020-11-05 Thread Michael Dewey
Dear John Your .bib file did not make it through the system as only a very limited set of attachments is supported. Try to make a minimal .bib file which triggers the problem by splitting your file in half, testing each half, repeat until you get the smallest file which triggers the error.

Re: [R] build a literature database

2020-11-04 Thread Michael Dewey
Dear John If you are doing a systematic review have you thought of investigating some of the packages in the CRAN Task View on MetaAnalysis like metagear or revtools? I have not used any of them but they claim to support that part of the process. Michael On 04/11/2020 09:22, John wrote:

Re: [R] Error: variable not found

2020-10-30 Thread Michael Dewey
Dear Hannah I think the problem is that attach() is not doing what you think it is. It does seem to make it easy to make mistakes. I would suggest switching to using with() instead or using the data = parameter to functions which support it. Michael On 30/10/2020 08:15, Hannah Van Impe

Re: [R] Need help in programming R on functional data

2020-10-29 Thread Michael Dewey
You have already asked this and people gave you a variety of answers. Just asking again without clarifying why those answers did not help you is not going to solve your problem. Tell us what you tried and why it failed might help. Michael On 29/10/2020 07:50, Ablaye Ngalaba wrote: Hello, I

Re: [R] Dear R experts, I have a question about meta-analysis

2020-10-28 Thread Michael Dewey
Dear Zhang There is a mailing list dedicated to meta-analysis in R where your question may get more attention. Before posting it would be a good idea to look at the archives as this issue comes up there repeatedly. https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis// For a link to

Re: [R] Stats help for dissertation project

2020-09-19 Thread Michael Dewey
Dear Raija This list is primarily for R programming questions so most of your post is off-topic here. If you are registered for a degree presumably someone is paying a fee to your institution and someone there is being paid to supervise your project so I would have thought they would be the

Re: [R] question including crossover trials in meta-analysis

2020-09-15 Thread Michael Dewey
Dear Vera In addition to what you already have you might like to know about the mailing list specifically dedicated to meta-analysis in R. https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis// You might like to search the archives first as this sort of issue does come up there. You

Re: [R] truth[(truth[, 1]=="G3" & truth[, 2]=="G2") | (truth[, 1]=="G2" & truth[, 2]=="G3"), 3]<-1

2020-09-06 Thread Michael Dewey
I am afraid this is completely unreadable because you posted in HTML ad this is a plain text list. Best to resend it having set your mailer to send plain text as HTML gets mangled here. Michael On 06/09/2020 10:58, Hesham A. AL-bukhaiti via R-help wrote:

Re: [R] hist from a list

2020-07-31 Thread Michael Dewey
Dear Pedro Some comments in-line On 30/07/2020 21:16, Pedro páramo wrote: Hi all, I attach my code, the think is I want to make a bar plot the last variable called "bwchist" so the X axis are "Accion" and the y axis are "reval" values. I have prove class(bwchist) and says dataframe but

Re: [R] Fortune nomination .... Re: Looping through a dataframe

2020-07-25 Thread Michael Dewey
Dear Chris Just send it to the maintainer Achim Zeileis to whom I have cc'ed this in case he is not reading the list at the moment. Michael On 25/07/2020 06:21, Chris Evans wrote: I really don't want to put too many Emails to the list but I had the same reaction to that lovely line ... so

Re: [R] Question about citing R packages in an academic paper

2020-07-23 Thread Michael Dewey
Dear Medhi It is good that you are going to cite the packages properly. I do not think it matters too much whether they are in the article itself or in the supplementary material, the important thing is that they are there. You do not have any legal obligation to cite them, as far as I know

Re: [R] Reading help functions

2020-07-23 Thread Michael Dewey
Dear Pedro If you prefer reading in Spanish which from your name may be the case then if you go to CRAN, click on Documentation | Contributed you will find some documentation in Spanish including several aimed at beginners. There are some in Portuguese as well if I have mis-interpreted your

Re: [R] how to calculate odd ratios with R?

2020-07-06 Thread Michael Dewey
Dear Luigi You could try the epitools package which gives a large number of ways of doing this. I would have thought that using Wald intervals for the log odds ration was not optimal with small frequencies. Michael On 06/07/2020 14:01, Luigi Marongiu wrote: Hello, Is it possible to

Re: [R] how to summarize results from studies?

2020-07-01 Thread Michael Dewey
Dear Frederik There is also a mailing list dedicated to meta-analysis in R https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis// Michael On 01/07/2020 16:40, Marc Schwartz via R-help wrote: Hi, It sounds like you will want to engage in a meta-analysis. There is a CRAN task view

Re: [R] Error message in meta-analysis package Metafor-weights =""

2020-06-23 Thread Michael Dewey
The two packages both define a function forest. When you load the second package R will have told you that it was masking the definition of forest from the first package. If you had loaded them in the other order it would have masked the other one. In fact it masked seven functions in total

Re: [R] comparing variances/distributions

2020-06-18 Thread Michael Dewey
Dear Ana This really depends on your scientific question. The two techniques you have shown do different things and there must be many more which could be applied. Michael On 17/06/2020 20:57, Ana Marija wrote: Hello, I have p values from two distributions, Pold and Pnew head(m) CHR

Re: [R] Convergence in Monte Carlo Simulation

2020-06-14 Thread Michael Dewey
, "Michael Dewey" wrote: I am not 100% clear what your code is doing as it gets a bit wangled as you posted in HTML but here are a couple of thoughts. You need to set the seed outside any loops so it happens once and for all. I would test after trycatc

Re: [R] Convergence in Monte Carlo Simulation

2020-06-14 Thread Michael Dewey
I am not 100% clear what your code is doing as it gets a bit wangled as you posted in HTML but here are a couple of thoughts. You need to set the seed outside any loops so it happens once and for all. I would test after trycatch and keep a separate count of failures and successes as the

Re: [R] Error in ordinal mixed effects

2020-06-09 Thread Michael Dewey
Dear Susana Without your dat it is hard to say (and it would have helped to know where mixor() comes from) but this almost always means that ne of your parameters to the call is not what you thought it was so trying str(res) might be enlightening. Also I do not see anywhere in your example

Re: [R] a question of etiquette

2020-06-01 Thread Michael Dewey
You might get better answers on the list dedicated to package development r-pkg-devel This may have already been discussed there so a quick look at the archive might also help you. On 01/06/2020 17:34, Adelchi Azzalini wrote: The new version of a package which I maintain will include a new

Re: [R] text annotation on Manhattn plot in qqman

2020-05-20 Thread Michael Dewey
ot;P",cex = 0.5,annotatePval = 0.0001) but I am getting this error: Error in textxy(topSNPs$pos, -log10(topSNPs$P), offset = 0.625, labs = topSNPs$SNP, : formal argument "cex" matched by multiple actual arguments Do you by any chance know how to do this? Cheers Ana On Wed, May 20, 2

Re: [R] text annotation on Manhattn plot in qqman

2020-05-20 Thread Michael Dewey
ing parameters. You may also need to play with the sep =, and collapse = parameters to paste() to get the precise layout you want. Michael On 19/05/2020 17:21, Ana Marija wrote: Hi Michael, can you please send me code how that would be done? Thanks Ana On Tue, May 19, 2020 at 11:18 AM Michael Dewey

Re: [R] text annotation on Manhattn plot in qqman

2020-05-19 Thread Michael Dewey
Dear Ana Perhaps paste together SNP and GENE using paste() and then supply that as the snp parameter. Michael On 19/05/2020 17:12, Ana Marija wrote: Hello, I am making manhattan plot with: library(qqman) manhattan(a, chr="CHR", bp="BP", snp="SNP", p="P",annotatePval = 0.0001) and I would

Re: [R] NA command in a 'for' loop

2020-04-20 Thread Michael Dewey
Just a thought Helen but is x being treated as a real and what you think are zero and are printed as zero are in fact some very small number? If so you need to alter your test appropriately. Michael On 20/04/2020 17:25, Helen Sawaya wrote: Thank you for your reply. I tried d[] <- lapply(d,

Re: [R] how to get all strings in a data frame that start with a particular string

2020-04-10 Thread Michael Dewey
Dear Ana Would it not be possible to use grep instead of grepl and get the values using the value = TRUE parameter? Michael On 10/04/2020 17:15, Ana Marija wrote: Hello, Hello, I have a data frame (tot) with about 2000 columns. How can I extract from it all strings that start with E14? I

Re: [R] Same results but different functions ?

2020-03-23 Thread Michael Dewey
The documentation suggests that the rlm method for a formula does not have psi as a parameter. Perhaps try using the method for a matrix x and a vector y. Michael On 23/03/2020 12:39, varin sacha via R-help wrote: Dear R-experts, The rlm command in the MASS package command implements

Re: [R] how do I add text lables on QQ plot

2020-03-12 Thread Michael Dewey
Dear Ana You can specify the first three parameters to text() as vectors so it is all done in one call. That may or may not answer your question. Michael On 12/03/2020 14:08, Ana Marija wrote: HI David, thank you for getting back to me. Is there is a way for qq() to pick up text label

Re: [R] R CMD check report error not finding function

2020-03-05 Thread Michael Dewey
Dear Ruben I do not think your export pattern matches .foo so that may be the problem. Michael On 05/03/2020 09:10, Ruben wrote: Dear R-help listers, I am creating a package for a friend using his scripts and I run into a problem when doing the check. Inside a macro created with

Re: [R] testing my package : unstated dependency to self in package tests

2020-02-16 Thread Michael Dewey
When something similar happened to me I found it went away when I added Suggests: to the DESCRIPTION file. Whether this will work for you I have no idea. Michael On 16/02/2020 11:03, Servet Ahmet Çizmeli wrote: I am updating my CRAN package geoSpectral. I get the following Warning during R

Re: [R] Information about pairwise Wilcoxon Sum Rank Tests

2020-02-06 Thread Michael Dewey
Dear Gabriela Apologies if you have already tried this but: 1 - you can see the code which it uses by typing pairwise.wilcox.test at the command line. 2 - to find you more about the method of adjustment you need to follow the documentation for p.adjust 3 - you may also need to look into

Re: [R] R package for meta-analysis from z scores

2020-01-20 Thread Michael Dewey
Dear James Your question really boil down to whether you can estimate tau^2, the between study variance of the effect sizes, if you only have p-values. As far as I can see the answer has to be no. Michael On 16/01/2020 13:10, james poweraid wrote: Hello, I have a set of Z scores from an

Re: [R] "In sqrt(VS) : NaNs produced"

2020-01-20 Thread Michael Dewey
Your script and data were stripped so we are none the wiser I am afraid. You need to embed the script in your post and give a minimal data-set which exhibits the problem using dput() and embed that in the post too. Michael On 19/01/2020 08:25, Atul Saini wrote: Hello R, I am

Re: [R] all the MAE metric values are missing (Error message)

2019-12-23 Thread Michael Dewey
What Jim is alluding to is that sometimes in the process of reading in data a small typo can mean that what was intended to be a numeric variable is read in as a factor. So he was suggesting that you double check that this has not happened to you. Michael On 23/12/2019 11:45, Neha gupta

Re: [R] Orthogonal polynomials used by R

2019-11-28 Thread Michael Dewey
Dear Ashim As John said your two examples give the same model to within rounding error so it is not clear what you see the problem as being. You can always remove some of the correlation by subtracting out a large constant from x before you use poly() on it. Michael On 28/11/2019 16:02,

Re: [R] [r] How to obtain the correlation coefficients between two variables in the R object

2019-11-25 Thread Michael Dewey
To get multiple plots on the same physical page try ?layout ot look through the various option under ?par I do not know enough about DCC models to answer how you get correlations. Michael On 25/11/2019 10:34, Subhamitra Patra wrote: Dear R-users, I am estimating the DCC-model by using the

Re: [R] QQ plot

2019-11-13 Thread Michael Dewey
Dear Ana As others have commented this is getting a bit off-topic but here are some hints. It is helpful to distinguish two sorts of plot: archival plots and impact plots. If you want to have an impact plot which gives you a picture but possibly at the cost of completeness and accuracy

Re: [R] how to calculate multiple meta p values

2019-10-30 Thread Michael Dewey
umz(x[2:4], weights = x[5:7]) : Must have at least two valid p values Please advise, Ana On Sun, Oct 27, 2019 at 9:49 AM Michael Dewey wrote: Dear Ana There must be several ways of doing this but see below for an idea with comments in-line. On 26/10/2019 00:31, Ana Marija wrote: Hello, I would li

Re: [R] how to calculate multiple meta p values

2019-10-27 Thread Michael Dewey
Dear Ana There must be several ways of doing this but see below for an idea with comments in-line. On 26/10/2019 00:31, Ana Marija wrote: Hello, I would like to use this package metap to calculate multiple o values I have my data frame with 3 p values head(tt) RSG

Re: [R] negative vector length when merging data frames

2019-10-24 Thread Michael Dewey
Dear Ana Since this appears to be genetics data have you thought of looking at Bioconductor for help? I do not use genetic data-sets but people there must use big files every day three times before breakfast. Michael On 24/10/2019 00:33, Ana Marija wrote: thanks but I would need solution

Re: [R] Creating a before-and-after variable in R

2019-10-03 Thread Michael Dewey
Dear Faradj I am afraid your post is unreadable since this is a plain text list and you sent in HTML. Michael On 03/10/2019 12:17, Faradj Koliev wrote: Dear R-users, I need an urgent help with the following: I have a country-year data covering the period 1982 - 2013. I want to assess how

Re: [R] I wrote the document for extended usage of my pyramid package

2019-08-29 Thread Michael Dewey
Dear Minato That is a nice idea, but why not make it a vignette of your package on CRAN so it is immediately accessible to anyone using the package? Michael On 29/08/2019 07:15, Minato Nakazawa wrote: Dear R-users, I have developed the package pyramid, which has been available from CRAN,

Re: [R] parsing files with "\" character

2019-08-27 Thread Michael Dewey
Dear April Can you show us an example of what you are trying to do and how it fails? There are rules about backspaces but I find that if one backspace does not work try two, three, four until it works. It would be better to understand the rules but life is short. Michael On 27/08/2019

Re: [R] Code modification for post-hoc power

2019-08-27 Thread Michael Dewey
ncluding as a parameter something which was not of the type which the program expected. I would suggest as a first step not using attach() at all and instead using data=geefile in the call to gee or investigating the with() command if it does not allow a data parameter. On 27/08/2019 12:25, Michael

Re: [R] Code modification for post-hoc power

2019-08-27 Thread Michael Dewey
with a confusing error message sent earlier I don't understand. Best, Anne -Message d'origine- De : Michael Dewey [mailto:li...@dewey.myzen.co.uk] Envoyé : lundi, 26 août 2019 18:29 À : Marc Schwartz ; CHATTON Anne Cc : R-help Objet : Re: [R] Code modification for post-hoc power Dear Anne

Re: [R] Code modification for post-hoc power

2019-08-26 Thread Michael Dewey
Dear Anne In addition to Marc's comments if you are forced to do this then, assuming your package computes sample size from power then just feed it a range of powers and find the one for which it calculates the sample size you had. There is a more elegant way to do this using uniroot but

Re: [R] Modelos neural netwok

2019-07-28 Thread Michael Dewey
Hay una lista en castellano https://stat.ethz.ch/mailman/listinfo/r-help-es Michael On 27/07/2019 18:41, Agustín Alonso Rodriguez wrote: Estimados miembros de la Lista: Me dirijo a vosotros para que, si pod�is, ilumin�is mi confusi�n. Estoy trabajando con los modelos neural netwok del

Re: [R] How to combine the results of different survey?

2019-06-04 Thread Michael Dewey
Dear Damien Comments in line On 03/06/2019 15:50, Damien Zurkinden wrote: I have results of different surveys (years) for different species that I called NAMA. I want to combining the results of those different surveys to produce one overall result by using meta-analysis I tried to

Re: [R] Problema estrazione file in formato "large STFDF" in formato "ASCII" compatibile con GIS

2019-05-22 Thread Michael Dewey
Dear Agostino I am afraid this list is an English language one but I suggest that as well as trying to translate your message you post on the list which specialises in this sort of thing https://stat.ethz.ch/mailman/listinfo/r-sig-geo where they may be better able to help you. Michael On

Re: [R] Survuval Anaysis

2019-05-02 Thread Michael Dewey
Without more details it is hard to answer but it is suspicious that it is dropping one of your predictors and the standard errors of the other are very large. This suggests you should investigate the joint distribution of your predictors and the events. Michael On 02/05/2019 13:37, Haddison

Re: [R] R 3.4.4 is released

2019-04-18 Thread Michael Dewey
Dear Stephen Questions about RStudio ae best asked in their help forums but I would definitely install the latest version of R and RStudio and do update.packages before asking Michael On 18/04/2019 13:19, Stephen Muldoon wrote: Hi, I am new to R studio. If my R studio continually asks to

Re: [R] lm fails on some large input

2019-04-18 Thread Michael Dewey
Perhaps subtract 1506705766 from y? Saying some other software does it well implies you know what the _correct_ answer is here but I would question what that means with this sort of data-set. On 17/04/2019 07:26, Dingyuan Wang wrote: Hi, This input doesn't have any interesting properties

Re: [R] combining data.frames with is.na & match (), two questions

2019-04-18 Thread Michael Dewey
Dear Drake See in-line comments On 18/04/2019 00:24, Drake Gossi wrote: Hello everyone, I'm working through this book, *Humanities Data in R* (Arnold & Tilton), and I'm just having trouble understanding this maneuver. In sum, I'm trying to combine data in two different data.frames. This

Re: [R] Unable to load built In Datasets From Stat2Data package.

2019-04-15 Thread Michael Dewey
Did you mean to do data(ICU) and then use ICU and so on On 15/04/2019 16:06, Anaanthan Pillai wrote: Hi, I’ve been trying to load datasets from Stat2Data, but couldn’t proceed to do so. I can install the package though, but the datasets are not being loaded. Is there any error that I do?

Re: [R] Problem with forest plot in 'meta' after update and byvar command

2019-03-22 Thread Michael Dewey
Dear Antonello There is a mailing list dedicated to meta-analysis and R. People from the team behind the package you mention do post there. https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis// Remember to register first please. Michael On 21/03/2019 22:04, Antonello Preti wrote:

Re: [R] plot.xmean.ordinaly vs plot() in package "rms"

2019-03-20 Thread Michael Dewey
Dear Kim See inline On 19/03/2019 22:29, Kim Jacobsen wrote: Mailing list now included (apologies, first time I post anything so not quite sure how it works). You are quite right, it was a typo. I meant to write that plot.xmean.ordinaly(). So please let me correct my last statement: the

Re: [R] Installing package into...

2019-03-12 Thread Michael Dewey
Dear Steven If you use install,packages() from within R there is a lib argument. I do not use RStudio so not sure how it works there. Michael On 12/03/2019 15:05, Steven Yen wrote: I install package using either the command line or Tools -> Install packages... (in RStudio) and get a

Re: [R] add points to lattice cloud plot (3D scatter)

2019-02-28 Thread Michael Dewey
Not sure whether this helps but try library(lattice) ?llines Note that is indeed a double ll at the start Michael On 28/02/2019 10:39, Luigi Marongiu wrote: Dear all, is it possible to add points to a lattice cloud plot (3D scatter)? I can plot the main data, but what if I wanted to add

Re: [R] Remove cases with -Inf from a data frame

2019-02-16 Thread Michael Dewey
Dear Abou Depends on exact details of your variables but ?is.finite Gives you the basic tool. On 16/02/2019 15:05, AbouEl-Makarim Aboueissa wrote: Dear All: good morning I have a log-transformed data frame with some *-Inf* data values. *my question: *how to remove all rows with *-Inf*

Re: [R] very slow code execution

2019-02-07 Thread Michael Dewey
umn in the table is AICc[,1] to store i and second column for j and so on. Maybe i am mistaken and this won't give me what i want, the code been running for 5 hours now. So i am waiting On Wed, Feb 6, 2019 at 4:59 PM Michael Dewey <mailto:li...@dewey.myzen.co.uk>> wrote: This is n

Re: [R] very slow code execution

2019-02-06 Thread Michael Dewey
This is not an answer to your speed problem but are your assignments to AICc[,1] and so on doing what you hope they are doing? Michael On 06/02/2019 12:03, salah maadawy wrote: i am a beginner regarding R but i am trying to do a simple thing, but it is taking too much time and i am asking if

Re: [R] A priori contrast for binomial GLM

2019-01-18 Thread Michael Dewey
Dear Rula That is really a statistical question not one for this list but the answer is that the fact that they are all zero for that category explains it. Search on-line for separation for more details. Michael On 18/01/2019 09:56, Rula Domínguez wrote: Hello to everyone, after much

Re: [R] Warning message: NAs introduced by coercion

2019-01-08 Thread Michael Dewey
Dear Meriam Your csv file did not come through as attachments are stripped unless of certain types and you post is very hard to read since you are posting in HTML. Try renaming the file to .txt and set your mailer to send plain text then people may be able to help you better. Michael

Re: [R] Mailinglist

2019-01-06 Thread Michael Dewey
Dear Rachel Not sure if this is going to help but if it is a csv file then read.csv() is your friend. Read the help first in case you need to specify what is being used for the decimal point and the separator as if it is from the Netherlands they may not be the default settings. michael On

Re: [R] Problem with Kruskal–Wallis test

2018-12-23 Thread Michael Dewey
answer. So, I'm not an expert in R and statistics, how can I create this interval of confidence of groups? Thanks Gc Il giorno sab 22 dic 2018, 13:34 Michael Dewey <mailto:li...@dewey.myzen.co.uk>> ha scritto: Dear Giuseppe If I understand you correctly you have a very large sa

Re: [R] Glitch in Kruskal-Wallis test?

2018-12-22 Thread Michael Dewey
Dear Jenny What exactly do you think you are testing here? You are telling K-W you have seven groups each with a single value which is not the usual situation for K-W. Michael On 22/12/2018 04:58, Jenny Liu wrote: Hi everyone, I have been running a K-W test with the attached data,

Re: [R] Problem with Kruskal–Wallis test

2018-12-22 Thread Michael Dewey
Dear Giuseppe If I understand you correctly you have a very large sample size so it is not surprising that you get very small p-values. Eevn a scientifically uninteresting difference can become statistically significant with large samples. You probably need to define a metric for meaningful

  1   2   3   4   5   6   >