similar to: Plot problem

Hi all: A short question: For example, a=[3,4,6,2,3], obviously the 3rd entry of the array has the maxium value, what I want is index of the maxium value: 3. is there a neat expression to get this index? Thank you! Best, Feng [[alternative HTML version deleted]]

Hello list, I am trying to apply the paired t.test between diseased and not diseased patients to identify genes that are more expressed in the one situation under the other. In order to retrieve the genes that are more expressed in the positive disease state I do: p.values<-c() for(i in 1:length(Significant[,1])){ p.values[i]<-try(t.test(positive[i,],negative[i,],alternative

Suppose there are 4 variables d is a function of a , b and c I want to know how a, b and c change will make d change It will be straightforward to see it if we can graph the d surface if d is only a function of a and b, I can use 'persp' to see the surface of d. I can easily see at what values of a and b, d will get the maxium or minium or multiple modes, etc But for 4 dimention graph,

## Question1: How to define IV with interaction alone, without main effects? ## Question2: Should Type III ANOVA in package car be independent of the factor level order? ## data from http://www.otago.ac.nz/sas/stat/chap30/sect52.htm drug <- c(t(t(rep(1,3)))%*%t(1:4)); disease <- c(t(t(1:3)) %*% t(rep(1,4))); y <- t(matrix(c( 42 ,44 ,36 ,13 ,19 ,22 ,33 ,NA ,26 ,NA ,33 ,21 ,31 ,-3 ,NA

I would be grateful for advice about the following problem. It's not directly R-related, but I'm hoping that R will help me analyse the following data. I have a table which indicates the progression of a certain age-related disease. At a certain point in time, a population was sampled; and I have measurements for the age of each individual, and their disease stage. (Disease stage is an

I'm having trouble figuring out how to format Date variables when used as axis labels in graphs. The particular case here is an attempt to re-create Nightingale's coxcomb graph with ggplot2, where I'd like the months to be labeled as "Mar 1885", "Apr 1885", using a date format of "%b %Y" applied to label the dates, or really anything other than

Dear R experts, Suppose I have an data frame likes this: > example <- data.frame(age=c(1,2,3, 4,5,6), height=c(100,110,120,130,140,150), disease=c(TRUE, TRUE, TRUE, FALSE, FALSE, FALSE)) > example age height disease 1 1 100 TRUE 2 2 110 TRUE 3 3 120 TRUE 4 4 130 FALSE 5 5 140 FALSE 6 6 150 FALSE Is there anyway to compare the age and

Hi all I am wandering how to control /proc/sys/fs/file-max ls it automatically? and what is for? and the maxium? Thank you --------------------------------- Ask a question on any topic and get answers from real people. Go to Yahoo! Answers. -------------- next part -------------- An HTML attachment was scrubbed... URL:

> > There is also Citrix XenServer for maxium > comfort. It offers GUI tools under Windows for managing VMs. The > drawback is, there are a lot of parameters you cannot set in the non > paid version. Just to throw it out there, there's also XCP-NG which is a rebranded (w/o Citrix) version of XenServer without all of the proprietary bits. The same people also produce

Hi, This code is trying to produce a series of graphics files, with plots of male and female disease rates by age, one plot per disease. The dataframe contains a variable 'Age' and a set of variables called 'Male_CVD, Female_CVD,Male_RTA,Female_RTA, and so on. For each disease, I want to pull out the column of data containing the word 'Male' and plot this against age, and then

Hi I got an error message using datadist() from Design package: > library(Design,T) > dd <- datadist(beta.final) > options(datadist="dd") > lrm(Disease ~ gsct+apcct+rarct, x=TRUE, y=TRUE) Error in eval(expr, envir, enclos) : object "Disease" not found All variables inclduing response variable "Disease" are in the data frame

Hi Chin Yi, If you are trying to correlate "Health" with "Disease", i.e. cydf<-read.table(text="OTU ID Health Disease Bacterial 1 0.29 0.34 Bacterial 2 0.25 0.07 Bacterial 3 0.06 0.06 Bacterial 4 0.07 0.09 Bacterial 5 0.02 0.05", header=TRUE) print(cor(cydf$Health,cydf$Disease)) [1] 0.7103517 If you are getting that error, it probably means that

I have been trying to open data that I have saved in an excel spread sheet. I saved it as a csv. Then I tried using the read.csv command. However, everytime I do this-- diseasedat<-read.csv("M:/sloan/R/disease/disease.csv", sep=, header = TRUE, fill= TRUE)-- I get an error message: Error in file(file, "r") : unable to open connection In addition: Warning message:

hello, I asked about this simulation a few days ago, but still i can't get what i need. I have 2 units: icu and regular. from icu I want to take 200 observations from binomial distribution, when probability for disease is: p=0.6. from regular I want to take 300 observation with the same probability: p=0.6 . the distribution to detect disease when disease occurred- *for someone from icu* - is:

Hi R help list I'm simulating logistic regression data with a specified odds ratio (beta) and have a problem/unexpected behaviour that occurs. The datasets includes a lognormal exposure and diseased and healthy subjects. Here is my loop: ors <- vector() for(i in 1:200){ # First, I create a vector with a lognormally distributed exposure: n <- 10000 # number of study subjects

On 2017-07-05 11:56, Jim Lemon wrote: > Hi Chin Yi, > If you are trying to correlate "Health" with "Disease", i.e. > > cydf<-read.table(text="OTU ID Health Disease > Bacterial 1 0.29 0.34 > Bacterial 2 0.25 0.07 > Bacterial 3 0.06 0.06 > Bacterial 4 0.07 0.09 > Bacterial 5 0.02 0.05", > header=TRUE) >

Dear UseRs, Suppose I have data regarding smoking habits of a prospective cohort and wish to determine the risk ratio of colorectal cancer in the smokers compared to the non-smokers. What do I do at the end of the study with people who die of heart disease? Can I just censor them exactly the same as people who become uncontactable or who die in a plane crash? If not, why not? I'm thinking

Dear all, I have a stat question that may not be related to R, but I would like to have your advice. I have just read a medical paper in which the authors report the 1-p (where p is the cumulative survival probability from the Kaplan Meier curve) as incidence of disease. Specifically, the study followed ~12000 women on drug A and ~20000 women on drug B for 12 months. During that period

hello R users, I didn't find a solution for a special problem. I have two dataframes. dataframe1: X value row col ID 1 8.973498062 5512625 3460000 1 2 11.656658570 5501625 3464000 2 3 11.121777570 5495625 3473000 3 4 9.310465964 5508625 3477000 4 5 8.883483845 5515625 3496000 5 dataframe2: X value

Hello R gurus, I have a data set from which i have to extract the gender and age matched rows from controls and disease group disease<-paste(rep(c('y','n'),11)) gender<-paste(rep(c('m','f'),11)) mcp<-rnorm(700,1400) age<-rnorm(32,34) dat<-data.frame(disease=disease,sex=gender,Dr_age=age[1:22],MCP=mcp[1:22]) I have other categorical variables also to