[R] converting to data.frame
dear allI have a data set with three types (Tree, Sapling, Seedling). I have estimated the correlation values. now i need to bring all the correlation values in a table like the one i have shown in attached file with R codes.could you please give me idea on this problem thanking you MSNepal _ Hotmail: Trusted email with Microsoft’s powerful SPAM protection. http://clk.atdmt.com/GBL/go/177141664/direct/01/read.table() colnames(c.df) [1] "id" "code""long""lat" "elev""age_obs" "year" [8] "Est_age" "Ht_cm" "BD_mm" "CDA_cm" "CDB_cm" "growth" "type" ### Correlation estimation cor(age_obs,Ht_cm) cor(age_obs[type=='SD'],Ht_cm[type=='SD']) #seedling cor(age_obs[type=='SPB'],Ht_cm[type=='SPB'])#sapling cor(age_obs[type=='T'],Ht_cm[type=='T'])#tree cor(age_obs,BD_mm) cor(age_obs[type=='SD'],BD_mm[type=='SD']) cor(age_obs[type=='SPB'],BD_mm[type=='SPB']) cor(age_obs[type=='T'],BD_mm[type=='T']) cor(age_obs,CDA_cm) cor(age_obs[type=='SD'],CDA_cm[type=='SD']) cor(age_obs[type=='SPB'],CDA_cm[type=='SPB']) cor(age_obs[type=='T'],CDA_cm[type=='T']) cor(age_obs,CDB_cm) cor(age_obs[type=='SD'],CDB_cm[type=='SD']) cor(age_obs[type=='SPB'],CDB_cm[type=='SPB']) cor(age_obs[type=='T'],CDB_cm[type=='T']) cor(Ht_cm,BD_mm) cor(Ht_cm[type=='SD'],BD_mm[type=='SD']) cor(Ht_cm[type=='SPB'],BD_mm[type=='SPB']) cor(Ht_cm[type=='T'],BD_mm[type=='T']) cor(Ht_cm,CDA_cm) cor(Ht_cm[type=='SD'],CDA_cm[type=='SD']) cor(Ht_cm[type=='SPB'],CDA_cm[type=='SPB']) cor(Ht_cm[type=='T'],CDA_cm[type=='T']) cor(Ht_cm,CDB_cm) cor(Ht_cm[type=='SD'],CDB_cm[type=='SD']) cor(Ht_cm[type=='SPB'],CDB_cm[type=='SPB']) cor(Ht_cm[type=='T'],CDB_cm[type=='T']) cor(BD_mm,CDA_cm) cor(BD_mm[type=='SD'],CDA_cm[type=='SD']) cor(BD_mm[type=='SPB'],CDA_cm[type=='SPB']) cor(BD_mm[type=='T'],CDA_cm[type=='T']) cor(BD_mm,CDB_cm) cor(BD_mm[type=='SD'],CDB_cm[type=='SD']) cor(BD_mm[type=='SPB'],CDB_cm[type=='SPB']) cor(BD_mm[type=='T'],CDB_cm[type=='T']) cor(CDA_cm,CDB_cm) cor(CDA_cm[type=='SD'],CDB_cm[type=='SD']) cor(CDA_cm[type=='SPB'],CDB_cm[type=='SPB']) cor(CDA_cm[type=='T'],CDB_cm[type=='T']) cor(age_obs,Est_age,use="pairwise.complete.obs") cor(age_obs[type=='SD'],Est_age[type=='SD'],use="pairwise.complete.obs") cor(age_obs[type=='SPB'],Est_age[type=='SPB'],use="pairwise.complete.obs") cor(age_obs[type=='T'],Est_age[type=='T'],use="pairwise.complete.obs") ### table looks like below one parameters combinedseedlingsapling tree age_obs Vs Ht_cm age_obs Vs BD_mm age_obs Vs CDA_cm age_obs Vs CDB_cm Ht_cm Vs BD_mm Ht_cm Vs CDA_cm Ht_cm Vs CDB_cm BD_mm Vs CDA_cm BD_mm Vs CDB_cm CDA_cm Vs CDB_cm age_obs Vs est_age __ R-help@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] converting to data.frame
> dear all I have a data set with three types (Tree, Sapling, Seedling). I have estimated the correlation values. now i need to bring all the correlation values in a table like the one i have shown in attached file with R codes.could you please give me idea on this problem > thanking you MS Nepal _ Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. http://clk.atdmt.com/GBL/go/171222985/direct/01/read.table() colnames(c.df) [1] "id" "code""long""lat" "elev""age_obs" "year" [8] "Est_age" "Ht_cm" "BD_mm" "CDA_cm" "CDB_cm" "growth" "type" ### Correlation estimation cor(age_obs,Ht_cm) cor(age_obs[type=='SD'],Ht_cm[type=='SD']) #seedling cor(age_obs[type=='SPB'],Ht_cm[type=='SPB'])#sapling cor(age_obs[type=='T'],Ht_cm[type=='T'])#tree cor(age_obs,BD_mm) cor(age_obs[type=='SD'],BD_mm[type=='SD']) cor(age_obs[type=='SPB'],BD_mm[type=='SPB']) cor(age_obs[type=='T'],BD_mm[type=='T']) cor(age_obs,CDA_cm) cor(age_obs[type=='SD'],CDA_cm[type=='SD']) cor(age_obs[type=='SPB'],CDA_cm[type=='SPB']) cor(age_obs[type=='T'],CDA_cm[type=='T']) cor(age_obs,CDB_cm) cor(age_obs[type=='SD'],CDB_cm[type=='SD']) cor(age_obs[type=='SPB'],CDB_cm[type=='SPB']) cor(age_obs[type=='T'],CDB_cm[type=='T']) cor(Ht_cm,BD_mm) cor(Ht_cm[type=='SD'],BD_mm[type=='SD']) cor(Ht_cm[type=='SPB'],BD_mm[type=='SPB']) cor(Ht_cm[type=='T'],BD_mm[type=='T']) cor(Ht_cm,CDA_cm) cor(Ht_cm[type=='SD'],CDA_cm[type=='SD']) cor(Ht_cm[type=='SPB'],CDA_cm[type=='SPB']) cor(Ht_cm[type=='T'],CDA_cm[type=='T']) cor(Ht_cm,CDB_cm) cor(Ht_cm[type=='SD'],CDB_cm[type=='SD']) cor(Ht_cm[type=='SPB'],CDB_cm[type=='SPB']) cor(Ht_cm[type=='T'],CDB_cm[type=='T']) cor(BD_mm,CDA_cm) cor(BD_mm[type=='SD'],CDA_cm[type=='SD']) cor(BD_mm[type=='SPB'],CDA_cm[type=='SPB']) cor(BD_mm[type=='T'],CDA_cm[type=='T']) cor(BD_mm,CDB_cm) cor(BD_mm[type=='SD'],CDB_cm[type=='SD']) cor(BD_mm[type=='SPB'],CDB_cm[type=='SPB']) cor(BD_mm[type=='T'],CDB_cm[type=='T']) cor(CDA_cm,CDB_cm) cor(CDA_cm[type=='SD'],CDB_cm[type=='SD']) cor(CDA_cm[type=='SPB'],CDB_cm[type=='SPB']) cor(CDA_cm[type=='T'],CDB_cm[type=='T']) cor(age_obs,Est_age,use="pairwise.complete.obs") cor(age_obs[type=='SD'],Est_age[type=='SD'],use="pairwise.complete.obs") cor(age_obs[type=='SPB'],Est_age[type=='SPB'],use="pairwise.complete.obs") cor(age_obs[type=='T'],Est_age[type=='T'],use="pairwise.complete.obs") ### table looks like below one parameters combinedseedlingsapling tree age_obs Vs Ht_cm age_obs Vs BD_mm age_obs Vs CDA_cm age_obs Vs CDB_cm Ht_cm Vs BD_mm Ht_cm Vs CDA_cm Ht_cm Vs CDB_cm BD_mm Vs CDA_cm BD_mm Vs CDB_cm CDA_cm Vs CDB_cm age_obs Vs est_age __ R-help@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] slope calculation
Dear all I am new R user, and trying to learn more. I am doing linear regression analysis in R with my data. I am trying to find the way to calculate the slope value (coefficient of x) to degree of slope. Please give me idea on this. Thanking you in anticipation Warm regardMS _ [[alternative HTML version deleted]] __ R-help@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] defining number of samples
dear all i am trying to perform t-test (t.test) in R here i am putting the commands and answer i got > t.test(Ht_cm[from_treeline=='above'][type=='SD'],Ht_cm[from_treeline=='below'][type=='SD']) Welch Two Sample t-test data: Ht_cm[from_treeline == "above"][type == "SD"] and Ht_cm[from_treeline == "below"][type == "SD"] t = 1.4181, df = 159.203, p-value = 0.1581alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -1.547018 9.426110 sample estimates:mean of x mean of y 21.34627 17.40672 > t.test(Ht_cm[from_treeline=='above'][type=='SPB'],Ht_cm[from_treeline=='below'][type=='SPB'])Error > in t.test.default(Ht_cm[from_treeline == "above"][type == "SPB"], : not > enough 'x' observations> its ok with the first command, but it says 'not enough 'x' observations', in the later case i have 31 samples. how can i do t.test with my day for latter case, will it work by defining the number of samples, if yes, how can i define it. thanking you in anticipation regard MS Nepal _ 24727::T:WLMTAGL:ON:WL:en-US:WWL_WIN_pcscout:102009 [[alternative HTML version deleted]] __ R-help@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] Sandard deviation calculation
Dear all I am trying to learn R I was trying to calculate standard deviation here are the commands and outputs in R > sd(Ht_cm[type=='SD'])[1] 3.283605> sd(Ht_cm[from_treeline=='above'])[1] > 16.83289> sd(Ht_cm[type=='SD'][from_treeline=='above'])[1] NA> the problem is that, i could not understand why the third command could not give any result. please help me out thanking you Regard MS Nepal _ [[elided Hotmail spam]] D24727::T:WLMTAGL:ON:WL:en-US:WWL_WIN_myidea:102009 [[alternative HTML version deleted]] __ R-help@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] x11(title=' ')
Dear all I was trying to put a title for my plot, but i got this result, > x11(width=10,height=5,title="seedling");par(mfrow=c(1,2))Error in x11(width = > 10, height = 5, title = "seedling") : unused argument(s) (title = > "seedling")> The title is not taking How can i give a title for the plot (where in need to make 2 plots within a window of x11() ) thanking you in anticipation Warm regard MS Nepal _ WLMTAGL:ON:WL:en-US:WWL_WIN_evergreen3:102009 [[alternative HTML version deleted]] __ R-help@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.
Re: [R] x11(title=' ')
Dear all sorry that i am late to reply for the question about my OS, I am using MS Windows XP 2003 OS thanking you Warm RegardMSNepal > Date: Wed, 28 Oct 2009 22:20:50 +1100 > From: j...@bitwrit.com.au > CC: r-help@r-project.org > Subject: Re: [R] x11(title=' ') > > On 10/28/2009 09:51 PM, Prof Brian Ripley wrote: > > ... > > Are you perchance x11() on Windows, when you should be using > > windows()? The posting guide asked you to tell us your OS, amongst > > other things > I say, my dear professor, if this "MS" chap who has used a Microsoft > email client and a Microsoft email service turns out to be using Linux, > he's covered his tracks pretty well. > > Jim > > __ > R-help@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. _ WLMTAGL:ON:WL:en-US:WWL_WIN_evergreen2:102009 [[alternative HTML version deleted]] __ R-help@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] z-test with NAs
Dear all I am learning R I am doing Z-test with package 'BSDA' here is my syntax and result: > sdgr.ztest<-z.test(growth[type=='SD'& from_treeline=='above'], + > growth[type=='SD'& from_treeline=='below'],alternative = "two.sided", + mu > = 0, sigma.x =(sd(growth[type=='SD'& from_treeline=='above'],na.rm = T)), + > sigma.y =(sd(growth[type=='SD'& from_treeline=='below'],na.rm = T)), > conf.level = 0.95)> sdgr.ztest # Z-test for Seedling Growth above and below > treeline Two-sample z-Test data: growth[type == "SD" & from_treeline == "above"] and growth[type == "SD" & from_treeline == "below"] z = NA, p-value = NAalternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: NA NA sample estimates:mean of x mean of yNANA I don't understand where my syntax went wrong. Similar syntax works for sapling case, but not for seedling case. is it due to the NAs in my data set for seedlings?, as the package says, NAs are allowed but will be removed during analysis, i don't think it is the problem. i am attaching the data file too. and one more information, I am using MS Windows 2003 OS. Please guide me, how can i go ahead. Thanking you in anticipation Warm regard MS Madan _ Bing brings you maps, menus, and reviews organized in one place. a=TEXT_MFESRP_Local_MapsMenu_Resturants_1x1id,code,growth,from_treeline,type 1,SD1,1.73,below,SD 2,SD2,1.12,above,SD 3,SD3,1.6,above,SD 4,SD4,1.66,above,SD 5,SD5,1.65,below,SD 6,SD6,1.43,above,SD 7,SD7,,above,SD 9,SD9,1.1,below,SD 10,SD10,1.44,above,SD 11,SD11,1.48,above,SD 12,SD12,2.25,above,SD 13,SD13,,above,SD 14,SD14,1.5,above,SD 15,SD15,,above,SD 16,SD16,1.06,above,SD 17,SD17,1.5,above,SD 18,SD18,4.17,below,SD 19,SD19,2.02,below,SD 20,SD20,1.38,above,SD 21,SD21,2,below,SD 22,SD22,1.41,above,SD 23,SD23,1.9,above,SD 24,SD24,1.7,above,SD 25,SD25,2.08,above,SD 26,SD26,1.92,below,SD 27,SD27,1.55,below,SD 28,SD28,1.48,below,SD 29,SD29,1.33,below,SD 30,SD30,1.5,above,SD 31,SD31,1.6,below,SD 32,SD32,1.24,below,SD 33,SD33,1.56,above,SD 34,SD34,1.58,below,SD 35,SD35,1.43,above,SD 36,SD36,1.3,above,SD 37,SD37,2.12,below,SD 38,SD38,1.63,below,SD 39,SD39,1.3,above,SD 40,SD40,1.16,below,SD 41,SD41,1.13,below,SD 42,SD42,1.43,below,SD 43,SD43,1.07,below,SD 44,SD44,1.25,below,SD 45,SD45,1.18,below,SD 46,SD46,1.17,above,SD 47,SD47,1,below,SD 48,SD48,1.43,above,SD 49,SD49,1.18,above,SD 50,SD50,1.17,below,SD 51,SD51,1.78,below,SD 52,SD52,1.67,above,SD 53,SD53,1,above,SD 54,SD54,1.8,below,SD 55,SD55,0.8,above,SD 56,SD56,0.95,above,SD 57,SD57,1.8,below,SD 58,SD58,,above,SD 59,SD59,,below,SD 60,SD60,,below,SD 61,SD61,,below,SD 62,SD62,,below,SD 63,SD63,,below,SD 64,SD64,,below,SD 65,SD65,0.93,below,SD 66,SD66,1.35,below,SD 67,SD67,1.33,below,SD 68,SD68,1.09,below,SD 69,SD69,1.41,below,SD 70,SD70,1.3,below,SD 71,SD71,1.2,below,SD 72,SD72,1.27,above,SD 73,SD73,0.7,below,SD 74,SD74,1.05,below,SD 75,SD75,1.68,below,SD 76,SD76,1.57,below,SD 77,SD77,1.22,above,SD 78,SD78,1.46,below,SD 79,SD79,2.32,below,SD 80,SD80,2.7,below,SD 81,SD81,1.57,above,SD 82,SD82,1.53,below,SD 83,SD83,1.63,below,SD 84,SD84,1.23,above,SD 85,SD85,1.77,below,SD 86,SD86,1.53,below,SD 87,SD87,1.39,below,SD 88,SD88,1.18,above,SD 89,SD89,1.85,below,SD 90,SD90,1.82,below,SD 91,SD91,,below,SD 92,SD92,1.74,below,SD 93,SD93,1.76,below,SD 94,SD94,1.16,below,SD 95,SD95,,below,SD 96,SD96,1.2,below,SD 97,SD97,1.42,below,SD 98,SD98,0.95,below,SD 99,SD99,2.43,below,SD 100,SD100,,below,SD 101,SD101,1.48,below,SD 102,SD102,1.33,below,SD 103,SD103,,below,SD 104,SD104,0.94,below,SD 105,SD105,,below,SD 106,SD106,1.15,below,SD 107,SD107,,below,SD 108,SD108,1.87,below,SD 109,SD109,1.23,below,SD 110,SD110,1.83,below,SD 111,SD111,1.47,below,SD 112,SD112,,below,SD 113,SD113,1.28,below,SD 114,SD114,,below,SD 115,SD115,1.57,below,SD 116,SD116,,below,SD 117,SD117,1.75,below,SD 118,SD118,1.8,below,SD 119,SD119,1.49,below,SD 120,SD120,1.54,below,SD 1,SPB1,2.19,above,SPB 2,SPB2,1.66,above,SPB 3,SPB3,1.78,above,SPB 4,SPB4,1.73,above,SPB 5,SPB5,8.63,above,SPB 6,SPB6,1.96,below,SPB 7,SPB7,5.17,above,SPB 8,SPB8,3.5,above,SPB 9,SPB9,2.78,above,SPB 10,SPB10,5.5,above,SPB 11,SPB11,2.75,below,SPB 12,SPB12,3,above,SPB 13,SPB13,4.71,above,SPB 14,SPB14,4.86,above,SPB 15,SPB15,4.63,below,SPB 16,SPB16,1.92,below,SPB 17,SPB17,3,below,SPB 18,SPB18,1.48,below,SPB 19,SPB19,2.35,below,SPB 20,SPB20,2.9,below,SPB 21,SPB21,2.35,below,SPB 22,SPB22,2.11,above,SPB 23,SPB23,1.92,below,SPB 24,SPB24,19.38,below,SPB 25,SPB25,2.33,above,SPB 27,SPB27,3.17,below,SPB 28,SPB28,1.73,below,SPB 29,SPB29,2.31,below,SPB 30,SPB30,1.9,below,SPB 31,SPB31,2.23,below,SPB 32,SPB32,2.04,below,SPB 33,SPB33,2.38,above,SPB 34,SPB34,4.9,above,SPB 35,SPB35,5.96,below,SPB 36,SPB36,1.77,below,SPB 37,SPB37,2.64,above,SPB 38,SPB38,4.01,above,SPB 39,SPB39,1.94,below,SPB 40,SPB40,1.8,above,SPB 41,SPB41,14.14,below,SPB 42,SPB
[R] optima in unimode
Dear all i could not estimate the optima value or range value in unimodal plot in glm please help me out thanking you regard madan _ Hotmail® has a new way to see what's up with your friends. orial_WhatsNew1_052009 [[alternative HTML version deleted]] __ R-help@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] montly mean temp plot
Dear all i got a problem in monthly mean temperature. here i am attaching the data set as well as the plot i got with the following command plot(month,type='n') plot(month,X1999) this command gave the plot where the month names are in alphabetic order, i want the plot in monthly sequence could you please suggest me how can i solve my problem? thanking you regard madan _ Lauren found her dream laptop. Find the PC that’s right for you. month 196719681969197019711972197319741975 1976197719781979198019811982198319841985 1986198719881989199019911992199319941995 1996199719981999200020012002200320042005 20062007 Janaury 13.312.912.711.914.513.618 DNA 11.5 11.912.911.413.712.512.314.511.5DNA 13.9 14.3DNA 14.812.716.213 13.513.315.412.7 14.212.714.415 14.614.213.914.713.613.4 18.214.6 Frebuary15 15.617.313.814.414.516.3DNA 14.315.414.914 14.614.116.513.214.2DNA 15.115.9DNA 16.715.814.317.613.615.915.2 14.316.113.616.420 14.517.617.215.315.9 17.519.113.7 March 16.819.119.118.818.620.219.7DNA 19.5 18.919.716.718.918.217.817.719.3DNA 21.1 18.9DNA 19.118.817 19.521.817.719.419.7 20.319 17.721.719.119.719.518.421.620 20 18.8 April DNA 22.421.921.618.521.523.7DNA 23.6 21.420.120.521.923.719.721.620.421.522.8 21.520.922.522.821.222.224 21.322.423.3 22.919.321.525 22.722.521.622.621.922.7 21.722.6 May 23.923.322.623.120.823.121.5DNA 23.6 21.220.421.924 22.521.823.821.522 22.7 21.624.423.623.722.524 22.423.124.125.6 24.623.323.724.223.823 22.623.623.623.4 23.324 June24.722.123 22.122.423.622.2DNA 23.4 22 22.823 23.422.423.423.324.423.525 23.922.923.724 25 23.924.924.324.623.8 23.724.325.824.124.124.223.924.123.924.7 24.223.6 JulyDNA 24 22.722.322.423 22.9DNA 22.1 22.622.922.422.922.822.823.523.922.422.6 23.122.723.923.223.924.523.824.325 23.6 23.224.523.924 24.524.623.824 23.223.9 24.623.5 August 24.722.622.422.821.822.822.6DNA 22.8 21.322.723.222.622.623.623.623.623.925.1 24 22.323.723.924.124 23.823.724.624.2 23.824.123.323.323.924.323.824.224.524.3 24.524.1 Spetemberr 23.823.121.822 21.921.521.5DNA 20.921.222.121.521.522.122.622.523.121.5 22.122.323 23.523 23.623.623.623.123.6 23.323.623.224.223.722.823.823 23.623 24.423.523 October 21.119.319.219.820.520.219.4DNA 20.4 19.318.819.719 19.920.919.7DNA 21.819.8 20.420.622.322 20.822.120.821.921.521.8 21.619.823.421.722.122 21.122.120.621.6 22.221.9 November15.917.116.516 16.717.816.3DNA 15.516.416.615.717.717.317.316.417.316.9 16.716.918.318.117.119.317 17.518 17.4 18.318.618.120 18.618.418.918.418.417 18.218.1DNA December15.714.513.915.614 17.213.1DNA 13.415 12.615.313.115.314 14.215 DNA 15.8DNA 16.215.414.815.814.714.316.115.1 14.715.814.416.515.415.315 15.314.915.3 16.415.716 __ R-help@r-project.org mailing list ht
[R] threat analysis
dear all i am in great need of threat analysis technique, could anyone give me information about this method in R, which package does the Threat Analysis (in animals) thanking you Ms.com _ Hotmail® has a new way to see what's up with your friends. orial_WhatsNew1_052009 [[alternative HTML version deleted]] __ R-help@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] upgrade R version
dear all contributors i am new R user i want to upgrade my R version from 2.7.1 to 2.8.1 i tried myself, but i could not do it as R is different from other software could anyone guide me please warm regard madan student University of Bergen Windows Live: Keep your life in sync. Check it out. Windows Live Hotmail®: Chat. Store. Share. Do more with mail. See how it works. _ [[alternative HTML version deleted]] __ R-help@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.