Re: [R] 1-pnorm values in a table

From: John Kane <jrkrideau_at_yahoo.ca>
Date: Thu, 06 Mar 2008 10:39:56 -0500 (EST)


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


> 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.
>



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. Received on Thu 06 Mar 2008 - 15:48:01 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 06 Mar 2008 - 16:30:20 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive