similar to: pheno package (PR#10674)

Dear R-people, I tried to use a command seqMK(x) of the pheno package. However, the installation of the package fails with following error message: Error in library(pkg, character.only = TRUE) : 'pheno' is not a valid package -- installed < 2.0.0? What is the matter? Can anybody help me? I installed the latest version of R. Cheers, Anja --

Dear R users, This command works (calling a programm -called whap- with file specifiers etc.): >system('cmd /c "c:\\pheno\\whap --file c:\\pheno\\smri --alt 1 --perm 500"', intern=TRUE) Now I need to call it from a loop to replace the "1" by different number, however I get lost using the quotes: I tried numerous versions of: >i<-1

Hello, When programming it is better to use dat[["variable"]] than dat$variable. So your code could be pfas.pheno[[cat.var]] <- NA pfas.pheno[[cat.var]][pfas.pheno[,i] <= quantile(pfas.pheno[,i],0.25, na.rm =T)] <- 0 etc. Untested. Hope this helps, Rui Barradas On 4/19/2018 7:20 PM, Ding, Yuan Chun wrote: > Hi All, > > I want to create a categorical variable,

Hi All, I want to create a categorical variable, cat.pfoa, in the file of pfas.pheno (a data frame) based on log2pfoa values. I can do it using the following code. pfas.pheno <-within(pfas.pheno, {cat.pfoa<-NA cat.pfoa[pfas.pheno$log2pfoa <=quantile(pfas.pheno$log2pfoa,0.25, na.rm =T)]<-0 cat.pfoa[pfas.pheno$log2pfoa >=quantile(pfas.pheno$log2pfoa,0.75, na.rm =T)]<-2

> On Apr 19, 2018, at 11:20 AM, Ding, Yuan Chun <ycding at coh.org> wrote: > > Hi All, > > I want to create a categorical variable, cat.pfoa, in the file of pfas.pheno (a data frame) based on log2pfoa values. I can do it using the following code. > > pfas.pheno <-within(pfas.pheno, {cat.pfoa<-NA > cat.pfoa[pfas.pheno$log2pfoa

> On Apr 19, 2018, at 1:22 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > >> On Apr 19, 2018, at 11:20 AM, Ding, Yuan Chun <ycding at coh.org> wrote: >> >> Hi All, >> >> I want to create a categorical variable, cat.pfoa, in the file of pfas.pheno (a data frame) based on log2pfoa values. I can do it using the following code.

I am trying to find a convenient way to control line colors when printing from a for loop using the lines command. Right now I have solved this by creating a colors vector that is refered to in the loop with index. However, the colors choosen here are just 1,2,3,4,5... I would like to get colors from the col = rainbow(x) that you can use in plot() and set the to be my number of lines (I think

Dear List Members, I have encountered two problems when using the step function to select models. To better illustrate the problems, attached is an R image which includes the objects needed to run the code attached. lm.data.frame have factor variables with 3 levels. The following run shows the first problem. AICs (* and **) are different. I noticed that the Df for rs13482096:rs13483699 is 4,

The function heatmap() allows to specify row/column labels via the options labRow/labCol. From the code of heatmap(), I understand that when no labels are specified, the row/column labels (or indices) of the input matrix are taken as labels and re-ordered together with the rows and columns of the matrix before plotting, whereas labels supplied via labRow/labCol are plotted in the original order.

Full_Name: Anja von Heydebreck Version: 1.9.1 OS: Linux Submission from: (NULL) (155.250.128.25) I obtained a negative p-value from the fisher.test() function: > fisher.test(matrix(c(14576,3023,89,68),2), alternative="g")$p.value [1] -8.426593e-13 With R 1.8.1, I got a p-value of -6.239231e-12 for this example. Anja von Heydebreck

Full_Name: Anja von Heydebreck Version: 1.3.0 OS: Alpha Unix Submission from: (NULL) (141.14.19.61) Hi, I repeatedly obtained meaningless results from the function 'cutree' in the 'mva' package, when the argument 'h' was greater or equal to the maximum height occuring: > library('mva') > y [,1] [,2] [,3] [,4] [1,] 0 1 -1 1 [2,] 0 -1

Hi, ? I am very new to R. So, pardon my dumb question. I was trying to write my own function to run a different model (perform an ordered logistic regression) using the example in website http://pngu.mgh.harvard.edu/~purcell/plink/rfunc.shtml But R returns a error `R Error in eval(expr, envir, enclos) : object 's' not found' when I run it. What am I doing wrong here? Here's

Dear R-help list, Apologies. I am trying to convert one table to another. It feels that it should be a very straightforward answer with a single (or two) commands with the right extensions, but I really can't figure this out right now. I have several hundred pheno factors actually, so manually doing this line by line is not an option. My original table is approximately like this : ID

>>>>>>>>>>>>>>><<<<<<<<<<<<<<<< Pheno-Mint upsmon # /sbin/upsdrvctl start Network UPS Tools - UPS driver controller 2.6.1 Network UPS Tools - Generic HID driver 0.35 (2.6.1) USB communication driver 0.31 This PowerCOM device (0d9f:0004) is not (or perhaps not yet) supported by usbhid-ups. Please make sure

hi all, I have a matrix and named each row and column as like below... a<-matrix(c(seq(3,45,3),seq(10,6,-1)),4,5,byrow=F) > col<-c("peter","david","richrd","vincent","selva") > rows<-c("julius","caeser","anja","maya") > dimnames(a)<-list(rows,col) > a peter david

Hi, I am running the 64bit version of R on a AIX RS6000 cluster on which is also running the ' SQL*Plus: Release 9.2.0.2.0' in the 'ORACLE-environment Version 64 bit AIX 5.1'. I have problems setting up an Oralce database connection. I can instantiate a driver using drv <- Oracle(). But when I try to connect to my database using dbConnect (dbConnect(drv,

Hallo, I'd like to write a GUI (first choice with GTK+). I've surfed through the R- an Omegahat-Pages, because I'd like to use RGTK2, GTK 2.10.11 in combination with glade on Windows XP (perhaps later Unix, Mac). I've found a lot of different information. Because of the information I'm not sure, if this combination is running on Windows XP and I'm unsure how it works. Is

Dear R-Team! I am a beginner in R and use the WISP library for a project. I am using mark-recapture to estimate abundance. I would like to know if it is possible to only survey a part of the created area (with the created population) and not all of it. I am trying to quantify bias introduced by non random sampling design (i.e. the area is 24x100 but I only wanna search 24x40) Thank you already in

Hello does anyone know how to set the confidence level within the seqMK() function in the pheno package. It seems to be set automatically at 0.05 and there seems to be neither an input function to set a different level, as there are with some changepoint functions, nor an output to give one the p-value of the changepoint identification which you get with functions ie func$p_value etc ... Thanks

Dear everyone, Your code is slow and you are interested in performance optimization for scientific software? Figuring out where’s the bottlenecks guided by performance evaluation tools. If you are interested in porting your code to a HPC hardware platform and architecture, than OpenACC as a user-driven directive-based performance-portable parallel programming model might be a solution. What are