Re: [R] glm.fit: fitted probabilities numerically 0 or 1 occurred?

Dear ted, thanks for your help. Now, everything is more clear. I read something about linear separation you mentioned, and my data set is very suitable for this problem. But, there is a confusing question for me; I can not controll the process adequately because of usage of bootstrap. So, this

Re: [R] glm.fit: fitted probabilities numerically 0 or 1 occurred?

Dear Ellison, Many thanks for your reply. The information you typed is clear and now I know what to do. Your suggestion about finding some coffee while running simulation is so good =) Regards - Best regards Ufuk -- View this message in context:

[R] glm.fit: fitted probabilities numerically 0 or 1 occurred?

Hi all, I am doing bootstrap with logistic regression by using glm function, and I get the errors; glm.fit: fitted probabilities numerically 0 or 1 occurred and glm.fit: algorithm did not converge I have read some things about this issue in the mailing list. I can guess what was the problem. My

[R] Robust Statistics for Outlier Detection

Hi Dear All, Can someone give me a suggestion about which robust statistics are most appropriate for outlier detection in linear models, and is available with R ? Thanks for any idea. -- View this message in context:

[R] Why unique(sample) decreases the performance ?

Hi, I' am interested in differences between sample's result when samples consist of full elements and consist of only distinct elements. When sample consist of full elements it take about 120 sec., but when consist of only distinct elements it take about 4.5 or 5 times more sec. I expected that

Re: [R] Why unique(sample) decreases the performance ?

Sorry, the last part of code does not work when uniqu() are used, the true version; e - rnorm(n=50, mean=0, sd=sqrt(0.5625)) x0 - c(rep(1,50)) x1 - rnorm(n=50,mean=2,sd=1) x2 - rnorm(n=50,mean=2,sd=1) x3 - rnorm(n=50,mean=2,sd=1) x4 - rnorm(n=50,mean=2,sd=1) y - 1+ 2*x1+4*x2+3*x3+2*x4+e x2[1] =

Re: [R] A question about list

Many thanks to all -- View this message in context: http://r.789695.n4.nabble.com/A-question-about-list-tp3385711p3387410.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] A question about list

Hi dear all, It may be a simple question, i have a list output with different number of elements as following; [[1]] [1] 0.86801402 -0.82974691 0.3974 -0.98566707 -4.96576856 -1.32056754 [7] -5.54093319 -0.07600462 -1.34457280 -1.04080125 1.62843297 -0.20473912 [13] 0.30659907

[R] selection of a subset from a loop

Hi dear all, The code like this; e - rnorm(n=50, mean=0, sd=sqrt(0.5625)) x0 - c(rep(1,50)) x1 - rnorm(n=50,mean=2,sd=1) x2 - rnorm(n=50,mean=2,sd=1) x3 - rnorm(n=50,mean=2,sd=1) x4 - rnorm(n=50,mean=2,sd=1) y - 1+ 2*x1+4*x2+3*x3+2*x4+e x2[1] = 10 #influential observarion y[1] = 10

[R] Problem about for loop

Hi everyones, my function like; e - rnorm(n=50, mean=0, sd=sqrt(0.5625)) x0 - c(rep(1,50)) x1 - rnorm(n=50,mean=2,sd=1) x2 - rnorm(n=50,mean=2,sd=1) x3 - rnorm(n=50,mean=2,sd=1) x4 - rnorm(n=50,mean=2,sd=1) y - 1+ 2*x1+4*x2+3*x3+2*x4+e x2[1] = 10 #influential observarion y[1] = 10

[R] selection statistics from function

Hi, My code: e - rnorm(n=50, mean=0, sd=sqrt(0.5625)) x0 - c(rep(1,50)) x1 - rnorm(n=50,mean=2,sd=1) x2 - rnorm(n=50,mean=2,sd=1) x3 - rnorm(n=50,mean=2,sd=1) x4 - rnorm(n=50,mean=2,sd=1) y - 1+ 2*x1+4*x2+3*x3+2*x4+e x2[1] = 10 #influential observarion y[1] = 10 #influential observarion

[R] Computing and Finding

Hi dear all, i am triying to do jackknife-after bootstrap for detection of influential observation. my data and resamples are following ; e - rnorm(n=50, mean=0, sd=sqrt(0.5625)) x0 - c(rep(1,50)) x1 - rnorm(n=50,mean=2,sd=1) x2 - rnorm(n=50,mean=2,sd=1) x3 - rnorm(n=50,mean=2,sd=1) x4 -

Re: [R] selection of outputs from the function

Thanks so much Joshua Wiley you are right your function better than mine :) -- View this message in context: http://r.789695.n4.nabble.com/selection-of-outputs-from-the-function-tp3163361p3163853.html Sent from the R help mailing list archive at Nabble.com.

[R] selection of outputs from the function

Hi Dear All, This is a function which contains Covariance Ratio and Likelihood Distance values (CVRi, LDi). i want to compute the all row's values, that is run this function for nrow(X) times. The X and Y matrices are;

Re: [R] selection of outputs from the function

Thank so much David ! -- View this message in context: http://r.789695.n4.nabble.com/selection-of-outputs-from-the-function-tp3163361p3163421.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] regression

Hi dear all, suppose that s is a statistic code; i have a matrix (x) which has 7 columns (1=x1,2=x23=x3,4=x4,5=x5,6=x6 and7=y) and has 20 rows. i want to do linear reggression like reg-lm(x[,7]~1+x[,1]+x[,2]+...+x[,6]) but i want to do delete i th row for nrows times and create regression

Re: [R] regression

Thank you so much i did with your idea.. thank you :) -- View this message in context: http://r.789695.n4.nabble.com/regression-tp3161328p3161524.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] regression

sorry not read i want to say write -- View this message in context: http://r.789695.n4.nabble.com/regression-tp3161328p3161551.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] regression

thank you so much, and i have a question too if i read some statistic for example cook-weisberg statistic or welsh-kuh distance and i say stat-welsh.kuh than i put this statistic in your idea, can i get the statistics each times? -- View this message in context:

Re: [R] delete-d jackknife

Thank you so much :) -- View this message in context: http://r.789695.n4.nabble.com/delete-d-jackknife-tp3058335p3059364.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] delete-d jackknife

Hi dear all, Can aynone help me about delete-d jackknife usually normal jackknife code for my data is: n - nrow(data) y - data$y z - data$z theta.hat - mean(y) / mean(z) print (theta.hat) theta.jack - numeric(n) for (i in 1:n) theta.jack[i] - mean(y[-i]) / mean(z[-i]) bias - (n - 1) *

[R] Counting

Hi dear all, i have a data (data.frame) which contain y and x coloumn(i.e. y x 1 0.58545723 0.15113102 2 0.02769361 -0.02172165 3 1.00927527 -1.80072610 4 0.56504053 -1.12236685 5 0.58332337 -1.24263981 6 -1.70257274 0.46238255 7 -0.88501561 0.89484429 8

Re: [R] Counting

thank you very much for your idea, if i write code as; my data name is data. samples-function(data,num){ resamples-lapply(1:num,function(i) sample(data,n,replace=TRUE)) list(resamples=resamples)} n=10 data-rnorm(n=10,mean=5,sd=2) data[1]=100 obj-samples(data,1000) i generate 1000 sample, i did

Re: [R] Counting

Thank you so much -- View this message in context: http://r.789695.n4.nabble.com/Counting-tp3045756p3045918.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] Counting

Thank you so much -- View this message in context: http://r.789695.n4.nabble.com/Counting-tp3045756p3045917.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

Re: [R] jackknife-after-bootstrap

I am always trying but i could not do it. Are there any example about this -- View this message in context: http://r.789695.n4.nabble.com/Re-jackknife-after-bootstrap-tp3043213p3043398.html Sent from the R help mailing list archive at Nabble.com. __