- R-br - DINF UFPR

Localizar palavras e substituir-as por outras
by Alisson Lucrecio 18 Set '12

18 Set '12

Caros colegas da r-br, Boa tarde. Gostaria de sabe qual função localizar palavras e substitui-as por outras de interesse? Obrigado. Alisson Lucrécio da Costa

3 2

Erro em solve.default(v$varcov, xmat) GeoR likfit
by Rodrigo Sant'Ana 18 Set '12

18 Set '12

Caros, estou tentando ajustar um modelo gaussiano aos de biomassa de peixe (kg), e estou caindo em um problema que não consigo identificar o motivo! Abaixo segue o CMR, bem como o Erro resultante do procedimento. Obrigado. dados <- structure(list(coords = structure(c(-42.2064355778, -42.3032008959, -41.5452059035, -42.1580529187, -41.1420170777, -42.9724943467, -41.7306727634, -43.2224714187, -42.601560627, -43.0773234414, -43.0692596649, -42.61768818, -42.6499432861, -42.7708999338, -42.6338157331, -43.0128132293, -42.2951371194, -42.319328449, -41.7387365399, -43.0692596649, -42.3032008959, -42.3596473316, -42.4402850967, -43.4643847142, -43.6579153505, -43.7304893392, -43.8272546573, -42.7547723808, -44.2143159301, -43.7224255627, -43.117642324, -43.4724484907, -44.3272088013, -43.7466168922, -43.4160020551, -44.149805718, -43.8595097634, -44.6175047559, -44.1981883771, -44.2143159301, -44.6497598619, -45.0690762407, -44.9884384756, -44.5368669907, -44.8836093809, -44.5852496498, -44.7707165096, -45.1497140059, -45.5448390551, -44.948119593, -45.1497140059, -45.6254768203, -46.1334947407, -45.0610124642, -45.4722650665, -45.6254768203, -45.2787344301, -45.5932217142, -45.1900328884, -45.9480278809, -45.9641554339, -46.2141325059, -45.5125839491, -46.4963646839, -46.6253851082, -45.4077548544, -45.6335405968, -46.6495764377, -46.6737677673, -46.6253851082, -45.6496681498, -46.1415585173, -45.9641554339, -45.746433468, -46.5528111195, -45.99641054, -47.0688928165, -46.0447931991, -46.1334947407, -46.4318544718, -46.3834718127, -46.52861979, -46.6737677673, -47.4478903128, -46.4318544718, -47.2624234529, -47.4398265362, -47.310806112, -46.9882550514, -46.3270253771, -46.810851968, -46.8592346271, -46.8269795211, -47.6817398317, -47.1898494643, -47.4962729718, -47.4801454188, -46.9721274983, -47.6978673847, -47.0043826044, -47.6575485022, -47.7059311612, -47.8913980211, -47.4156352067, -47.843015362, -47.6736760552, -47.4720816423, -47.1575943582, -47.06082904, -47.2221045703, -47.6252933961, -47.7381862673, -47.3188698885, -47.1253392521, -47.3027423355, -47.1092116991, -47.8026964794, -47.9317169037, -47.6333571726, -47.5365918544, -47.576910737, -47.0527652635, -47.3591887711, -46.9801912749, -47.7220587143, -47.2059770173, -47.5688469605, -48.0849286575, -47.6575485022, -47.9800995628, -47.560783184, -48.1978215287, -47.9559082332, -48.4800537067, -48.6332654605, -48.6493930135, -48.5445639188, -48.4961812598, -48.3671608355, -48.2381404113, -47.9962271158, -48.8106685439, -48.7139032257, -47.9155893506, -47.9317169037, -48.5606914719, -48.4074797181, -48.9477527446, -47.9720357863, -48.5203725893, -48.7461583317, -48.0607373279, -49.0848369454, -48.8671149795, -49.0929007219, -48.3913521651, -48.9638802976, -49.1977298166, -48.5606914719, -49.2461124757, -48.3268419529, -48.7300307787, -48.1736301991, -49.407388006, -48.4477986007, -49.0525818393, -48.5284363658, -49.7138115136, -48.7622858848, -49.5041533242, -49.6815564075, -49.0929007219, -49.3590053469, -49.3509415704, -49.6654288545, -49.8831508204, -49.8428319378, -49.8347681613, -49.9637885855, -49.8105768318, -50.1492554454, -49.8670232674, -50.0363625742, -50.455678953, -50.3669774113, -50.3024671992, -50.4072962939, -50.3589136348, -50.5121253886, -50.6492095894, -50.3831049643, -50.9153142144, -50.6814646954, -51.1491637333, -50.7540386841, -51.4959061235, -51.4797785704, -51.3991408053, -51.6007352182, -51.5281612295, -23.1939349125, -23.2357031667, -23.26911777, -23.3359469767, -23.36936158, -23.3777152309, -23.4027761834, -23.4111298342, -23.4194834851, -23.4278371359, -23.4361907867, -23.4528980884, -23.4612517392, -23.4779590409, -23.5280809459, -23.5280809459, -23.5614955492, -23.5614955492, -23.5698492001, -23.5782028509, -23.6366784067, -23.7118612643, -23.7619831693, -23.7619831693, -23.7619831693, -23.8121050743, -23.8288123759, -23.8538733284, -23.8622269793, -23.8705806301, -23.9039952335, -23.9039952335, -23.9541171385, -24.079421901, -24.1378974568, -24.1462511077, -24.2715558702, -24.2966168227, -24.3634460294, -24.3885069819, -24.5472263477, -24.6808847611, -24.7393603169, -24.7477139677, -24.7727749202, -24.7811285711, -24.8479577778, -24.9649088894, -25.1152746045, -25.1570428586, -25.3241158753, -25.3241158753, -25.340823177, -25.3658841295, -25.390945082, -25.3992987328, -25.4160060345, -25.4577742887, -25.4744815904, -25.5078961937, -25.5078961937, -25.5078961937, -25.5747254004, -25.6499082579, -25.6499082579, -25.6666155595, -25.6666155595, -25.6833228612, -25.7334447662, -25.8169812746, -25.8503958779, -25.8503958779, -25.8587495287, -25.8671031796, -25.9339323863, -25.9506396879, -26.0675907996, -26.0759444504, -26.0759444504, -26.2263101655, -26.2430174671, -26.301493023, -26.4435050872, -26.5604561988, -26.5855171513, -26.6439927072, -26.7191755647, -26.7609438189, -26.8611876289, -26.9614314389, -26.9781387406, -27.0199069947, -27.0616752489, -27.0700288998, -27.0867362014, -27.0950898523, -27.1452117573, -27.1702727098, -27.2621628689, -27.362406679, -27.4542968381, -27.5211260448, -27.5294796957, -27.5796016007, -27.5963089023, -27.613016204, -27.6213698548, -27.6297235057, -27.6714917598, -27.6714917598, -27.8051501732, -27.8719793799, -28.0223450949, -28.0641133491, -28.0975279524, -28.2228327149, -28.2980155724, -28.3147228741, -28.3230765249, -28.3899057316, -28.3982593824, -28.4567349383, -28.4817958908, -28.5319177958, -28.5653323991, -28.6572225583, -28.6572225583, -28.6822835108, -28.77417367, -28.7825273208, -28.8493565275, -28.924539385, -28.9579539884, -28.9579539884, -29.1918562117, -29.3004536726, -29.3255146251, -29.3422219267, -29.4424657368, -29.4675266893, -29.5092949434, -29.5761241501, -29.5928314518, -29.6095387535, -29.6513070076, -29.684721611, -29.7264898651, -29.7766117701, -29.8267336752, -29.8434409768, -29.8434409768, -30.035574946, -30.035574946, -30.0439285969, -30.0689895494, -30.2110016135, -30.2444162169, -30.3195990744, -30.4365501861, -30.4449038369, -30.4950257419, -30.5284403453, -30.5451476469, -30.5702085994, -30.5952695519, -30.6119768536, -30.6286841553, -30.6537451078, -30.6704524094, -30.7205743145, -30.7289279653, -30.845879077, -30.8792936803, -31.0630739987, -31.0714276495, -31.1299032053, -31.3303908254, -31.4306346354, -31.4389882862, -31.5141711437, -31.6561832079, -31.6812441604, -31.6812441604, -32.5834384506, -32.6753286097, -32.7672187689, -32.8256943248, -32.9175844839, -33.0763038498, -33.1264257548, -33.2350232157, -33.2767914698, -33.3853889307, -33.3853889307, -33.5608155982, -33.7195349641, -33.8615470282, -33.8699006791, -33.8866079807, -33.9701444891), .Dim = c(200L, 2L)), data = c(70.347798763, 76.2305279637, 56.7935202123, 91.2747172018, 36.6574714598, 74.5676848522, 80.45544927, 63.7985345319, 98.6062126882, 76.1440959745, 78.0342424622, 103.211817252, 103.469894054, 101.543807018, 113.430501289, 96.7287932817, 119.081813263, 119.357653249, 91.6231191698, 101.038405309, 123.946425703, 127.264845138, 129.378943162, 101.547939554, 90.634442506, 95.7125681057, 94.798698159, 136.088211184, 83.1631121458, 105.907150692, 134.683969664, 120.236176625, 93.2955838471, 135.734205354, 134.793068333, 146.056344048, 156.071027198, 118.33813258, 176.165459485, 178.101033469, 141.733992676, 68.1513901758, 85.906096787, 160.683155231, 106.672886988, 153.732695652, 126.282924957, 71.171101282, 46.3063932604, 94.8066154966, 79.6077933159, 66.3022759053, 53.8001871156, 81.8971879578, 75.2725941278, 76.3895591688, 78.1481792228, 84.0088538824, 79.2339837253, 95.3285594814, 94.9048569376, 79.7371402431, 94.8361107261, 75.6734008995, 56.8583003566, 93.9377590667, 114.925013457, 58.7902731855, 63.3651704313, 87.7623028428, 138.243923048, 174.14679455, 176.924595279, 155.265108276, 130.599056347, 197.328254597, 43.1264520732, 221.751423068, 226.374023855, 231.22155796, 242.436720324, 225.68515641, 211.182339909, 30.7901884494, 268.635407904, 66.9770777548, 43.4651532814, 69.6985834637, 162.069565781, 201.17261197, 206.563852105, 197.412358817, 200.569088904, 36.5304326403, 134.230440546, 68.9447508778, 77.3092256209, 179.991789488, 46.7833022819, 178.837329648, 67.6041455343, 61.735642779, 30.9986850158, 136.367462528, 40.2420355711, 74.1549949452, 125.706398802, 182.134045089, 183.669620178, 179.176911697, 96.3292516587, 71.0485302772, 177.743449233, 187.3998893, 179.464856475, 177.962224878, 72.9460880655, 52.2011730864, 109.173183978, 129.8967523, 121.836259287, 147.025253713, 152.031441102, 124.646779658, 100.076387829, 133.011240036, 120.939828964, 64.3940362035, 114.996266127, 89.786038347, 118.008543994, 84.0282040216, 119.041263465, 44.0234944515, 41.8771712674, 48.2651512851, 73.4117997844, 87.9741710827, 142.802105252, 186.596945575, 228.955878351, 39.727571142, 59.5842097784, 237.196336949, 244.661132149, 111.885460631, 172.978503244, 31.9606458285, 262.185785391, 142.37968626, 71.9893685188, 261.41293353, 32.6187049615, 61.3911365603, 33.9234837497, 194.293576816, 59.9894516509, 40.4781868063, 136.131434893, 43.8955202545, 166.161274956, 101.623813429, 160.809302111, 39.3742142225, 134.084319986, 65.9772777256, 119.705366465, 27.1765413035, 90.2079815554, 41.3695883805, 32.9317985059, 61.9912390925, 53.1711001574, 53.8325000656, 52.2218112743, 49.1524441926, 62.6024730954, 67.1897383436, 66.0508521024, 69.425604381, 69.874832943, 73.9822629935, 75.6320803161, 98.9442750104, 108.632920447, 111.310932929, 117.076441302, 119.455819615, 125.059523339, 117.758307932, 117.561279774, 90.3341173828, 113.921889667, 67.6379513723, 99.9119129685, 58.2801014975, 63.4483498456, 66.4356707327, 57.6791519282, 62.3488331665), borders = structure(list(coords.x1 = c(-52.0676766304347, -52.0039483695653, -51.956152173913, -51.7968315217391, -51.6853070652174, -51.6853070652174, -51.5419184782609, -51.4144619565217, -51.2710733695652, -51.1914130434782, -51.0002282608695, -50.8090434782609, -50.7134510869565, -50.6178586956521, -50.5541304347826, -50.5541304347826, -50.5700624999999, -50.4744701086956, -50.4585380434783, -50.410741847826, -50.4426059782608, -50.4426059782608, -50.4585380434783, -50.3310815217391, -50.2195570652173, -50.0283722826087, -49.8531195652174, -49.6619347826086, -49.47075, -49.3751576086956, -49.2636331521738, -49.1202445652174, -48.9768559782608, -48.8971956521739, -48.8334673913043, -48.7060108695651, -48.6422826086957, -48.658214673913, -48.4670298913043, -48.3077092391304, -48.1643206521739, -48.0368641304347, -48.0687282608695, -48.0209320652174, -48.0368641304347, -47.9890679347826, -47.9253396739131, -47.8297472826087, -47.7341548913043, -47.7341548913043, -47.6066983695652, -47.4951739130435, -47.4633097826087, -47.4473777173913, -47.3039891304348, -47.2083967391304, -47.0331440217391, -46.9216195652173, -46.7304347826087, -46.5073858695652, -46.2524728260869, -45.9816277173912, -45.7426467391305, -45.4558695652173, -45.2965489130434, -45.1212961956522, -44.9619755434783, -44.8504510869565, -44.7070624999999, -44.5955380434783, -44.4202853260869, -44.0857119565217, -43.894527173913, -43.6555461956522, -43.3847010869565, -43.0341956521739, -42.795214673913, -42.6040298913043, -42.4128451086956, -42.221660326087, -42.0464076086956, -41.9348831521739, -41.7436983695652, -41.4728532608696, -41.2976005434782, -41.1542119565217, -41.0108233695652, -41.0426875, -41.1542119565217, -41.2816684782609, -41.329464673913, -41.4887853260869, -41.6640380434783, -41.8233586956522, -41.9986114130435, -42.2535244565217, -42.4447092391304, -42.5243695652173, -42.7792826086957, -42.9704673913043, -43.3050407608695, -43.4802934782608, -43.6396141304347, -43.8626630434782, -44.1016440217391, -44.2928288043478, -44.4840135869565, -44.627402173913, -44.7548586956522, -44.8982472826086, -45.0097717391304, -45.1372282608695, -45.2646847826086, -45.3602771739131, -45.4558695652173, -45.5992581521739, -45.5514619565217, -45.5992581521739, -45.7426467391305, -45.9816277173912, -46.1728125, -46.3162010869565, -46.5073858695652, -46.7145027173913, -46.8419592391304, -46.9375516304348, -46.9375516304348, -47.0331440217391, -47.0331440217391, -47.0012798913043, -46.9056875, -46.8578913043478, -46.8419592391304, -46.9056875, -46.985347826087, -47.0650081521739, -47.2243288043478, -47.4155135869565, -47.6066983695652, -47.7341548913043, -47.7819510869565, -47.8297472826087, -47.7819510869565, -47.8138152173913, -47.8616114130435, -47.9253396739131, -47.9253396739131, -48.0687282608695, -48.291777173913, -48.4829619565218, -48.7219429347826, -48.8175353260869, -48.9927880434782, -49.1680407608695, -49.3114293478261, -49.5344782608695, -49.6778668478261, -49.7734592391304, -49.9487119565217, -50.0921005434783, -50.1398967391304, -50.2673532608696, -50.2832853260869, -50.2832853260869, -50.3151494565217, -50.3629456521739, -50.5541304347826, -50.6815869565217, -50.8727717391304, -51.0002282608695, -51.2073451086957, -51.3029375, -51.3348016304348, -51.4463260869564, -51.3507336956521, -51.4941222826087, -51.5737826086956, -51.6215788043478, -51.7808994565218, -51.9083559782608, -52.0676766304347), coords.x2 = c(-34.511214673913, -34.3200298913043, -34.2244375, -34.1607092391304, -34.0651168478261, -33.9217282608695, -33.7464755434783, -33.6349510869565, -33.5074945652174, -33.3800380434783, -33.2685135869565, -33.0932608695652, -33.1410570652174, -33.0613967391304, -32.9020760869565, -32.7427554347826, -32.5675027173913, -32.471910326087, -32.2647934782608, -32.0736086956522, -31.8824239130434, -31.7390353260869, -31.5956467391304, -31.372597826087, -31.1495489130434, -30.9583641304347, -30.7193831521739, -30.5600625, -30.4644701086957, -30.2732853260869, -30.0980326086956, -29.9865081521739, -29.7953233695652, -29.6041385869565, -29.4288858695652, -29.3332934782608, -29.1739728260869, -29.0305842391305, -28.8393994565217, -28.648214673913, -28.4092336956522, -28.281777173913, -28.1224565217391, -28.0268641304347, -27.8834755434783, -27.7082228260869, -27.5966983695652, -27.3895815217391, -27.2302608695651, -27.0709402173912, -26.9275516304348, -26.7363668478261, -26.6248423913044, -26.5133179347825, -26.3858614130435, -26.1946766304347, -25.9875597826087, -25.8441711956521, -25.7007826086956, -25.5573940217391, -25.3662092391304, -25.2387527173913, -25.1272282608695, -24.9997717391304, -24.8563831521739, -24.6651983695652, -24.4580815217391, -24.250964673913, -24.1394402173913, -24.0597798913043, -23.8845271739131, -23.7570706521739, -23.6614782608695, -23.5340217391304, -23.4384293478261, -23.3428369565217, -23.2950407608695, -23.2153804347826, -23.1675842391305, -23.1038559782609, -23.0560597826087, -23.0560597826087, -23.1038559782609, -23.2153804347826, -23.2153804347826, -23.2153804347826, -23.1994483695652, -23.3109728260869, -23.4702934782608, -23.5658858695652, -23.6774103260869, -23.8048668478261, -23.8207989130435, -23.8207989130435, -23.7730027173913, -23.7730027173913, -23.7730027173913, -23.8685951086956, -23.9641875, -24.0597798913043, -24.2031684782609, -24.2828288043478, -24.3465570652173, -24.4262173913044, -24.5377418478261, -24.7129945652174, -24.8404510869565, -24.9519755434783, -25.143160326087, -25.3343451086956, -25.509597826087, -25.7007826086956, -25.796375, -25.8760353260869, -25.8441711956521, -25.9716277173912, -26.0353559782608, -26.1309483695651, -26.1787445652174, -26.5133179347825, -26.8000951086956, -27.0072119565218, -27.1346684782609, -27.182464673913, -27.2461929347825, -27.3417853260869, -27.4214456521739, -27.4851739130435, -27.5807663043478, -27.7560190217391, -27.915339673913, -28.0905923913044, -28.2499130434783, -28.3933016304348, -28.5844864130434, -28.7119429347826, -28.7119429347826, -28.807535326087, -28.9031277173913, -29.0943124999999, -29.3014293478261, -29.4766820652174, -29.7475271739131, -29.9865081521739, -30.2254891304348, -30.3688777173912, -30.5122663043478, -30.5122663043478, -30.6078586956521, -30.6875190217391, -30.7353152173913, -30.7353152173913, -30.7831114130435, -30.9105679347826, -31.0380244565217, -31.2292092391304, -31.4044619565217, -31.6115788043478, -31.8346277173912, -32.1054728260869, -32.3125896739131, -32.4878423913044, -32.6949592391304, -32.8861440217391, -33.125125, -33.3003777173913, -33.4915624999999, -33.6030869565217, -33.6827472826087, -33.6986793478261, -33.7942717391304, -33.7942717391304, -33.9057961956522, -33.9057961956522, -34.0173206521739, -34.0969809782609, -34.1607092391304, -34.2881657608695, -34.4474864130434, -34.4793505434782, -34.511214673913 )), .Names = c("coords.x1", "coords.x2"), class = "data.frame", row.names = c(NA, -181L))), .Names = c("coords", "data", "borders"), class = "geodata") ini.fit <- ini.fit <- expand.grid((2000:5000), (0:12)/12) mod.gau <- likfit(dados.geo, ini.cov.pars = ini.fit, cov.model = "gau", trend = "cte", fix.nugget = FALSE, fix.lambda = FALSE, lik.method = "ML") kappa not used for the gaussian correlation function likfit: searching for best initial value ... selected values: sigmasq phi tausq kappa lambda psiR initial.value "2952.00" " 0.83" " 0.00" " 0.50" " 1.00" " 1.00" status "est" "est" "est" "fix" "est" "fix" psiA initial.value " 0.00" status "fix" likelihood value: -319.704123974414 --------------------------------------------------------------- likfit: likelihood maximisation using the function optim. likfit: Use control() to pass additional arguments for the maximisation function. For further details see documentation for optim. likfit: It is highly advisable to run this function several times with different initial values for the parameters. likfit: WARNING: This step can be time demanding! --------------------------------------------------------------- trying another decomposition (svd) Erro em solve.default(v$varcov, xmat) : sistema é computacionalmente singular: condição recíproca número = 6.36899e-20 _______________ Rodrigo Sant'Ana Oceanógrafo

2 5

Recuperar coordenadas de locais de uma estrada a partir do Google Maps
by Augusto Ribas 18 Set '12

18 Set '12

Ola a todos, estou devolta aqui. Estou com o seguinte problema. Olhando os dados de atropelamentos de varias especies de bixinhos ao longo de uma estrada, eu tenho os seguintes dados: > head(dados) Especie km 1 Especie 5 27 2 Especie 3 32 3 Especie 3 64 4 Especie 4 40 5 Especie 1 51 6 Especie 6 57 A especie que foi atropelada, e qual o km da estrada. Eu gostaria de fazer um gráfico espacial, de onde estão acontecendo esses atropelamentos. Os dados em questão são da Rodovia ES-060. Eu fui no google maps e encontrei ela. Ai peguei as coordenadas de alguns pontos e importei no R, por exemplo: Lon Lat Local 1 -40.32012 -20.39675 Inicio 2 -40.32477 -20.40193 KM 12 3 -40.39341 -20.54320 KM 30 4 -40.51056 -20.63095 KM 49 5 -40.55583 -20.73349 KM 66 6 -40.55161 -20.73921 Fim Peguei esses valores assim, eu ia no google street view, ai olhava na estrada as plaquinha km 30, ai ia la e marcava um ponto. Olhava outra plaquinha e marcava outro ponto. Se você traçar uma rota do fim para o inicio, da pra ver a estrada inteira (se a rota for ao contrario, não da muito certo, ja que tem que voltar pra uma rotatória no trajeto.). Eu gostaria de saber se tem uma estrategia simples, usando essas informações de saber as coordenadas para cada km do 11 ao 67. Pra completar minha primeira planilha, dados, com coordenadas. Tipo importar a rota e marcar pontos usando estes de referencia, pra saber a coordenada do km 11, do km 12 e assim por diante. Dai eu veria num mapa onde estão os atropelamentos. Segue um exemplo dos dados e um plot do mapa da estrada. library(RgoogleMaps) #exemplo de dados especie<-paste("Especie",1:6) km<-seq(11,67,by=1) dados<-data.frame(Especie=NA,km=NA) for(i in 1:20) { dados[i,1]<-sample(especie,1) dados[i,2]<-sample(km,1) } head(dados) #Locais no google maps em um data.frame dados.pos<-structure(list(Lon = c(-40.3201206572012, -40.3247698942156, -40.3934139596311, -40.5105609306218, -40.5558289282114, -40.551610815962), Lat = c(-20.3967530595136, -20.4019278830141, -20.5432001637965, -20.6309481844606, -20.733489883797, -20.7392078853725), Local = c("Inicio", "KM 12", "KM 30", "KM 49", "KM 66", "Fim")), .Names = c("Lon", "Lat", "Local"), row.names = c(NA, -6L), class = "data.frame") dados.pos #Baixando o maps lat = dados.pos[c(1,6),2] lon = dados.pos[c(1,6),1] center = c(lat=mean(lat), lon=mean(lon)) zoom <- min(MaxZoom(range(lat), range(lon))) MyMap <- GetMap(center=center,zoom=zoom,maptype ="roadmap") #plot PlotOnStaticMap(MyMap) #pontos conhecidos TextOnStaticMap(MyMap,claudia.coord[,2],claudia.coord[,1],labels = claudia.coord[,3]) Desde ja agradeço a atenção Abraços a todos -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 2

Dúvida função no R
by Felipe Buchbinder 17 Set '12

17 Set '12

Prezados professores, Estou com um problema em uma função no R. A função é a seguinte: resp <- function(subject,item,trait) > { > > p.vec = NULL > > for(category in 1:k) # k = 5 > { > p = ogive ( theta[subject,trait] , a[item] , b[item,category] , D ) > # aqui é uma logit (D=1) ou probit (D=1.7) > p.vec = c(p.vec,p) > } > > #Generates a random number > > seed<-runif(1,0.2,0.8) > > #Selects an answer > > temp = which ( p.vec>seed ) > if( length (temp) > 0 ) { resp = max (temp) } else { resp = 1 } > resp > } Esta função está inserida em um loop e deveria me retornar um número inteiro entre 1 e 5. Da primeira vez, isto acontece. Mas a partir da segunda, ela começa a me retornar valores como 20, 16,... Alguma idéia de porque isso acontece??? Desde já, obrigado! Felipe Buchbinder

7 25

Comando for em cadeia
by Guilherme Heiden 17 Set '12

17 Set '12

Olá, eu talvez não tenha sido claro no assunto do e-mail, mas na minha concepção o problema é este: um comando for dentro de outro. Eu estou trabalhando com a seguinte matriz: [,1] [,2] [,3] [,4] [,5][1,] 20 5 11 7 49[2,] 18 9 10 2 45[3,] 11 35 30 15 7[4,] 10 3 7 4 26 E a partir desta eu gostaria de gerar outra matriz com a seguinte configuração: [,1] [,2] [,3] [,4] [,5][1,] dv11 dv12 dv13 dv14 dv15[2,] - dv22 dv23 dv24 dv25[3,] - - dv33 dv34 dv35[4,] - - - dv44 dv45[5,] - - - - dv55 Para achar os devidos valores, eu criei a seguinte função:read.csv2("C:/multivar.csv",header=T) multivar<-read.csv2("C:/multivar.csv",header=T) attach(multivar) ##Criar DV11: vetor1<-NULL for(i in 1:4) { vetor1[i]<-(multivar[i,1]-multivar[i,1])^2 } dv11<-(sum(vetor1)^0.5) rm(vetor1) ##Criar DV12: vetor1<-NULL for(i in 1:4) { vetor1[i]<-(multivar[i,1]-multivar[i,2])^2 } dv12<-(sum(vetor1)^0.5) rm(vetor1) E assim sucessivamente. Minha dúvida: Há a possibilidade de criar um comando "for" para os valores de dv serem criados automaticamente, ao invés de gerar um de cada vez, manualmente?

3 6

Selecionar colunas de interesse apartir de uma matriz...
by andrebvs 16 Set '12

16 Set '12

3 6

R + Banco de dados orientado a objetos (db4o)
by Roney Fraga 16 Set '12

16 Set '12

Caros, É possível trabalhar com o R usando banco de dados orientado a objetos (db4o)? Caso positivo, qual pacote que faz essa conexão? Atenciosamente -- Roney

2 1

não-convergência glm poisson
by Heloíse Pavanato 15 Set '12

15 Set '12

Caros, Gostaria de alguma sugestão sobre o que fazer quando um glm (Poisson) não converge. Erros: Mensagens de aviso: 1: glm.fit: algoritmo não convergiu 2: glm.fit: taxas ajustadas numericamente 0 ocorreu Obrigada. Heloise --- Oceanógrafa Heloise Pavanato Julião Pós-graduação em Oceanogra Biológica Instituto de Oceanografia Universidade Federal do Rio Grande

4 7

Usando table() com datas.
by Augusto Ribas 14 Set '12

14 Set '12

Boa noite a todos. Estou com uma duvida aqui sobre como fazer tabelas com dados de tempo. Bom eu tenho um dataframe que é uma lista de bixos atropelados na estrada. Cada linha é um registro, e tem a especie do bixinho que morreu e a data e outras coisas. Mais ou menos assim: > head(dados) Especie Data 1 Especie 1 2003-06-01 2 Especie 3 2001-06-01 3 Especie 1 2003-07-01 4 Especie 2 2004-02-01 5 Especie 3 2000-07-01 6 Especie 1 2002-06-01 Então tem varias especies, e varios anos de coleta. Eu gostaria de fazer uma tabela pra saber quantos individuos de cada especie morreram por mes, separando os anos. Eu fiz isso usando o comando table(), porem tem alguns meses que não teve registro para nenhuma especie. Dai a minha tabela fica incompleta. Eu gostaria que do periodo inicial dos registros até o final, tivesse todos os meses, mesmo que aparecendo zero para todas as especies. Segue um CMR pra exemplificar melhor o exemplo #gerando dados de exemplo: set.seed(5) datas<-seq(from =as.Date('1/1/2000',format='%d/%m/%Y'),to=as.Date('1/1/2005',format='%d/%m/%Y'),by='month') especies<-paste("Especie",1:3) dados<-data.frame(Especie=NA,Data=as.Date('1/1/2000',format='%d/%m/%Y')) for(i in 1:30) { dados[i,1]<-sample(especies,1) dados[i,2]<-sample(datas,1) } head(dados) #Fazendo a tabela table(format(dados$Data,'%m/%Y'),dados$Especie) > table(format(dados$Data,'%m/%Y'),dados$Especie) Especie 1 Especie 2 Especie 3 01/2001 1 0 0 01/2002 0 1 0 01/2003 0 1 0 02/2000 1 0 0 02/2001 1 0 0 02/2004 0 1 0 03/2001 0 0 1 03/2003 0 0 1 04/2000 0 1 0 04/2001 1 0 0 04/2004 0 1 0 05/2000 0 0 1 06/2001 0 0 1 06/2002 3 0 0 06/2003 1 0 0 07/2000 0 0 1 07/2003 1 0 0 07/2004 0 1 0 08/2003 0 0 1 09/2000 1 0 0 09/2002 0 0 1 10/2002 0 0 1 11/2002 1 1 0 11/2004 1 1 0 12/2001 0 0 1 12/2004 0 0 1 Por exemplo ali onde destaquei em vermelho, tinha que ter o mes 2 pro ano de 2002 e 2003 Nos dados eu vejo que o trabalho acabou em 2004, dezembro > max(dados$Data) [1] "2004-12-01" Mas começou em 2000, fevereiro. > min(dados$Data) [1] "2000-02-01" Ou seja ta faltando muitos meses. Tem como o table completar isso automaticamente? Eu entendo que ele não completou esses meses pois não tem um registro, e eles não entraram como um nivel de fator. Mas eu gostaria que todos os meses, mesmo sem registros tivessem ali. A principio eu estava tentando fazer um vetor com os meses/anos, ai fazer uma tabela com zeros e usar o resultado do table pra preencher. mas imagino que essa não deve ser a solução mais inteligente. Alguém tem alguma sugestão? Abraços Augusto Ribas -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 3

Problema com read.xls do pacote gdata
by Edwin Isensee 14 Set '12

14 Set '12

Prezados, Estou tentando ler uma planilha Excel simples, com uma única aba, usando o comando read.xls do pacote gdata e sempre recebo a seguinte mensagem de erro: Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 2 Alguém já passou por isso? Alguma dica? Grato, Edwin

4 4

Dúvida- resposta do Walmes
by Fernando A. Souza 13 Set '12

13 Set '12

Caro walmes, Ha algum tempo atras postei uma dúvida a respeito de testes de comparações multiplas para dados desbalanceados. E você então postou a seguinte resposta a qual reproduzo novamente abaixo e que foi de grande ajuda. Depois de entender tudo que você postou fiquei com a seguinte dúvida: Em seu exemplo no final da resposta, vc executa a função glht do pacote multicomp, baseado em seu contraste de interesse. A função acima utiliza o procedimento step-one (Default), o qual segundo o autor do livro do pacote multicomp é um procedimento conservador. Com base em sua abordagem como fica esta questão de escolher um teste menos conservador (ex. SNK) ou mais conservador (Ex.Tukey)?Não importa mais qual o tipo de teste utilizado, desde que ele controle a taxa de erro? grato pela atenção "Eu prefiro usar a multcomp porque lá eu tenho controle total das comparações que quero fazer. Eu aprendi a usá-la antes de pacotes mais dedicados estarem disponíveis. Basicamente, qual o problema das comparações múltiplas? É o fato de você testar muitas hipóteses e isso elevar rapidamente a taxa de erro tipo I global. Qual a solução? Lançar mão de métodos que mantenham essa taxa global próximo do valor nominal, na sua maioria 5%. Tudo o que se conhece sobre testes de comparações múltiplas são conceitualmente maneiras de fazer isso. Cada um lança mão de um princípio para fazê-lo e por isso possem alguma divergência. Alguns deles se baseiam no fato das estimativas terem mesma precisão e serem independentes, o que ocorre em todo experimento de fatores categóricos completamente cruzados e balanceado. Nesses casos a matriz de covariância das estimativas das médias é diagonal. Com fatores cruzados mas não balanceado (diferente número de repetições), será diagonal com precisões diferentes para cada nível do fator. Com cruzamento incompleto, chamado de perda de caselas, você terá elementos fora da diagonal não nulos. É importante que você entenda que médias amostrais só são estimadores admissíveis para médias (parâmetro) em raras situações. Deve usar médias obtidas matricialmente, que foi batizado por aí de médias ajustadas. Eu pessoalmente não vejo porque o uso da palavra ajustado se de fato é o procedimento correto, o caso geral para estimar uma média. As médias amostrais é que deviam ser chamadas de "médias possivelmente desajustadas/viesadas". Eu gosto da multcomp porque o que tá implementado lá resolve todos os casos e faz o que tem que fazer, controlar a taxa nominal de significância (conduza um estudo de simulação e verifique se for o caso). Os tipos de somas de quadrados em si não têm importância porque as incertezas sobre as médias estão na matriz de covariância das estimativas e são invariantes a soma de quadrado. Além do mais, essa coisa de soma de quadrado é particular da anova e comparações múltiplas de estimativas podem ser aplicadas em qualquer situação, seja em GLM e até mesmo para parâmetros em modelos de regressão não linear. Esqueça uso de letras seguindo as médias. Faça um gráfico dos contrastes do teu interesse e verifique se o valor nulo está contido no intervalo de confiança. Para ilustrar, rode o CMR #------------------------------------------------------------------------------------------ require(doBy) # caso 1 # experimento balanceado e completamente cruzados níveis de bloc e trat da <- expand.grid(trat=gl(3,1), bloc=gl(2,2)) da$y1 <- rnorm(nrow(da), 10, 0.1) str(da) m0 <- lm(y1~bloc+trat, da) ctr <- popMatrix(m0, effect="trat") dim(ctr) V0 <- vcov(m0) dim(V0) vm0 <- ctr%*%V0%*%t(ctr) vm0 # mesma precisão e sem correlação #------------------------------------------------------------------------------------------ # caso 2 # experimento desbalanceado e completamente cruzados níveis de bloc e trat da$y2 <- da$y1 da$y2[1:2] <- NA # trat 1 e 2 têm apenas 2 repetições (desbalanceamento) m1 <- lm(y2~bloc+trat, da) V1 <- vcov(m1) vm1 <- ctr%*%V1%*%t(ctr) vm1 # trat 1 e 2 tem maior variância mas fora da diagonal é zero! #------------------------------------------------------------------------------------------ # caso 3 # experimento desbalanceado e completamente cruzados níveis de bloc e trat da$y3 <- da$y1 da$y3[c(1,4,9,12)] <- NA # trat 1 não observado no bloc 1 (perda completa de combinação) m2 <- lm(y3~bloc+trat, da) V2 <- vcov(m2) vm2 <- ctr%*%V2%*%t(ctr) vm2 # diferença de precisão round(cov2cor(vm2),3) # agora tem presença de correção devido não ortogonalidade #------------------------------------------------------------------------------------------ # em qualquer caso, multcomp pode comparar as médias, escolha 1 deles e faça comb <- combn(1:nlevels(da$trat), 2) # 1 vs 2, 1 vs 3, 2 vs 3 mcomp <- ctr[comb[1,],]-ctr[comb[2,],] # matriz de contrastes 2 a 2 require(multcomp) g0 <- glht(m2, linfct=mcomp) summary(g0) plot(g0) # diferenças nas amplitudes refletem o número de repetições #------------------------------------------------------------------------------------------ À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br <mailto:walmes@ufpr.br> twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218

2 3

Auto-correlação residual utilizando a gnls
by Tales Fernandes 13 Set '12

13 Set '12

Olá Pessoal, gostaria de contar com a experiência de vocês. Estou com problemas no momento de estimar o parâmetro de auto-correlação residual utilizando a gnls. Fiz o teste Durbin-Watson e como conclusão os resíduos são dependentes. Mas quando estimo os parâmetros de acordo com a fórmula abaixo a função me retorna que a estimativa do parâmetro de autocorrelação, "Phi1", é zero. Minha variável independente é daf e a variável dependente é mfs2. reg3=gnls(mfs2~((alfa)/(1 + exp(beta - gama*daf))), start=c(alfa=1.5,beta=3,gama=0.015), weights=varExp(form=~daf),correlation=corAR1(form=~daf)) Li a documentação da "nlme" mas não consigo compreender a diferença de declarar a correlação de outras formas, como por exemplo: reg2=gnls(mfs2~((alfa)/(1 + exp(beta - gama*daf))), start=c(alfa=1.5,beta=3,gama=0.015), weights=varExp(form=~daf),correlation=corAR1(form=~1|daf)) reg3=gnls(mfs2~((alfa)/(1 + exp(beta - gama*daf))), start=c(alfa=1.5,beta=3,gama=0.015), weights=varExp(form=~daf),correlation=corAR1()) Se possível alguém poderia esclarecer qual eu deveria utilizar e qual a diferença entre estas formas. DETALHE: Apenas na reg3 o parâmetro Phi1 é diferente de zero, mas utilizando este método ele dá diferente de zero mesmo quando os resíduos são independentes. Grato. Tales Jesus Fernandes Doutorando em Estatística UFLA Universidade Federal de Lavras

3 5

[Dúvida] Pacote em R.
by Pedro Rafael 13 Set '12

13 Set '12

Pessoal a pessoa que cria ou mantém algum pacote no R ou qualquer outra linguagem/software terá esse trabalho reconhecido em concursos assim como uma publicação em revista? Existe algum reconhecimento de pontuação para esses trabalhos? Falo mais com respeitos à concursos acadêmicos. [ ], Pedro Rafael Diniz Marinho.

3 4

Re: [R-br] Ajuda com substituição de valores em um data frame
by Walmes Zeviani 12 Set '12

12 Set '12

Olá. É mais interessante que você envie suas dúvidas para a lista pois (1) você pode obter diversas soluções (2) você pode ser respondida em menos tempo (3) sua mensagem pode interessar a mais pessoas (4) ela fica no histórico da lista para consulta futura. Segue uma possível solução usando match(). da <- matrix(sample(c("a","c","g","t"), 20, replace=TRUE), ncol=5) da typeof(da) # matrix de character, converter seu data.frame para matriz com as.matrix() atual <- c("a","c","g","t") recod <- c(1,2,1,2) data.frame(atual, recod) indic <- match(da, atual) # matriz de indices db <- matrix(recod[indic], nrow=nrow(da)) db # matriz recodificada, as.data.frame(db) # pode ser convertida para data.frame com as.data.frame À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/9/11 <elisabpereira(a)gmail.com> > Olá, > > tudo bem? > Desculpe enviar um e-mail diretamente para você, mas já enviei minha > dúvida para o fórum do R e não consegui resolver meu problema. Vi que vc > ajuda muitas pessoas por aqui, e se não for pedir muito, também queria uma > ajudinha. > > Tenho um data frame de 21 linhas e 329 colunas. As células do data frame > são compostas pelos caracteres "a", "c", "g" ou "t". Quero fazer com que as > células preenchidas pelos caracteres "a" ou "g" sejam substituídas pelo > numeral 1, e as células com os caracteres "c" ou "t" sejam substituídas por > 2. > > Quando uso: > ifelse(dataframe=="a",1,0) > > Nesse caso, as células com "a" são substituídas pelo valor 1 e as demais > são substituídas por 0. > Tem alguma forma de fazer com q elas sejam substituídas por 1 e que o > restante permaneça como está? Desta forma eu poderia usar o ifelse quatro > vezes, uma para cada caracter. > Também tentei usar: > if (dataframe=="a") {1} else > if (dataframe=="g") {1} else > if (dataframe=="t") {2} else > if (dataframe=="c") {2} > Mas recebo o seguinte warning: "a condição tem comprimento > 1 e somente o > primeiro elemento será usado" > > Também tentei usar o recode, mas não consegui fazer funcionar. > > Muito obrigada pela ajuda!!! >

2 1

título Gráfico para dois graficos ou mais
by Edson Lira 12 Set '12

12 Set '12

Pessoal, tenho três gráficos em uma mesma janela. Vejam a rotina abaixo: y<-c(1,2,3,4,5) x<-c(2,3,4,5,6) par(mfrow=c(1,2)) plot(y,main="y") plot(x,main="x") tenho os dois gráficos lado a lado com seus títulos respectivos. Como faço para incluir um título para os dois, tipo: Valores de X e Y, centralizado entre os dois gráfics. Boa noite Edson Lira Estatístico Manaus-Amazonas

3 4

Combinar listas...
by andrebvs 12 Set '12

12 Set '12

2 3

Edição de figuras no R
by Ari Clecius 11 Set '12

11 Set '12

Prezados colegas, estou submetendo um artigo em uma revista que pede que as figuras sejam editáveis, mas o formato fornecido pelo R não permite, alguém teria alguma solução? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)87687786 (085)33669042

4 5

OT Fwd: Edital de Seleção Mestrado ENCE
by Leandro Marino 11 Set '12

11 Set '12

Enviado de um dispositivo móvel. Sent from a mobile device. ---------- Mensagem encaminhada ---------- De: "Suzana Cavenaghi" <suzana_cavenaghi(a)uol.com.br> Data: 10/09/2012 23:49 Assunto: Edital de Seleção Mestrado ENCE Para: <cpgence(a)ibge.gov.br> Colegas, Envio o Edital de Seleção do Mestrado em Estudos Populacionais e Pesquisas Sociais, para Turma 2013, da ENCE e peço colaboração para ampla divulgação em listas que possam ter interessados. Caso tenham dificuldades em abrir os arquivos anexos, estes estão disponíveis em http://www.ence.ibge.gov.br/web/ence/processo-seletivo-mestrado. Muito obrigada. Abraços. Suzana

1 0

Re: [R-br] Problemas com rgl!
by Ivan Bezerra Allaman 10 Set '12

10 Set '12

Boa tarde Oswando! Fiz aquilo que recomendastes, não entanto não surgiu nenhum efeito, veja: ivan-Inspiron-N5110 ivan # sudo apt-get build-dep r-cran-rgl Lendo listas de pacotes... Pronto Construindo árvore de dependências Lendo informação de estado... Pronto A escolher 'rgl' como pacote pacote de código fonte em vez de 'r-cran-rgl' 0 pacotes atualizados, 0 pacotes novos instalados, 0 a serem removidos e 0 não atualizados. Outra sugestão? Grato! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] Erro GVim+R
by Roney Fraga 09 Set '12

09 Set '12

Jakson, Acabei optando pelo plugin do R para o gedit. Obrigado novamente pela ajuda e pela elaboração tando do vim-r quanto do Rgedit. Atenciosamente -- Roney

3 2

Re: [R-br] Plot
by Ivan Bezerra Allaman 08 Set '12

08 Set '12

Enquanto o Walmes não posta algo, vai se contentando com isso. Dê uma olhada na ggplot2, pois é outra opção à graphics. y <- rnorm(20) x <- rep(c(0,10,20,30),rep(5,4)) medias <- tapply(y,x,mean) desv <- tapply(y,x,sd) plot(tapply(y,x,mean) ~ seq(0,30,10),ylim=c(-2,2)) segments(0,medias[1],0,medias[1]-desv[1]) segments(0,medias[1],0,medias[1]+desv[1]) segments(10,medias[2],10,medias[2]-desv[2]) segments(10,medias[2],10,medias[2]+desv[2]) segments(20,medias[3],20,medias[3]-desv[3]) segments(20,medias[3],20,medias[3]+desv[3]) segments(30,medias[4],30,medias[4]-desv[4]) segments(30,medias[4],30,medias[4]+desv[4]) \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Re: [R-br] Problemas com rgl!
by Ivan Bezerra Allaman 08 Set '12

08 Set '12

Já está instalado! mesa-common-dev; libegl1-mesa-dev; libgl1-mesa-dev; libglu1-mesa-dev; libopenvg1-mesa-dev; Outras não, libgels1-mesa-dev; liblgles2-mesa-dev; libglw1-mesa-dev; libosmesa6-dev Outra sugestão? Grato! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature} ________________________________ De: Ricardo Borges Gama Neto <ricardobgneto(a)gmail.com> Para: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>; Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Enviadas: Sábado, 8 de Setembro de 2012 16:51 Assunto: Re: [R-br] Problemas com rgl! Instale mesa-devel. Como você esta usando um derivado Debian deve ser algo como sudo apt-get install mesa-devel ou algo assim, ir de Synaptic é um caminho mais confiável. Abs, Prof Ricardo Borges Gama Neto Universidade Federal de Pernambuco Departamento de Ciência Política Programa de Pós-graduação em Ciência Política Enviado via iPhone Em 08/09/2012, às 15:30, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> escreveu: Boa tarde senhores! > > >Ao rodar um gráfico em 3D com a rgl percebi que estava muito lento, travando. Então descobri que não estava instalado o driver de vídeo da minha máquina. Ao instalar e tentar usar a rgl obtive erros. Então removi o pacote que estava dentro do usr blablab e tentei reinstalar o pacote sem sucesso com o seguinte erro: > > >* installing *source* package ‘rgl’ ... >** package ‘rgl’ successfully unpacked and MD5 sums checked >checking for gcc... gcc -std=gnu99 >checking whether the C compiler works... yes >checking for C compiler default output file name... a.out >checking for suffix of executables... >checking whether we are cross compiling... no >checking for suffix of object files... o >checking whether we are using the GNU C compiler... yes >checking whether gcc -std=gnu99 accepts -g... yes >checking for gcc -std=gnu99 option to accept ISO C89... none needed >checking how to run the C preprocessor... gcc -std=gnu99 -E >checking for gcc... (cached) gcc -std=gnu99 >checking whether we are using the GNU C compiler... (cached) yes >checking whether gcc -std=gnu99 accepts -g... (cached) yes >checking for gcc -std=gnu99 option to accept ISO C89... (cached) none needed >checking for libpng-config... yes >configure: using libpng-config >configure: using libpng dynamic linkage >checking for X... libraries , headers >checking GL/gl.h usability... yes >checking GL/gl.h presence... yes >checking for GL/gl.h... yes >checking GL/glu.h usability... yes >checking GL/glu.h presence... yes >checking for GL/glu.h... yes >checking for glEnd in -lGL... yes >checking for gluProject in -lGLU... yes >checking for freetype-config... yes >configure: using Freetype and FTGL >configure: creating ./config.status >config.status: creating src/Makevars > >. >. >. > >Error : .onLoad failed in loadNamespace() for 'rgl', details: > call: dyn.load(file, DLLpath = DLLpath, ...) > error: impossível carregar objeto compartilhado '/usr/local/lib/R/site-library/rgl/libs/rgl.so': > libGL.so.1: cannot open shared object file: No such file or directory >Erro: loading failed >Execução interrompida >ERROR: loading failed >* removing ‘/usr/local/lib/R/site-library/rgl’ > >The downloaded source packages are in > ‘/tmp/RtmpS6wZHV/downloaded_packages’ >Mensagens de aviso perdidas: >In install.packages("rgl") : > installation of package ‘rgl’ had non-zero exit status > > >Já pesquisei em vários lugares, instalar várias coisas e o erro persiste. Grato se alguém puder ajudar. Utilizo o mint 13 maya ubuntu. > > >Abraços! > > > >\begin{signatu ><<>>= >Prof. Dr. Ivan Bezerra Allaman >Universidade Estadual de Santa Cruz >Departamento de Ciências Exatas e Tecnológicas >Ilhéus/BA - Brasil >Fone: +55 73 3680-5596 >E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com >@ >\end{signature} _______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Problemas com rgl!
by Ivan Bezerra Allaman 08 Set '12

08 Set '12

Boa tarde senhores! Ao rodar um gráfico em 3D com a rgl percebi que estava muito lento, travando. Então descobri que não estava instalado o driver de vídeo da minha máquina. Ao instalar e tentar usar a rgl obtive erros. Então removi o pacote que estava dentro do usr blablab e tentei reinstalar o pacote sem sucesso com o seguinte erro: * installing *source* package ‘rgl’ ... ** package ‘rgl’ successfully unpacked and MD5 sums checked checking for gcc... gcc -std=gnu99 checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc -std=gnu99 accepts -g... yes checking for gcc -std=gnu99 option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -std=gnu99 -E checking for gcc... (cached) gcc -std=gnu99 checking whether we are using the GNU C compiler... (cached) yes checking whether gcc -std=gnu99 accepts -g... (cached) yes checking for gcc -std=gnu99 option to accept ISO C89... (cached) none needed checking for libpng-config... yes configure: using libpng-config configure: using libpng dynamic linkage checking for X... libraries , headers checking GL/gl.h usability... yes checking GL/gl.h presence... yes checking for GL/gl.h... yes checking GL/glu.h usability... yes checking GL/glu.h presence... yes checking for GL/glu.h... yes checking for glEnd in -lGL... yes checking for gluProject in -lGLU... yes checking for freetype-config... yes configure: using Freetype and FTGL configure: creating ./config.status config.status: creating src/Makevars . . . Error : .onLoad failed in loadNamespace() for 'rgl', details: call: dyn.load(file, DLLpath = DLLpath, ...) error: impossível carregar objeto compartilhado '/usr/local/lib/R/site-library/rgl/libs/rgl.so': libGL.so.1: cannot open shared object file: No such file or directory Erro: loading failed Execução interrompida ERROR: loading failed * removing ‘/usr/local/lib/R/site-library/rgl’ The downloaded source packages are in ‘/tmp/RtmpS6wZHV/downloaded_packages’ Mensagens de aviso perdidas: In install.packages("rgl") : installation of package ‘rgl’ had non-zero exit status Já pesquisei em vários lugares, instalar várias coisas e o erro persiste. Grato se alguém puder ajudar. Utilizo o mint 13 maya ubuntu. Abraços! \begin{signatu <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

rJava
by William Amorim 08 Set '12

08 Set '12

Oi, estou tentando usar o pacote rJava mas, quando tento carregá-lo, aparece esse erro: > require(rJava) Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: stop("No CurrentVersion entry in '", key, "'! Try re-installing Java and make sure R and Java have matching architectures.") erro: objeto 'key' não encontrado Alguém teria uma solução?

4 6

Mensagem de erro
by Renato Braga 07 Set '12

07 Set '12

Lista, o que significa o erro "Unrecognized record type 7, subtype 18 encountered in system file" na importação de uma base de dados do SPSS (.sav)? Obrigado,

2 1

Re: [R-br] Séries Temporais dividir dados em teste e treino
by Hugo Sena 06 Set '12

06 Set '12

Obrigado Alisson<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=user_nodes…> Já tentei aplicar e não consigo. Alguém pode ajudar como se faz neste caso, tipo (80% de teste 20% de treino). Será que posso e/ou é mais fácil se dividir os dados manualmente em 2 ficheiros? O exemplo poderá ser com estes dados. http://robjhyndman.com/tsdldata/data/cola.dat Já googlei e não consigo aplicar. Cumps.

2 1

Como substituir caracteres de um data frame por valores numéricos?
by Elisa Barreto 06 Set '12

06 Set '12

Bom dia, tenho um data frame composto por quatro caracteres "a","g","t" e "c". Nas linhas do data frame estão as espécies (sp1, sp2, sp3 ...) e nas colunas estão os sítios, que vão de 1 a 320 e podem conter qualquer um dos quatro caracteres. Gostaria de transformar todos os caracteres "a" e "g" do data frame em 1 e todos os caracteres "t" e "c" em 2. Como faço isso? Tentei usar o *transform*, mas não consegui. Obrigada, -- Elisa

2 1

Plot
by Stephan M. CARVALHO 05 Set '12

05 Set '12

Boa noite amigos do R, Estou precisando de uma ajuda para plotar uma figura! Na verdade é algo bem simples, pois preciso somente plotar meus dados médios e o desvio padrão. O grande problema é que quando uso o pacote gplots ele não mantem a proporção no eixo x, o que acaba distorcendo os dados. Abaixo tem uma tentativa minha usando a função "plot", mas ainda não consegui uma saída pra plotar a média + desvio padrão. Por acaso alguém tem uma solução???? Att, Stephan par(cex.lab=2,cex.axis=2,bty="n",yaxs="r",xaxs="r",lwd=2,family="",mar=c(5,5,2.5,2.5),mex=1.7) plot(M1~Conc,ylim=c(0,0.015),yaxt="n",xaxt="n",type="b",ylab="v(OD/min)",xlab="[AcSCh](mM)") y=seq(0,0.015,0.003) axis(2,y,cex.axis=1.6) x=c(0.0,0.006,0.009,0.012,0.018,0.024,0.051,0.099,0.201,0.399,0.0624,1.251,2.5,5,10,20,30,40,50,60,70,80) axis(1,x,cex.axis=1.5) Conc M1 M2 1 0.006 0.001878293 0.002188962 2 0.006 0.001979844 0.002324205 3 0.006 0.001923994 0.002458228 4 0.009 0.002777392 0.002851871 5 0.009 0.002666793 0.003053945 6 0.009 0.002961166 0.002890259 7 0.012 0.003439677 0.003595863 8 0.012 0.003258286 0.003884679 9 0.012 0.003643608 0.004078671 10 0.018 0.004100375 0.004927258 11 0.018 0.004246799 0.004782028 12 0.018 0.004502533 0.005023878 13 0.024 0.005001329 0.005659239 14 0.024 0.005206655 0.005692063 15 0.024 0.005172655 0.005850228 16 0.051 0.006123975 0.006956247 17 0.051 0.005811863 0.007316030 18 0.051 0.006167323 0.006911872 19 0.099 0.006962573 0.007377967 20 0.099 0.007139235 0.008086381 21 0.099 0.006480800 0.008198115 22 0.201 0.007123331 0.008603852 23 0.201 0.007168465 0.008660410 24 0.201 0.007407053 0.008852498 25 0.399 0.006926494 0.008316612 26 0.399 0.007173567 0.008623142 27 0.399 0.007018680 0.008502061 28 0.624 0.006472567 0.008146881 29 0.624 0.007144195 0.008449461 30 0.624 0.006987643 0.007874225 31 1.251 0.005747177 0.007065195 32 1.251 0.006040083 0.007266183 33 1.251 0.006227965 0.006856509 34 2.500 0.006262694 0.007632029 35 2.500 0.006196804 0.007620238 36 2.500 0.006325028 0.007639225 37 5.000 0.006178142 0.007060139 38 5.000 0.006097117 0.006880805 39 5.000 0.006209771 0.006781390 40 10.000 0.005537297 0.006051157 41 10.000 0.004966642 0.006025582 42 10.000 0.005145429 0.006104559 43 20.000 0.004470247 0.005349909 44 20.000 0.004525358 0.005015830 45 20.000 0.004317146 0.005267166 46 30.000 0.004269059 0.004942385 47 30.000 0.003840345 0.004537099 48 30.000 0.004240015 0.004720365 49 40.000 0.003817233 0.004359037 50 40.000 0.003796501 0.004416372 51 40.000 0.003692410 0.004425881 52 50.000 0.003174361 0.003873741 53 50.000 0.003279383 0.003953878 54 50.000 0.003280700 0.004253781 55 60.000 0.002868184 0.003442964 56 60.000 0.003055768 0.003253154 57 60.000 0.002782196 0.003396175 58 70.000 0.002111723 0.003052187 59 70.000 0.002271047 0.002656482 60 70.000 0.002223503 0.002681296 61 80.000 0.001827056 0.001720578 62 80.000 0.001311233 0.001930022 63 80.000 0.001532169 0.001639464 -- *************************************************** STEPHAN M. CARVALHO "Honeybee protection and ecotoxicology" Universidade Estadual Paulista Centro de Estudo de Insetos Sociais Av. 24A, 1515 Bela Vista 13.506-900 Rio Claro/SP Brazil smalfitano(a)uol.com.br msn: pefasmc(a)hotmail.com / Skype: smalfitano +55 19 3526-4284 (university) +55 19 3534-8523 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! ***************************************************

1 0

Intervalo de Confiança para Binomial Negativa...
by andrebvs 05 Set '12

05 Set '12

5 6

RgoogleMaps e rgdal
by Mauricio Cardeal 05 Set '12

05 Set '12

Olá caros mestRes ! Estou reenviando essa mensagem porque não tenho certeza de que ela chegou à lista, uma vez que não recebi no meu email a notificação. Estou tentando plotar alguns pontos (latitude e longitude) no mapa de Monte Santo, Bahia, para uma pesquisa em genética médica, através do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo, apesar de as coordenadas estarem de acordo com as de Monte Santo. O código não apresenta erro, mas o resultado não está certo. Uso o Ubuntu 12.04, o R 2.15.1 e o ESS/Emacs. O que estou fazendo errado ? Segue o CMR: #bibliotecas library(RgoogleMaps) library(rgdal) #dados de latitude bdlt <- c(-43.565,-43.564,-43.563,-43.565,-43.564) #dados de longitude bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474) dados <- as.data.frame(cbind(bdlt,bdlg)) names(dados) <- c("lat","lon") attach(dados,warn.conflicts=F) dados dados$size <- "tiny" dados$col <- "red" dados$char <- "" #get the bounding box: bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"]) #download the map: mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, markers = dados); #determine the max zoom, so that all points fit on the plot zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); plot.new() PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"], cex=1.5,pch=20,col=c('blue', 'green', 'red'),add=F); Agradeço a ajuda, Mauricio Cardeal Complexo Hospitalar Universitário Prof. Edgard Santos Pós-graduação Processos Interativos dos Órgãos e Sistemas Instituto de Ciências da Saúde Universidade Federal da Bahia

2 4

Re: [R-br] RgoogleMaps e rgdal
by Elias T. Krainski 04 Set '12

04 Set '12

Troquei os pacotes. É necessário o pacote plotGoogleMaps. Corrigindo: y <- c(-43.565,-43.564,-43.563,-43.565,-43.564) x <- c(-10.518,-10.466,-10.412,-10.516,-10.474) require(sp) xy <- SpatialPoints(cbind(xy)) require(plotGoogleMaps) plotGoogleMaps(xy) Elias T. Krainski >________________________________ > De: Mauro Sznelwar <sznelwar(a)uol.com.br> >Para: Elias T. Krainski <eliaskrainski(a)yahoo.com.br> >Enviadas: Terça-feira, 4 de Setembro de 2012 22:06 >Assunto: Re: [R-br] RgoogleMaps e rgdal > > > >Tentei rodar e veja o resultado: > >> plotGoogleMaps(x,y) >Erro: não foi possível encontrar a função "plotGoogleMaps" > >Veja se isso serve: >> >> >>y <- c(-43.565,-43.564,-43.563,-43.565,-43.564) >> >>x <- c(-10.518,-10.466,-10.412,-10.516,-10.474) >> >> >>require(sp) >>xy <- SpatialPoints(cbind(xy)) >> >> >> >>require(RgoogleMaps) >>plotGoogleMaps(xy) >> >> >>Att. >>Elias T. Krainski >> >> >>>________________________________ >>> De: Mauricio Cardeal <mcardeal2010(a)gmail.com> >>>Para: r-br(a)listas.c3sl.ufpr.br >>>Enviadas: Terça-feira, 4 de Setembro de 2012 15:36 >>>Assunto: [R-br] RgoogleMaps e rgdal >>> >>> >>>Olá caros mestRes ! Estou reenviando essa mensagem porque não tenho certeza de que ela chegou à lista, >>>uma vez que não recebi no meu email a notificação. >>> >>>Estou tentando plotar alguns pontos (latitude e longitude) no mapa de Monte Santo, Bahia, >>>para uma pesquisa em genética médica, através >>>do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo, >>>apesar de as coordenadas estarem de acordo com as de Monte Santo. >>>O código não apresenta erro, mas o resultado não está certo. >>> >>>Uso o Ubuntu 12.04, o R 2.15.1 e o ESS/Emacs. >>> >>>O que estou fazendo errado ? Segue o CMR: >>> >>>#bibliotecas >>>library(RgoogleMaps) >>>library(rgdal) >>> >>>#dados de latitude >>>bdlt <- c(-43.565,-43.564,-43.563,-43.565,-43.564) >>>#dados de longitude >>>bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474) >>>dados <- as.data.frame(cbind(bdlt,bdlg)) >>>names(dados) <- c("lat","lon") >>>attach(dados,warn.conflicts=F) >>>dados >>> >>>dados$size <- "tiny" >>>dados$col <- "red" >>>dados$char <- "" >>>#get the bounding box: >>>bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"]) >>> >>>#download the map: >>>mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, >>>markers = dados); >>> >>>#determine the max zoom, so that all points fit on the plot >>>zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); >>>plot.new() >>>PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"], >>> cex=1.5,pch=20,col=c('blue', 'green', 'red'),add=F); >>> >>>Agradeço a ajuda, >>> >>>Mauricio Cardeal >>>Complexo Hospitalar Universitário Prof. Edgard Santos >>>Pós-graduação Processos Interativos dos Órgãos e Sistemas >>>Instituto de Ciências da Saúde >>>Universidade Federal da Bahia >>> > >

1 0

Re: [R-br] Erro GVim+R
by Roney Fraga 04 Set '12

04 Set '12

Jakson, Tinha instalado o vim-r-plugin pelo yaourt do Arch Linux, que é uma instalação via terminal, simples e fácil. Pressupondo que o erro poderia ter origem na forma de instalação, resolvi seguir passo a passo o material elaborado por você e por Jose Claudio Faria, que por sinal está muito bem feito, mas o erro continua. Também tentei usar o vim pelo terminal, mas não resolveu. Talvez não exclui todos os arquivos da primeira instalação, basta substituir os arquivos em ~/.vim ou eles ficam salvos em outro lugar e devo deletá-los? Ao selecionar várias linhas e executá-las com '\ss' as vezes aparece no console do R: source("/tmp/r-plugin-roney/Rsource-9742-JD.R") Obrigado pela ajuda Roney

2 1

Dúvida distribuição gamma incompleta
by Leonardo Monteiro 04 Set '12

04 Set '12

Boa tarde, tenho uns dados de chuva, e quero ajustá-los pela distr. gamma incompleta ( que considera os valores zero também). Alguem saberia de alguma função (ou alternativa) para fazer isso?? obrigado, -- Leonardo Amaral Monteiro - Eng° Agrônomo Mestre em "Física do Ambiente Agrícola" - Agrometeorologia Doutorando do PPG em "Engenharia de Sistemas Agrícolas" Escola Superior de Agricultura "Luiz de Queiroz" - Universidade de São Paulo Piracicaba,SP. Fone: (19)34294123 - Ramal:236

3 2

Séries Temporais dividir dados em teste e treino
by Hugo Sena 04 Set '12

04 Set '12

Viva, Sou novato em R e datamining pelo que preciso muito de ajuda e agradeço a vossa paciência. Estou a fazer uma análise de séries temporais com rminer e vou precisar de dividir a série em 2 conjuntos: de treino (in-samples) e de teste (out-samples). O objetivo é fazer a análise nos dados de treino sem nunca olhar para os dados de teste. Como é feita esta divisão? É suposta ter 2 ficheiros distintos ou faço o read de todos os dados e separam-se de alguma forma no R e como?? Obrigado. Cumps

2 1

Re: [R-br] Somar vendas por semana e por mês
by Hugo Sena 04 Set '12

04 Set '12

Boas. Obrigado pela dica, penso que a função aggregate() resolve mas não consigo utilizar. Tenho um ficheiro CSV e o que pretendo é em cada linha somar as vendas totais por SEMANA. Tenho os dados seguintes mas apenas com as colunas DATA e QTD: *SEMANA* *DATA* *QTD* 45 05-11-2003 1 46 11-11-2003 1 46 11-11-2003 1 46 11-11-2003 1 46 11-11-2003 1 46 11-11-2003 1 46 11-11-2003 1 47 20-11-2003 1 48 26-11-2003 1 48 27-11-2003 1 E o que pretendo é somar as quantidades em cada dia de forma a ficar igualmente com as colunas DATA e QTD: *SEMANA* *DATA* *QTD* 45 05-11-2003 1 46 11-11-2003 6 47 20-11-2003 1 48 26-11-2003 2 Acrescentei a coluna SEMANA apenas para o exemplo porque é o que pretendo ter em consideração, até porque se não der para fazer de forma automática terei de fazer isto manualmente e somar para as mais de 5 000 linhas... Podem dar uma ajuda? E depois é suposto exportar esta agregação para um ficheiro? Como? Já procurei no google mas ainda não consegui fazer no meu caso em particular. Obrigado. Cumps.

4 14

Erro GVim+R
by Roney Fraga 03 Set '12

03 Set '12

Caros, Visto que o RStudio em meu computador as vezes apresenta problemas no cursor do mouse, ele fica umas 3 letras a frente de onde deveria estar, tanto no Linux quanto no Windows, estou migrando para o GVim-R-plugin. Achei a curva de aprendizado pequena e gostei muito do resultado, contudo, tenho 2 dúvidas. Primeira, as vezes ocorre um erro quando executo muitas linhas de comando de uma única vez, daí a janela do terminal é fechada, sempre utilizo os comandos '\r' ou '\l'. Segue a mensagem de erro: [warn] Epoll ADD(1) on fd 0 failed. Old events were 0; read change was 1 (add); write change was 0 (none): Operation not permitted [warn] Epo ll ADD(1) on fd 0 failed. Old events were 0; read change was 1 (add); write change was 0 (none): Operation not permitted failed to connect to server failed to connect to server failed to connect to server failed to connect to server failed to connect to server A segunda dúvida é quanto a possibilidade de mudar a tecla padrão dos comandos gvim+r que é '\' para 'AltGr' ou para a tecla que representa o botão direito do mouse. Meu sistema operacional é o Arch Linux 64 bits com R 2.15.1. Agradeço a ajuda Roney

2 1

Extrair elementos comuns
by clara luz 03 Set '12

03 Set '12

Oi gente! pergunta rápida e idiota: tenho 2 vetores: s<- 4467 4468 4469 4470 4471 4745 4746 4747 4748 4749 5023 5024 5025 5026 5027 5301 5302 5303 5304 5305 5579 5580 5581 5582 5583 sn<- 4746 4747 4748 5024 5025 5026 5302 5303 5304 preciso criar outro vetor, 'se' com os valores que tem em 's', mas que não tem em 'sn'. como eu extraio esses elementos de 's'?? se é 1 único valor até rola, mas com vários valores não sei como fazer.... valeu galera!! Clara Luz B. Sant'Anna Bióloga, Bacharel e Licenciada Mestranda em Ecologia - Unicamp Laboratório de Biogeografia da Conservação e Mudanças Climáticas- UFG (19) 9233.1062

4 4

[Dúvida] Reordenar as linhas de uma matriz de acordo com a sequência de uma coluna
by Carolina R. Caiado Gomes 03 Set '12

03 Set '12

Olá Pessoal, Minha dúvida é a seguinte: tenho uma matriz em que na primeira coluna tenho coordenadas de longitude, na segunda coluna tenho coordenadas de latitude e nas demais colunas tenho espécies, que podem estar presentes (valor 1) ou ausentes (valor 0) na célula (par de longitude+latitude), como no exemplo abaixo... a<-c(2,5,0,0,1) b<-c(2,1,0,1,1) c<-c(2,3,0,1,0) d<-c(2,2,1,1,0) e<-c(2,4,1,0,1) matriz<-rbind(a,b,c,d,e) colnames(matriz)<-c('longitude','latitude','sp1','sp2','sp3') Nesse caso, todas as minhas células estão na mesma longitude, mas em latitudes diferentes. Quero reordenar minha matriz de modo que na primeira linha fique a menor latitude (no ex:1) com sua longitude (no ex:2) e espécies correspondentes (no ex: 0,1,1) , na segunda linha fique a segunda menor latitude (no ex: 2) com sua longitude (no ex: 2) e espécies correspondentes (no ex: 1,1,0) e assim sucessivamente... Enfim, quero reordenar todas as linhas da minha matriz de acordo com uma ordem crescente dos valores da coluna de latitudes. Alguém sabe como fazer isso? (Segue exemplo de como quero que a matriz fique...) matriz_ordenada<-rbind(b,d,c,a,e) colnames(matriz_ordenada)<-c('longitude','latitude','sp1','sp2','sp3') Espero ter sido clara e que alguém possa me ajudar... Obrigada pela atenção, -- Carolina Ramos Caiado Gomes Mestranda em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos +55 62 8215 3595

2 1

Rstudio n Ubuntu 12.04
by Andre Oliveira 03 Set '12

03 Set '12

Baixei o RStudio aqui em meu netbook com Ubuntu 12.04. Quando clico para instalar vai abrir o gerenciador de pacotes, mas diz que não pode abrir o pacote. Já baixei duas vezes e nada. Alguma dica da possível causa? André Oliveira Souza

3 3

sobre Revolution R Enterprise
by clara luz 03 Set '12

03 Set '12

Olá pessoal!! Alguém conhece ou já trabalhou com o Revolution R Enterprise?? Se puder dar umas dicas.... aparentemente é bastante semelhante ao R clássico, mas não estou conseguindo nem importar uma tabela... rs Valeu!! Clara Luz B. Sant'Anna Bióloga, Bacharel e Licenciada Mestranda em Ecologia - Unicamp Laboratório de Biogeografia da Conservação e Mudanças Climáticas- UFG (19) 9233.1062

4 3

[OFF] Scientific Linux
by Pedro Rafael 02 Set '12

02 Set '12

Pessoal alguém já ouviu falar no Scientific Linux? (http://en.wikipedia.org/wiki/Scientific_Linux) Vocês sabem informar se essa distrbuição é uma boa distribuição Linux? [ ], Pedro Rafael Diniz Marinho.

2 2

Plotar pontos no mapa de Monte Santo - Bahia com as bibliotecas
by Mauricio Cardeal 02 Set '12

02 Set '12

Desculpem, não havia escrito o carregamento das bibliotecas Olá caros colegas ! Estou tentando plotar alguns pontos no mapa de Monte Santo, Bahia através do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo. Por favor, o que estou fazendo errrado. Segue o CMR: #bibliotecas library(RgoogleMaps) library(rgdal) #dados de latitude bdlt <- c(-43.565,-43.564,-43.563,-43. 565,-43.564) #dados de longitude bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474) dados <- as.data.frame(cbind(bdlt,bdlg)) names(dados) <- c("lat","lon") attach(dados,warn.conflicts=F) dados dados$size <- "tiny" dados$col <- "red" dados$char <- "" #get the bounding box: bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"]) #download the map: mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, markers = dados); #determine the max zoom, so that all points fit on the plot zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); plot.new() PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"], cex=1.5,pch=20,col=c('blue', 'green', 'red'), add=F); Agradeço a ajuda

1 0

Plotar pontos no mapa de Monte Santo - Bahia
by Mauricio Cardeal 02 Set '12

02 Set '12

Olá caros colegas ! Estou tentando plotar alguns pontos no mapa de Monte Santo, Bahia através do Rgoogle e rgdal, mas o mapa que aparece ao final não é o de Monte Santo. Por favor, o que estou fazendo errrado. Segue o CMR: #dados de latitude bdlt <- c(-43.565,-43.564,-43.563,-43.565,-43.564) #dados de longitude bdlg <- c(-10.518,-10.466,-10.412,-10.516,-10.474) dados <- as.data.frame(cbind(bdlt,bdlg)) names(dados) <- c("lat","lon") attach(dados,warn.conflicts=F) dados dados$size <- "tiny" dados$col <- "red" dados$char <- "" #get the bounding box: bb <- qbbox(lat = dados[,"lat"], lon = dados[,"lon"]) #download the map: mapa <- GetMap.bbox(bb$lonR, bb$latR, destfile = "dados.png", GRAYSCALE =F, markers = dados); #determine the max zoom, so that all points fit on the plot zoom <- min(MaxZoom(latrange=bb$latR,lonrange=bb$lonR)); plot.new() PlotOnStaticMap(mapa,lat = dados[,"lat"], lon = dados[,"lon"], cex=1.5,pch=20,col=c('blue', 'green', 'red'), add=F); Agradeço a ajuda Mauricio Cardeal UFBA

1 0

Tutorial do pacote Hmisc
by Carlos Andrade 01 Set '12

01 Set '12

Olá, pessoal da lista. Alguém referência de links de tutorias para o pacote Hmisc? -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

2 1

Duvida de estatistica - GLM binomial
by Augusto Ribas 01 Set '12

01 Set '12

Ola pessoal. Eu estava olhando aqui uns dados com um amigo e me surgiu uma duvida. O problema é o seguinte: Os dados são referentes a ocorrência de uma doença em passarinhos. Pegamos muitos passarinhos em redes em um lugar e testamos se ele esta doente ou não. Durante a captura, registramos varias caracteristicas dos passarinhos, entre elas o sexo e a idade. Ai gostaríamos de perguntar se essas variáveis influenciam na chance do passarinho estar doente ou não. O problema começa no seguinte, a idade é fácil de determinar, são 3 classes, jovem, sub-adulto e adulto e isso é feito sem muito erro, Agora o sexo, era pra ser macho ou femea, mas quando o passarinho é muito novo, é impossível saber o sexo dele, então a gente tava colocando indeterminado, pq visualmente não da pra saber o sexo do passarinho, e não tem outro jeito. Então eu tenho 2 variáveis, de 3 níveis cada uma e fiz um glm binomial. So que qd eu faço a analise o sexo da significativo (a classe indeterminado), sendo que os indeterminados tem menor chance de estar doente. Mas na verdade eu acho que o que interessa é a idade (Que faz sentido biológico ser significativo). Se vc olha o gráfico, vc ve o padrão é o mesmo no sexo e na idade. E a inabilidade de identificar o sexo nos passarinhos vem exatamente dele ser jovem demais, não desenvolveu ainda. Ai uma coisa confunde com a outra. Será que alguem tem aluma sugestão de o que fazer? Tem como representar o sexo como "Não sei", é um dos 2 níveis mas pra esse passarinho eu não sei? Eu não sei como fazer isso. É possível fazer isso pro sexo no glm no meu caso? Segue um exemplo para ilustrar mais ou menos o que eu to vendo: #Gerando dados de exemplo set.seed(666) n<-90 #eu acredito que a idade tem um efeito significativo #ja que por exemplo, o cara muito jovem vai ter tempo de exposição menor #a indivíduos doentes idade<-sample(c("Adulto","Sub-Adulto","Jovem"),n,replace=T) idade efeito<-ifelse(idade=="Jovem",1,0) efeito diagnostico<-rbinom(n,1,0.75+efeito*(-0.25)) #E o sexo do passarinho não tem efeito, mas os jovens sempre são inderterminados sexo<-NA for(i in 1:length(idade)) { if(idade[i]=="Jovem") { sexo[i]<-c("Indeterminado") } else{ sexo[i]<-sample(c("Macho","Femea"),1) } } dados<-data.frame(sexo,idade,diagnostico) #so completando os casos dados<-rbind(dados,c("Femea","Jovem",0)) dados<-rbind(dados,c("Indeterminado","Adulto",0)) dados<-rbind(dados,c("Macho","Jovem",0)) dados<-rbind(dados,c("Indeterminado","Sub-Adulto",0)) #amostras table(dados$sexo,dados$idade) head(dados) #Grafico #no final o que eu vejo é mais ou menos isso: par(mfrow=c(1,2),cex=0.6) barplot(table(diagnostico,sexo),ylim=c(0,40),legend.text =c("Parasitado","Não Parasitado")) barplot(table(diagnostico,idade),ylim=c(0,40),legend.text =c("Parasitado","Não Parasitado")) #O indeterminado e Jovem é mais meio a meio no numero de parasitados que as outras classes #mas eu vejo um resultado significativo pro sexo, ou jovem se mudar a ordem #isso esta errado modelo01<-glm(diagnostico~factor(idade)+factor(sexo),family="binomial") summary(modelo01) plogis(coef(modelo01)[1:2]) #eu acredito que cai naquele problema do aov, de experimento não balanceado e tal #se olhar individualmente tudo ok, mas não consigo olhar as 2 variaveis juntas. modelo02<-glm(diagnostico~factor(idade),family="binomial") summary(modelo02) plogis(coef(modelo02)[1:2]) modelo03<-glm(diagnostico~factor(sexo),family="binomial") summary(modelo03) plogis(coef(modelo03)[1:2]) Se alguem puder dar uma olhada e dar alguma sugestão, fiz um exemplo, mas os dados deixam um grafico exatamente como aquele, se alguem quiser olhar os dados originais eu posso mandar. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

6 16

Mais códigos HTML e função readlines
by Paulo Nogueira 01 Set '12

01 Set '12

Pessoal, não sei se a lista se presta a esclarecer dúvidas como essa minha, mas pelo menos estou usando R para resolver os problemas. rs Ultimamente tenho trabalhado bastante com determinadas buscas na web para coletar dados e estruturá-los. Porém, como meus conceitos de web se restringem ao conhecimento de usuário, e não de desenvolvedor, acabo tendo alguns problemas. Vou exemplificar uma das minhas dificuldades com um problema prático: Utilizar readlines() e os pacores RCurl, XML, etc, para tratar o código fonte é algo que já faço até que bem. Mas problemas surgem quando o conteúdo que eu preciso extrair de informação não vem no código fonte? Vejam esse exemplo: http://www.skyscanner.com.br/passagens-aereas/saoa/slz/121002/130409/tarifa… Após carregar toda a busca, o site SkyScanner exibe os preços encontrados para o trecho São Paulo - São Luis. [1] No navegador Firefox ou Chrome, ao clicar em "Exibir Código Fonte", não aparecem as informações do vôos, apenas o código que estrutura a página. [2] Ao clicar em "Inspecionar Elemento" eu consigo analisar o código fonte do meu interesse. [3] Ao utilizar readlines passando o a url acima o R importa o código fonte [1] e não o [2], com as informações de interesse. Imagino que isso tenha a ver com as funções de javascript e coisas desse tipo... Alguém sabe contornar essa situação com o R? Se não, o que devo estudar para agilizar a busca dessa solução? Estou meio perdido nessa parte e agradeço ajudas de qualquer tipo. rs Peço desculpas se por acaso isso for muito off-topic. Abraços, Paulo Nogueira

4 5

Ajuda com o Gedit
by Luciano F. Sgarbi 01 Set '12

01 Set '12

Olá pessoal, Alguém ai usa o gedit no R, com o sistema operacional Windows 7? Estou tentando faze-lo funcionar mas não tem jeito... Se alguém der um help fico muito agradecido. Abraços Luciano -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

2 2

Bug do Rcmdr?
by Emerson Cotta Bodevan 31 Ago '12

31 Ago '12

Prezados, recorro aos senhores para mais uma dúvida. No momento estou utilizando o R para windows... Instalei o R e depois o Rcmdr. Após a instalação do R criei um atalho na barra de ferramentas. Estou com o seguinte problema: - ao iniciar o R deste atalho e carregar o Rcmdr pelo comando (require(Rcmdr)) aparece no menu *Dados *a opção de *importar arquivos *e lá a opção de *arquivos excel, access, etc*. - ao iniciar o R por um atalho em um subdiretório de um projeto esta opção *arquivos excel, access, etc*. não aparece. Alguém sabe como resolver? Agradeço qualquer ajuda. Abs, Emerson. -- "A educação não muda o mundo. A educação muda as pessoas. As pessoas mudam o mundo." * Paulo Freire.*

2 2

[Off-topic] Corretor ortográfico pt_br Kile 2.1
by Andre Oliveira 31 Ago '12

31 Ago '12

Pessoal como posso colocar corretor ortográfico para o Kile 2.1 no Ubuntu 12.04? obrigado André Oliveira Souza

1 0

Manual sobre o pacote tcltk
by Ary Casado 31 Ago '12

31 Ago '12

Pessoal Estou desenvolvendo um pequeno sistema com interfaces gráficas para o meu trabalho de dissertação. Escolhi o pacote tcltk para esse trabalho porque achei muito fácil de lidar com ele, porém quanto mais escrevo e crio forms mais aparecem funções, eventos, propriedades e métodos. Alguém conhece algum manual bem didático sobre o tcltk???? O site www.sciviews.org/rgui já me ajudou bastante com vários exemplos, porém não explica tudo. Arykerne Nascimento Mestrando Proten/UFPE

1 0