##. In this script I will remove Non T cells from PBMC

1. load libraries

2. Load Seurat Object

Cell type Distribution to check clusters

Cell type Distribution barplot

3. filter cells just keep CD4 T cells

Cell type Distribution to check clusters

4. filter B cells from L4

Cell type Distribution to check clusters

5. Save the Seurat object as an Robj file

LS0tCnRpdGxlOiAidXNlIEFubm90YXRlZCBSb2JqIGluY2x1ZGluZyBQQk1DMTB4IHRvIHJlbW92ZSBOb25DRDRUY2VsbHMgZnJvbSBDb250cm9sIGFuZCBCIGNlbGxzIGZyb20gTDQgYW5kIGtlZXAganVzdCBDRDRUY2VsbHMiCmF1dGhvcjogTmFzaXIgTWFobW9vZCBBYmJhc2kKZGF0ZTogImByIFN5cy5EYXRlKClgIgpvdXRwdXQ6CiAgIyBwZGZfZG9jdW1lbnQ6IGRlZmF1bHQKICAjIHdvcmRfZG9jdW1lbnQ6IGRlZmF1bHQKICAjIGh0bWxfZG9jdW1lbnQ6IGRlZmF1bHQKICAjcm1kZm9ybWF0czo6cmVhZHRoZWRvd24KICBodG1sX25vdGVib29rOgogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICAgIHRvY19jb2xsYXBzZWQ6IHRydWUKLS0tCgoKCiMjLiBJbiB0aGlzIHNjcmlwdCBJIHdpbGwgcmVtb3ZlIE5vbiBUIGNlbGxzIGZyb20gUEJNQyAKCiMgMS4gbG9hZCBsaWJyYXJpZXMKYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9CgpsaWJyYXJ5KFNldXJhdCkKbGlicmFyeShTZXVyYXRPYmplY3QpCmxpYnJhcnkoU2V1cmF0RGF0YSkKbGlicmFyeShwYXRjaHdvcmspCgpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KFJDb2xvckJyZXdlcikKbGlicmFyeShtYWdyaXR0cikKbGlicmFyeShkYnBseXIpCmxpYnJhcnkocm1hcmtkb3duKQpsaWJyYXJ5KGtuaXRyKQpsaWJyYXJ5KHRpbnl0ZXgpCiNBemltdXRoIEFubm90YXRpb24gbGlicmFyaWVzCmxpYnJhcnkoQXppbXV0aCkKI1Byb2plY1RpbHMgQW5ub3RhdGlvbiBsaWJyYXJpZXMKbGlicmFyeShTVEFDQVMpCmxpYnJhcnkoUHJvamVjVElMcykKI3NpbmdsZVIgQW5ub3RhdGlvbiBsaWJyYXJpZXMKbGlicmFyeShTaW5nbGVSKQpsaWJyYXJ5KGNlbGxkZXgpCmxpYnJhcnkoU2luZ2xlQ2VsbEV4cGVyaW1lbnQpCgpsaWJyYXJ5KGNsdXN0cmVlKQoKYGBgCgoKIyAyLiBMb2FkIFNldXJhdCBPYmplY3QgCmBgYHtyIGxvYWRfc2V1cmF0fQoKI0xvYWQgU2V1cmF0IE9iamVjdCBtZXJnZWQgZnJvbSBjZWxsIGxpbmVzIGFuZCBhIGNvbnRyb2woUEJNQykgYWZ0ZXIgZmlsdHJhdGlvbgpsb2FkKCIuLi8wLVJfT2JqZWN0cy9BbGxfU2FtcGxlc19NZXJnZWRfd2l0aF8xMHhfQXppdG11dGhfQW5ub3RhdGVkLnJvYmoiKQoKCkFsbF9zYW1wbGVzX01lcmdlZAoKCmBgYAoKIyMgIENlbGwgdHlwZSBEaXN0cmlidXRpb24gdG8gY2hlY2sgY2x1c3RlcnMKYGBge3IgRGlzdHJpYnV0aW9uMSwgZmlnLmhlaWdodD0xMCwgZmlnLndpZHRoPTEwfQojIFdlIGNhbiBhcHBseSBpdCBsYXRlciBvbiBSIG9iaiB0byBnZXQgdGhlc2UgdGFibGVzIHRvIGNvbXBhcmUgaXQgdG8gZGVjaWRlIGFib3V0IHJlc29sdXRpb24uCgojIEF6aW11dGggbDEKamFuaXRvcjo6dGFieWwoQWxsX3NhbXBsZXNfTWVyZ2VkQG1ldGEuZGF0YSwgcHJlZGljdGVkLmNlbGx0eXBlLmwxLCBjZWxsX2xpbmUpCgoKIyBBemltdXRoIGwyCmphbml0b3I6OnRhYnlsKEFsbF9zYW1wbGVzX01lcmdlZEBtZXRhLmRhdGEsIHByZWRpY3RlZC5jZWxsdHlwZS5sMiwgY2VsbF9saW5lKQoKIyBBemltdXRoIGwzCmphbml0b3I6OnRhYnlsKEFsbF9zYW1wbGVzX01lcmdlZEBtZXRhLmRhdGEsIHByZWRpY3RlZC5jZWxsdHlwZS5sMywgY2VsbF9saW5lKQoKYGBgCgojIyAgQ2VsbCB0eXBlIERpc3RyaWJ1dGlvbiBiYXJwbG90CmBgYHtyIERpc3RyaWJ1dGlvbjIsIGZpZy5oZWlnaHQ9NiwgZmlnLndpZHRoPTEwfQoKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KFJDb2xvckJyZXdlcikgIAoKIyBBc3N1bWluZyB5b3UgaGF2ZSAxMCBkaWZmZXJlbnQgY2VsbCBsaW5lcywgZ2VuZXJhdGluZyBhIGNvbG9yIHBhbGV0dGUgd2l0aCAxMCBjb2xvcnMKY2VsbF9saW5lX2NvbG9ycyA8LSBicmV3ZXIucGFsKDEwLCAiU2V0MyIpCgojIEFzc3VtaW5nIEFsbF9zYW1wbGVzX01lcmdlZCRjZWxsX2xpbmUgaXMgYSBmYWN0b3Igb3IgY2hhcmFjdGVyIHZlY3RvciBjb250YWluaW5nIGNlbGwgbGluZSBuYW1lcwpkYXRhIDwtIGFzLmRhdGEuZnJhbWUodGFibGUoQWxsX3NhbXBsZXNfTWVyZ2VkJGNlbGxfbGluZSkpCmNvbG5hbWVzKGRhdGEpIDwtIGMoImNlbGxfbGluZSIsICJuVU1JIikgICMgQ2hhbmdlIGNvbHVtbiBuYW1lIHRvIG5VTUkKCm5jZWxscyA8LSBnZ3Bsb3QoZGF0YSwgYWVzKHggPSBjZWxsX2xpbmUsIHkgPSBuVU1JLCBmaWxsID0gY2VsbF9saW5lKSkgKyAKICBnZW9tX2NvbCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGdlb21fdGV4dChhZXMobGFiZWwgPSBuVU1JKSwgCiAgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjkpLCAKICAgICAgICAgICAgdmp1c3QgPSAtMC4yNSkgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IGNlbGxfbGluZV9jb2xvcnMpICsgCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgaGp1c3QgPSAxKSwKICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KGhqdXN0ID0gMC41KSkgKyAgIyBBZGp1c3QgdGhlIHRpdGxlIHBvc2l0aW9uCiAgZ2d0aXRsZSgiRmlsdGVyZWQgY2VsbHMgcGVyIHNhbXBsZSIpICsKICB4bGFiKCJDZWxsIGxpbmVzIikgKyAgIyBBZGp1c3QgeC1heGlzIGxhYmVsCiAgeWxhYigiRnJlcXVlbmN5IikgICAgIyBBZGp1c3QgeS1heGlzIGxhYmVsCgpwcmludChuY2VsbHMpCgoKYGBgCgojIDMuIGZpbHRlciBjZWxscyBqdXN0IGtlZXAgQ0Q0IFQgY2VsbHMKYGBge3IgRGlzdHJpYnV0aW9uMywgZmlnLmhlaWdodD02LCBmaWcud2lkdGg9MTB9CgojIFNldCBpZGVudGl0eSB0byBjZWxsX2xpbmUgCklkZW50cyhBbGxfc2FtcGxlc19NZXJnZWQpIDwtICJjZWxsX2xpbmUiCgojIElkZW50aWZ5IENENCBUIGNlbGxzIGluIFBCTUMgYW5kIFBCTUNfMTB4CmNkNF9jZWxsc19wYm1jIDwtIFdoaWNoQ2VsbHMoQWxsX3NhbXBsZXNfTWVyZ2VkLCBleHByZXNzaW9uID0gCiAgZ3JlcGwoIl5DRDQiLCBBbGxfc2FtcGxlc19NZXJnZWQkcHJlZGljdGVkLmNlbGx0eXBlLmwxKSAmCiAgZ3JlcGwoIl5DRDQiLCBBbGxfc2FtcGxlc19NZXJnZWQkcHJlZGljdGVkLmNlbGx0eXBlLmwyKSAmCiAgZ3JlcGwoIl5DRDQiLCBBbGxfc2FtcGxlc19NZXJnZWQkcHJlZGljdGVkLmNlbGx0eXBlLmwzKSAmCiAgKGNlbGxfbGluZSAlaW4lIGMoIlBCTUMiLCAiUEJNQ18xMHgiKSkKKQoKIyBTZWxlY3QgYWxsIGNlbGxzIGZyb20gb3RoZXIgY2VsbCBsaW5lcwpvdGhlcl9jZWxsc19zZXVyYXQgPC0gc3Vic2V0KEFsbF9zYW1wbGVzX01lcmdlZCwgc3Vic2V0ID0gY2VsbF9saW5lICVpbiUgYygiTDEiLCAiTDIiLCAiTDMiLCAiTDQiLCAiTDUiLCAiTDYiLCAiTDciKSkKCiMgR2V0IHRoZSBjZWxsIG5hbWVzIGZyb20gdGhlIHN1YnNldApvdGhlcl9jZWxscyA8LSBjb2xuYW1lcyhvdGhlcl9jZWxsc19zZXVyYXQpCgojIENvbWJpbmUgdGhlIGNlbGwgbGlzdHMKY2VsbHNfdG9fa2VlcCA8LSBjKGNkNF9jZWxsc19wYm1jLCBvdGhlcl9jZWxscykKCiMgQ3JlYXRlIHRoZSBmaW5hbCBmaWx0ZXJlZCBTZXVyYXQgb2JqZWN0CmZpbHRlcmVkX3NldXJhdCA8LSBzdWJzZXQoQWxsX3NhbXBsZXNfTWVyZ2VkLCBjZWxscyA9IGNlbGxzX3RvX2tlZXApCgoKCmxpYnJhcnkoZ2dwbG90MikKbGlicmFyeShSQ29sb3JCcmV3ZXIpICAKCiMgQXNzdW1pbmcgeW91IGhhdmUgMTAgZGlmZmVyZW50IGNlbGwgbGluZXMsIGdlbmVyYXRpbmcgYSBjb2xvciBwYWxldHRlIHdpdGggMTAgY29sb3JzCmNlbGxfbGluZV9jb2xvcnMgPC0gYnJld2VyLnBhbCgxMCwgIlNldDMiKQoKIyBBc3N1bWluZyBBbGxfc2FtcGxlc19NZXJnZWQkY2VsbF9saW5lIGlzIGEgZmFjdG9yIG9yIGNoYXJhY3RlciB2ZWN0b3IgY29udGFpbmluZyBjZWxsIGxpbmUgbmFtZXMKZGF0YSA8LSBhcy5kYXRhLmZyYW1lKHRhYmxlKGZpbHRlcmVkX3NldXJhdCRjZWxsX2xpbmUpKQpjb2xuYW1lcyhkYXRhKSA8LSBjKCJjZWxsX2xpbmUiLCAiblVNSSIpICAjIENoYW5nZSBjb2x1bW4gbmFtZSB0byBuVU1JCgpuY2VsbHMgPC0gZ2dwbG90KGRhdGEsIGFlcyh4ID0gY2VsbF9saW5lLCB5ID0gblVNSSwgZmlsbCA9IGNlbGxfbGluZSkpICsgCiAgZ2VvbV9jb2woKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICBnZW9tX3RleHQoYWVzKGxhYmVsID0gblVNSSksIAogICAgICAgICAgICBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC45KSwgCiAgICAgICAgICAgIHZqdXN0ID0gLTAuMjUpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjZWxsX2xpbmVfY29sb3JzKSArIAogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gNDUsIGhqdXN0ID0gMSksCiAgICAgICAgcGxvdC50aXRsZSA9IGVsZW1lbnRfdGV4dChoanVzdCA9IDAuNSkpICsgICMgQWRqdXN0IHRoZSB0aXRsZSBwb3NpdGlvbgogIGdndGl0bGUoIkZpbHRlcmVkIGNlbGxzIHBlciBzYW1wbGUiKSArCiAgeGxhYigiQ2VsbCBsaW5lcyIpICsgICMgQWRqdXN0IHgtYXhpcyBsYWJlbAogIHlsYWIoIkZyZXF1ZW5jeSIpICAgICMgQWRqdXN0IHktYXhpcyBsYWJlbAoKcHJpbnQobmNlbGxzKQoKCmBgYAoKCgoKIyMgQ2VsbCB0eXBlIERpc3RyaWJ1dGlvbiB0byBjaGVjayBjbHVzdGVycwpgYGB7ciBEaXN0cmlidXRpb240LCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KCiMgV2UgY2FuIGFwcGx5IGl0IGxhdGVyIG9uIFIgb2JqIHRvIGdldCB0aGVzZSB0YWJsZXMgdG8gY29tcGFyZSBpdCB0byBkZWNpZGUgYWJvdXQgcmVzb2x1dGlvbi4KCgpqYW5pdG9yOjp0YWJ5bChmaWx0ZXJlZF9zZXVyYXRAbWV0YS5kYXRhLCBwcmVkaWN0ZWQuY2VsbHR5cGUubDEsIGNlbGxfbGluZSkKCgpqYW5pdG9yOjp0YWJ5bChmaWx0ZXJlZF9zZXVyYXRAbWV0YS5kYXRhLCBwcmVkaWN0ZWQuY2VsbHR5cGUubDIsIGNlbGxfbGluZSkKCgpqYW5pdG9yOjp0YWJ5bChmaWx0ZXJlZF9zZXVyYXRAbWV0YS5kYXRhLCBwcmVkaWN0ZWQuY2VsbHR5cGUubDMsIGNlbGxfbGluZSkKCmBgYAoKCiMgNC4gZmlsdGVyIEIgY2VsbHMgZnJvbSBMNApgYGB7ciBEaXN0cmlidXRpb241LCBmaWcuaGVpZ2h0PTYsIGZpZy53aWR0aD0xMH0KCiMgU2V0IGlkZW50aXR5IHRvIGNlbGxfbGluZQpJZGVudHMoZmlsdGVyZWRfc2V1cmF0KSA8LSAiY2VsbF9saW5lIgoKCiMgTG9hZCBuZWNlc3NhcnkgbGlicmFyaWVzCmxpYnJhcnkoU2V1cmF0KQpsaWJyYXJ5KGdncGxvdDIpCmxpYnJhcnkoUkNvbG9yQnJld2VyKQoKIyBJZGVudGlmeSBCIGNlbGxzIGluIEw0IHRvIGV4Y2x1ZGUKYl9jZWxsc19sNCA8LSBXaGljaENlbGxzKGZpbHRlcmVkX3NldXJhdCwgZXhwcmVzc2lvbiA9IAogIGdyZXBsKCJeQiIsIGZpbHRlcmVkX3NldXJhdCRwcmVkaWN0ZWQuY2VsbHR5cGUubDEpICYKICBncmVwbCgiXkIiLCBmaWx0ZXJlZF9zZXVyYXQkcHJlZGljdGVkLmNlbGx0eXBlLmwyKSAmCiAgZ3JlcGwoIl5CIiwgZmlsdGVyZWRfc2V1cmF0JHByZWRpY3RlZC5jZWxsdHlwZS5sMykgJgogIGNlbGxfbGluZSA9PSAiTDQiCikKCiMgSWRlbnRpZnkgY2VsbHMgdG8ga2VlcCAoZXhjbHVkaW5nIEIgY2VsbHMgaW4gTDQpCmNlbGxzX3RvX2tlZXAgPC0gc2V0ZGlmZihDZWxscyhmaWx0ZXJlZF9zZXVyYXQpLCBiX2NlbGxzX2w0KQoKIyBTdWJzZXQgdGhlIFNldXJhdCBvYmplY3Qgd2l0aCBzZWxlY3RlZCBjZWxscwpmaWx0ZXJlZF9zZXVyYXQgPC0gc3Vic2V0KGZpbHRlcmVkX3NldXJhdCwgY2VsbHMgPSBjZWxsc190b19rZWVwKQoKIyBQbG90IHRoZSBmaWx0ZXJlZCBkYXRhCgojIERlZmluZSBjZWxsIGxpbmUgY29sb3JzCmNlbGxfbGluZV9jb2xvcnMgPC0gYnJld2VyLnBhbChsZW5ndGgodW5pcXVlKGZpbHRlcmVkX3NldXJhdCRjZWxsX2xpbmUpKSwgIlNldDMiKQoKIyBHZW5lcmF0ZSBkYXRhIGZyYW1lIGZvciBwbG90dGluZwpkYXRhIDwtIGFzLmRhdGEuZnJhbWUodGFibGUoZmlsdGVyZWRfc2V1cmF0JGNlbGxfbGluZSkpCmNvbG5hbWVzKGRhdGEpIDwtIGMoImNlbGxfbGluZSIsICJuVU1JIikKCiMgQ3JlYXRlIHRoZSBiYXIgcGxvdApuY2VsbHMgPC0gZ2dwbG90KGRhdGEsIGFlcyh4ID0gY2VsbF9saW5lLCB5ID0gblVNSSwgZmlsbCA9IGNlbGxfbGluZSkpICsKICBnZW9tX2NvbCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGdlb21fdGV4dChhZXMobGFiZWwgPSBuVU1JKSwgcG9zaXRpb24gPSBwb3NpdGlvbl9kb2RnZSh3aWR0aCA9IDAuOSksIHZqdXN0ID0gLTAuMjUpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjZWxsX2xpbmVfY29sb3JzKSArCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA0NSwgaGp1c3QgPSAxKSwKICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KGhqdXN0ID0gMC41KSkgKwogIGdndGl0bGUoIkZpbHRlcmVkIENlbGxzIHBlciBTYW1wbGUiKSArCiAgeGxhYigiQ2VsbCBMaW5lcyIpICsKICB5bGFiKCJGcmVxdWVuY3kiKQoKcHJpbnQobmNlbGxzKQoKYGBgCgojIyBDZWxsIHR5cGUgRGlzdHJpYnV0aW9uIHRvIGNoZWNrIGNsdXN0ZXJzCmBgYHtyIERpc3RyaWJ1dGlvbjYsIGZpZy5oZWlnaHQ9MTAsIGZpZy53aWR0aD0xMH0KIyBXZSBjYW4gYXBwbHkgaXQgbGF0ZXIgb24gUiBvYmogdG8gZ2V0IHRoZXNlIHRhYmxlcyB0byBjb21wYXJlIGl0IHRvIGRlY2lkZSBhYm91dCByZXNvbHV0aW9uLgoKIyBBemltdXRoIGwxCmphbml0b3I6OnRhYnlsKGZpbHRlcmVkX3NldXJhdEBtZXRhLmRhdGEsIHByZWRpY3RlZC5jZWxsdHlwZS5sMSwgY2VsbF9saW5lKQoKIyBBemltdXRoIGwyCmphbml0b3I6OnRhYnlsKGZpbHRlcmVkX3NldXJhdEBtZXRhLmRhdGEsIHByZWRpY3RlZC5jZWxsdHlwZS5sMiwgY2VsbF9saW5lKQoKIyBBemltdXRoIGwzCmphbml0b3I6OnRhYnlsKGZpbHRlcmVkX3NldXJhdEBtZXRhLmRhdGEsIHByZWRpY3RlZC5jZWxsdHlwZS5sMywgY2VsbF9saW5lKQoKCmBgYAoKCgoKCiMgNS4gU2F2ZSB0aGUgU2V1cmF0IG9iamVjdCBhcyBhbiBSb2JqIGZpbGUKYGBge3Igc2F2ZVJPQkosIGVjaG89RkFMU0V9CgpzYXZlKGZpbHRlcmVkX3NldXJhdCwgZmlsZSA9ICIuLi8wLVJvYmovU1NfQ0Q0X1RjZWxsc19QQk1DMTB4X2ZpbmFsX2Zvcl9TQ1RfQXppbXV0aEFubm90YXRpb25fYW5kX0ludGVncmF0aW9uLnJvYmoiKQoKCmBgYAoKCgoKCgoK