search for: yulei

Hi, there. Can I ask how to set up missing values in R? Suppose I want to assign the missing value to the elements in vector which is greater than zero like this: x<-c(1,3,-1,0,4); after the missing value assignment, x becomes (NA,NA,-1,0,NA). Thanks! Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi, there. I have two questions about using R to create boxplots. 1. The function boxplot() plots the outliers. How can I label the exact values arount these outlier points? Does R have an option allow me to do that? 2. How can I put two boxplots in one x-y axis? Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 276 Grove St. Apt 3 Newton, MA 02466 617-796-7834(H) 617-432-3428(O) 617-432-3435(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

...tool to check if a multivariate sample is from multivariate normal distribution? If there is one, is there any function doing it in R? 2. Is there any function of testing if two multivariate distribution are same, i.e. the multivariate extension of Kolomogrov-Smirnov test? Thanks for your help. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

...X_2_3, .. X_2_10,..., X_10_1, X_10_2, ... X_10_10. It looks like I need to use 2 indexes I and J so I is looped from 1 to 10, and J is looped from 1 to 10. But I don't know how to automatically produce X with these combination of indexes. Should I use paste function? Can someone help? Thanks. Yulei [[alternative HTML version deleted]]

Hi, there. Does anybody know how to plot a smooth density plot for some data simulated from certain distribution? Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi, there. Suppose I have a bivarariate data matrix y1 and y2. I want to plot a 3-D picture of the estimated density f(y1, y2) against y1 and y2? How can I do that? Do I use persp() or density()? Thanks for your help. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi, there. Suppose I already have sensitivities and specificities. What is the quick R-function to calculate AUC for the ROC plot? There seem to be many R functions to calculate AUC. Thanks. Yulei [[alternative HTML version deleted]]

Hi, there. Is there any function in R to plot the probability-probability plot (PP plot)? Suppose I am testing some data against normal. Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi, there. Does anyone know the R function for calculating the cdf of bivariate normal distribution function? Thanks. Yulei [[alternative HTML version deleted]]

Hi, there: Could someone tell me a simple function of plot ROC curve and calculate AUC in R? My setting is very simple, a column of the true binary response and another column of predicted probabilities. Thanks! Yulei [[alternative HTML version deleted]]

Hi, there. Suppose I have a bivariate data set y1 and y2. Can anybody tell me how to estimate the conditional density of f(y1|y2) and vice versa? Thanks. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

Hi, there. How do I calculate the cross-product in the form of \sum_{i=1}^{n}X_{i}^{t} \Sigma X_{i} using R code without using do loop? X_{i} is the covariate matrix for subject I, \Sigma is the covariance matrix. Thanks for your help. Yulei [[alternative HTML version deleted]]

Hi, there: Does anyone know how to specify the ranges in the axises when I make scatter plots using pairs()? In the general plot function, I can use xlim and ylim option. But how can I do this if I use pairs()? Thanks. Yulei [[alternative HTML version deleted]]

...5990 -0.19470 3 FL -0.47900 -0.62930 -0.32130 I would like to have a plot the 95% CI (characterized by the mean, lower, and upper bound, and the lines linking them) for each state, with x-axis is the parameter estimate, and y-axis is the state. What kind of plot functions can I use? Thanks. Yulei [[alternative HTML version deleted]]

...nt(s) (summary.only = FALSE) The program is running: > schools.sim <- bugs(data, inits, parameters, model.file, + n.chains=3, n.iter=5000, + bugs.directory="c:/Program Files/WinBUGS14/", + working.directory=NULL, clearWD=TRUE) Thanks for your help! Yulei [[alternative HTML version deleted]]

...in r[i1] - r[-length(r):-(length(r) - lag + 1)] : non-numeric argument to binary operator I also try to do simple code with it, and getting error message as well, such as >for(j in 1:length(x)){w<-x/1} >Error in x/1 : non-numeric argument to binary operator Thanks for your help! Yulei [[alternative HTML version deleted]]

...ve) alias analysis algorithm implemented in LLVM. (I use LLVM 3.9, hope to switch to 4.0 soon.) As far as I know, none of the built-in analysis (basicAA, globals-modref, andersAA, etc.) is intended to be flow-sensitive. So I searched and came across these two 1. https://github.com/unsw-corg/SVF by Yulei Sui (for LLVM 3.8) 2. http://www.cs.ucsb.edu/~benh/research/downloads.html by Ben Hardekopf (for LLVM 2.5) Are there other implementations which use the AA-Interface? Giving a little context, I need some functionality in LLVM that answers the following question. For a given argument of a call ins...

.... Since the Gibbs chains are independent, so does the draw of the random effects for the data matrix, I am just wondering if there exists faster command in R to do above operation. I happen to see function sapply(), will it be faster than my double for loops? Your help will be greatly appreciated. Yulei $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ Yulei He 1586 Murfin Ave. Apt 37 Ann Arbor, MI 48105-3135 yuleih at umich.edu 734-647-0305(H) 734-763-0421(O) 734-763-0427(O) 734-764-8263(fax) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

...The reason why I chose not to use Multiple Imputation(MI) is because according to [http://blogs.iq.harvard.edu/sss/archives/2008/09/a_handy_trick_f.shtml] "MI packages assume the Multivariate Normal Distribution which may not hold for certain types of categorical and binary data. Yucel Recai, Yulei He, and Alan Zaslavsky point out in their May 2008 article in The American Statistician, naive rounding MI imputations can bias estimates, particularly when the underlying data are asymmetric or multimodal." However instead of using Yucel Recai, Yulei He, and Alan Zaslavsky's rounding st...

...didn't set the GATEWAY in ifcfg-eno1.5, so I added the GATEWAY=192.168.81.254(My local area network gateway), it didn't work. Could anybody teach me how to create a VLAN interface in CentOS7, or tell me what' s wrong with my configuration, more detailed more better. Thanks Li yulei