[R] bnlearn and cpquery

2017-07-12 Thread Ross Chapman
Hi all I have built a Bayesian network using discrete data using the bnlearn package. When I try to run the cpquery function on this data it returns NaN for some some cases. Running the cpquery in debug mode for such a case (n=10^5, method="lw") creates the following output:

[R] Quadratic function with interaction terms for the PLS fitting model?

2017-07-12 Thread Ng, Kelvin Sai-cheong
Dear all, I am using the pls package of R to perform partial least square on a set of multivariate data. Instead of fitting a linear model, I want to fit my data with a quadratic function with interaction terms. But I am not sure how. I will use an example to illustrate my problem: Following

Re: [R] about plotting a special case

2017-07-12 Thread Jim Lemon
values A and B have different ranges. Therefore you would want horizontal error bars for the B mean. The dispersion function doesn't do horizontal error bars, so you can use plotCI also in plotrix. Here is an example with two databases: DF2<-read.table(text="A B C 62 22 54 69 24 55

Re: [R] about plotting a special case

2017-07-12 Thread lily li
Thanks, Jim. The code works, but I don't understand why you use q1090 <- quantile(DF1$B, probs=c()), rather than DF1$A? Also, how to add a legend for both points DF1 and DF2? On Wed, Jul 12, 2017 at 8:25 PM, Jim Lemon wrote: > Hi lily, > Here is the first plot: > >

Re: [R] about plotting a special case

2017-07-12 Thread Jim Lemon
Hi lily, Here is the first plot: plot(DF1$A,DF1$B,pch=19,col="red") meanA<-mean(DF1$A) meanB<-mean(DF1$B) points(meanA,meanB,pch=18,col="red") q1090<-quantile(DF1$B,probs=c(0.1,0.9)) library(plotrix) dispersion(meanA,meanB,q1090[2],q1090[1], intervals=FALSE,col="red") The same code will work

[R] about plotting a special case

2017-07-12 Thread lily li
Hi R users, I have a question about plotting. There is the dataframe below, while each row represents a record. If I want to plot on a A-B plot, i.e., x-axis represents A, while y-axis represents B values. However, I want to plot the mean value from records 1-10 as one point, while the 10th and

Re: [R] How to make a figure plotting p-values by range of different adjustment values?

2017-07-12 Thread Jim Lemon
Hi Kirsten, Perhaps this will help: set.seed(3) kmdf<-data.frame(group=rep(1:4,each=20), prop=c(runif(20,0.25,1),runif(20,0.2,0.92), runif(20,0.15,0.84),runif(20,0.1,0.77))) km.glm<-glm(prop~group,kmdf,family=quasibinomial(link="logit")) summary(km.glm) pval<-0.00845 padjs<-NA npadj<-1 # assume

[R] Question on Simultaneous Equations & Forecasting

2017-07-12 Thread OseiBonsu, Frances
Hello, I have estimated a simultaneous equation model (similar to Klein's model) in R using the system.fit package. I have an identity equation, along with three other equations. Do you know how to explicitly identify the identity equation in R? I am also trying to forecast the dependent

[R] metRology package

2017-07-12 Thread mark . hogue
I'm having trouble with a simple application with metRology. I need to estimate the uncertainty of the density thickness of seven sheets of film. This is calculated from measurements of mass, length and width of rectangular samples of film. It's not too hard to calculate the whole thing with a

[R] How to make a figure plotting p-values by range of different adjustment values?

2017-07-12 Thread Kirsten Morehouse
Hi all, Thank you for taking the time to read my message. I'm trying to make a figure that plots p-values by a range of different adjustment values. (Using the **logit** function in package **car**) My Statistical analyses were conducted on probability estimates ranging from 0% to 100%. As it's

Re: [R-es] consulta por varTypes de sqlSave en RODBC

2017-07-12 Thread Javier Marcuzzi
Estimado José Ramirez Costa Yo tengo un archivo donde guarde lo que realizaba cuándo aprendía, habría que ver si actualmente es así o hay cambios, por lo que usted dice creo que usa sql server, si no es así busco el archivo para otras bases de datos. Si es sqlserver se podría realizar con la

Re: [R] Extracting sentences with combinations of target words/terms from cancer patient text medical records

2017-07-12 Thread Robert McGehee
Hi Paul, Sounds like you have your answer, but for fun I thought I'd try solving your problem using only a regular expression query and base R. I believe this works: > txt <- "Patient had stage IV breast cancer. Nothing matches this sentence. > Metastatic and breast match this sentence. French

[R-es] consulta por varTypes de sqlSave en RODBC

2017-07-12 Thread Jose Ramirez Costa
Buen dia compañerxs de R, les traigo una pregunta bastante especifica. En el trabajo me estan requiriendo el uso de SQL, asi q estoy aprendiendo a usarlo desde R, con el paquete RODBC. Mi problema puntual es que al crear una tabla en el servidor SQL a partir de un data frame en R (usando la

Re: [R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread John Kane via R-help
What is cbplatte? Where do we find it? On Wednesday, July 12, 2017, 6:20:38 AM EDT, Shanu Singh wrote: Please solve the problem. on which R Version works cbplatte? Regards Shanu     [[alternative HTML version deleted]]

Re: [R] Extracting sentences with combinations of target words/terms from cancer patient text medical records

2017-07-12 Thread Paul Miller via R-help
Hi Bert, Thanks for your reply. It appears that I didn't replace the variable name "sampletxt" with the argument "x" in my function. I've corrected that and now my code seems to be working fine. Paul From: Bert Gunter Cc: R-help

Re: [R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread Rui Barradas
Hello, Google is your friend. The only reference to cbPalette I've found is in http://www.cookbook-r.com/Graphs/Colors_%28ggplot2%29/ and it is defined as: A colorblind-friendly palette These are color-blind-friendly palettes, one with gray, and one with black. To use with ggplot2, it is

Re: [R] Importing Big data to R

2017-07-12 Thread Peter Dalgaard
> On 12 Jul 2017, at 08:44 , Mangalani Peter Makananisa > wrote: > > Dear R-Gurus, > > I am trying to read in data with 24,349,113 rows to R-3.3.3 (64 bit) and have > used the library "data.table" and It managed to read 23,347,070 rows and > the remainder was

Re: [R] Importing Big data to R

2017-07-12 Thread jim holtman
A little more information would be useful. Why did it stop? Was there an error message? Can you show the commands/console log of what you did. Provide information on how much memory your computer has on it. When the operation completed, how much memory was used. An important aspect is how many

Re: [R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread Uwe Ligges
On 12.07.2017 12:35, Shanu Singh wrote: Hi John, Thanks for reply ! Sorry It was my mistake. It is *cbPalette. *Basically its a R package. Come on. - Where do you want to get thsi from? - What have you tried to do so? Please re-read the posting guide that helps to ask good questions.

Re: [R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread Shanu Singh
Thanks for reply! sorry my my mistake It is *cbPalette* *Shanu* On Wed, Jul 12, 2017 at 4:04 PM, Uwe Ligges wrote: > > > On 12.07.2017 12:03, Shanu Singh wrote: > >> Please solve the problem. on which R Version works cbplatte? >> > > What is cbPlatte or

Re: [R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread Shanu Singh
Hi John, Thanks for reply ! Sorry It was my mistake. It is *cbPalette. *Basically its a R package. Shanu On Wed, Jul 12, 2017 at 3:57 PM, John Kane wrote: > What is cbplatte? Where do we find it? > > > > On Wednesday, July 12, 2017, 6:20:38 AM EDT, Shanu Singh < >

Re: [R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread Uwe Ligges
On 12.07.2017 12:03, Shanu Singh wrote: Please solve the problem. on which R Version works cbplatte? What is cbPlatte or cbplatte (case matters). Best, Uwe Ligges Regards Shanu [[alternative HTML version deleted]] __

Re: [R] devtools::check "hangs"

2017-07-12 Thread Uwe Ligges
On 12.07.2017 10:31, Sigbert Klinke wrote: Hi, checking my package with check from devtools leads to the output Functions with \usage entries need to have the appropriate \alias entries, and all their arguments documented. The \usage entries must correspond to syntactically valid R code. See

[R] Problem for installing cbPlatte on R Version 3.4.1

2017-07-12 Thread Shanu Singh
Please solve the problem. on which R Version works cbplatte? Regards Shanu [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] How to formulate quadratic function with interaction terms for the PLS fitting model?

2017-07-12 Thread Ng, Kelvin Sai-cheong
Dear all, I am using the pls package of R to perform partial least square on a set of multivariate data. Instead of fitting a linear model, I want to fit my data with a quadratic function with interaction terms. But I am not sure how. I will use an example to illustrate my problem: Following

[R] Importing Big data to R

2017-07-12 Thread Mangalani Peter Makananisa
Dear R-Gurus, I am trying to read in data with 24,349,113 rows to R-3.3.3 (64 bit) and have used the library "data.table" and It managed to read 23,347,070 rows and the remainder was 2,043 rows only. Could you please advise me as to which library/R-commands is suitable to read the full

[R] devtools::check "hangs"

2017-07-12 Thread Sigbert Klinke
Hi, checking my package with check from devtools leads to the output Functions with \usage entries need to have the appropriate \alias entries, and all their arguments documented. The \usage entries must correspond to syntactically valid R code. See chapter ‘Writing R documentation files’ in

Re: [R] submitting R scripts with command_line_arguments to PBS HPC clusters

2017-07-12 Thread Bogdan Tanasa
Dear Peter, that is very very helpful, many thanks for your suggestions ;) ! On Tue, Jul 11, 2017 at 11:34 PM, Anthoni, Peter (IMK) < peter.anth...@kit.edu> wrote: > Hi, > > The problem is most likely, you need to call a R CMD BATCH with your > arguments and the R-script inside of a shell script

Re: [R] submitting R scripts with command_line_arguments to PBS HPC clusters

2017-07-12 Thread Anthoni, Peter (IMK)
Hi, The problem is most likely, you need to call a R CMD BATCH with your arguments and the R-script inside of a shell script that you submit to your qsub. Unfortunately we don't use qsub anymore so can't test it, but it should be as follows: R-script eg. test.R: > ##First read in the arguments