An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/r-help/attachments/20080306/e003d5bb/attachment.pl
1-pnorm values in a table
2 messages · Jabez Wilson, John Kane
This is a bit ugly but I think it works.
================================================
myf <- function(x) 1-pnorm(x,mean(x), sd(x))
results <- apply(test, 2, myf)
mymeans <- apply(test, 2, mean); mymeans
for (i in 1:length(test)){
test[,i][test[,1]>=mymeans[i]] <- NA
}
results[is.na(test)] <- NA
====================================================
--- Jabez Wilson <jabezwuk at yahoo.co.uk> wrote:
Hi,
I've read in a csv file (test.csv) which gives me
the following table:
Hin1 Hin2 Hin3 Hin4 Hin5
Hin6
HAI1 9534.83 4001.74 157.16 3736.93 484.60 59.25
HAI2 13272.48 1519.88 36.35 33.64 46.68 82.11
HAI3 12587.71 5686.94 656.62 572.29 351.60 136.91
HAI4 15240.81 10031.57 426.73 275.29 561.30 302.38
HAI5 15878.32 10517.14 18.93 22.00 16.91 21.17
I would like to find a way of finding the 1-pnorm
of each value in the table based on the mean and sd
of the data only in the column in which the value
lies. I can do it using a for loop, but would like
to know if it can be done using e.g. apply or
something similar, so that the whole table is
printed out with the 1-pnorm values.
1-pnorm(test[,1],mean([,1]), sd([,1])) gives me
the values for col1 only, but that's as far as I've
got.
tia
p.s. I know I'm asking a lot, but ideally, I'd
like to print out the table with those 1-pnorm
values only if they are in the right hand tail (i.e.
= mean) and if not nothing or NA be written.
---------------------------------
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.